• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    瑪湖凹陷風(fēng)城組頁巖油地質(zhì)工程甜點地震預(yù)測方法及應(yīng)用

    2022-11-21 11:10:52于江龍陳剛吳俊軍李維楊森唐廷明
    新疆石油地質(zhì) 2022年6期
    關(guān)鍵詞:白云質(zhì)巖相儲集層

    于江龍,陳剛,吳俊軍,李維,楊森,唐廷明

    (中國石油 新疆油田分公司 勘探開發(fā)研究院,新疆 克拉瑪依 834000)

    近年來,隨著非常規(guī)油氣勘探的不斷深入,頁巖油已成為非常規(guī)油氣勘探的重點領(lǐng)域。目前,準噶爾盆地吉木薩爾凹陷二疊系蘆草溝組頁巖油勘探已獲得進展,瑪湖凹陷下二疊統(tǒng)風(fēng)城組頁巖油正成為準噶爾盆地陸相頁巖油勘探新的潛力領(lǐng)域[1]?,敽枷蒿L(fēng)南7 井、風(fēng)南14 井、瑪頁1 井等井均在風(fēng)城組獲得工業(yè)油流,展現(xiàn)了準噶爾盆地下二疊統(tǒng)風(fēng)城組陸相頁巖油豐富的資源潛力。

    甜點預(yù)測對頁巖油勘探開發(fā)至關(guān)重要,所謂甜點,即在致密儲集層背景下,經(jīng)壓裂改造后能夠獲得較高工業(yè)產(chǎn)能的優(yōu)質(zhì)儲集層發(fā)育區(qū)[2]。近年來,頁巖油甜點的主控地質(zhì)因素研究得到越來越多的關(guān)注。在主控因素分析基礎(chǔ)上,頁巖油甜點的評價與優(yōu)選成為生產(chǎn)中需要解決的關(guān)鍵問題。頁巖油甜點的評價與優(yōu)選主要采用地質(zhì)評價、地震預(yù)測和地震地質(zhì)一體化。頁巖油甜點的地質(zhì)評價主要以地質(zhì)資料為依據(jù),從烴源巖品質(zhì)、儲集層品質(zhì)、脆性指數(shù)等方面開展評價工作,如利用烴源巖特征、巖性、物性、含油性、電性、脆性和地應(yīng)力各向異性進行頁巖油甜點評價的七性評價方法,是目前中國頁巖油評價的主要方法[3];頁巖油甜點的地震預(yù)測主要依據(jù)地震信息,通過地震巖石物理分析優(yōu)選出地震彈性參數(shù),再運用地震彈性參數(shù)反演,獲得頁巖油甜點分布的預(yù)測結(jié)果[4];地震地質(zhì)一體化的頁巖油甜點評價方法以地震信息為依據(jù),以地質(zhì)信息為約束條件,建立頁巖油甜點評價與優(yōu)選標準[5]。

    針對瑪湖凹陷風(fēng)城組甜點預(yù)測,前人已做過很多工作,如對風(fēng)城組頁巖油發(fā)育的主控因素進行分析總結(jié),認為風(fēng)城組發(fā)育優(yōu)質(zhì)烴源巖,其巖性、物性及裂縫發(fā)育程度共同決定了頁巖油甜點分布,聯(lián)合地震多屬性、波阻抗反演、隨機模擬反演等技術(shù)預(yù)測了風(fēng)三段優(yōu)勢甜點發(fā)育區(qū)[6];利用OVT 域偏移數(shù)據(jù)進行風(fēng)三段白云質(zhì)巖儲集層預(yù)測,認為甜點主要發(fā)育于儲集層厚度大、脆性強及高壓的疊合區(qū)域,運用不同方位角疊前偏移數(shù)據(jù)預(yù)測白云質(zhì)巖儲集層厚度及脆性,利用有效應(yīng)力與縱橫波速度關(guān)系預(yù)測異常壓力區(qū)[7];針對瑪湖凹陷風(fēng)二段地震振幅較弱且連續(xù)性差導(dǎo)致甜點預(yù)測異常的現(xiàn)象,采用基于屬性建模的地震正反演聯(lián)合方法,證實了存在2 個甜點,并利用分頻均方根振幅屬性預(yù)測了甜點的發(fā)育范圍[8]。前人傾向于地質(zhì)甜點預(yù)測,對工程甜點預(yù)測還需進一步探索。

    針對準噶爾盆地瑪湖凹陷風(fēng)城組頁巖油儲集層,首先明確了儲集層地質(zhì)甜點與工程甜點的特征,總結(jié)了甜點發(fā)育的主控因素,優(yōu)選巖性、脆性指數(shù)及水平主應(yīng)力差作為地質(zhì)工程甜點預(yù)測的主控因素,利用疊前同時反演技術(shù),結(jié)合巖相概率分析、Rickman 脆性指數(shù)法及組合彈簧模型,建立了一套以疊前同時反演為基礎(chǔ)的地質(zhì)工程甜點預(yù)測流程,預(yù)測了瑪湖凹陷風(fēng)三段底部地質(zhì)工程甜點有利區(qū)。

    1 區(qū)域地質(zhì)概況

    瑪湖凹陷位于準噶爾盆地西北緣,面積近5 000 km2,東南側(cè)分布有英西凸起、夏鹽凸起和達巴松凸起,西南接中拐凸起,西側(cè)為克百斷裂帶,北部為烏夏斷裂帶(圖1)。研究區(qū)位于瑪湖凹陷北部,東南鄰近瑪湖凹陷沉積中心,處于烏夏斷裂帶與瑪湖凹陷斜坡區(qū),目的層下二疊統(tǒng)風(fēng)城組埋深為4 000~8 000 m。發(fā)育一套沉積于堿湖環(huán)境,以半封閉深湖—半深湖沉積為主,受機械、化學(xué)、火山等沉積作用的混積巖,巖石類型復(fù)雜多樣,總沉積厚度為800~1 800 m。風(fēng)城組是準噶爾盆地瑪湖凹陷的主力烴源巖層,自下而上劃分為風(fēng)一段、風(fēng)二段和風(fēng)三段,烴源巖全層段發(fā)育,風(fēng)一段、風(fēng)三段發(fā)育粉砂巖類和白云質(zhì)巖類甜點,風(fēng)二段發(fā)育白云質(zhì)巖類甜點[9],風(fēng)一段局部發(fā)育火山巖。

    圖1 研究區(qū)構(gòu)造位置Fig.1.Structural location of the study area

    2 甜點特征

    瑪湖凹陷風(fēng)城組儲集層巖性復(fù)雜,白云質(zhì)巖、碎屑巖、火山巖等多種巖石類型共存,主要為陸相碎屑巖、爆發(fā)相火山巖和湖盆內(nèi)化學(xué)沉積碳酸鹽巖的細?;旆e巖。風(fēng)城組沉積中心在研究區(qū)西部,為深湖—半深湖相沉積環(huán)境,水體安靜,碳酸鹽巖類較發(fā)育,是非常規(guī)頁巖儲集層的集中發(fā)育區(qū)。研究區(qū)東部為扇三角洲沉積,是常規(guī)碎屑巖儲集層的集中發(fā)育區(qū)。由于裂隙式火山噴發(fā),風(fēng)城組底部局部發(fā)育爆發(fā)相火山巖。通過巖心觀察、巖石薄片鑒定等方法,將研究區(qū)風(fēng)城組儲集層巖性劃分為火山巖、砂礫巖、白云質(zhì)粉砂巖、白云質(zhì)細粉砂巖和泥巖。

    風(fēng)城組儲集層致密,孔隙度為0.10%~13.60%,平均為4.61%;滲透率為0.010~13.800 mD,平均為0.134 mD,為典型的頁巖油儲集層特征,頁巖油的開發(fā)很大程度上依賴于大規(guī)模的壓裂改造[10]。巖石脆性是頁巖油儲集層壓裂改造需要考慮的重要巖石力學(xué)參數(shù)之一,脆性越好,形成的裂縫越復(fù)雜,從而能更有效地改造儲集層[11-12]。水平主應(yīng)力差是影響壓裂裂縫形態(tài)的重要因素,水平主應(yīng)力差越小,壓裂后越易形成網(wǎng)狀裂縫,有利于水平井壓裂改造[13]。

    綜合以上分析,頁巖油儲集層物性普遍較差,白云質(zhì)粉砂巖含油性最好,普遍以油浸和油斑為主,因此,白云質(zhì)粉砂巖的分布范圍,可認為是地質(zhì)甜點的分布范圍。在確定地質(zhì)甜點分布范圍的基礎(chǔ)上,為有效改造儲集層,需進行脆性指數(shù)和水平主應(yīng)力差的預(yù)測,優(yōu)選脆性指數(shù)高且水平主應(yīng)力差小的區(qū)域,也就是工程甜點的范圍。對研究區(qū)頁巖油儲集層地質(zhì)工程甜點進行預(yù)測,以瑪頁1 井為例,目的層位于風(fēng)三段底部,深度為4 593—4 617 m,厚度為24 m,在地震剖面上位于強波谷的下半部分與波峰的上半部分,強波谷同相軸橫向分布穩(wěn)定,連續(xù)性好,可識別且可連續(xù)追蹤。

    3 甜點預(yù)測方法

    3.1 疊前同時反演技術(shù)

    疊前同時反演是將疊前共反射點道集按照遠、中、近偏移距進行部分疊加,通過求解Zoeppritz 方程的近似公式,得到縱波速度、橫波速度、波阻抗等參數(shù),間接得到楊氏模量、泊松比等參數(shù)[14]。為了避免疊前反演中密度不準確的問題,分別基于Zoeppritz方程的Aki-Richards近似式及Fatti近似式作2次疊前同時反演,Aki-Richards近似式:

    3.2 地質(zhì)甜點預(yù)測方法

    3.2.1 巖相對比分析

    測井數(shù)據(jù)頻率可達幾萬赫茲,測井解釋巖相分辨率高,地震數(shù)據(jù)的頻率(10~80 Hz)限制了地震預(yù)測巖相的尺度。地震預(yù)測巖相是測井解釋巖相的綜合反映,在一定深度范圍內(nèi)若某種測井解釋巖相累計厚度最大或出現(xiàn)概率最高,地震預(yù)測巖相則表現(xiàn)為該巖相(圖2)。研究區(qū)瑪頁1 井在4 593~4 617 m,測井解釋巖相包含白云質(zhì)粉砂巖、白云質(zhì)細粉砂巖和砂礫巖,其中,白云質(zhì)細粉砂巖累計厚度最大,地震預(yù)測巖相表現(xiàn)為白云質(zhì)細粉砂巖。

    圖2 研究區(qū)測井解釋巖相與地震預(yù)測巖相對比Fig.2.Lithofacies derived from logging interpretation and seismic prediction in the study area

    3.2.2 巖石物理分析

    通過鉆井資料巖石物理分析,對比縱橫波速度比、泊松比、楊氏模量、體積模量和縱波阻抗對瑪湖凹陷風(fēng)城組巖性的識別效果,其中,縱橫波速度比與縱波阻抗(圖3a)和泊松比與縱波阻抗(圖3b)交會圖對巖性的識別效果較好,楊氏模量與縱波阻抗(圖3c)和體積模量與縱波阻抗(圖3d)交會圖對巖性的識別效果略差,考慮到由縱橫波速度比得到泊松比誤差的累加與運算量的增加,優(yōu)選縱橫波速度比和縱波阻抗作為識別巖性的敏感參數(shù)。因此,可以通過縱波阻抗與縱橫波速度比交會識別風(fēng)城組優(yōu)勢巖性白云質(zhì)粉砂巖,這也體現(xiàn)了疊前同時反演在研究區(qū)風(fēng)城組頁巖油儲集層巖性識別中的適用性。

    圖3 研究區(qū)巖性識別交會圖Fig.3.Cross-plots of lithology identification in the study area

    3.2.3 巖相概率分析

    巖相概率分析技術(shù)是基于疊前同時反演獲得的多個彈性參數(shù)體,在鉆井資料的基礎(chǔ)上,利用概率密度函數(shù)結(jié)合彈性參數(shù)反演結(jié)果,采用貝葉斯判定公式實現(xiàn)從彈性參數(shù)體轉(zhuǎn)化為含地質(zhì)信息的概率數(shù)據(jù)體,是將反演的彈性參數(shù)轉(zhuǎn)化為地質(zhì)結(jié)果的主要手段之一。針對瑪湖凹陷風(fēng)城組頁巖油儲集層,考慮巖石物理分析結(jié)果,優(yōu)選縱波阻抗和縱橫波速度比作為巖相概率分析的彈性參數(shù)體。利用測井解釋巖相確定各個巖相先驗概率分布(圖4a),通過調(diào)整概率密度函數(shù)(圖4b),使不同巖相得到區(qū)分,白云質(zhì)粉砂巖與白云質(zhì)細粉砂巖在交會圖中部分重疊,調(diào)整巖性概率分析范圍,使白云質(zhì)粉砂巖與白云質(zhì)細粉砂巖得到區(qū)分。白云質(zhì)粉砂巖縱波阻抗為14 000~16 000 g/cm3·m/s 的概率較高,平均縱波阻抗為14 940 g/cm3·m/s,標準差為317;白云質(zhì)細粉砂巖縱波阻抗為13 000~15 000 g/cm3·m/s的概率較高,平均縱波阻抗為13 930 g/cm3·m/s,標準差為246(圖4c)。白云質(zhì)粉砂巖縱橫波速度比為1.60~1.80的概率較高,平均為1.70,標準差為0.038;白云質(zhì)細粉砂巖縱橫波速度比為1.65~1.84 的概率較高,平均為1.75,標準差為0.036(圖4d)。因此,兩者雖然存在部分重疊,但主體部分可以識別,滿足巖相預(yù)測的要求。概率密度函數(shù)調(diào)整后,結(jié)合彈性參數(shù)反演結(jié)果,采用貝葉斯判定公式,對研究區(qū)風(fēng)城組巖相進行定性解釋,刻畫平面優(yōu)勢巖相分布范圍,實現(xiàn)地質(zhì)甜點的地震預(yù)測。

    圖4 研究區(qū)巖相概率分布Fig.4.Lithofacies probability distribution in the study area

    3.3 工程甜點預(yù)測方法

    3.3.1 各向異性分析

    陸相頁巖油儲集層由于沉積環(huán)境等因素的影響,其縱向上巖性變化相對頻繁,表現(xiàn)出明顯的各向異性。受頻率范圍影響,地震數(shù)據(jù)對頁巖油儲集層各向異性的反映受到限制,從正交偶極聲波測井斯通利波數(shù)據(jù)中分別提取水平方向與垂直方向的縱波速度和橫波速度,并將其濾波至80 Hz,使其與地震數(shù)據(jù)頻率接近,水平方向與垂直方向的縱波速度基本重合,橫波速度也基本重合(圖5),因此,在地震數(shù)據(jù)頻率范圍內(nèi)頁巖油儲集層的各向異性很弱,可以近似看作彈性各向同性介質(zhì)。

    圖5 研究區(qū)縱橫波速度各向異性分析Fig.5.Analysis on anisotropy of P-and S-wave velocities in the study area

    3.3.2 脆性指數(shù)預(yù)測

    頁巖儲集層與常規(guī)碎屑巖儲集層相比,具有低孔低滲特征,覆壓基質(zhì)滲透率不大于0.1 mD,單井一般無自然產(chǎn)能或自然產(chǎn)能低于工業(yè)油流下限[15],最有效的增產(chǎn)方法是通過人工壓裂造縫來溝通天然裂縫和基質(zhì)孔隙。巖石脆性是儲集層改造需要考慮的重要巖石力學(xué)參數(shù)之一,反映巖石在一定條件下形成裂縫的能力,脆性越強,形成的裂縫越復(fù)雜[16],從而有效改善頁巖油儲集層低滲情況,獲得較高的單井產(chǎn)量。材料斷裂或破壞前表現(xiàn)出極少或沒有塑性形變的特征為脆性[17]。

    脆性指數(shù)的預(yù)測主要有3 種方法[18-20]:①利用頁巖中脆性礦物含量評價頁巖脆性的方法,應(yīng)用石英含量脆性指數(shù)公式求得頁巖的脆性指數(shù);②采用巖石力學(xué)實驗方法,通過應(yīng)力—應(yīng)變曲線直接評價巖石脆性程度,如三軸應(yīng)力實驗;③利用地球物理及測井方法求取楊氏模量、泊松比等彈性力學(xué)參數(shù),進而通過Rickman 脆性指數(shù)公式求得頁巖的脆性指數(shù)。比較這3 種方法,巖石力學(xué)實驗費用昂貴且實驗周期長,所得結(jié)果不連續(xù),在實際應(yīng)用中有較大局限;脆性礦物含量評價頁巖脆性方法簡單易操作,但由于頁巖礦物的多樣性,導(dǎo)致用礦物含量來評價頁巖脆性精度低。楊氏模量、泊松比等彈性力學(xué)參數(shù)更能反映巖石的脆性,利用地震與測井方法可快速預(yù)測研究區(qū),尤其是井間的脆性變化情況。針對瑪湖凹陷風(fēng)城組頁巖油儲集層,基于疊前同時反演得到的縱橫波速度比、縱波阻抗等數(shù)據(jù)體,利用巖石物理關(guān)系運算得到楊氏模量、泊松比等彈性參數(shù)體,再運用Rickman脆性指數(shù)公式計算研究區(qū)風(fēng)城組脆性指數(shù)體,從而實現(xiàn)研究區(qū)脆性指數(shù)的地震預(yù)測。

    Rickman脆性指數(shù)公式:

    3.3.3 水平主應(yīng)力差預(yù)測

    前人提出了有物理基礎(chǔ)的確定最小主應(yīng)力的方法,該方法與泊松比有關(guān)[21],并假設(shè)在沉積中,只考慮上覆壓力的影響,地層只在垂向上產(chǎn)生應(yīng)變,水平方向應(yīng)變?yōu)?;有學(xué)者對該方法作了改進,假設(shè)地下巖層的地應(yīng)力主要為上覆壓力和水平方向構(gòu)造應(yīng)力,且水平方向構(gòu)造應(yīng)力與上覆壓力呈正比,考慮了構(gòu)造應(yīng)力的影響,但沒有考慮巖性對地應(yīng)力的影響,適用于構(gòu)造平緩區(qū)域[22-26];對于構(gòu)造運動比較劇烈的區(qū)域,水平主應(yīng)力很大部分來源于地質(zhì)構(gòu)造運動產(chǎn)生的構(gòu)造應(yīng)力,不同性質(zhì)的地層由于其抵抗外力的變形特點不同,因而其承受的構(gòu)造應(yīng)力也不相同,根據(jù)組合彈簧的構(gòu)造運動模型推導(dǎo)的分層地應(yīng)力計算模型,即組合彈簧模型[27]:

    由于瑪湖凹陷風(fēng)城組經(jīng)歷多期強烈構(gòu)造運動,選擇組合彈簧模型作為地應(yīng)力預(yù)測方法。上覆壓力是由上覆巖層重量導(dǎo)致的應(yīng)力,分析瑪北地區(qū)多口井的密度和縱波阻抗發(fā)現(xiàn),兩者相關(guān)性好,由全層系縱波阻抗反演體可轉(zhuǎn)換得到密度體,密度積分獲得上覆壓力;孔隙壓力由Eaton 公式獲得[28]。靜態(tài)泊松比和靜態(tài)楊氏模量由三軸應(yīng)力實驗獲得的線性擬合公式,通過疊前同時反演得到的動態(tài)泊松比和動態(tài)楊氏模量轉(zhuǎn)換得到,其中,靜態(tài)楊氏模量與動態(tài)楊氏模量具有顯著的線性關(guān)系(圖6a),靜態(tài)泊松比與動態(tài)泊松比的相關(guān)性略差(圖6b)。利用聲發(fā)射Kaiser 效應(yīng)實驗測得的水平主應(yīng)力計算得到研究區(qū)風(fēng)城組在最小水平主應(yīng)力和最大水平主應(yīng)力方向上的應(yīng)變系數(shù)分別為0.584 81和0.825 34。

    圖6 研究區(qū)動靜態(tài)彈性參數(shù)轉(zhuǎn)換關(guān)系Fig.6.Relationships between dynamic and static elastic parameters in the study area

    利用組合彈簧模型預(yù)測風(fēng)城組儲集層最大水平主應(yīng)力、最小水平主應(yīng)力及水平主應(yīng)力差,與脆性指數(shù)結(jié)合,實現(xiàn)研究區(qū)工程甜點的地震預(yù)測。根據(jù)以上分析結(jié)果,建立了瑪湖凹陷風(fēng)城組頁巖油儲集層地質(zhì)工程甜點預(yù)測流程(圖7)。

    圖7 研究區(qū)地質(zhì)工程甜點預(yù)測流程Fig.7.Process of predicting geological and engineering sweet spots in the study area

    4 應(yīng)用效果分析

    4.1 地質(zhì)甜點預(yù)測

    研究區(qū)風(fēng)三段底部白云質(zhì)粉砂巖分布范圍較廣,主要位于西部(圖8a),從西向東,巖性整體由白云質(zhì)粉砂巖過渡為白云質(zhì)細粉砂巖、砂礫巖,碳酸鹽巖含量逐漸減少,陸源碎屑含量逐漸增多。沉積環(huán)境控制著碳酸鹽巖的發(fā)育規(guī)模,風(fēng)城組沉積中心位于研究區(qū)西部,沉積時水體較深且安靜,早期火山噴發(fā)釋放的大量Mg2+富集,為白云石化提供了良好的環(huán)境,因此,研究區(qū)西部細粒白云質(zhì)巖類發(fā)育,東部逐漸遠離風(fēng)城組沉積中心,水體變淺,沉積環(huán)境從深湖相過渡到扇三角洲相,陸源碎屑供給增多,白云質(zhì)巖類含量減少,以碎屑巖沉積為主,巖性預(yù)測結(jié)果與地質(zhì)認識相吻合。A井處地震預(yù)測巖相與過A井巖相概率分析預(yù)測剖面吻合較好(圖8b),證實了巖相概率分析預(yù)測巖相的準確性。研究區(qū)西部白云質(zhì)粉砂巖連片發(fā)育,為地質(zhì)甜點有利區(qū),面積共135.5 km2。

    圖8 研究區(qū)風(fēng)三段底部巖相預(yù)測結(jié)果Fig.8.Prediction of lithofacies at the bottom of Feng 3 member in the study area

    4.2 工程甜點預(yù)測

    研究區(qū)西部風(fēng)三段底部整體脆性指數(shù)較高,大于0.60(圖9a),且高脆性指數(shù)范圍與白云質(zhì)巖類發(fā)育范圍對應(yīng)性好,白云質(zhì)含量的增加有助于提高頁巖油儲集層的脆性。將白云質(zhì)粉砂巖發(fā)育范圍疊加在脆性指數(shù)預(yù)測平面分布圖上,區(qū)域①脆性最好,脆性指數(shù)為0.75~0.90;區(qū)域②脆性較好,脆性指數(shù)為0.70~0.80;區(qū)域③脆性相對略差,脆性指數(shù)為0.60~0.75。將A 井預(yù)測脆性指數(shù)與由巖石三軸抗壓測試獲得的脆性指數(shù)對比(圖9b),兩者一致性好,證實了脆性指數(shù)預(yù)測的準確性。

    圖9 研究區(qū)風(fēng)三段底部脆性指數(shù)預(yù)測結(jié)果Fig.9.Prediction of brittleness index at the bottom of Feng 3 member in the study area

    根據(jù)上述地應(yīng)力預(yù)測方法,首先利用密度積分得到上覆壓力,利用Eaton公式計算孔隙壓力,基于疊前同時反演間接得到靜態(tài)楊氏模量、靜態(tài)泊松比等彈性參數(shù)體,在此基礎(chǔ)上,由(4)式和(5)式計算得到最小水平主應(yīng)力、最大水平主應(yīng)力及水平主應(yīng)力差。最小水平主應(yīng)力和最大水平主應(yīng)力縱向上與深度有一定的相關(guān)性,深度越深,最小水平主應(yīng)力和最大水平主應(yīng)力越大(圖10),橫向上水平主應(yīng)力與地層構(gòu)造形態(tài)也存在一定的相關(guān)性。A 井處水平主應(yīng)力為測井計算并經(jīng)過Kaiser 實驗實測值校正后的結(jié)果,與最小水平主應(yīng)力和最大水平主應(yīng)力預(yù)測剖面吻合性好,證實了最小水平主應(yīng)力和最大水平主應(yīng)力預(yù)測的準確性。

    圖10 研究區(qū)風(fēng)城組地應(yīng)力預(yù)測剖面(剖面位置見圖8a)Fig.10.Profile of predicted in-situ stress for the Fengcheng formation in the study area(profile location is shown in Fig.8a)

    研究區(qū)風(fēng)三段底部最小水平主應(yīng)力分布呈臺階狀(圖11a),第一臺階位于構(gòu)造高部位,深度為4 100~4 900 m,最小水平主應(yīng)力相對較小,小于90.0 MPa;第二臺階深度為5 100~6 000 m,最小水平主應(yīng)力中等,為95.0~115.0 MPa;第三臺階為構(gòu)造低部位,深度為6 900~8 100 m,最小水平主應(yīng)力相對較大,大于120.0 MPa,最小水平主應(yīng)力的分布符合研究區(qū)風(fēng)三段底部的構(gòu)造趨勢,最大水平主應(yīng)力分布與最小水平主應(yīng)力相似。目前,研究區(qū)內(nèi)鉆遇風(fēng)城組的井大部分部署在第一臺階,深度偏小,最小水平主應(yīng)力相對較小。研究區(qū)風(fēng)三段底部整體水平主應(yīng)力差較小,為5.5~10.5 MPa,研究區(qū)西部水平主應(yīng)力差相對偏高,東部水平主應(yīng)力差相對偏低,水平主應(yīng)力差的分布與構(gòu)造形態(tài)關(guān)系不明顯,將白云質(zhì)粉砂巖發(fā)育范圍疊加在水平主應(yīng)力差預(yù)測平面圖上,區(qū)域①和區(qū)域②水平主應(yīng)力差相對偏大,為8.0~10.0 MPa;區(qū)域③水平主應(yīng)力差相對較小,為7.0~9.0 MPa(圖11b)。

    圖11 研究區(qū)風(fēng)三段底部地應(yīng)力預(yù)測立體圖Fig.11.3D view of predicted in-situ stress at the bottom of Feng 3 member in the study area

    在地質(zhì)工程甜點預(yù)測的基礎(chǔ)上,在A 井西北方向風(fēng)三段底部部署了水平井B 井,水平段長1 000 m,已完成壓裂,壓裂段共28 級,用液量為30 119 m3,加砂量為1 464 m3,共排液5 466 m3,返排率為18.1%,試油最高日產(chǎn)油108.00 m3。統(tǒng)計每級壓裂段的地面瞬時停泵壓力,與水平井軌跡處預(yù)測最小水平主應(yīng)力對比(圖12),兩者變化趨勢相似,相差約45.0 MPa。工程上最小水平主應(yīng)力上限等效于地面瞬時停泵壓力與液柱壓力之和,水平井風(fēng)三段底部平均深度為4 555 m,此時液柱壓力約為46.0 MPa,與兩者差值吻合,進一步驗證了地應(yīng)力預(yù)測的準確性。

    圖12 研究區(qū)水平井B井風(fēng)三段底部最小水平主應(yīng)力與瞬時停泵壓力對比Fig.12.Comparison of the minimum horizontal principal stress and the instantaneous pump-off pressure at the bottom of Feng 3 member in horizontal well B in the study area

    研究區(qū)C井在風(fēng)三段壓裂,用液強度為26.18 m3/m,加砂強度為2.14 m3/m,試油獲日產(chǎn)油2.03 m3;D 井在風(fēng)三段壓裂,用液強度為48.70 m3/m,加砂強度為2.95 m3/m,試油獲日產(chǎn)油16.93 m3。D 井位于白云質(zhì)粉砂巖發(fā)育區(qū),含油性好,C 井位于白云質(zhì)細粉砂巖發(fā)育區(qū),含油性略差;D 井脆性指數(shù)為0.80,高于C 井脆性指數(shù)0.65;D井水平主應(yīng)力差為9.5 MPa,大于C井水平主應(yīng)力差8.5 MPa。結(jié)合巖性、脆性指數(shù)及水平主應(yīng)力差預(yù)測結(jié)果,區(qū)域①白云質(zhì)粉砂巖發(fā)育,脆性最好,水平主應(yīng)力差相對偏大;區(qū)域②白云質(zhì)粉砂巖發(fā)育,脆性中等,水平主應(yīng)力差相對偏大;區(qū)域③白云質(zhì)粉砂巖發(fā)育,脆性略差,水平主應(yīng)力差相對偏小。

    5 結(jié)論

    (1)建立了一套以疊前同時反演為基礎(chǔ)的地質(zhì)工程甜點預(yù)測流程,并預(yù)測了瑪湖凹陷風(fēng)三段底部地質(zhì)工程甜點有利區(qū)。

    (2)巖相概率分析的巖相預(yù)測結(jié)果與地質(zhì)認識及測井巖相解釋結(jié)果吻合,實現(xiàn)了瑪湖凹陷風(fēng)城組巖相定性預(yù)測,刻畫了優(yōu)勢巖性白云質(zhì)粉砂巖的分布范圍,預(yù)測了地質(zhì)甜點有利區(qū),為利用疊前地震數(shù)據(jù)進行巖相預(yù)測提供了支撐。

    (3)Rickman 脆性指數(shù)法可基于疊前同時反演快速預(yù)測研究區(qū)脆性指數(shù),預(yù)測結(jié)果與巖石三軸抗壓測試結(jié)果吻合,從而刻畫了脆性發(fā)育區(qū),為預(yù)測工程甜點有利區(qū)提供支持。組合彈簧模型地應(yīng)力預(yù)測方法可實現(xiàn)地震頻率范圍內(nèi)的各向同性介質(zhì)地應(yīng)力預(yù)測,預(yù)測結(jié)果與Kaiser 實驗實測值吻合,其與脆性指數(shù)結(jié)合,在地質(zhì)甜點發(fā)育區(qū)優(yōu)選工程甜點分布范圍,優(yōu)選儲集層壓裂改造有利區(qū),為頁巖油勘探開發(fā)部署提供指導(dǎo)。

    (4)隨著鉆探井位的不斷增加,將微地震等技術(shù)應(yīng)用于監(jiān)測壓裂效果,下一步將細化研究脆性、地應(yīng)力等對儲集層壓裂改造的影響,深化地質(zhì)工程甜點的認識,更好地支撐頁巖油勘探開發(fā)。

    符號注釋

    BI——脆性指數(shù),無量綱;

    E——動態(tài)楊氏模量,GPa;

    Emax——最大動態(tài)楊氏模量,GPa;

    Emin——最小動態(tài)楊氏模量,GPa;

    Es——靜態(tài)楊氏模量,GPa;

    Ip——縱波阻抗,g/cm3·m/s;

    Is——橫波阻抗,g/cm3·m/s;

    pp——孔隙壓力,MPa;

    pv——上覆壓力,MPa;

    RppA(θ)——Aki-Richards 近似式推導(dǎo)的縱波反射系數(shù),無量綱;

    RppF(θ)——Fatti近似式推導(dǎo)的縱波反射系數(shù),無量綱;

    vp——縱波速度,m/s;

    vs——橫波速度,m/s;

    α——Biot系數(shù),無量綱,經(jīng)驗常數(shù)為0.8;

    ΔIp——介質(zhì)分界面兩側(cè)縱波阻抗差,g/cm3·m/s;

    ΔIs——介質(zhì)分界面兩側(cè)橫波阻抗差,g/cm3·m/s;

    Δvp——介質(zhì)分界面兩側(cè)縱波速度差,m/s;

    Δvs——介質(zhì)分界面兩側(cè)橫波速度差,m/s;

    Δρ——介質(zhì)分界面兩側(cè)密度差,g/cm3;

    θ——平均入射角,(°);

    ν——動態(tài)泊松比,無量綱;

    νmax——最大動態(tài)泊松比,無量綱;

    νmin——最小動態(tài)泊松比,無量綱;

    νs——靜態(tài)泊松比,無量綱;

    ξh——在最小水平主應(yīng)力方向上的應(yīng)變系數(shù);

    ξH——在最大水平主應(yīng)力方向上的應(yīng)變系數(shù);

    ρ——密度,g/cm3;

    σh——最小水平主應(yīng)力,MPa;

    σH——最大水平主應(yīng)力,MPa。

    猜你喜歡
    白云質(zhì)巖相儲集層
    準噶爾盆地瑪湖凹陷二疊系風(fēng)城組白云質(zhì)巖儲集層白云石成因
    渤中34-9油田古近系火山巖巖相特征與分布預(yù)測
    銀-額盆地哈日凹陷蘇紅圖組復(fù)雜巖類分析
    地下水(2018年6期)2018-12-14 01:07:10
    方解石顆粒尺寸對白云質(zhì)灰?guī)r中堿白云石反應(yīng)的影響
    硅酸鹽通報(2016年9期)2016-11-10 10:35:53
    川中震旦系燈影組儲集層形成及演化研究
    花崗巖儲集層隨鉆評價方法及應(yīng)用
    四川盆地普光氣田須家河組四段儲集層特征
    麥蓋提斜坡東部構(gòu)造帶奧陶系巖相、成巖作用帶與儲層成因
    多源多相復(fù)雜儲集層預(yù)測——以春風(fēng)油田新近系沙灣組一段1砂組為例
    斷塊油氣田(2014年5期)2014-03-11 15:33:44
    喀斯特區(qū)白云質(zhì)砂石山天然次生灌木林植物多樣性研究
    国产精品亚洲一级av第二区| 久久久色成人| 51午夜福利影视在线观看| 午夜两性在线视频| 成熟少妇高潮喷水视频| 亚洲性夜色夜夜综合| 大型黄色视频在线免费观看| 久久精品国产亚洲av香蕉五月| 麻豆一二三区av精品| 色av中文字幕| 国产亚洲av嫩草精品影院| 欧美日韩乱码在线| a级毛片a级免费在线| 长腿黑丝高跟| 国产一区二区激情短视频| 有码 亚洲区| 午夜福利免费观看在线| 九九热线精品视视频播放| 亚洲av电影在线进入| 99国产综合亚洲精品| 亚洲欧美日韩卡通动漫| 一夜夜www| 老汉色av国产亚洲站长工具| 亚洲专区中文字幕在线| 国产精品一区二区免费欧美| 91在线精品国自产拍蜜月 | 村上凉子中文字幕在线| 老熟妇乱子伦视频在线观看| 亚洲人成网站高清观看| 99久国产av精品| 精品一区二区三区视频在线 | 午夜视频国产福利| 夜夜躁狠狠躁天天躁| 人妻久久中文字幕网| 国产亚洲精品久久久com| 偷拍熟女少妇极品色| 亚洲av成人av| 欧美av亚洲av综合av国产av| 国产av在哪里看| 日本与韩国留学比较| 美女免费视频网站| 热99在线观看视频| 国产精品一区二区三区四区久久| 国产国拍精品亚洲av在线观看 | 亚洲精品影视一区二区三区av| 美女cb高潮喷水在线观看| 99国产精品一区二区蜜桃av| 国产亚洲欧美在线一区二区| 黄片小视频在线播放| 亚洲,欧美精品.| 长腿黑丝高跟| 制服丝袜大香蕉在线| 国产精品自产拍在线观看55亚洲| 国产黄a三级三级三级人| 女生性感内裤真人,穿戴方法视频| a级毛片a级免费在线| 女人十人毛片免费观看3o分钟| 91久久精品国产一区二区成人 | 亚洲第一电影网av| 亚洲国产精品久久男人天堂| 淫妇啪啪啪对白视频| 国产色爽女视频免费观看| 日日摸夜夜添夜夜添小说| 国产免费男女视频| 国产精品1区2区在线观看.| 国产极品精品免费视频能看的| 日韩免费av在线播放| 此物有八面人人有两片| 18禁黄网站禁片免费观看直播| 亚洲va日本ⅴa欧美va伊人久久| 亚洲无线观看免费| 老司机午夜十八禁免费视频| 亚洲av电影在线进入| 90打野战视频偷拍视频| 好男人电影高清在线观看| 波多野结衣巨乳人妻| 亚洲av成人精品一区久久| 欧美日韩亚洲国产一区二区在线观看| 在线免费观看不下载黄p国产 | 在线观看av片永久免费下载| 最好的美女福利视频网| 色噜噜av男人的天堂激情| 免费电影在线观看免费观看| 国产一区二区在线av高清观看| 久久精品综合一区二区三区| 乱人视频在线观看| 两个人视频免费观看高清| 夜夜爽天天搞| 欧美一区二区亚洲| 欧美午夜高清在线| 亚洲成人免费电影在线观看| 亚洲国产中文字幕在线视频| 欧美xxxx黑人xx丫x性爽| 长腿黑丝高跟| 国产蜜桃级精品一区二区三区| 日本五十路高清| 午夜福利高清视频| 午夜福利在线观看免费完整高清在 | 国产真实伦视频高清在线观看 | 色综合站精品国产| 色吧在线观看| 搡老岳熟女国产| 欧美bdsm另类| 床上黄色一级片| 女人十人毛片免费观看3o分钟| 国产精华一区二区三区| 国内精品久久久久久久电影| 日本在线视频免费播放| 免费一级毛片在线播放高清视频| 亚洲欧美日韩东京热| 动漫黄色视频在线观看| 全区人妻精品视频| h日本视频在线播放| 成人高潮视频无遮挡免费网站| 欧美乱色亚洲激情| 欧美最新免费一区二区三区 | 女人被狂操c到高潮| 在线看三级毛片| 成熟少妇高潮喷水视频| 免费人成在线观看视频色| 国产成+人综合+亚洲专区| 小说图片视频综合网站| 成人国产一区最新在线观看| 村上凉子中文字幕在线| 丰满的人妻完整版| 欧美性猛交╳xxx乱大交人| 岛国在线免费视频观看| 免费看十八禁软件| 亚洲人成伊人成综合网2020| 精品99又大又爽又粗少妇毛片 | 亚洲成人免费电影在线观看| 国产毛片a区久久久久| 一区二区三区国产精品乱码| 久久久国产成人精品二区| 国产单亲对白刺激| 欧美xxxx黑人xx丫x性爽| 90打野战视频偷拍视频| 岛国在线免费视频观看| 长腿黑丝高跟| 午夜福利免费观看在线| 欧美日韩乱码在线| 国产亚洲精品久久久com| 亚洲av不卡在线观看| 亚洲成av人片在线播放无| 国产亚洲精品久久久久久毛片| 精品人妻1区二区| 午夜免费观看网址| 国产精品美女特级片免费视频播放器| 国产欧美日韩一区二区精品| 久久久久久久精品吃奶| 最近视频中文字幕2019在线8| 美女高潮喷水抽搐中文字幕| 色吧在线观看| 国产欧美日韩精品亚洲av| 老汉色av国产亚洲站长工具| 国产一区二区亚洲精品在线观看| 国产精品日韩av在线免费观看| 丰满乱子伦码专区| 制服丝袜大香蕉在线| 淫秽高清视频在线观看| 欧美性猛交黑人性爽| 我要搜黄色片| 18禁美女被吸乳视频| 女人十人毛片免费观看3o分钟| 日本 欧美在线| 久9热在线精品视频| 亚洲精品久久国产高清桃花| 内地一区二区视频在线| 两性午夜刺激爽爽歪歪视频在线观看| h日本视频在线播放| 亚洲自拍偷在线| 90打野战视频偷拍视频| 成人av在线播放网站| 亚洲不卡免费看| 国产免费男女视频| 日本一本二区三区精品| 久久性视频一级片| 在线看三级毛片| 久久久精品大字幕| h日本视频在线播放| 色在线成人网| 高清日韩中文字幕在线| 男人和女人高潮做爰伦理| 色精品久久人妻99蜜桃| 亚洲人成电影免费在线| www国产在线视频色| 夜夜躁狠狠躁天天躁| 免费av不卡在线播放| 日本五十路高清| 麻豆国产97在线/欧美| 一区二区三区免费毛片| 两人在一起打扑克的视频| 免费在线观看日本一区| 国语自产精品视频在线第100页| 国产高清视频在线播放一区| 一个人看的www免费观看视频| www日本在线高清视频| 免费无遮挡裸体视频| 小说图片视频综合网站| 日本免费一区二区三区高清不卡| 国产高清视频在线播放一区| 男人舔奶头视频| 日韩成人在线观看一区二区三区| 久久精品国产亚洲av香蕉五月| www日本在线高清视频| 日韩人妻高清精品专区| 久久精品夜夜夜夜夜久久蜜豆| 国产伦精品一区二区三区视频9 | 国产av在哪里看| 97人妻精品一区二区三区麻豆| 男女床上黄色一级片免费看| 亚洲av成人av| 男女之事视频高清在线观看| 久久久久国内视频| 1024手机看黄色片| 又爽又黄无遮挡网站| 精品国内亚洲2022精品成人| 999久久久精品免费观看国产| 久久久久久久久大av| 丝袜美腿在线中文| 女人被狂操c到高潮| 9191精品国产免费久久| 老熟妇仑乱视频hdxx| 国产欧美日韩精品亚洲av| 男人的好看免费观看在线视频| 日本在线视频免费播放| 波野结衣二区三区在线 | 国产高潮美女av| 午夜精品一区二区三区免费看| 国产精品国产高清国产av| 日韩欧美国产一区二区入口| 亚洲熟妇中文字幕五十中出| 亚洲最大成人手机在线| 欧美成人一区二区免费高清观看| 日本黄色片子视频| 国产淫片久久久久久久久 | 亚洲精品在线观看二区| 久久久国产精品麻豆| 国产亚洲欧美在线一区二区| 真人一进一出gif抽搐免费| 高清毛片免费观看视频网站| 亚洲不卡免费看| 国产又黄又爽又无遮挡在线| 熟妇人妻久久中文字幕3abv| 国产高潮美女av| 美女黄网站色视频| 中亚洲国语对白在线视频| 老司机福利观看| 亚洲va日本ⅴa欧美va伊人久久| 人人妻人人澡欧美一区二区| 丰满乱子伦码专区| 51国产日韩欧美| 精品一区二区三区av网在线观看| 亚洲电影在线观看av| 黄片小视频在线播放| 亚洲最大成人手机在线| 看片在线看免费视频| 熟女电影av网| 日韩欧美三级三区| 亚洲国产精品成人综合色| 欧美日韩亚洲国产一区二区在线观看| 欧美不卡视频在线免费观看| 色吧在线观看| 欧美日韩乱码在线| 麻豆国产av国片精品| 在线播放国产精品三级| 夜夜躁狠狠躁天天躁| 国产精品99久久久久久久久| 国产午夜福利久久久久久| 免费看a级黄色片| 久久精品国产亚洲av香蕉五月| or卡值多少钱| 一本一本综合久久| 精品午夜福利视频在线观看一区| 看片在线看免费视频| 免费无遮挡裸体视频| 大型黄色视频在线免费观看| 日韩欧美一区二区三区在线观看| 黄色视频,在线免费观看| 亚洲七黄色美女视频| 久99久视频精品免费| 最近最新中文字幕大全电影3| 国产探花极品一区二区| 亚洲一区二区三区色噜噜| av欧美777| 久久久久国内视频| 真实男女啪啪啪动态图| 免费在线观看影片大全网站| 哪里可以看免费的av片| 在线观看免费午夜福利视频| 91九色精品人成在线观看| 在线视频色国产色| 欧美xxxx黑人xx丫x性爽| 久久久精品大字幕| 又粗又爽又猛毛片免费看| 欧美一区二区亚洲| 在线观看日韩欧美| 国产成+人综合+亚洲专区| 国产成人系列免费观看| 国产精品爽爽va在线观看网站| 日韩亚洲欧美综合| 观看免费一级毛片| 又紧又爽又黄一区二区| 国产私拍福利视频在线观看| 亚洲av免费高清在线观看| 啦啦啦韩国在线观看视频| 小说图片视频综合网站| 黄片小视频在线播放| 欧美成狂野欧美在线观看| xxx96com| 国产精品野战在线观看| 一区二区三区激情视频| 色吧在线观看| 成人性生交大片免费视频hd| 看黄色毛片网站| 91av网一区二区| 性色avwww在线观看| www.熟女人妻精品国产| 欧美三级亚洲精品| 免费电影在线观看免费观看| av女优亚洲男人天堂| 99久久99久久久精品蜜桃| 亚洲美女视频黄频| 欧美成狂野欧美在线观看| 人妻丰满熟妇av一区二区三区| av在线蜜桃| av在线天堂中文字幕| 校园春色视频在线观看| 亚洲精品一卡2卡三卡4卡5卡| 黄片小视频在线播放| 少妇人妻精品综合一区二区 | 欧美大码av| 国产亚洲精品一区二区www| 婷婷六月久久综合丁香| 国产精品日韩av在线免费观看| 97超级碰碰碰精品色视频在线观看| 国内久久婷婷六月综合欲色啪| 亚洲国产精品合色在线| 国产爱豆传媒在线观看| 九九热线精品视视频播放| 欧美成狂野欧美在线观看| 国产精品电影一区二区三区| 欧美又色又爽又黄视频| 欧美性猛交╳xxx乱大交人| 精品久久久久久成人av| 中文亚洲av片在线观看爽| 久久精品影院6| 51国产日韩欧美| 天堂影院成人在线观看| 九九在线视频观看精品| www.色视频.com| 性色av乱码一区二区三区2| 国产精品永久免费网站| 国产一级毛片七仙女欲春2| 久久久久久久午夜电影| 欧美一级毛片孕妇| 手机成人av网站| 国产精品久久久久久精品电影| tocl精华| 成人特级黄色片久久久久久久| 国产真人三级小视频在线观看| 老鸭窝网址在线观看| 国产精品嫩草影院av在线观看 | 亚洲精品乱码久久久v下载方式 | 19禁男女啪啪无遮挡网站| 午夜精品久久久久久毛片777| 免费在线观看日本一区| 人妻夜夜爽99麻豆av| 哪里可以看免费的av片| 国产美女午夜福利| 欧美最黄视频在线播放免费| 国产乱人视频| 国产精品自产拍在线观看55亚洲| 久久午夜亚洲精品久久| a级一级毛片免费在线观看| 两个人的视频大全免费| 午夜福利免费观看在线| 日本黄色片子视频| 日本熟妇午夜| 少妇丰满av| 久99久视频精品免费| 成人高潮视频无遮挡免费网站| 88av欧美| 性欧美人与动物交配| 国产亚洲精品久久久久久毛片| 亚洲av电影在线进入| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 欧美中文日本在线观看视频| 夜夜看夜夜爽夜夜摸| 亚洲av一区综合| 亚洲欧美日韩高清在线视频| 97人妻精品一区二区三区麻豆| 午夜福利视频1000在线观看| 久久人妻av系列| 精品一区二区三区视频在线观看免费| 人人妻,人人澡人人爽秒播| 免费大片18禁| 国产高清videossex| 人人妻,人人澡人人爽秒播| 欧美激情在线99| 怎么达到女性高潮| 日韩精品青青久久久久久| 露出奶头的视频| 欧美一级毛片孕妇| 波野结衣二区三区在线 | 日韩av在线大香蕉| 日韩大尺度精品在线看网址| 亚洲欧美日韩无卡精品| 亚洲最大成人手机在线| 特级一级黄色大片| 五月玫瑰六月丁香| 成熟少妇高潮喷水视频| 欧美一区二区精品小视频在线| 老鸭窝网址在线观看| 99精品在免费线老司机午夜| 午夜福利欧美成人| 欧美日韩福利视频一区二区| 99久久九九国产精品国产免费| 国产精品日韩av在线免费观看| 久久精品国产99精品国产亚洲性色| 小说图片视频综合网站| 精品免费久久久久久久清纯| 性欧美人与动物交配| 国产成+人综合+亚洲专区| 深夜精品福利| 国产精品嫩草影院av在线观看 | 久久久久久人人人人人| 有码 亚洲区| 99精品久久久久人妻精品| 丁香六月欧美| 国产乱人视频| 午夜福利视频1000在线观看| 欧美激情久久久久久爽电影| 黄色视频,在线免费观看| 久久久精品大字幕| 日韩免费av在线播放| 深夜精品福利| 女警被强在线播放| 亚洲一区二区三区不卡视频| 国内揄拍国产精品人妻在线| 在线观看日韩欧美| 日本成人三级电影网站| 舔av片在线| 一本精品99久久精品77| 女人高潮潮喷娇喘18禁视频| 日本一二三区视频观看| 欧美+日韩+精品| 久久精品影院6| 亚洲国产精品久久男人天堂| 老司机午夜十八禁免费视频| 中文字幕高清在线视频| 热99在线观看视频| 欧美色视频一区免费| 亚洲国产日韩欧美精品在线观看 | 一级作爱视频免费观看| 内射极品少妇av片p| 亚洲精品色激情综合| 欧美3d第一页| 国产极品精品免费视频能看的| 亚洲av一区综合| 99热这里只有是精品50| 无遮挡黄片免费观看| 成年人黄色毛片网站| 亚洲国产欧美网| 高清日韩中文字幕在线| 美女 人体艺术 gogo| 欧美不卡视频在线免费观看| 国产三级黄色录像| 午夜福利成人在线免费观看| 久久久久久九九精品二区国产| 免费人成视频x8x8入口观看| 国产精品 欧美亚洲| 一个人看视频在线观看www免费 | 成人18禁在线播放| av中文乱码字幕在线| 亚洲国产精品成人综合色| 久久伊人香网站| 色在线成人网| 99久久无色码亚洲精品果冻| 精品人妻一区二区三区麻豆 | 午夜福利成人在线免费观看| 人人妻人人澡欧美一区二区| 久久久国产成人免费| 国产精品1区2区在线观看.| 一二三四社区在线视频社区8| 久久精品国产99精品国产亚洲性色| 精品一区二区三区视频在线观看免费| 精品久久久久久久毛片微露脸| 日韩欧美国产在线观看| 国产一区在线观看成人免费| 国产成人av教育| 色哟哟哟哟哟哟| 久久性视频一级片| 激情在线观看视频在线高清| 国产一区二区激情短视频| 最新美女视频免费是黄的| 国产不卡一卡二| 成年免费大片在线观看| 中文字幕av成人在线电影| 十八禁人妻一区二区| 婷婷亚洲欧美| 日本免费a在线| 看黄色毛片网站| 十八禁网站免费在线| 日韩欧美国产在线观看| 搡老岳熟女国产| 国产97色在线日韩免费| 国内少妇人妻偷人精品xxx网站| 美女cb高潮喷水在线观看| 久久国产精品影院| 亚洲欧美日韩卡通动漫| 成人永久免费在线观看视频| 亚洲无线在线观看| 啦啦啦韩国在线观看视频| 18禁国产床啪视频网站| 国产精品 欧美亚洲| 男人舔女人下体高潮全视频| 色综合欧美亚洲国产小说| 欧美一区二区精品小视频在线| 狠狠狠狠99中文字幕| 欧美+日韩+精品| 中文字幕人妻丝袜一区二区| 亚洲一区高清亚洲精品| 日本五十路高清| 给我免费播放毛片高清在线观看| 一本综合久久免费| 免费看美女性在线毛片视频| 精品一区二区三区av网在线观看| 欧美av亚洲av综合av国产av| 桃色一区二区三区在线观看| 亚洲电影在线观看av| 男女之事视频高清在线观看| 在线观看舔阴道视频| 亚洲激情在线av| 精品日产1卡2卡| 精品久久久久久,| e午夜精品久久久久久久| 日韩免费av在线播放| 麻豆成人午夜福利视频| 欧美日韩黄片免| 成年人黄色毛片网站| 成年女人看的毛片在线观看| 一本久久中文字幕| 九九在线视频观看精品| 日韩高清综合在线| 欧美区成人在线视频| 久久久久九九精品影院| 在线观看免费视频日本深夜| 久久天躁狠狠躁夜夜2o2o| h日本视频在线播放| 亚洲国产精品成人综合色| 激情在线观看视频在线高清| 99国产精品一区二区蜜桃av| 欧美绝顶高潮抽搐喷水| 亚洲av一区综合| 9191精品国产免费久久| 国产精品爽爽va在线观看网站| 香蕉久久夜色| 在线观看美女被高潮喷水网站 | 美女cb高潮喷水在线观看| 一个人看的www免费观看视频| 最近最新中文字幕大全免费视频| 国产激情偷乱视频一区二区| 日本精品一区二区三区蜜桃| 精品人妻一区二区三区麻豆 | 9191精品国产免费久久| 欧美性猛交黑人性爽| 精品久久久久久久人妻蜜臀av| 国产av不卡久久| 午夜精品一区二区三区免费看| 深夜精品福利| 天堂网av新在线| 亚洲av免费在线观看| 国产在线精品亚洲第一网站| 免费av毛片视频| 国产在线精品亚洲第一网站| 国产精品一区二区三区四区久久| 脱女人内裤的视频| 国产三级黄色录像| 搡老熟女国产l中国老女人| 最近视频中文字幕2019在线8| 久久久久久久久中文| 久久精品综合一区二区三区| 欧美最黄视频在线播放免费| 天堂动漫精品| 啦啦啦观看免费观看视频高清| 欧美黄色片欧美黄色片| 男女那种视频在线观看| 宅男免费午夜| 久久伊人香网站| 久久这里只有精品中国| 人妻久久中文字幕网| 美女免费视频网站| 欧美日韩亚洲国产一区二区在线观看| 日本 av在线| 久久人妻av系列| 亚洲18禁久久av| 免费一级毛片在线播放高清视频| 国产精品爽爽va在线观看网站| 内射极品少妇av片p| 91久久精品国产一区二区成人 | 国产成人系列免费观看| 一卡2卡三卡四卡精品乱码亚洲| 丰满人妻一区二区三区视频av | 久久久久久人人人人人| 亚洲人成伊人成综合网2020| 在线视频色国产色| 国产精品永久免费网站| 国产国拍精品亚洲av在线观看 | 身体一侧抽搐| 亚洲成人久久爱视频| 久久精品人妻少妇| av片东京热男人的天堂| 最好的美女福利视频网| 人妻久久中文字幕网|