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

    天然氣水合物降壓開(kāi)采過(guò)程中儲(chǔ)層應(yīng)力規(guī)律分析

    2020-02-07 05:15:14張俊斌何玉發(fā)
    海洋科學(xué) 2020年11期
    關(guān)鍵詞:水合物壓差滲流

    金 顥, 張俊斌, 何玉發(fā)

    天然氣水合物降壓開(kāi)采過(guò)程中儲(chǔ)層應(yīng)力規(guī)律分析

    金 顥, 張俊斌, 何玉發(fā)

    (中海石油(中國(guó))有限公司深圳分公司深水工程技術(shù)中心, 廣東 深圳 518000)

    為了研究天然氣水合物降壓開(kāi)采過(guò)程的儲(chǔ)層應(yīng)力及其穩(wěn)定性, 運(yùn)用線性多孔彈性力學(xué)和巖石力學(xué)知識(shí), 考慮水合物儲(chǔ)層原始應(yīng)力、孔隙壓力、滲流附加應(yīng)力及降壓開(kāi)采水合物過(guò)程中水合物飽和度的變化, 建立了降壓開(kāi)采天然氣水合物儲(chǔ)層的力學(xué)模型, 結(jié)合墨西哥灣某處水合物藏的基本參數(shù), 對(duì)降壓開(kāi)采水合物儲(chǔ)層應(yīng)力變化和開(kāi)采過(guò)程的儲(chǔ)層穩(wěn)定性進(jìn)行研究。結(jié)果表明: 井底壓力是影響水合物儲(chǔ)層應(yīng)力變化的關(guān)鍵因素之一; 滲流附加應(yīng)力在一定程度上減小了儲(chǔ)層的應(yīng)力; 水合物分解儲(chǔ)層應(yīng)力發(fā)生變化, 儲(chǔ)層應(yīng)力在井壁處的波動(dòng)最大, 井壁處是整個(gè)儲(chǔ)層所受軸向偏應(yīng)力最大的位置, 因此井壁處是優(yōu)先發(fā)生剪切破壞的位置; 為了儲(chǔ)層的穩(wěn)定性, 降壓開(kāi)采水合物生產(chǎn)壓差應(yīng)小于2.19 MPa。

    降壓開(kāi)采; 天然氣水合物; 滲流附加應(yīng)力; 水合物飽和度; 臨界生產(chǎn)壓差; 儲(chǔ)層穩(wěn)定性

    目前, 天然氣被認(rèn)為是21世紀(jì)的一種未來(lái)優(yōu)質(zhì)、高效、潔凈的能源, 天然氣的蘊(yùn)藏量為現(xiàn)有地球化石燃料總碳量的2倍左右[1]。然而, 天然氣能否真正成為未來(lái)的可用能源, 關(guān)鍵在于能否有效地對(duì)天然氣水合物進(jìn)行合理的開(kāi)采[2]。

    天然氣水合物廣泛分布于海底沉積層和極地凍土帶。通常天然氣水合物以固態(tài)形式充填于多孔介質(zhì)儲(chǔ)層中, 天然氣水合物成藏相比于常規(guī)儲(chǔ)層條件要復(fù)雜, 并且不具備流動(dòng)性。因此, 根據(jù)天然氣水合物的存在條件和基本性質(zhì)決定了天然氣水合物的開(kāi)采方式[3]。天然氣水合物理論上的開(kāi)采思路是: 首先通過(guò)一定的方式打破天然氣水合物儲(chǔ)層的平衡性, 促使水合物在儲(chǔ)層中發(fā)生分解, 進(jìn)而開(kāi)采出分解的天然氣。目前, 天然氣水合物廣泛的開(kāi)采方法有: 注化學(xué)劑法、加熱法和降壓法等三種[4]。注化學(xué)劑法較加熱法作用緩慢, 且二者使用的費(fèi)用昂貴; 降壓法是通過(guò)降低儲(chǔ)層壓力使其低于水合物相平衡壓力引起儲(chǔ)層中水合物分解的方法, 降壓法開(kāi)采天然氣水合物由于其不需要昂貴的連續(xù)激發(fā), 所以被認(rèn)為是最具有商業(yè)前景的開(kāi)采方法[5]。

    天然氣水合物在降壓開(kāi)采過(guò)程中, 由于降低井底壓力, 打破了水合物儲(chǔ)層的穩(wěn)定性, 促使水合物發(fā)生分解。受水合物分解的影響, 分解區(qū)水合物儲(chǔ)層膠結(jié)性變差、力學(xué)強(qiáng)度大幅度降低, 進(jìn)而近井區(qū)儲(chǔ)層呈現(xiàn)出弱膠結(jié)、低強(qiáng)度、高孔滲的特點(diǎn), 類(lèi)似于常規(guī)疏松砂巖[6]。由于水合物降壓分解引起水合物儲(chǔ)層物性參數(shù)及巖石力學(xué)性質(zhì)等發(fā)生一系列變化, 如水合物飽和度減小、儲(chǔ)層內(nèi)聚力降低、儲(chǔ)層孔隙度增大、滲透率增大等, 導(dǎo)致水合物儲(chǔ)層應(yīng)力發(fā)生改變, 致使水合物儲(chǔ)層穩(wěn)定性降低, 達(dá)到或超過(guò)本身固有的剪切強(qiáng)度, 使近井區(qū)水合物儲(chǔ)層發(fā)生剪切破壞, 從而影響天然氣水合物安全有效的開(kāi)采。

    由于天然生成的天然氣水合物巖樣取芯困難, 目前, 國(guó)內(nèi)外學(xué)者對(duì)天然氣水合物的研究大多是在實(shí)驗(yàn)室內(nèi)完成的。因此, 基于前人的研究, 通過(guò)分析水合物儲(chǔ)層的穩(wěn)定性, 本文在采用降壓法開(kāi)采天然氣水合物過(guò)程中考慮儲(chǔ)層滲流附加應(yīng)力的影響, 同時(shí)利用疊加原理建立了天然氣水合物降壓開(kāi)采過(guò)程中儲(chǔ)層穩(wěn)定性的力學(xué)模型, 進(jìn)行了降壓開(kāi)采水合物過(guò)程中儲(chǔ)層應(yīng)力的模擬計(jì)算, 分析天然氣水合物在降壓開(kāi)采過(guò)程中儲(chǔ)層應(yīng)力變化特征及儲(chǔ)層的穩(wěn)定性, 為有效開(kāi)采天然氣水合物提供一定的理論基礎(chǔ)。

    1 水合物儲(chǔ)層模型建立

    通常天然氣水合物存在于疏松砂巖沉積層中, 埋藏深度淺。在降壓開(kāi)采天然氣水合物過(guò)程中, 可將整個(gè)水合物儲(chǔ)層劃分為兩個(gè)區(qū)域, 即水合物分解區(qū)、儲(chǔ)層穩(wěn)定區(qū), 如圖1所示[4, 7]。降壓開(kāi)采天然氣水合物時(shí), 當(dāng)水合物儲(chǔ)層內(nèi)壓力低于水合物相平衡時(shí), 打破了水合物儲(chǔ)層的平衡條件, 促使水合物發(fā)生分解。天然氣水合物的分解將會(huì)導(dǎo)致儲(chǔ)層的穩(wěn)定性發(fā)生改變, 而水合物儲(chǔ)層的應(yīng)力在一定程度上影響儲(chǔ)層的穩(wěn)定性。因此, 水合物儲(chǔ)層應(yīng)力的變化成為影響水合物儲(chǔ)層穩(wěn)定性的重要因素, 而井底壓力是分析水合物儲(chǔ)層應(yīng)力變化的關(guān)鍵因素。

    圖1 水合物儲(chǔ)層區(qū)域劃分示意圖

    在研究降壓開(kāi)采天然氣水合物儲(chǔ)層應(yīng)力變化時(shí), 結(jié)合國(guó)內(nèi)外學(xué)者建立的有關(guān)水合物飽和度、儲(chǔ)層孔隙度、內(nèi)聚力及儲(chǔ)層應(yīng)力等相關(guān)模型之間的關(guān)系, 構(gòu)建水合物儲(chǔ)層力學(xué)模型。由于降壓開(kāi)采天然氣水合物的過(guò)程是非等溫滲流過(guò)程, 因此該模型在研究過(guò)程中做了如下簡(jiǎn)化假設(shè): (1) 天然氣水合物儲(chǔ)層為均質(zhì)、各向同性的線彈性體。(2) 水合物儲(chǔ)層僅考慮三相三組分。三相: 氣、水、水合物; 三組分: 氣、水、水合物。其中, 水合物為固相, 氣體僅含甲烷, 且不重新生成新的水合物。(3) 降壓開(kāi)采水合物過(guò)程中, 水合物儲(chǔ)層的溫度變化忽略不計(jì), 本文注重考慮井底壓力的變化影響水合物儲(chǔ)層應(yīng)力改變的問(wèn)題。(4) 僅考慮氣、水兩相流, 且兩相滲流符合廣義Darcy定律。(5) 井為垂直井。

    1.1 水合物儲(chǔ)層穩(wěn)定性模型建立

    在降壓開(kāi)采水合物前, 水合物儲(chǔ)層受上覆巖層壓力、最大水平主應(yīng)力和最小水平主應(yīng)力, 水合物儲(chǔ)層處于應(yīng)力平衡狀態(tài)。當(dāng)水合物儲(chǔ)層被鉆開(kāi)時(shí), 近井區(qū)水合物儲(chǔ)層應(yīng)力平衡被打破, 儲(chǔ)層應(yīng)力發(fā)生變化, 引起水合物儲(chǔ)層應(yīng)力的重新分布, 水合物儲(chǔ)層的穩(wěn)定性發(fā)生改變。因此, 井周?chē)膽?yīng)力分布可表示為[8-9]:

    劉玉石等[10]通過(guò)運(yùn)用有限元方法分析了井周儲(chǔ)層穩(wěn)態(tài)達(dá)西滲流時(shí)儲(chǔ)層井壁圍巖應(yīng)力的變化, 從而得出井壁處的徑向應(yīng)力和周向應(yīng)力與線彈性理論計(jì)算的結(jié)果存在很大的差別。水合物儲(chǔ)層為多孔介質(zhì)儲(chǔ)層, 考慮到水合物儲(chǔ)層的滲透性, 同時(shí)又由于儲(chǔ)層流體在向孔隙中的徑向流動(dòng)過(guò)程中會(huì)對(duì)儲(chǔ)層周?chē)a(chǎn)生應(yīng)力, 因此滲流附加應(yīng)力會(huì)對(duì)儲(chǔ)層產(chǎn)生一定的影響, 其相關(guān)表達(dá)式為[11]:

    應(yīng)用線性疊加原理, 由公式(1)和(2)得到天然氣水合物儲(chǔ)層近井區(qū)任意一點(diǎn)(≥1)處的應(yīng)力分布式如下:

    (1) 初始條件

    天然氣水合物開(kāi)采前儲(chǔ)層壓力滿(mǎn)足:

    (2) 近井區(qū)水合物儲(chǔ)層初始分量

    1.2 水合物儲(chǔ)層物性參數(shù)

    (1) 水合物儲(chǔ)層孔隙壓力的分布

    開(kāi)采前, 整個(gè)水合物儲(chǔ)層處于穩(wěn)定的狀態(tài), 儲(chǔ)層中各點(diǎn)壓力均等于原始地層應(yīng)力。假設(shè)天然氣水合物儲(chǔ)層為一圓形區(qū)域, 在降壓開(kāi)采水合物過(guò)程中, 圓形封閉儲(chǔ)層中任意一點(diǎn)時(shí)刻的壓力滿(mǎn)足彈性不穩(wěn)定滲流邊界條件, 其控制方程可表示為[12-13]:

    整理方程(8), 得到降壓過(guò)程中水合物儲(chǔ)層孔隙壓力的分布函數(shù):

    當(dāng)=1時(shí), 得到井底壓力w隨時(shí)間變化的關(guān)系式為:

    式中:為水合物儲(chǔ)層厚度, m;為通過(guò)井孔的砂巖流量, cm3/s。

    通過(guò)式(9)及式(10)可得:

    式(11)可化為:

    (2) 水合物飽和度隨時(shí)間變化

    降壓開(kāi)采水合物過(guò)程中水合物分解動(dòng)力學(xué)方程可表示為:

    根據(jù)降壓法開(kāi)采中Kim-Bishnoi模型, 天然氣水合物分解速率可表示為:

    水興百業(yè)旺。結(jié)合飲水安全工程建設(shè),各地的鄉(xiāng)村振興建設(shè)也搞得有聲有色。有了水,農(nóng)戶(hù)就把改廁、改灶、改路結(jié)合起來(lái),很多農(nóng)村的面貌煥然一新。

    分解區(qū)水合物儲(chǔ)層內(nèi)氣、水、水合物之間的關(guān)系式:

    H+W+g= 1, (16)

    其中,H=HH/H, 式中:H為分解水合物的體積;H為水合物物質(zhì)的量, mol;H為水合物密度, 0.918 g/cm3。

    多孔介質(zhì)孔隙體積為:

    其中,s= 2πd。

    (3) 水合物儲(chǔ)層孔隙度

    隨著水合物儲(chǔ)層埋深、狀態(tài)分布等特點(diǎn), 水合物儲(chǔ)層物性參數(shù)發(fā)生變化。隨著埋深的增加, 地層沉積壓實(shí)作用增大, 水合物儲(chǔ)層孔隙度減小。多孔介質(zhì)儲(chǔ)層孔隙度與有效應(yīng)力關(guān)系式為[15-16]:

    Nazridoust等[17]假設(shè)天然氣水合物儲(chǔ)層孔隙度均勻分布且忽略了溫度的影響, 得出儲(chǔ)層孔隙度與水合物飽和度之間的關(guān)系式為:

    由式(20)和式(21), 可得到孔隙度變化動(dòng)態(tài)模型:

    式中:為儲(chǔ)層有效應(yīng)力, MPa。

    (4) 水合物儲(chǔ)層內(nèi)聚力及有效應(yīng)力

    在降壓開(kāi)采水合物過(guò)程中, 由于水合物的分解導(dǎo)致儲(chǔ)層膠結(jié)性變差, 儲(chǔ)層內(nèi)聚力減小。隨著水合物的分解, 水合物儲(chǔ)層的內(nèi)聚力降低, 儲(chǔ)層內(nèi)聚力與孔隙度之間的數(shù)學(xué)表達(dá)式如下[18]:

    式中:為水合物分解后儲(chǔ)層內(nèi)聚力, MPa;0為初始水合物儲(chǔ)層內(nèi)聚力, MPa;為假設(shè)為1.2;為水合物儲(chǔ)層孔隙度增量。

    應(yīng)用Terzaghi有效應(yīng)力原理, 設(shè)儲(chǔ)層巖石應(yīng)力受壓為負(fù), 受拉為正[19], 其表達(dá)式為:

    式中:為有效應(yīng)力, MPa;v為儲(chǔ)層原始垂向應(yīng)力, MPa;為Biot系數(shù);為Biot有效應(yīng)力參數(shù)。

    (5) 天然氣水合物相平衡

    天然氣水合物廣泛分布于極地凍土帶和海底沉積層中, 水合物在儲(chǔ)層中存在穩(wěn)定帶即水合物在儲(chǔ)層中特定溫度和壓力條件下的一個(gè)范圍, 在這個(gè)范圍內(nèi)水合物達(dá)到相平衡, 甲烷水合物相溫度和壓力平衡采用Makogon[20]模型計(jì)算, 得出表達(dá)式為:

    lge=(–0) +(–0)2+, (25)

    式中:e、0為分別是水合物分解的平衡壓力(MPa)和平衡溫度(K);= 0.034 2K–1,= 0.000 5K–2,= 6.480 4。

    由于降壓開(kāi)采水合物導(dǎo)致的變形幾乎是注熱的10倍[21], 因此本文在研究降壓開(kāi)采水合物過(guò)程中, 假設(shè)水合物儲(chǔ)層內(nèi)的溫度變化忽略不計(jì)。

    (6) 儲(chǔ)層穩(wěn)定性判據(jù)準(zhǔn)則

    在研究水合物儲(chǔ)層的穩(wěn)定性時(shí), 運(yùn)用Drucker- Prager準(zhǔn)則進(jìn)行判斷:

    即,

    降壓開(kāi)采水合物過(guò)程中生產(chǎn)壓差的計(jì)算公式表示為:

    Δ=0–w, (28)

    式中: Δ為水合物儲(chǔ)層臨界生產(chǎn)壓差, MPa;0為原始水合物儲(chǔ)層的孔隙壓力, MPa;w為井底壓力, MPa。

    將式(3)、(29)代入(28)化簡(jiǎn), 最終得到降壓開(kāi)采水合物過(guò)程中儲(chǔ)層穩(wěn)定性指數(shù)與生產(chǎn)壓差之間的關(guān)系式為:

    2 降壓開(kāi)采水合物儲(chǔ)層應(yīng)力變化分析

    2.1 模擬采用的基本參數(shù)

    模型計(jì)算所選用的參數(shù)采用墨西哥灣某處甲烷水合物藏的部分基本參數(shù)[18]來(lái)分析降壓開(kāi)采水合物過(guò)程中近井區(qū)儲(chǔ)層應(yīng)力的變化特征。

    墨西哥灣某處水合物藏的部分基本參數(shù): 水合物儲(chǔ)層溫度288 K; 儲(chǔ)層初始?jí)毫?6.9 MPa; 井孔半徑0.15 m; 儲(chǔ)層垂向原始主應(yīng)力21.8 MPa; 儲(chǔ)層最大水平原始主應(yīng)力20.45 MPa; 儲(chǔ)層最小水平原始主應(yīng)力19.70 MPa; 儲(chǔ)層厚度30 m; Biot有效應(yīng)力參數(shù)0.6; Biot系數(shù)1.0; 水合物飽和度0.5; 儲(chǔ)層孔隙度0.4; 泊松比0.4; 儲(chǔ)層初始滲透率0.2 μm2; 多孔介質(zhì)壓縮系數(shù)10–9Pa–1; 水的壓縮系數(shù)4.6×10–10Pa–2; 儲(chǔ)層內(nèi)聚力1.8 MPa; 內(nèi)摩擦角30°。其中, 泄油半徑100 m。

    2.2 水合物儲(chǔ)層物性參數(shù)及力學(xué)性質(zhì)變化分析

    圖2、圖3分別給出了生產(chǎn)壓差3 MPa, 不同生產(chǎn)時(shí)間下水合物飽和度分布、儲(chǔ)層滲透率分布曲線圖。降壓開(kāi)采水合物過(guò)程中, 初始水合物儲(chǔ)層飽和度受已存自由氣體的影響。

    圖2 水合物飽和度分布曲線

    從圖2中可以看出, 由降壓開(kāi)采引起的水合物飽和度變化主要發(fā)生在近井區(qū), 在降壓開(kāi)采水合物初期(20 d), 靠近井壁處的水合物飽和度減小的相對(duì)緩慢, 這主要是因?yàn)樵衔飪?chǔ)層內(nèi)存在的自由氣體使得儲(chǔ)層內(nèi)的壓力升高, 在一定程度上阻礙了水合物的分解速率, 使得水合物的分解速率減慢; 在較遠(yuǎn)處, 水合物飽和度趨于穩(wěn)定。隨著時(shí)間不斷增加(60 d), 儲(chǔ)層內(nèi)的自由氣體逐漸減少, 水合物飽和度減小的速度逐漸加快, 這與唐良廣等[22]研究的水合物飽和度變化規(guī)律相符, 進(jìn)一步證明了該水合物飽和度變化模型的準(zhǔn)確性。

    圖3 水合物儲(chǔ)層滲透率分布曲線

    分析圖3可以看出, 井壁處的儲(chǔ)層滲透率增大的最快, 遠(yuǎn)離井壁儲(chǔ)層滲透率逐漸恢復(fù)到原始狀態(tài)。在水合物分解區(qū), 特別是井壁處水合物分解速率最快, 水合物飽和度減小的最快, 根據(jù)式(19)可知, 水合物飽和度和儲(chǔ)層滲透率呈反比。因此, 生產(chǎn)時(shí)間越長(zhǎng), 水合物飽和度變得越小, 儲(chǔ)層滲透率就越大。

    2.3 儲(chǔ)層滲流附加應(yīng)力的影響

    水合物儲(chǔ)層為多孔介質(zhì)儲(chǔ)層, 通過(guò)模擬計(jì)算得到在生產(chǎn)壓差3 MPa, 開(kāi)井生產(chǎn)100 d后, 研究?jī)?chǔ)層應(yīng)力過(guò)程中考慮滲流附加應(yīng)力與不考慮滲流附加應(yīng)力兩種情況下對(duì)水合物儲(chǔ)層徑向應(yīng)力、周向應(yīng)力進(jìn)行了對(duì)比, 如圖4所示。

    圖4 存在滲流附加應(yīng)力與無(wú)滲流附加應(yīng)力兩種情況下儲(chǔ)層應(yīng)力的對(duì)比

    從圖4中可以得出, 考慮滲流附加應(yīng)力時(shí)的儲(chǔ)層徑向、周向應(yīng)力小于不考慮滲流附加應(yīng)力時(shí)的儲(chǔ)層徑向、周向應(yīng)力。通過(guò)對(duì)式(2)分析可知, 水合物降壓分解過(guò)程中, 飽和度逐漸減小, 儲(chǔ)層有效孔隙度逐漸增大, 這說(shuō)明滲流附加應(yīng)力降低了儲(chǔ)層的穩(wěn)定性。

    在井底壓力13.9 MPa, 井眼半徑0.15 m的條件下, 分析考慮近井區(qū)水合物儲(chǔ)層滲流附加應(yīng)力和不考慮滲流附加應(yīng)力時(shí), 儲(chǔ)層所受最大軸向偏應(yīng)力分布曲線如圖5所示。結(jié)合圖4, 由圖5不難得出, 考慮滲流附加應(yīng)力時(shí)儲(chǔ)層所受最大軸向偏應(yīng)力小于不考慮滲流附加應(yīng)力時(shí)儲(chǔ)層所受最大軸向偏應(yīng)力。因此, 在建立水合物儲(chǔ)層的力學(xué)模型時(shí), 必須考慮滲流附加應(yīng)力對(duì)水合物儲(chǔ)層的影響, 以減小預(yù)測(cè)結(jié)果的誤差。

    圖5 儲(chǔ)層最大軸向偏應(yīng)力分布

    2.4 水合物儲(chǔ)層應(yīng)力分析

    通過(guò)降低井底壓力, 使水合物儲(chǔ)層壓力低于其平衡壓力, 促使水合物發(fā)生分解生成天然氣和水。井底壓力變化會(huì)改變一定區(qū)域的原始儲(chǔ)層應(yīng)力, 造成分解區(qū)水合物儲(chǔ)層發(fā)生應(yīng)力集中, 從而降低水合物儲(chǔ)層的穩(wěn)定性。在降壓開(kāi)采水合物過(guò)程中水合物儲(chǔ)層徑向、周向、垂向應(yīng)力表現(xiàn)出不同的變化。圖6顯示了在生產(chǎn)壓差3 MPa, 開(kāi)采時(shí)間100 d后, 塑性區(qū)水合物儲(chǔ)層應(yīng)力的分布曲線。

    通過(guò)式(3)不難看出, 水合物儲(chǔ)層應(yīng)力是角度的函數(shù), 儲(chǔ)層周向、垂向應(yīng)力隨角度呈周期性變化, 儲(chǔ)層徑向應(yīng)力不發(fā)生變化。近井區(qū)水合物儲(chǔ)層中最大、最小應(yīng)力之間的差值隨著距井孔距離的增大逐漸減小, 在井周上某點(diǎn)與水平最大主應(yīng)力方向的夾角90°或270°方位處, 水合物儲(chǔ)層的最大、最小應(yīng)力之間的差值最大, 即近井區(qū)水合物儲(chǔ)層所受到的擠壓應(yīng)力最大。根據(jù)Drucker-Prager巖石破壞準(zhǔn)則, 儲(chǔ)層周向、徑向應(yīng)力之間的差值越大, 儲(chǔ)層越容易發(fā)生剪切破壞。因此, 這兩個(gè)方位處也是儲(chǔ)層最容易發(fā)生剪切破壞的位置, 儲(chǔ)層發(fā)生剪切破壞一般是從井壁處開(kāi)始的。

    圖6 儲(chǔ)層應(yīng)力隨角度的變化曲線

    Drucker-Prager準(zhǔn)則認(rèn)為儲(chǔ)層發(fā)生的破壞主要是剪切破壞, 根據(jù)式(29)進(jìn)一步說(shuō)明儲(chǔ)層破壞的臨界條件。保持其他參數(shù)不變, 水合物儲(chǔ)層穩(wěn)定性隨生產(chǎn)壓差變化的關(guān)系圖如圖7所示。由圖7可知, 水合物儲(chǔ)層穩(wěn)定性指數(shù)隨生產(chǎn)壓差的增大而減小。這是因?yàn)榫讐毫υ叫? 即生產(chǎn)壓差越大, 儲(chǔ)層有效應(yīng)力增大越快, 造成儲(chǔ)層失穩(wěn)的可能性就越大; 隨著生產(chǎn)壓差的增大, 水合物分解速率加快, 儲(chǔ)層物性變化越顯著。

    圖7 儲(chǔ)層穩(wěn)定性隨生產(chǎn)壓差變化關(guān)系圖

    通過(guò)對(duì)圖7曲線進(jìn)行線性回歸得到儲(chǔ)層穩(wěn)定性指數(shù)與生產(chǎn)壓差之間的關(guān)系式(30)所示:

    = 0.247 3Δ2– 2.803 7Δ+ 5.092 9, (30)

    式中:為儲(chǔ)層穩(wěn)定性指數(shù), MPa; Δ為生產(chǎn)壓差, MPa。

    通過(guò)解式(30), 得到水合物儲(chǔ)層臨界生產(chǎn)壓差為2.19 MPa。當(dāng)生產(chǎn)壓差小于2.19 MPa, 水合物儲(chǔ)層保持穩(wěn)定; 相反, 當(dāng)生產(chǎn)壓差大于2.19 MPa, 水合物儲(chǔ)層發(fā)生剪切強(qiáng)度破壞, 有可能會(huì)引起儲(chǔ)層出砂。因此, 降壓開(kāi)采水合物過(guò)程中, 應(yīng)保持井底壓力大于14.71 MPa。

    綜合以上分析可知, 在降壓開(kāi)采水合物過(guò)程中, 水合物分解區(qū)是整個(gè)儲(chǔ)層穩(wěn)定性較差的區(qū)域。降低井底壓力: 一方面, 促進(jìn)了水合物的分解, 儲(chǔ)層內(nèi)聚力減小, 儲(chǔ)層各應(yīng)力發(fā)生改變; 另一方面, 不合理的降壓會(huì)造成儲(chǔ)層發(fā)生剪切破壞, 進(jìn)而影響水合物的有效開(kāi)采。因此, 在降壓開(kāi)采水合物過(guò)程中應(yīng)保證合理的井底壓力。

    3 結(jié)論

    1) 降壓開(kāi)采水合物過(guò)程中, 由于水合物發(fā)生分解, 分解區(qū)儲(chǔ)層呈現(xiàn)出弱膠結(jié)、低強(qiáng)度、高孔滲的特點(diǎn), 在不同生產(chǎn)壓差下, 近井區(qū)儲(chǔ)層應(yīng)力發(fā)生了明顯的改變, 距井孔較遠(yuǎn)處儲(chǔ)層應(yīng)力不再發(fā)生變化而趨于原始儲(chǔ)層應(yīng)力。說(shuō)明近井區(qū)儲(chǔ)層出現(xiàn)應(yīng)力集中的現(xiàn)象, 井底壓力是影響儲(chǔ)層穩(wěn)定性的重要因素。

    2) 生產(chǎn)壓差越大, 塑性區(qū)水合物儲(chǔ)層所受軸向偏應(yīng)力就越大, 井壁處所受軸向偏應(yīng)力最大, 是整個(gè)水合物儲(chǔ)層優(yōu)先破壞的位置。其中, 生產(chǎn)壓差2.19 MPa作為降壓開(kāi)采水合物過(guò)程中儲(chǔ)層發(fā)生剪切破壞的臨界壓差。

    3) 有待進(jìn)一步開(kāi)展對(duì)天然氣水合物降壓開(kāi)采過(guò)程中儲(chǔ)層應(yīng)力變化分析, 確定水合物儲(chǔ)層發(fā)生剪切破壞的生產(chǎn)壓差的合理臨界點(diǎn), 為今后有效開(kāi)采天然氣水合物奠定一定的基礎(chǔ)。

    [1] Makogon Y F, Holditch S A, Makogon T Y. Natural gas- hydrates: A potential energy source for the 21st Century[J]. Journal of Petroleum Science & Engineering, 2007, 56(1-3): 14-31.

    [2] Moridis G J, Collett T S, Dallimore S R, et al. Numerical studies of gas production from several CH4, hydrate zones at the Mallik site, Mackenzie Delta, Canada[J]. Journal of Petroleum Science & Engineering, 2004, 43(3-4): 219-238.

    [3] Sloan E D. Clathrate hydrates of natural gases[J]. Journal of the American Chemical Society, 1998, 120(43): 11212.

    [4] Ji C, Ahmadi G, Smith D H. Natural gas production from hydrate decomposition by depressurization[J]. Chemical Engineering Science, 2001, 56(20): 5801-5814.

    [5] Wang Yi, Li Xiaosen, Xu Wenyue, et al. Experimental investigation into factors influencing methane hydrate formation and a novel method for hydrate formation in porous media[J]. Energy and Fuel, 2013, 27(7): 3751- 3757.

    [6] 喻西崇, 吳應(yīng)湘, 安維杰, 等. 開(kāi)采地層中的天然氣水合物的數(shù)學(xué)模型[J]. 天然氣工業(yè), 2004, 24(1): 63-67. Yu Xichong, Wu Yingxiang, An Weijie, et al. Mathema-tical model to recover gas hydrate from formations[J]. Natural Gas Industry, 2004, 24(1): 63-67.

    [7] Ahmadi G, Ji C, Smith D H. Numerical solution for natural gas production from methane hydrate dissociation[J]. Journal of Petroleum Science & Engineering, 2004, 41(4): 269-285.

    [8] Aadnoy B S, Chenevert M E. Stability of highly inclined boreholes[J]. SPE Drilling Engineering, 1987, 2: 364-374.

    [9] Fjar E. Petroleum Related Rock Mechanics[M]. Amsterdam: Elsevier, 2008.

    [10] 劉玉石, 黃克累. 孔隙流體對(duì)井眼穩(wěn)定的影響[J]. 石油鉆探技術(shù), 1995, 23(3): 4-6. Liu Yushi, Huang Kelei. Effects of pore fluids on borehole stability[J]. Petroleum Drilling Techniques, 1995, 23(3): 4-6.

    [11] 閻鐵, 李士斌. 深部井眼巖石力學(xué)理論與應(yīng)用[M]. 北京: 石油工業(yè)出版社, 2002. Yan Tie, Li Shibin. Deep Borehole Rock Mechanics Theory and Applications[M]. Beijing: Petroleum Industry Press, 2002.

    [12] 張建國(guó), 杜殿發(fā), 侯建, 等. 油氣層滲流力學(xué)[M]. 東營(yíng): 中國(guó)石油大學(xué)出版社, 2010. Zhang Jianguo, Du Dianfa, Hou Jian, et al. The Oil and Gas Flow through Porous Media[M]. Dongying: China University of Petroleum Press, 2010.

    [13] Jacob B. Dynamics of Fluids in Porous Media[M]. New York: American Elsevier Pub. Co., 1972.

    [14] Masuda Y, Kurihara M. A filed scale simulation study on gas productivity of formation containing gas hydrates[J]. Proceedings of the Fourth International Conference on Gas Hydrates, 2002, 15(6): 622-627.

    [15] 張新紅, 秦積舜. 低滲巖心物性參數(shù)與應(yīng)力關(guān)系的試驗(yàn)研究[J]. 中國(guó)石油大學(xué)學(xué)報(bào): 自然科學(xué)版, 2001, 25(4): 56-57. Zhang Xinhong, Qin Jishun. Relation of Ground Stress with the Properties of Low Permeable Core[J]. Journal of China University of Petroleum (Edition of Natural Science), 2001, 25(4): 56-57.

    [16] 沈海超. 天然氣水合物藏降壓開(kāi)采流固耦合數(shù)值模擬研究[D]. 青島: 中國(guó)石油大學(xué), 2009. Shen Haichao. Fluid-solid coupling numerical simulation on natural gas production from hydrate reservoirs by depressurization[D]. Qingdao: China University of Petroleum, 2009.

    [17] Nazridoust K, Ahmadi G. Computational modeling of methane hydrate dissociation in a sandstone core[J]. Che-mical Engineering Science, 2007, 62(22): 6155-6177.

    [18] Freij-Ayoub R, Tan C, Clennell B, et al. A wellbore stability model for hydrate bearing sediments[J]. Journal of Petroleum Science & Engineering, 2007, 57(1): 209-220.

    [19] 冉啟全, 顧小蕓. 油藏滲流與應(yīng)力耦合分析[J]. 巖土工程學(xué)報(bào), 1998, 20(2): 69-73. Ran Qiquan, Gu Xiaoyun. Coupling analysis of multiphase flow and stress for oil reservoir[J]. Chinese Journal of Geotechnical Engineering, 1998, 20(2): 69-73.

    [20] Makogon Y F. Hydrates of hydrocarbons[J]. Penn Well Books, 1997, 610(52): 50-62.

    [21] Kimoto S, Oka F, Fushita T, et al. A chemo-thermo- mechanically coupled numerical simulation of the subsurface ground deformations due to methane hydrate dissociation[J]. Computers & Geotechnics, 2007, 34(4): 216-228.

    [22] 唐良廣, 馮自平, 李小森, 等. 天然氣水合物降壓開(kāi)采的數(shù)值模擬研究[C]//中國(guó)科學(xué)技術(shù)協(xié)會(huì).中國(guó)工程熱物理學(xué)會(huì)多相流學(xué)術(shù)會(huì)議. 北京: 中國(guó)科學(xué)技術(shù)協(xié)會(huì), 2006: 691-698. Tang Lianguang, Feng Ziping, Li Xiaoshen, et al. Numerical Simulation on Natural Gas Production from Hydrate Reservoirs by Depressurization[C]// Science and Technology Conference of Chinese. Multiphase Flow Academic Conference of Chinese Society of Engineering Thermophysics. Beijing: Science and Technology Conference of Chinese, 2006: 691-698.

    Geomechanical stress law analysis of hydrate reservoir for gas hydrate production by depressurization

    JIN Hao, ZHANG Jun-bin, HE Yu-fa

    (Deepwater Engineering Technology Center of CNOOC (China) Co., Ltd., Shenzhen 518000, China)

    A model combining linear porous elasticity and rock mechanics knowledge is developed to analyze the geomechanical stress of the hydrate reservoir and assess the wellbore stability during depressurization. The proposed model considers the original reservoir stress, pore pressure, seepage additional stress, and depressurization of hydrate mining during the process of hydrate saturation change. This model is applied to the case of the Gulf of Mexico, where the hydrate reservoir basic parameters have been publicly published and depressurization of the hydrate mining reservoir has been simulated and analyzed. Results show that the effect of bottom hole pressure is the dominant factor affecting the formation of reservoir stress. Moreover, reservoir stress significantly affects the wall. Therefore, the wall is the location that is most prone to shear failure. To assess the wellbore stability during the process of depressurization, the production pressure difference of the hydrate reservoir must be <2.19 MPa.

    depressurization; gas hydrate reservoir; seepage additional stress; hydrate saturation; critical production pressure differential; wellbore stability

    Oct. 31, 2019

    TE319

    A

    1000-3096(2020)11-0001-09

    10.11759/hykx20191031002

    2019-10-31;

    2019-12-10

    國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目 (2016YFC0304000);國(guó)家重點(diǎn)研發(fā)計(jì)劃課題 (2016YFC03044007);國(guó)家重點(diǎn)研發(fā)計(jì)劃課題 (2018YFC0310205)

    [The National Key R&D Program of China, No. 2016YFC0304000; The National Key R&D Program of China, No. 2016YFC03044007; The National Key R&D Program of China, No. 2018YFC0310205]

    金顥(1984-),男,廣東深圳人,高級(jí)工程師,主要從事鉆井、完井、測(cè)試、水合物等領(lǐng)域相關(guān)內(nèi)容研究,電話:15899755760,E-mail:jinhao@cnooc.com.cn

    (本文編輯: 劉珊珊)

    猜你喜歡
    水合物壓差滲流
    氣井用水合物自生熱解堵劑解堵效果數(shù)值模擬
    在線處理脫丁烷塔壓差高
    熱水吞吐開(kāi)采水合物藏?cái)?shù)值模擬研究
    天然氣水合物保壓轉(zhuǎn)移的壓力特性
    我國(guó)海域天然氣水合物試采成功
    二維氣體動(dòng)力學(xué)中壓差方程的特征分解和簡(jiǎn)單波
    基于改進(jìn)經(jīng)濟(jì)壓差算法的主站AVC的研究
    淺析濾紙透氣度與初始?jí)翰畹年P(guān)系
    簡(jiǎn)述滲流作用引起的土體破壞及防治措施
    河南科技(2014年12期)2014-02-27 14:10:26
    關(guān)于渠道滲流計(jì)算方法的選用
    河南科技(2014年11期)2014-02-27 14:09:48
    成年女人毛片免费观看观看9| 亚洲欧美清纯卡通| 久久中文看片网| 男女做爰动态图高潮gif福利片| 亚洲美女黄片视频| 男人舔女人下体高潮全视频| 免费电影在线观看免费观看| 久久精品国产亚洲av香蕉五月| 热99在线观看视频| avwww免费| 国产不卡一卡二| 国产 一区精品| 久久精品国产自在天天线| 热99在线观看视频| 99热网站在线观看| 人人妻人人澡欧美一区二区| 久久精品人妻少妇| 999久久久精品免费观看国产| 女同久久另类99精品国产91| 亚洲经典国产精华液单| 亚洲色图av天堂| 伦理电影大哥的女人| 国产毛片a区久久久久| 亚洲色图av天堂| 999久久久精品免费观看国产| 亚洲电影在线观看av| 久久精品91蜜桃| 亚洲一区二区三区色噜噜| 99久久久亚洲精品蜜臀av| 大型黄色视频在线免费观看| 日日啪夜夜撸| 久久人人精品亚洲av| 精品久久久久久久久亚洲 | 日韩欧美精品v在线| 22中文网久久字幕| 男女啪啪激烈高潮av片| av中文乱码字幕在线| 亚洲精品一区av在线观看| av国产免费在线观看| 久久精品国产自在天天线| 久久婷婷人人爽人人干人人爱| 亚洲av.av天堂| 亚洲18禁久久av| 亚洲av免费在线观看| 国产亚洲91精品色在线| 久久热精品热| 日韩中文字幕欧美一区二区| 啦啦啦韩国在线观看视频| 日韩一区二区视频免费看| 国产在线男女| 欧美激情国产日韩精品一区| 日韩一区二区视频免费看| 小蜜桃在线观看免费完整版高清| 韩国av一区二区三区四区| 免费av观看视频| 中文资源天堂在线| 亚洲欧美日韩东京热| .国产精品久久| 成年女人看的毛片在线观看| 国产乱人视频| 亚洲精品日韩av片在线观看| 免费av毛片视频| 91麻豆精品激情在线观看国产| 十八禁网站免费在线| 中文字幕熟女人妻在线| 偷拍熟女少妇极品色| 日日干狠狠操夜夜爽| 成人无遮挡网站| 一进一出抽搐动态| 成年女人毛片免费观看观看9| 国产精品永久免费网站| 午夜精品一区二区三区免费看| 日韩一区二区视频免费看| 精品福利观看| 国产免费男女视频| 少妇高潮的动态图| 成人特级黄色片久久久久久久| 国产乱人视频| 色在线成人网| 国产主播在线观看一区二区| 丰满的人妻完整版| 精品久久久久久久久av| 最近视频中文字幕2019在线8| 女人十人毛片免费观看3o分钟| 亚洲最大成人手机在线| 啦啦啦啦在线视频资源| 国产又黄又爽又无遮挡在线| 亚洲午夜理论影院| 欧美绝顶高潮抽搐喷水| 不卡一级毛片| 成人三级黄色视频| 亚洲一级一片aⅴ在线观看| 一级毛片久久久久久久久女| 久久久久国内视频| 狂野欧美白嫩少妇大欣赏| 97超级碰碰碰精品色视频在线观看| 国产精品久久电影中文字幕| 联通29元200g的流量卡| 日本熟妇午夜| 午夜福利在线观看免费完整高清在 | 成人特级黄色片久久久久久久| 亚洲 国产 在线| 草草在线视频免费看| 久久久久久大精品| 精品免费久久久久久久清纯| 一本久久中文字幕| 别揉我奶头~嗯~啊~动态视频| 欧美另类亚洲清纯唯美| 欧美最黄视频在线播放免费| 国产伦在线观看视频一区| 国产精品亚洲一级av第二区| 亚洲乱码一区二区免费版| 免费av毛片视频| 超碰av人人做人人爽久久| 日本一二三区视频观看| 免费看日本二区| 淫妇啪啪啪对白视频| 黄色欧美视频在线观看| 一个人看视频在线观看www免费| 在线天堂最新版资源| 午夜福利在线在线| 亚洲综合色惰| 乱人视频在线观看| 午夜福利视频1000在线观看| 国产日本99.免费观看| 亚洲在线自拍视频| 国产成人一区二区在线| 99久久久亚洲精品蜜臀av| 99久久精品热视频| a级毛片a级免费在线| 精品久久久久久久末码| aaaaa片日本免费| 久久国产精品人妻蜜桃| 亚洲性久久影院| 成年女人看的毛片在线观看| 久久久午夜欧美精品| 美女高潮喷水抽搐中文字幕| 高清在线国产一区| 美女被艹到高潮喷水动态| 亚洲中文字幕日韩| 婷婷色综合大香蕉| 国产精品日韩av在线免费观看| 免费人成视频x8x8入口观看| 国产麻豆成人av免费视频| 岛国在线免费视频观看| 我的老师免费观看完整版| 麻豆av噜噜一区二区三区| 国产精品无大码| 精品国内亚洲2022精品成人| 91久久精品国产一区二区三区| 99久久久亚洲精品蜜臀av| 日本 av在线| 99久久精品热视频| 99九九线精品视频在线观看视频| 欧美xxxx黑人xx丫x性爽| 国产视频内射| 99热精品在线国产| 亚洲人成网站在线播放欧美日韩| 一级a爱片免费观看的视频| 变态另类丝袜制服| 精品久久久久久久末码| 最新在线观看一区二区三区| 国产av一区在线观看免费| 久久久久免费精品人妻一区二区| 国产熟女欧美一区二区| 一级a爱片免费观看的视频| 国国产精品蜜臀av免费| 国模一区二区三区四区视频| 色视频www国产| 国产精品一区二区免费欧美| 国产毛片a区久久久久| 简卡轻食公司| 别揉我奶头 嗯啊视频| 亚洲无线观看免费| 全区人妻精品视频| 一级a爱片免费观看的视频| 欧美日韩中文字幕国产精品一区二区三区| 高清毛片免费观看视频网站| 成人二区视频| 国产激情偷乱视频一区二区| 中文亚洲av片在线观看爽| 精品一区二区免费观看| 免费不卡的大黄色大毛片视频在线观看 | 99在线视频只有这里精品首页| 欧美高清性xxxxhd video| 99热精品在线国产| 国产精品1区2区在线观看.| 久久人人精品亚洲av| 女生性感内裤真人,穿戴方法视频| 性插视频无遮挡在线免费观看| 欧美xxxx性猛交bbbb| 天堂网av新在线| 国产毛片a区久久久久| 午夜影院日韩av| 亚洲美女视频黄频| 日本免费a在线| 午夜激情福利司机影院| 精品久久久噜噜| 国产老妇女一区| 亚洲欧美日韩无卡精品| 88av欧美| 成人无遮挡网站| 蜜桃久久精品国产亚洲av| 老司机深夜福利视频在线观看| 国产免费一级a男人的天堂| 国产精品一及| 日韩av在线大香蕉| 国产v大片淫在线免费观看| 国产亚洲精品av在线| 波多野结衣高清无吗| 亚洲内射少妇av| 国产欧美日韩一区二区精品| 国产在视频线在精品| h日本视频在线播放| 黄片wwwwww| 69人妻影院| 亚洲欧美精品综合久久99| 人妻制服诱惑在线中文字幕| 午夜免费成人在线视频| 欧美性猛交黑人性爽| 99在线视频只有这里精品首页| 亚洲在线观看片| 亚洲电影在线观看av| 国产男人的电影天堂91| 国产精品嫩草影院av在线观看 | 国产精品98久久久久久宅男小说| 美女免费视频网站| 深爱激情五月婷婷| 在线免费观看不下载黄p国产 | 国产精品,欧美在线| 男女视频在线观看网站免费| 特大巨黑吊av在线直播| 亚洲va日本ⅴa欧美va伊人久久| 国产三级中文精品| 麻豆国产97在线/欧美| 1024手机看黄色片| 观看免费一级毛片| 国产老妇女一区| 色吧在线观看| videossex国产| 99热这里只有精品一区| 熟妇人妻久久中文字幕3abv| 免费看a级黄色片| 乱人视频在线观看| x7x7x7水蜜桃| 赤兔流量卡办理| 国产女主播在线喷水免费视频网站 | 真人做人爱边吃奶动态| 日韩,欧美,国产一区二区三区 | 日本精品一区二区三区蜜桃| 亚洲va在线va天堂va国产| 又爽又黄a免费视频| 精品不卡国产一区二区三区| 久久久成人免费电影| 99riav亚洲国产免费| 亚洲男人的天堂狠狠| 看免费成人av毛片| 国产不卡一卡二| 色播亚洲综合网| 亚洲不卡免费看| 少妇熟女aⅴ在线视频| 3wmmmm亚洲av在线观看| 极品教师在线免费播放| 91久久精品国产一区二区成人| 97超级碰碰碰精品色视频在线观看| 免费看日本二区| 免费av观看视频| 2021天堂中文幕一二区在线观| 中文字幕熟女人妻在线| 18禁黄网站禁片午夜丰满| 亚洲中文日韩欧美视频| 国产男靠女视频免费网站| 欧美性猛交╳xxx乱大交人| 99热这里只有是精品在线观看| 内射极品少妇av片p| 国产不卡一卡二| 成人无遮挡网站| 国产在线精品亚洲第一网站| 日韩欧美 国产精品| 国国产精品蜜臀av免费| 麻豆av噜噜一区二区三区| 久久久久久久久久成人| 精品久久久久久,| av视频在线观看入口| 中文字幕免费在线视频6| 看十八女毛片水多多多| 国产精品久久久久久精品电影| 日韩高清综合在线| 日韩欧美国产在线观看| av.在线天堂| 天堂√8在线中文| 亚洲成人久久爱视频| 91在线观看av| 亚洲成人精品中文字幕电影| 一级毛片久久久久久久久女| 久久草成人影院| 亚洲国产精品合色在线| 少妇被粗大猛烈的视频| 欧美黑人巨大hd| 联通29元200g的流量卡| 欧美日韩综合久久久久久 | 国产精品女同一区二区软件 | 一级a爱片免费观看的视频| 亚洲黑人精品在线| 国产精品久久电影中文字幕| 国产精品亚洲一级av第二区| 亚洲18禁久久av| 少妇人妻精品综合一区二区 | 3wmmmm亚洲av在线观看| 欧美xxxx黑人xx丫x性爽| 国产色爽女视频免费观看| 两个人的视频大全免费| 哪里可以看免费的av片| 亚洲人成伊人成综合网2020| 中文资源天堂在线| 国产三级在线视频| 一个人免费在线观看电影| 亚洲国产精品久久男人天堂| 色综合色国产| 国产精品国产三级国产av玫瑰| 在现免费观看毛片| 麻豆一二三区av精品| 人人妻人人澡欧美一区二区| 欧美成人性av电影在线观看| 成人亚洲精品av一区二区| 欧美区成人在线视频| 波多野结衣高清作品| bbb黄色大片| 亚洲va日本ⅴa欧美va伊人久久| 色噜噜av男人的天堂激情| 一区二区三区免费毛片| 亚洲天堂国产精品一区在线| 欧美日韩黄片免| or卡值多少钱| 女的被弄到高潮叫床怎么办 | 91久久精品国产一区二区成人| 国产大屁股一区二区在线视频| 国产高潮美女av| 听说在线观看完整版免费高清| 久久久国产成人精品二区| 精品不卡国产一区二区三区| 有码 亚洲区| 成人鲁丝片一二三区免费| 三级男女做爰猛烈吃奶摸视频| 国产欧美日韩精品亚洲av| xxxwww97欧美| 网址你懂的国产日韩在线| 亚洲av成人av| 免费看日本二区| 久久精品91蜜桃| 中文字幕熟女人妻在线| 91精品国产九色| 国产免费一级a男人的天堂| 一级毛片久久久久久久久女| 高清在线国产一区| 精品一区二区三区人妻视频| 亚洲最大成人av| 精品久久久久久成人av| 国产精品一区二区性色av| 欧美日韩乱码在线| av黄色大香蕉| 成人av一区二区三区在线看| 日本爱情动作片www.在线观看 | 成人一区二区视频在线观看| 亚洲 国产 在线| 校园人妻丝袜中文字幕| 最近最新中文字幕大全电影3| 又粗又爽又猛毛片免费看| 女人十人毛片免费观看3o分钟| 男女那种视频在线观看| 神马国产精品三级电影在线观看| 免费人成视频x8x8入口观看| 成人一区二区视频在线观看| 少妇人妻一区二区三区视频| 欧美xxxx黑人xx丫x性爽| 久久久久久久久中文| 精品久久久久久,| 美女cb高潮喷水在线观看| 真人做人爱边吃奶动态| 欧美激情在线99| 一本精品99久久精品77| 国产精品久久久久久亚洲av鲁大| 亚洲七黄色美女视频| 欧美一级a爱片免费观看看| 亚洲美女黄片视频| 国产精品久久视频播放| 日韩一区二区视频免费看| 国产高清激情床上av| 一本一本综合久久| 特级一级黄色大片| 亚洲熟妇中文字幕五十中出| 久久草成人影院| 免费观看在线日韩| 校园人妻丝袜中文字幕| 国产色爽女视频免费观看| av天堂在线播放| 一个人观看的视频www高清免费观看| 国产精品电影一区二区三区| 91av网一区二区| 亚洲最大成人手机在线| 联通29元200g的流量卡| 国产精品一区二区性色av| 在线观看美女被高潮喷水网站| aaaaa片日本免费| av在线亚洲专区| 国产精品一及| 成人二区视频| 国产极品精品免费视频能看的| 我要搜黄色片| 亚洲欧美日韩高清在线视频| 亚洲精品成人久久久久久| 国内精品宾馆在线| 伊人久久精品亚洲午夜| 精品乱码久久久久久99久播| 麻豆av噜噜一区二区三区| 日本-黄色视频高清免费观看| 日本免费a在线| 久久久精品大字幕| 精品久久久久久久久亚洲 | 亚洲第一区二区三区不卡| 欧美性感艳星| 国产高清三级在线| 97人妻精品一区二区三区麻豆| 99久久精品国产国产毛片| 日韩欧美在线乱码| 99久久中文字幕三级久久日本| 久久99热6这里只有精品| 成年女人看的毛片在线观看| 久久久久久久精品吃奶| 中文字幕高清在线视频| 色吧在线观看| 国产v大片淫在线免费观看| 又粗又爽又猛毛片免费看| 欧美成人一区二区免费高清观看| 国产精品人妻久久久久久| 精品午夜福利在线看| 九九爱精品视频在线观看| 久久人人精品亚洲av| 亚洲成人久久爱视频| 亚洲人成网站在线播放欧美日韩| 深爱激情五月婷婷| 黄色欧美视频在线观看| 在线观看一区二区三区| 国产精品美女特级片免费视频播放器| 中国美白少妇内射xxxbb| 一进一出好大好爽视频| 啦啦啦啦在线视频资源| 日本黄大片高清| 一个人观看的视频www高清免费观看| 国产在线男女| 麻豆成人av在线观看| 天堂√8在线中文| 日本在线视频免费播放| 国产男靠女视频免费网站| 最近在线观看免费完整版| 国语自产精品视频在线第100页| 欧美日韩中文字幕国产精品一区二区三区| 一本精品99久久精品77| 精品人妻视频免费看| 在线免费观看不下载黄p国产 | 91在线精品国自产拍蜜月| 小蜜桃在线观看免费完整版高清| 国产精品,欧美在线| 亚洲欧美精品综合久久99| 欧美又色又爽又黄视频| 美女cb高潮喷水在线观看| 亚洲经典国产精华液单| 免费观看人在逋| 波多野结衣巨乳人妻| 亚洲av成人av| 真实男女啪啪啪动态图| 日本 av在线| 又爽又黄无遮挡网站| 午夜福利在线在线| 能在线免费观看的黄片| 不卡一级毛片| 久久九九热精品免费| 婷婷精品国产亚洲av| 观看美女的网站| 国产精品人妻久久久久久| 久久久久久国产a免费观看| 久久久精品大字幕| 99热只有精品国产| 亚洲图色成人| 欧美日韩精品成人综合77777| 很黄的视频免费| 又粗又爽又猛毛片免费看| 亚洲欧美日韩高清专用| 国产精品亚洲一级av第二区| 日韩av在线大香蕉| 欧美激情国产日韩精品一区| 日韩强制内射视频| 日日摸夜夜添夜夜添小说| 国产欧美日韩精品亚洲av| 久久久久久久亚洲中文字幕| 日韩欧美一区二区三区在线观看| 成熟少妇高潮喷水视频| 村上凉子中文字幕在线| xxxwww97欧美| 亚洲精品在线观看二区| 久99久视频精品免费| 国产午夜福利久久久久久| 日韩欧美精品免费久久| 亚洲色图av天堂| 69av精品久久久久久| 特级一级黄色大片| 亚洲性夜色夜夜综合| 国产激情偷乱视频一区二区| 亚洲内射少妇av| 人妻丰满熟妇av一区二区三区| 99热只有精品国产| 最近最新中文字幕大全电影3| 日日撸夜夜添| 老师上课跳d突然被开到最大视频| 日韩中字成人| 老司机午夜福利在线观看视频| 久久人人精品亚洲av| 国产视频内射| 人人妻人人澡欧美一区二区| 日本撒尿小便嘘嘘汇集6| 国产三级中文精品| 三级毛片av免费| 日韩精品青青久久久久久| 能在线免费观看的黄片| 一区二区三区高清视频在线| 国产精品嫩草影院av在线观看 | 深夜精品福利| 国产aⅴ精品一区二区三区波| 一进一出好大好爽视频| 国产av一区在线观看免费| 日韩欧美一区二区三区在线观看| 免费看av在线观看网站| 人人妻人人澡欧美一区二区| 一级黄片播放器| 乱系列少妇在线播放| 全区人妻精品视频| 男插女下体视频免费在线播放| 女同久久另类99精品国产91| 最后的刺客免费高清国语| 精品不卡国产一区二区三区| 桃红色精品国产亚洲av| 欧美区成人在线视频| videossex国产| 国产精品久久久久久久电影| 毛片女人毛片| 国产女主播在线喷水免费视频网站 | 亚洲精华国产精华液的使用体验 | 亚洲精品一卡2卡三卡4卡5卡| 中亚洲国语对白在线视频| 亚洲av二区三区四区| 欧美成人a在线观看| 麻豆成人av在线观看| 嫁个100分男人电影在线观看| 国产熟女欧美一区二区| 给我免费播放毛片高清在线观看| 成年女人看的毛片在线观看| 国产欧美日韩一区二区精品| 欧美日韩黄片免| 能在线免费观看的黄片| 久久精品人妻少妇| 热99re8久久精品国产| 波多野结衣高清无吗| 很黄的视频免费| 国产午夜精品久久久久久一区二区三区 | www日本黄色视频网| 九九爱精品视频在线观看| 一进一出好大好爽视频| 国产精品一区二区免费欧美| 国产亚洲91精品色在线| 国产精品野战在线观看| 露出奶头的视频| 欧美激情国产日韩精品一区| 99国产精品一区二区蜜桃av| 老熟妇乱子伦视频在线观看| 日韩高清综合在线| 动漫黄色视频在线观看| 午夜爱爱视频在线播放| 午夜免费成人在线视频| 极品教师在线视频| 日韩精品有码人妻一区| 国产在视频线在精品| 十八禁网站免费在线| 午夜免费男女啪啪视频观看 | 亚洲七黄色美女视频| 精品福利观看| a在线观看视频网站| 99久久久亚洲精品蜜臀av| 18禁黄网站禁片午夜丰满| 国产精品无大码| 内地一区二区视频在线| 高清日韩中文字幕在线| 成人美女网站在线观看视频| 欧美一区二区国产精品久久精品| 中文资源天堂在线| 岛国在线免费视频观看| 又紧又爽又黄一区二区| 国产淫片久久久久久久久| 婷婷精品国产亚洲av在线| 又爽又黄无遮挡网站| 久久婷婷人人爽人人干人人爱| 精品一区二区三区av网在线观看| a级毛片免费高清观看在线播放| 久久久久免费精品人妻一区二区| 国产高清有码在线观看视频| 国产老妇女一区| 欧美色视频一区免费| 国产三级中文精品| 精品午夜福利视频在线观看一区| 99热网站在线观看| a级一级毛片免费在线观看| 天美传媒精品一区二区| 国产精品人妻久久久久久| 在线观看美女被高潮喷水网站| 精品福利观看| 麻豆国产av国片精品|