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

    考慮支撐劑壓實和嵌入作用的滑脫效應及滲流機制

    2021-11-03 07:03:38成巧耘李波波李建華
    煤田地質與勘探 2021年5期
    關鍵詞:支撐劑滲流滲透率

    成巧耘,李波波,2,3,李建華,高 政,王 斌

    考慮支撐劑壓實和嵌入作用的滑脫效應及滲流機制

    成巧耘1,李波波1,2,3,李建華1,高 政1,王 斌1

    (1. 貴州大學 礦業(yè)學院,貴州 貴陽 550025;2. 貴州大學 喀斯特地區(qū)優(yōu)勢礦產(chǎn)資源高效利用國家地方聯(lián)合工程實驗室,貴州 貴陽 550025;3. 貴州大學 貴州省非金屬礦產(chǎn)資源綜合利用重點實驗室,貴州 貴陽 550025)

    水力壓裂技術是煤礦瓦斯災害防治與煤層氣開采的關鍵技術之一,在實施水力壓裂過程中,支撐劑的嵌入往往會誘發(fā)煤儲層裂隙寬度的一系列變化。其中,滑脫效應的強度與滲透率的變化主要由裂隙寬度決定。因此,支撐劑嵌入將影響水力壓裂技術的有效性。為探究水力壓裂背景下氣體的滑脫效應與煤的滲流規(guī)律,采用赫茲接觸理論量化支撐劑的嵌入深度,并構建支撐劑與有效應力綜合作用的氣體滑脫系數(shù)計算方程與滲透率模型。結果表明:在不同瓦斯壓力下,煤的滲透率隨有效應力的增大先減小后趨于平緩;恒定有效應力條件下,瓦斯壓力越低,滲透率相對越高;且鋪置多層砂的增透效果相對鋪置單層砂的增透效果更佳;兩種鋪置條件下,滑脫因子在不同瓦斯壓力下呈相同的變化趨勢,均隨支撐劑嵌入深度的增大而增大;不同形態(tài)裂縫的滑脫因子均隨著有效應力的增大而增大,其中球形裂縫的滑脫因子最大,圓柱形次之,狹縫形最小。同時,不同形態(tài)裂縫煤的滲透率均隨有效應力的增加而減小,而3種形態(tài)裂縫滲透率之間的大小關系與滑脫因子大小關系一致;考慮到有效應力與支撐劑對裂縫寬度的貢獻,構建了考慮支撐劑和有效應力綜合作用的裂隙滲透率模型,并通過公開發(fā)布的試驗數(shù)據(jù)驗證其合理性。研究結果將有助于水力壓裂技術在煤礦瓦斯災害防治與煤層氣開采中的進一步應用。

    煤;支撐劑嵌入;滲透率;裂縫形態(tài);滑脫效應

    經(jīng)過十余年的商業(yè)化發(fā)展,我國煤層氣產(chǎn)業(yè)已經(jīng)初具規(guī)模[1]。但煤層是一種致密多孔的非常規(guī)儲層,滲透率極低,導致開采作業(yè)過程中總存在難以預測的安全隱患。為提高煤層透氣性,通常采用水力壓裂技術來輔助煤層氣抽采,該技術實施過程中將向裂縫中充入支撐劑使裂隙進一步發(fā)育,從而提高滲透率[2]。但在有效應力影響下,煤內部孔裂隙被壓縮,支撐劑普遍嵌入煤中使增透效果減弱[3]。同時,在這種緊密多孔的介質中,氣體滑脫效應的影響也不應該被忽視[4-5],且煤本身的吸附性及孔裂隙形態(tài)的多樣性將使?jié)B流演化規(guī)律十分復雜。而煤的滲流規(guī)律對于煤層氣開發(fā)與煤炭安全開采至關重要[6]。因此,研究有效應力作用下不同物理形態(tài)裂縫中煤的滑脫行為和滲流機理十分必要。

    在復雜的儲層環(huán)境中,有效應力對滲透率的影響通常起主導作用[7],祝捷等[8]通過滲流試驗發(fā)現(xiàn)有效應力是使煤樣變形的重要因素,在有效應力增大過程中,煤的裂隙空間被壓縮,導致滲流通道變窄而影響煤層透氣性[9]。同時,滑脫效應強度與裂隙寬度的變化相關,對于低滲煤層,即使在10 MPa的高壓下也不能忽略滑脫效應對滲透率的影響[10]。當孔隙特征長度與氣體分子平均自由程十分接近時,滑脫效應將使氣測滲透率大于絕對滲透率[11]。此外,由于煤是一種多孔介質,主要通過范德華分子作用力來吸附氣體分子[12],煤基質吸附氣體分子后將造成基質膨脹[13],其膨脹變形大小可用煤的表面能進行量化[14],該膨脹變形將使煤的裂隙寬度減小,進而降低滲透率。由于煤中有不同物理形態(tài)的裂縫,這些裂縫是煤層氣的主要儲存空間和運移通道。大量研究表明,裂縫物理形態(tài)將在一定程度上影響瓦斯?jié)B流規(guī)律[15],吸附實驗發(fā)現(xiàn)煤中多種形態(tài)的裂縫對滲透率影響不盡相同[16]。

    水力壓裂過程中,在有效應力影響下,支撐劑將逐漸被壓密、隨后嵌入煤中[17],M. Masowski等[18]利用成像方法使支撐劑嵌入可視化,發(fā)現(xiàn)支撐劑嵌入是降低裂縫寬度和滲透率的嚴重問題。為探明支撐劑嵌入作用下煤的滲透特性,Zheng Wenbo等[19]對4種不同類型的支撐劑進行改性粉碎性試驗,推導出了考慮支撐劑嵌入的經(jīng)驗方程,發(fā)現(xiàn)支撐劑嵌入深度隨應力的增大而增大,使裂縫寬度相對減小[20]。Zhao Junlong等[21]討論了有效應力與滑脫效應對滲透率的影響,發(fā)現(xiàn)有效應力增大,滑脫效應將在一定程度上增大。滑脫效應的強弱由滑脫因子表示,其大小受儲層性質及氣體性質的影響[22]。Wang Gongda等[23]通過裂隙寬度的變化建立了滑脫因子計算模型,對應力影響下滑脫效應的變化規(guī)律進行分析。綜上可知,煤的裂隙寬度變化將對滲透特性與滑脫效應強弱影響復雜。

    綜上,國內外眾多學者對煤的滲透率影響因素開展了大量研究,建立的數(shù)學模型可量化計算應力與支撐劑嵌入對滲透率的貢獻率。但考慮應力與支撐劑壓嵌綜合作用的滑脫效應與滲透率演化規(guī)律的研究卻鮮見報道。同時,針對存在支撐劑嵌入影響下不同形態(tài)裂縫的滑脫效應變化規(guī)律及滲流機理還需進一步闡明。因此,筆者基于前期研究認識[24],以水力壓裂為背景,進一步探討支撐劑壓實和嵌入作用對滑脫效應的影響及其對滲透率的貢獻情況,構建支撐劑與有效應力耦合下煤的滑脫系數(shù)計算方程和滲透率模型,并通過公開發(fā)表的數(shù)據(jù)對滲透率模型進行驗證。并分析討論不同形態(tài)裂縫中氣體的滑脫效應及煤的滲流機理,以期為水力壓裂技術的有效實施和瓦斯災害防治提供理論支持。

    1 滲透率計算模型

    1.1 考慮支撐劑壓實和嵌入作用的滲流機制

    煤吸附氣體后,將使煤基質產(chǎn)生膨脹變形,此時煤基質寬度的變化量Δs可表示為[25]:

    式中:0為初始基質寬度,μm;s為吸附誘導的膨脹應變。

    通常,吸附誘導的膨脹應變可用變形的Langmuir方程[26]來表示:

    式中:max為最大膨脹應變,%;L為Langmuir壓力,MPa;為瓦斯壓力,MPa。

    基于立方模型的幾何特征[27],初始基質寬度與初始裂隙寬度之間的關系為:

    式中:0為初始孔隙率,%;0為初始裂隙寬度,μm。

    煤吸附氣體后,煤基質產(chǎn)生的膨脹變形量將影響裂隙寬度的大小。通常,可引入修正因子來描述氣體吸附膨脹作用對裂縫寬度變化的貢獻率[28]:

    式中:Δs為吸附膨脹作用下裂隙寬度的變化量,μm;為吸附變形修正因子,在0~1范圍內取值。

    將式(1)—式(3)代入式(4)中,可得吸附作用下裂隙寬度的變化量:

    同時,有效應力會直接壓縮煤的裂隙空間,改變裂隙壓縮性大小并使裂縫寬度發(fā)生變化。此時裂隙體積模量可用裂隙壓縮性系數(shù)f表示,基于f=1/f,有效應力作用下裂隙寬度的變化量為:

    式中:Δe為有效應力作用下裂隙寬度的變化量,μm;e為有效應力引起的斷裂應變;f為煤的裂隙壓縮性系數(shù),MPa–1;Δ為應力變化量,MPa;Δ為瓦斯壓力變化量,MPa;為Biot系數(shù),通常取1。

    當煤中壓入支撐劑后,裂隙寬度將隨支撐劑的嵌入進一步發(fā)生變化,基于赫茲接觸理論[29],其嵌入深度可量化為:

    式中:為嵌入深度,μm;為支撐劑半徑,且=02[30],μm;1、2分別為支撐劑、煤的泊松比;1、2分別為支撐劑、煤的彈性模量,MPa。

    煤是一個裂隙與孔隙并存的復雜網(wǎng)絡體,如圖1所示。假設煤的裂隙為狹縫形,且氣體僅在軸方向上流動。在水力壓裂過程中,有效應力、吸附膨脹作用下裂隙變形量為:

    式中:Δw為裂隙寬度的變化量,μm。

    同時,煤裂隙寬度與滲透率之間的關系密不可分,裂隙寬度的改變將影響滲透率的大小。假設煤各向同性,則滲透率可表示為[31-32]:

    式中:0為初始滲透率,μm2。

    基于式(5)—式(10)可進一步得到有效應力與支撐劑綜合作用的滲透率模型:

    1.2 考慮支撐劑壓實和嵌入作用的氣體滑脫效應

    煤具有多孔致密的特點,當流動在煤中氣體分子的平均自由路徑接近納米尺度的孔徑時,滑脫效應將影響煤的滲透特性,使其表觀滲透率大于絕對滲透率,兩者之間的關系可用Klinkenberg方程[33]表示:

    式中:abs為絕對滲透率,μm2;為滑脫因子,MPa。

    相關研究表明,滑脫因子決定滑脫效應影響的強弱,其大小受儲層性質和氣體性質的影響[22]:

    式中:為裂隙寬度,μm;為氣體動力學黏度,Pa·s;為通用氣體常數(shù),J/(mol·K);為溫度,;g為氣體的摩爾質量,kg/mol;為常數(shù),通常取0.9。

    一般情況下,=2e[34],e為有效裂隙半徑。由于本文將裂隙形態(tài)假設為狹縫形,則對于狹縫形裂縫,有效裂隙半徑可根據(jù)狹縫形裂縫的幾何體積進行計算:

    式中:sli為狹縫形裂縫體積,μm3;e,slisli、sli分別為狹縫形裂縫的有效裂隙半徑、寬度和長度,μm。

    狹縫形裂縫在有效應力、吸附膨脹作用、支撐劑嵌入影響下,與SI Leilei等[34]的研究類似,其裂縫幾何體積為:

    式中:為應力引起的變形;為吸附膨脹誘導的變形;為支撐劑嵌入引起的變形。

    根據(jù)式(15)可得有效應力、吸附膨脹作用、支撐劑嵌入影響下的狹縫體積表達式:

    式中:0,sli為狹縫形裂縫初始體積,μm3,且0,sli=20slisli;0為初始裂隙半徑,μm。

    將式(16)代入式(14)中,可得狹縫形裂縫的有效裂隙半徑:

    因此,狹縫形裂縫的有效裂隙寬度為:

    式中:sli為狹縫形裂縫的有效裂隙寬度,μm。

    將sli值代入式(13)可得支撐劑與有效應力綜合作用的氣體滑脫系數(shù)計算方程:

    將式(19)代入式(12)可得考慮支撐劑壓嵌與滑脫效應綜合作用的滲透率模型:

    2 模型驗證與分析

    2.1 考慮支撐劑壓實與嵌入作用的滲透率演化規(guī)律

    Wu Yanting等[35]通過4組對比試驗對滑脫效應與支撐劑影響下的裂隙滲透率演化規(guī)律進行研究,該試驗與本文所建立的滲透率模型具有相同的約束條件與所求目標,故選取其中條件1(未鋪置支撐劑)、條件3(鋪置多層砂)、條件4(鋪置單層砂)的實驗數(shù)據(jù)對考慮支撐劑與有效應力綜合作用下的滲透率模型進行驗證。其中,將未鋪置支撐劑條件下滲透率模型中的嵌入深度考慮為0,經(jīng)計算后得到實測值與模型值之間的關系曲線(圖2)。驗證過程中,3種鋪置條件下滲透率模型擬合參數(shù)取值情況及平均絕對偏差百分比(AAD%)見表1,模型引用參數(shù)見表2。

    表1 滲透率模型擬合參數(shù)

    注:AAD*數(shù)據(jù)來自Wu Yanting等[35]。

    表2 滲透率模型引用參數(shù)

    表1中擬合參數(shù)的取值范圍為0

    由圖2可知,3種鋪置條件下的滲流模型曲線與實測數(shù)據(jù)點都有較好的匹配性,能較好反映有效應力作用下的滲流規(guī)律。對比未鋪置支撐劑與鋪置支撐劑的滲流數(shù)據(jù),發(fā)現(xiàn)鋪置支撐劑的滲透率顯著大于未鋪置支撐劑下的滲透率,且從圖2b、圖2c可以看出,除瓦斯壓力在0.54 MPa條件下的滲透率外,其他瓦斯壓力條件下鋪置多層砂的滲透率顯著大于鋪置單層砂的滲透率。在不同瓦斯壓力下,滲透率實測值與模型計算值都隨有效應力的增大而減小,以圖2c中瓦斯壓力為0.54 MPa下的滲流曲線為例,有效應力從1 MPa增大到5 MPa的過程中,相鄰有效應力的變化使?jié)B透率分別減小15.49× 10–3、14.61×10–3、13.74×10–3、12.90×10–3μm2,不難看出有效應力的繼續(xù)增大將使?jié)B透率減小趨勢變緩。究其原因可知,有效應力增加的初期,裂隙應力敏感性較強,具體表現(xiàn)為裂隙寬度快速減小,滲透率變化明顯,而有效應力繼續(xù)增大,應力敏感性降低,使?jié)B透率減小趨勢變緩[40]。

    圖2 不同鋪置條件下滲透率實測值與模型值對比關系曲線

    為探究不同瓦斯壓力作用下3種鋪置條件滲透率的演化規(guī)律,繪制不同有效應力下的滲透率與瓦斯壓力的變化關系曲線(圖3):在3種鋪置條件下,滲透率與瓦斯壓力間有明顯的變化關系。對比三種鋪置條件下滲透率的變化趨勢,整體上都隨瓦斯壓力的增加呈減小趨勢,但鋪置條件不同,在各個瓦斯壓力階段的變化程度也不盡相同。且在相同有效應力下,瓦斯壓力越小,滲透率越大,可能是由于滑脫效應在低瓦斯壓力下起到了顯著提高滲透率的作用[41]。但滑脫效應在此時并不會對滲透率的變化起主導作用,因此,在有效應力影響下,無論瓦斯壓力有多小,滲透率變化曲線也不會呈上升趨勢。

    2.2 支撐劑壓嵌下的氣體滑脫效應

    為了進一步探究有效應力作用下支撐劑嵌入對滑脫效應的影響規(guī)律,通過式(19)計算鋪置多層砂與鋪置單層砂下支撐劑嵌入與滑脫因子的大小關系,如圖4所示,由圖中可知:在2種鋪置條件下,滑脫因子都有相同的變化趨勢,都隨支撐劑嵌入深度的增大而增大。究其原因可知,滑脫因子與有效裂隙寬度呈負相關,因此,當嵌入深度逐漸增大時,有效裂隙寬度逐漸減小,使滑脫因子在支撐劑嵌入過程中呈逐漸增大的趨勢。同時,從圖中還可以看出,不同瓦斯壓力下的值不盡相同,但在鋪置單層砂條件下,不同壓力下值大小關系顯著。且嵌入深度在0.002~0.003 mm時,滑脫因子隨瓦斯壓力的增大而減小,但隨著嵌入深度越來越大,滑脫因子與瓦斯壓力間的變化規(guī)律逐漸復雜。在鋪置多層砂條件下,滑脫因子與瓦斯壓力間更沒有單一的增加或減少關系。究其原因可知,當嵌入深度較小時,裂隙寬度變化不大,瓦斯壓力越小,滑脫效應越明顯,但嵌入深度逐漸增大后,裂隙變形不規(guī)律,使滑脫因子隨瓦斯壓力變化也不規(guī)律。同理,在鋪置多層砂條件下,孔隙內部變化復雜,反而在瓦斯壓力相對較大時滑脫因子較大。

    圖4 不同鋪置條件下支撐劑嵌入對氣體滑脫特性的影響

    3 討論

    3.1 滑脫效應對滲透率的影響

    為探究滑脫效應對滲透率的影響機制,分別利用式(11)與式(20)計算鋪置單層砂時,不同瓦斯壓力條件下考慮滑脫效應與不考慮滑脫效應作用下的煤滲透率,如圖5所示。

    無論在低氣壓還是中低氣壓條件下,滑脫效應對滲透率的影響是顯著的,考慮滑脫效應影響的滲透率普遍大于不考慮滑脫效應影響的滲透率。對比圖5a與圖5b不難看出,低氣壓下滑脫效應作用顯著,而隨著瓦斯壓力的增大,滑脫效應越不明顯,這與前人研究結果[41]一致。

    3.2 不同形態(tài)裂縫煤的滲透率演化規(guī)律

    煤儲層具有復雜的裂隙系統(tǒng),P. Eisenklam[42]和Nie Baisheng等[43]將煤孔隙分為圓柱形、狹縫形、球形等。裂縫形態(tài)的差異使有效孔隙半徑各不相同,由于有效裂隙寬度受有效孔隙半徑的影響,導致不同形態(tài)裂縫下的滑脫效應及其對滲透率的影響具有差異性[44]。因此,探討不同物理形態(tài)裂縫中氣體的滑脫行為和滲流機理具有實際意義。圖6為狹縫形、圓柱形、球形三種形態(tài)裂縫的幾何簡圖。

    目前,針對滲透率的研究,通常將裂隙假設為狹縫形[34]。實際上,煤中裂縫形態(tài)還包括圓柱形和球形等。同理,不同形態(tài)裂縫下的滑脫效應如下:

    其中,圓柱形和球形裂隙中的有效裂隙寬度可表示為:

    此時,圓柱形和球形裂縫中的滲透率模型為:

    由式(20)、式(25)、式(26)可知,狹縫形、圓柱形、球形3種形態(tài)裂縫的滲透率形式上一致,不同的是3種形態(tài)裂縫的滑脫因子大小不同,導致3種孔的滲透率存在差異。將狹縫形孔中的各種參數(shù)取值代入式(25)—式(26),可得不同形態(tài)裂縫煤的滲透率。由于不同瓦斯壓力下、不同鋪置條件下的滲透率變化趨勢一致,故選取鋪置單層砂條件下,瓦斯壓力為0.54、0.90、1.76 MPa的3條曲線為例,3種孔的滲透率具體大小關系如圖7所示。

    由圖7可知:在不同瓦斯壓力下,3種形態(tài)裂縫的滲透率都隨有效應力的增大而減小,且在0.54 MPa下的滲透率最大,1.76 MPa下的滲透率最小,且瓦斯壓力越大,不同形態(tài)裂縫滲透率的值相差越小??赡苁堑蜌鈮簵l件下,滑脫效應對滲透率起到了積極作用。相同瓦斯壓力下,3種形態(tài)裂縫的滲透率趨勢上比較接近,其大小關系為:球形>圓柱形>狹縫形。

    圖7 有效應力增大過程中不同形態(tài)裂縫的滲透率變化曲線

    3.3 裂縫形態(tài)對氣體滑脫效應的影響

    由式(13)可知,滑脫因子的大小取決于裂隙寬度,而狹縫形、圓柱形、球形3種形態(tài)裂縫的有效孔隙半徑不同,加之又考慮支撐劑嵌入對裂縫寬度造成的影響,使3種形態(tài)裂縫的有效裂隙寬度不盡相同。將狹縫形孔條件下的各種參數(shù)代入另外2種形態(tài)裂縫中可計算出瓦斯壓力為0.54 MPa下3種形態(tài)裂縫的有效裂隙寬度的大小,得到嵌入深度與有效裂隙寬度之間的關系(圖8)。

    圖8 支撐壓嵌作用下裂隙寬度變化趨勢線

    由圖8可知,隨著支撐劑嵌入深度的增加,3種形態(tài)裂縫的有效裂隙寬度逐漸減小,且同一嵌入深度下,有效裂隙寬度的大小關系為:狹縫形>圓柱形>球形,說明在相同條件下,狹縫形裂縫更發(fā)育,圓柱形次之,球形發(fā)育程度最低。究其原因可知:有效裂隙寬度大小為有效孔隙半徑與支撐劑嵌入深度之差,在同一有效應力下,嵌入深度不變,有效裂隙半徑越大,有效裂隙寬度就越大,由式(18)、式(23)—式(24)可知,3種形態(tài)裂縫的有效裂隙半徑大小關系為:狹縫形>圓柱形>球形,故3種裂隙寬度的大小關系與有效孔隙半徑的大小關系一致。

    由于不同形態(tài)裂縫的有效裂隙寬度不同,在有效應力作用下產(chǎn)生的滑脫效應也存在差異,圖9為不同形態(tài)裂縫的滑脫因子在有效應力作用下的變化趨勢線。由圖中可知,有效應力作用下,3種形態(tài)裂縫的滑脫因子呈普遍上升趨勢,且有效應力大于3 MPa后,滑脫因子的增大趨勢明顯的減緩。究其原因可知,有效應力作用下,煤內部裂隙被壓縮,導致氣體滲流通道變窄,裂隙寬度越來越接近氣體分子的平均自由程[11,41],使氣體分子的平均自由程與流動通道特征維數(shù)的比值變化較大,滑脫因子增加趨勢就越明顯,而隨著有效應力繼續(xù)增大,裂隙寬度變化較小,使滑脫因子增大趨勢減緩。而3種形態(tài)裂縫的滑脫因子在相同有效應力條件下存在以下大小關系:球形>圓柱形>狹縫形,與有效孔隙寬度大小關系恰好相反。

    圖9 有效應力增大過程中不同形態(tài)裂縫滑脫因子的變化趨勢

    4 結論

    a.支撐劑增透效果受鋪置層數(shù)的影響,在不同鋪置條件下,鋪置多層砂的滲透率整體大于鋪置單層砂的滲透率;滲透率隨著有效應力的增大先減小后趨于平緩;恒定有效應力條件下,瓦斯壓力越低,氣體分子的滑脫效應越明顯,滲透率相對越高。

    b.根據(jù)赫茲接觸理論計算得出的嵌入深度除了與支撐劑和煤的力學性質、支撐劑半徑有關外,還與有效應力的大小有關,當其他參數(shù)為恒定值時,支撐劑嵌入深度與有效應力呈正相關趨勢,且支撐劑嵌入深度越大,有效裂隙寬度越小,滲透率越低。

    c.有效應力作用下,不同形態(tài)裂縫的滑脫因子均隨著有效應力的增大而增大,其中球形的滑脫因子最大,圓柱形次之,狹縫形最小。同時,不同形態(tài)裂縫煤的滲透率都隨有效應力的增加而減小,而3種形態(tài)裂縫滲透率之間的大小關系與滑脫因子大小相關。

    d.本文建立的有效應力與支撐劑綜合作用的滲透率模型與實測值吻合度較高,能在合理的誤差精度內反映支撐劑嵌入條件下,滑脫效應、吸附膨脹耦合作用下滲透率的演化規(guī)律。研究結果可用于進一步了解水力壓裂使用過程中煤的滲流率演化規(guī)律,為壓裂過程中煤層氣的高效開采與礦井瓦斯災害防治提供理論基礎。

    [1] 鄒才能,楊智,黃士鵬,等. 煤系天然氣的資源類型、形成分布與發(fā)展前景[J]. 石油勘探與開發(fā),2019,46(3):433–442.

    ZOU Caineng,YANG Zhi,HUANG Shipeng,et al. Resource types,formation,distribution and prospects of coal–measure gas[J]. Petroleum Exploration and Development,2019,46(3):433–442.

    [2] ZHANG Rui,NING Zhengfu,YANG Feng,et al. Impacts of nanopore structure and elastic properties on stress–dependent permeability of gas shales[J]. Journal of Natural Gas Science and Engineering,2015,26:1663–1672.

    [3] CHEN Dong,YE Zhihui,PAN Zhejun,et al. A permeability model for the hydraulic fracture filled with proppant packs under combined effect of compaction and embedment[J]. Journal of Petroleum Science and Engineering,2017,149:428–435.

    [4] Busch A,Schweinar K,Kampman N,et al. Determining the porosity of mud rocks using methodological pluralism[J]. Geological Society London Special Publications,2017,454:15–38.

    [5] ZHAO Huawei,ZHAO Tianyi,NING Zhengfu,et al. Petrophysical characterization of tight oil sandstones by microscale X–ray computed tomography[J]. Marine and Petroleum Geology,2019,102:604–614.

    [6] LYU Qiao,LONG Xinping,RANJITH P G,et al. Experimental investigation on the mechanical behaviours of a low-clay shale under water-based fluids[J]. Engineering Geology,2018,233:124–138.

    [7] 李祥春,黃濤,陳小龍,等. 煤基質變形影響下含瓦斯煤滲透率動態(tài)變化規(guī)律[J]. 天然氣工業(yè),2020,40(1):83–87.

    LI Xiangchun,HUANG Tao,CHEN Xiaolong,et al. Dynamic change laws of the permeability of coal containing gas under the effect of coal matrix deformation[J]. Natural Gas Industry,2020,40(1):83–87.

    [8] 祝捷,唐俊,王琪,等. 受載煤樣滲透率與應變的關聯(lián)性研究[J]. 煤炭學報,2019,44(S2):566–573.

    ZHU Jie,TANG Jun,WANG Qi,et al.Correlation between permeability and strain of coal samples under loading[J]. Journal of China Coal Society,2019,44(S2):566–573.

    [9] 祝捷,王學,于鵬程,等. 有效應力對煤樣變形和滲透性的影響研究[J]. 巖石力學與工程學報,2017,36(9):2213–2219.

    ZHU Jie,WANG Xue,YU Pengcheng,et al. Effects of the effective stress on deformation and permeability of coal[J]. Chinese Journal of Rock Mechanics and Engineering,2017,36(9):2213–2219.

    [10] FINK R,KROOSS B M,GENSTERBLUM Y,et al. Apparent permeability of gas shales Superposition of fluid dynamic and poro-elastic effects[J]. Fuel,2017,199:532–550.

    [11] FENG Yang,HE Zheng,LYU Bin,et al. Experimental investigation about gas transport in tight shales:An improved relationship between gas slippage and petrophysical properties[J]. Energy Fuels,2021,35:3937–3950.

    [12] HAO Dingyi,ZHANG Lei,LI Mingxue,et al. Experimental study of the moisture content influence on CH4adsorption and deformation characteristics of cylindrical bituminous coal core[J]. Adsorption Science & Technology,2018,36(7/8):1512–1537.

    [13] 馬衍坤,王恩元,李忠輝,等. 煤體瓦斯吸附滲流過程及聲發(fā)射特性實驗研究[J]. 煤炭學報,2012,37(4):641–647.

    MA Yankun,WANG Enyuan,LI Zhonghui,et al. Methane sorption and seepage in coal and characteristics of acoustic emission[J]. Journal of China Coal Society,2012,37(4):641–647.

    [14] PAN Zhejun,CONNELL L D. Modelling of anisotropic coal swelling and its impact on permeability behaviour for primary and enhanced coalbed methane recovery[J]. International Journal of Coal Geology,2011,85(3–4):257–267.

    [15] LI Jing,CHEN Zhangxing,WU Keliu,et al. Effect of water saturation on gas slippage in circular and angular pores[J]. AIChE Journal,2018,64:3529–3541.

    [16] XU Hao,TANG Dazhen,MATHEWS,et al. Evaluation of coal macrolithotypes distribution by geophysical logging data in the Hancheng Block,Eastern Margin,Ordos Basin,China[J]. International Journal of Coal Geology,2016,165:265–277.

    [17] 黃炳香,李浩澤,程慶迎,等. 煤層壓裂裂縫內支撐劑的壓嵌特性[J]. 天然氣工業(yè),2019,39(4):48–54.

    HUANG Bingxiang,LI Haoze,CHENG Qingying,et al. Compaction and embedment characteristics of proppant in hydraulic fractures of coal seams[J]. Natural Gas Industry,2019,39(4):48–54.

    [18] MASLOWSKI M,LABUS M. Preliminary studies on the proppant embedment in baltic basin shale rock[J]. Rock Mechanics and Rock Engineering,2021,54:2233–2248.

    [19] ZHENG Wenbo,SILVA S C,TANNANT D D. Crushing characteristics of four different proppants and implications for fracture conductivity[J]. Journal of Natural Gas Science and Engineering,2018,53:125–138.

    [20] 郭建春,盧聰,趙金洲,等. 支撐劑嵌入程度的實驗研究[J]. 煤炭學報,2008,33(6):661–664.

    GUO Jianchun,LU Cong,ZHAO Jinzhou,et al. Experimental research on proppant embedment[J]. Journal of China Coal Society,2008,33(6):661–664.

    [21] ZHAO Junlong,TANG Dazhen,LIN Wenji,et al. Permeability dynamic variation under the action of stress in the medium and high rank coal reservoir[J]. Journal of Natural Gas Science and Engineering,2015,26:1030–1041.

    [22] SI Leilei,LI Zenghua,YANG Yongliang. Evolution characteristics of gas permeability under multiple factors[J]. Transport in Porous Media,2019,127:415–432.

    [23] WANG Gongda,REN Ting,WANG Kai,et al. Improved apparent permeability models of gas flow in coal with Klinkenberg effect[J]. Fuel,2014,128:53–61.

    [24] 李波波,李建華,楊康,等. 考慮支撐劑嵌入作用的煤巖滲透率模型[J]. 中國安全科學學報,2019,29(11):14–19.

    LI Bobo,LI Jianhua,YANG Kang,et al. Research on coal permeability model considering proppant embedding[J]. China Safety Science Journal,2019,29(11):14–19.

    [25] LIU Huihua,RUTQVIST J. A new coal–permeability model:Internal swelling stress and fracture–matrix interaction[J]. Transport Porous Media,2010,82:157–171.

    [26] ROBERTSON E P. Measurement and modeling of sorption induced strain and permeability changes in coal[M]. Department of Energy,USA,2005.

    [27] ZHOU Yinbo,LI Zenghua,YANG Yongliang,et al. Improved porosity and permeability models with coal matrix block deformation effect[J]. Rock Mechanics and Rock Engineering,2016,49:3687–3697.

    [28] CHEN Zhongwei,LIU Jishan,PAN Zhejun,et al. Influence of the effective stress coefficient and sorption-induced strain on the evolution of coal permeability:Model development and analysis[J]. International Journal of Greenhouse Gas Control,2012,8(5):101–110.

    [29] HERTZ H. Ueber die Berührung fester elastischer K?rper[J]. Journal für die Reine und Angewandte Mathematik,1882,92:156–171.

    [30] KLINKENBWRG L J. The permeability of porous media to liquids and gases[J]. API Drilling and Production Practice,1941:200–213.

    [31] PEREA M S A,RANJITH P G,CHOI S K. Coal cleat permeability for gas movement under triaxial,non–zero lateral strain condition:A theoretical and experimental study[J]. Fuel,2013,109:389–399.

    [32] 榮騰龍,周宏偉,王路軍,等. 三向應力條件下煤體滲透率演化模型研究[J]. 煤炭學報,2018,43(7):1930–1937.

    RONG Tenglong,ZHOU Hongwei,WANG Lujun,et al. Coal permeability model for gas movement under the three-dimensional stress[J]. Journal of China Coal Society,2018,43(7):1930–1937.

    [33] TIAN Ye,YU Xiangyu,LI Jun,et al. Scaling law for slip flow of gases in nanoporous media from nanofluidics,rocks,and pore-scale simulations[J]. Fuel,2019,236:1065–1077.

    [34] SI Leilei,LI Zenghua,YANG Yongliang. Influence of the pore geometry structure on the evolution of gas permeability[J]. Transport in Porous Media,2018,123:321–339.

    [35] WU Yanting,PAN Zhejun,ZHANG Dingyu,et al. Experimental study of permeability behaviour for proppant supported coal fracture[J]. Journal of Natural Gas Science and Engineering,2018,51:18–26.

    [36] GAO Chao,XIE Lingzhi,XIE Heping,et al. Estimation of the equivalent elastic modulus in shale formation:Theoretical model and experiment[J]. Journal of Petroleum Science and Engineering,2017,151:468–479.

    [37] 申衛(wèi)兵,張保平. 不同煤階煤巖力學參數(shù)測試[J]. 巖石力學與工程學報,2000,19(增刊1):860–862.

    SHEN Weibing,ZHANG Baoping. Testing study on mechanical parameters of coal[J]. Chinese Journal of Rock Mechanics and Engineering,2000,19(Sup.1):860–862.

    [38] ZHANG Xiangxiang,WANG JianGuo,GAO Feng,et al. Impact of water and nitrogen fracturing fluids on fracturing initiation pressure and flow pattern in anisotropic shale reservoirs[J]. Computers and Geotechnics,2017,81:59–76.

    [39] GASH B W. Measurement of rock properties in coal for coalbed methane production[J]. AAPG Bulletin,1991,75(8):1395–1403.

    [40] 劉永茜,張浪,李浩蕩,等. 含水率對煤層氣滲流的影響[J]. 煤炭學報,2014,39(9):1840–1844.

    LIU Yongqian,ZHANG Lang,LI Haodang,et al. Effect of moisture content on CBM seepage[J]. Journal of China Coal Society,2014,39(9):1840–1844.

    [41] ZHANG Shaojie,SANG Qian,DONG Mingzhe. Experimental study of pressure sensitivity in shale rocks:Effects of pore shape and gas slippage[J]. Journal of Natural Gas Science and Engineering,2021,89:103885.

    [42] EISENKLAM P. The structure and properties of porous materials[J]. Chemical Engineering Science,1960,12(1):71.(未找到)

    [43] NIE Baisheng,LIU Xianfeng,YANG Longlong,et al. Pore structure characterization of different rank coals using gas adsorption and scanning electron microscopy[J]. Fuel,2015,158:908–917.

    [44] WANG Xianglong,PAN Zhejun,WANG Kai,et al. Characterizing the shape,size,and distribution heterogeneity of pore-fractures in high rank coal based on X-ray CT image analysis and mercury intrusion porosimetry[J]. Fuel,2020,282:118754.

    Slippage effect and the seepage mechanism considering the compaction and embedding action of proppant

    CHENG Qiaoyun1, LI Bobo1,2,3, LI Jianhua1, GAO Zheng1, WANG Bin1

    (1. College of Mining, Guizhou University, Guiyang 550025, China; 2. National & Local Joint Laboratory of Engineering for Effective Utilization of Regional Mineral Resources from Karst Areas, Guizhou University, Guiyang 550025, China; 3. Guizhou Key Laboratory of Comprehensive Utilization of Non-Metallic Mineral Resources, Guizhou University, Guiyang 550025, China)

    Hydraulic fracturing technology is an effective method to improve the efficiency of coalbed methane drainage, but the universal embedding of proppant under stress will change the width of the fractures, which will affect the gas slippage effect and the permeability mechanism. In order to explore the gas slippage effect and coal seepage law under the background of hydraulic fracturing, the Hertzian contact theory is used to quantify the embedding depth of proppant, and the gas slippage coefficient calculation equation and permeability model of the combined effect of proppant and effective stress are constructed. The results show that under different gas pressures, the permeability of coal first decreases with the increase of effective stress and then tends to be flat. Under the constant effective stress conditions, the lower the gas pressure, the higher the permeability; multiple layers of sand are laid. The anti-reflective effect of slab is better than that of single-layer sand. Under the two paving conditions, the slippage factorshows the same changing trend under different gas pressures, and both increase with the increase of the proppant embedding depth. The slippage factors of different forms of cracks all increase with the increase of effective stress, of which the slippage factor of spherical cracks is the largest, followed by the cylindrical cracks, and slit shape is the smallest. At the same time, the permeability of different forms of fractured coal decreases with the increase of effective stress, and the relationship between the permeability of the three forms of fractures is consistent with the relationship between the size of the slippage factor. Taking into account the effective stress and the proppant contribution to the fracture width , a fracture permeability model considering the comprehensive effects of proppant and effective stress was constructed, and its rationality was verified through publicly released test data. The conclusions obtained will help the further application of hydraulic fracturing technology in coalbed methane drainage.

    coal; proppant embedding; permeability; pore shape; slippage effect

    移動閱讀

    語音講解

    X936

    A

    1001-1986(2021)05-0088-10

    2021-04-02;

    2021-05-11

    國家自然科學基金項目(52064007,51804085);貴州省科學技術基金項目(黔科合基礎-ZK〔2021〕重點052)

    成巧耘,1996年生,女,重慶萬州人,碩士研究生,從事巖石力學、礦山安全與災害防治方面的研究工作. E-mail:654877 982@qq.com

    李波波,1985年生,男,貴州修文人,博士,教授,博士生導師,從事巖石力學、礦山安全與災害防治方面的教學與研究工作. E-mail:bbli@gzu.edu.cn

    成巧耘,李波波,李建華,等. 考慮支撐劑壓實和嵌入作用的滑脫效應及滲流機制[J]. 煤田地質與勘探,2021,49(5):88–97. doi: 10.3969/j.issn.1001-1986.2021.05.010

    CHENG Qiaoyun,LI Bobo,LI Jianhua,et al. Slippage effect and the seepage mechanism considering the compaction and embedding action of proppant[J]. Coal Geology & Exploration,2021,49(5):88–97. doi: 10.3969/j.issn.1001-1986. 2021.05.010

    (責任編輯 范章群 郭東瓊)

    猜你喜歡
    支撐劑滲流滲透率
    水平井多簇支撐劑分布數(shù)值模擬
    粉煤灰陶粒石油壓裂支撐劑的制備與表征
    陶瓷學報(2021年1期)2021-04-13 01:33:48
    中煤階煤層氣井排采階段劃分及滲透率變化
    不同滲透率巖芯孔徑分布與可動流體研究
    SAGD井微壓裂儲層滲透率變化規(guī)律研究
    高滲透率風電并網(wǎng)對電力系統(tǒng)失步振蕩的影響
    油氣壓裂支撐劑應用進展及發(fā)展趨勢
    簡述滲流作用引起的土體破壞及防治措施
    河南科技(2014年12期)2014-02-27 14:10:26
    關于渠道滲流計算方法的選用
    河南科技(2014年11期)2014-02-27 14:09:48
    尾礦壩滲流計算及排滲設計
    金屬礦山(2013年6期)2013-03-11 16:54:05
    亚洲精品乱码久久久久久按摩| 国产成人a∨麻豆精品| 男人爽女人下面视频在线观看| 伊人久久国产一区二区| 精品酒店卫生间| 一区二区三区精品91| 国产亚洲欧美精品永久| 大香蕉97超碰在线| tube8黄色片| 欧美最新免费一区二区三区| 日本午夜av视频| 精品第一国产精品| 欧美成人午夜精品| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲av综合色区一区| 欧美精品一区二区免费开放| 街头女战士在线观看网站| 亚洲精品久久久久久婷婷小说| 男女免费视频国产| 国产成人精品无人区| 免费女性裸体啪啪无遮挡网站| 久久久久久久久久成人| 久久精品国产a三级三级三级| 免费黄网站久久成人精品| 亚洲伊人色综图| 老女人水多毛片| 两个人免费观看高清视频| 国产成人一区二区在线| 日本免费在线观看一区| 久久免费观看电影| 久久久国产一区二区| 国产精品熟女久久久久浪| 男女午夜视频在线观看 | 国产男女内射视频| 人人妻人人添人人爽欧美一区卜| 性色avwww在线观看| 青春草国产在线视频| 亚洲少妇的诱惑av| 国产淫语在线视频| 人人妻人人爽人人添夜夜欢视频| 国产在线免费精品| 狠狠婷婷综合久久久久久88av| 视频在线观看一区二区三区| 欧美日韩精品成人综合77777| 日韩三级伦理在线观看| 国产1区2区3区精品| 欧美成人精品欧美一级黄| 欧美+日韩+精品| 久久99精品国语久久久| 久久人人爽av亚洲精品天堂| 国产精品国产三级国产专区5o| 少妇的逼水好多| 日本-黄色视频高清免费观看| 极品人妻少妇av视频| 精品少妇内射三级| 亚洲第一av免费看| 大话2 男鬼变身卡| 九色成人免费人妻av| 老司机影院成人| tube8黄色片| 精品国产一区二区三区久久久樱花| 国产国语露脸激情在线看| 日日爽夜夜爽网站| 一个人免费看片子| 欧美精品人与动牲交sv欧美| 有码 亚洲区| 中文乱码字字幕精品一区二区三区| 国产免费又黄又爽又色| 午夜日本视频在线| 成人漫画全彩无遮挡| 麻豆精品久久久久久蜜桃| kizo精华| 美女主播在线视频| 好男人视频免费观看在线| 亚洲图色成人| 你懂的网址亚洲精品在线观看| 桃花免费在线播放| 黄色视频在线播放观看不卡| 两个人免费观看高清视频| 母亲3免费完整高清在线观看 | av天堂久久9| 午夜91福利影院| 久久久久久人妻| 看十八女毛片水多多多| 黑人高潮一二区| 亚洲综合精品二区| 考比视频在线观看| 国产极品天堂在线| 国产成人免费无遮挡视频| 亚洲美女搞黄在线观看| 人人妻人人澡人人爽人人夜夜| 少妇人妻 视频| 好男人视频免费观看在线| 精品亚洲成国产av| 日韩av不卡免费在线播放| 欧美另类一区| 精品99又大又爽又粗少妇毛片| 国产精品一区二区在线观看99| 亚洲av.av天堂| 日本欧美国产在线视频| 蜜桃在线观看..| 国产一区二区三区av在线| 久久99热这里只频精品6学生| 飞空精品影院首页| 久久女婷五月综合色啪小说| 99热网站在线观看| 最近的中文字幕免费完整| 亚洲精品一二三| 一区二区三区精品91| 久久青草综合色| 男人舔女人的私密视频| 亚洲精品aⅴ在线观看| 2018国产大陆天天弄谢| 亚洲精品中文字幕在线视频| 岛国毛片在线播放| xxx大片免费视频| 97精品久久久久久久久久精品| 宅男免费午夜| 精品一区二区三区四区五区乱码 | 午夜福利,免费看| 日本av免费视频播放| 少妇猛男粗大的猛烈进出视频| 99久久人妻综合| 中文字幕另类日韩欧美亚洲嫩草| 高清av免费在线| 欧美成人午夜精品| 国产成人精品福利久久| 国产精品免费大片| 女人被躁到高潮嗷嗷叫费观| 成人亚洲欧美一区二区av| 国产欧美另类精品又又久久亚洲欧美| a级毛色黄片| 欧美日本中文国产一区发布| 极品人妻少妇av视频| 免费日韩欧美在线观看| 亚洲精品美女久久久久99蜜臀 | 又黄又爽又刺激的免费视频.| 边亲边吃奶的免费视频| 91在线精品国自产拍蜜月| 亚洲少妇的诱惑av| 超色免费av| 一区二区三区四区激情视频| 国产日韩欧美亚洲二区| 成人手机av| 久久久久久人妻| 久久女婷五月综合色啪小说| 啦啦啦在线观看免费高清www| 少妇 在线观看| 欧美激情国产日韩精品一区| 国产精品秋霞免费鲁丝片| 人人妻人人添人人爽欧美一区卜| 国产片内射在线| 久久99精品国语久久久| 十分钟在线观看高清视频www| 亚洲av中文av极速乱| 男人操女人黄网站| 极品人妻少妇av视频| 亚洲国产欧美在线一区| 久久人人爽av亚洲精品天堂| 中文字幕av电影在线播放| 亚洲人与动物交配视频| 美女中出高潮动态图| 久久久久久久久久人人人人人人| 国产免费福利视频在线观看| 如日韩欧美国产精品一区二区三区| 久久 成人 亚洲| 51国产日韩欧美| 欧美精品高潮呻吟av久久| 免费观看在线日韩| 久久精品国产自在天天线| 欧美国产精品一级二级三级| 韩国av在线不卡| 国产精品国产av在线观看| 交换朋友夫妻互换小说| 多毛熟女@视频| av女优亚洲男人天堂| 午夜影院在线不卡| 天天操日日干夜夜撸| 国产亚洲最大av| 欧美成人午夜精品| 热99久久久久精品小说推荐| 两性夫妻黄色片 | 2021少妇久久久久久久久久久| 精品国产一区二区三区四区第35| 天天躁夜夜躁狠狠久久av| 在线观看国产h片| a级毛色黄片| 伦理电影大哥的女人| av免费在线看不卡| 少妇猛男粗大的猛烈进出视频| 亚洲精品国产av蜜桃| 人体艺术视频欧美日本| 黑人高潮一二区| 国产黄色视频一区二区在线观看| 久久久欧美国产精品| 一区在线观看完整版| 草草在线视频免费看| 国产av国产精品国产| 飞空精品影院首页| 看免费成人av毛片| 99热网站在线观看| a级毛片黄视频| 欧美激情 高清一区二区三区| 成人黄色视频免费在线看| 丝袜美足系列| 精品人妻一区二区三区麻豆| 90打野战视频偷拍视频| 免费人妻精品一区二区三区视频| √禁漫天堂资源中文www| 曰老女人黄片| 精品视频人人做人人爽| 久久精品国产亚洲av涩爱| 视频在线观看一区二区三区| 一边亲一边摸免费视频| 曰老女人黄片| 日韩欧美一区视频在线观看| 人成视频在线观看免费观看| 国产激情久久老熟女| 黄色视频在线播放观看不卡| 精品人妻在线不人妻| 亚洲精品自拍成人| 丰满乱子伦码专区| 日韩 亚洲 欧美在线| 午夜免费观看性视频| 亚洲激情五月婷婷啪啪| 精品午夜福利在线看| 国产亚洲一区二区精品| 精品福利永久在线观看| 久久韩国三级中文字幕| 蜜桃在线观看..| 十分钟在线观看高清视频www| 性色av一级| 亚洲欧美一区二区三区黑人 | 亚洲成人一二三区av| 男女免费视频国产| 精品国产国语对白av| 国产精品不卡视频一区二区| 一本—道久久a久久精品蜜桃钙片| 国产精品一二三区在线看| 久久人妻熟女aⅴ| 日韩成人av中文字幕在线观看| 巨乳人妻的诱惑在线观看| av福利片在线| 久久99热6这里只有精品| 中文字幕人妻丝袜制服| 国精品久久久久久国模美| 日韩电影二区| 黑人高潮一二区| 在线亚洲精品国产二区图片欧美| 亚洲av中文av极速乱| 夜夜骑夜夜射夜夜干| 99国产精品免费福利视频| 自拍欧美九色日韩亚洲蝌蚪91| 中文字幕免费在线视频6| 亚洲精品乱码久久久久久按摩| 精品久久久精品久久久| 看十八女毛片水多多多| 一区二区三区四区激情视频| 国产精品人妻久久久久久| 色吧在线观看| 久久人人97超碰香蕉20202| 国产极品天堂在线| 51国产日韩欧美| 香蕉丝袜av| 久久99热这里只频精品6学生| 日韩欧美一区视频在线观看| 成人国产av品久久久| 久久 成人 亚洲| 午夜免费男女啪啪视频观看| 精品午夜福利在线看| 久久人人97超碰香蕉20202| 亚洲精品,欧美精品| 久久久久网色| 日本免费在线观看一区| 一二三四中文在线观看免费高清| 久久av网站| 啦啦啦在线观看免费高清www| 亚洲美女搞黄在线观看| 欧美bdsm另类| 久久久久久人人人人人| 日本av免费视频播放| 熟女人妻精品中文字幕| 欧美精品一区二区大全| 免费日韩欧美在线观看| 69精品国产乱码久久久| 亚洲欧美清纯卡通| 90打野战视频偷拍视频| 午夜福利视频精品| 成年人免费黄色播放视频| 亚洲精品国产色婷婷电影| 国产免费又黄又爽又色| 亚洲国产欧美在线一区| 国产日韩欧美在线精品| 中国三级夫妇交换| 亚洲精品一区蜜桃| 人人妻人人澡人人看| a 毛片基地| 精品国产一区二区三区久久久樱花| 亚洲 欧美一区二区三区| 国产午夜精品一二区理论片| 中文精品一卡2卡3卡4更新| 99久久人妻综合| 免费观看a级毛片全部| 亚洲五月色婷婷综合| 男女免费视频国产| 精品第一国产精品| 日日啪夜夜爽| 午夜免费男女啪啪视频观看| 亚洲,欧美精品.| 国产一区二区三区av在线| 久久午夜福利片| 亚洲少妇的诱惑av| 汤姆久久久久久久影院中文字幕| 国产成人免费无遮挡视频| 婷婷色av中文字幕| 青春草视频在线免费观看| 18在线观看网站| 久久国产亚洲av麻豆专区| 少妇精品久久久久久久| 91午夜精品亚洲一区二区三区| 久久精品久久精品一区二区三区| 乱人伦中国视频| 纵有疾风起免费观看全集完整版| 精品第一国产精品| 亚洲伊人色综图| 久久久亚洲精品成人影院| 精品少妇久久久久久888优播| av在线app专区| 视频在线观看一区二区三区| 一二三四中文在线观看免费高清| 一级爰片在线观看| 大香蕉久久成人网| 国产精品久久久久久av不卡| 日韩精品免费视频一区二区三区 | 欧美精品国产亚洲| 亚洲精品中文字幕在线视频| 国产精品偷伦视频观看了| 美女内射精品一级片tv| 国产免费又黄又爽又色| freevideosex欧美| 国产片内射在线| 亚洲五月色婷婷综合| 亚洲av电影在线进入| 又粗又硬又长又爽又黄的视频| 精品一区二区三区视频在线| 国产一区二区三区综合在线观看 | 老女人水多毛片| 亚洲精品,欧美精品| 又黄又粗又硬又大视频| 日韩在线高清观看一区二区三区| 精品久久久久久电影网| 九色亚洲精品在线播放| 久久久久国产网址| 日韩av免费高清视频| 欧美 日韩 精品 国产| av网站免费在线观看视频| xxxhd国产人妻xxx| 国产成人精品一,二区| 男人操女人黄网站| www日本在线高清视频| 捣出白浆h1v1| 免费av中文字幕在线| 久久精品国产a三级三级三级| 日韩av在线免费看完整版不卡| 亚洲高清免费不卡视频| 男人添女人高潮全过程视频| 成人手机av| 久久韩国三级中文字幕| 成年女人在线观看亚洲视频| 欧美日韩精品成人综合77777| 美女国产高潮福利片在线看| 看免费成人av毛片| 久久久国产一区二区| 国产精品蜜桃在线观看| 天天躁夜夜躁狠狠久久av| 99久久综合免费| 久久毛片免费看一区二区三区| 精品久久久久久电影网| 欧美成人午夜精品| 日本-黄色视频高清免费观看| 日本爱情动作片www.在线观看| 男人添女人高潮全过程视频| 国产国语露脸激情在线看| 在线免费观看不下载黄p国产| 99国产综合亚洲精品| 9色porny在线观看| 精品一区二区三卡| 国产成人免费观看mmmm| 在线观看一区二区三区激情| 亚洲综合色网址| 成年动漫av网址| 多毛熟女@视频| 久久久欧美国产精品| 国产精品一区二区在线不卡| 国产伦理片在线播放av一区| 久久青草综合色| 国产在线免费精品| 亚洲国产最新在线播放| 国产福利在线免费观看视频| 最近的中文字幕免费完整| 18在线观看网站| 国产69精品久久久久777片| 成人手机av| 亚洲国产欧美日韩在线播放| 久久精品久久久久久噜噜老黄| 中文字幕免费在线视频6| 草草在线视频免费看| videos熟女内射| 一本—道久久a久久精品蜜桃钙片| 一级爰片在线观看| 黑人猛操日本美女一级片| 亚洲精品一二三| 综合色丁香网| 亚洲人成77777在线视频| 最近手机中文字幕大全| 母亲3免费完整高清在线观看 | 久久久久久人人人人人| 美女国产视频在线观看| 久久狼人影院| 国产高清国产精品国产三级| 人妻一区二区av| 色视频在线一区二区三区| 精品酒店卫生间| 午夜91福利影院| 天堂中文最新版在线下载| www日本在线高清视频| a级毛片在线看网站| 建设人人有责人人尽责人人享有的| 如日韩欧美国产精品一区二区三区| 在线亚洲精品国产二区图片欧美| 桃花免费在线播放| 美女中出高潮动态图| 国产亚洲最大av| 国产成人午夜福利电影在线观看| 国产亚洲欧美精品永久| 咕卡用的链子| av女优亚洲男人天堂| 免费在线观看完整版高清| 亚洲伊人久久精品综合| 飞空精品影院首页| 一级,二级,三级黄色视频| 欧美国产精品va在线观看不卡| 成人亚洲欧美一区二区av| 国产成人免费无遮挡视频| 久久毛片免费看一区二区三区| 亚洲精品乱久久久久久| 精品99又大又爽又粗少妇毛片| 午夜日本视频在线| 久久国内精品自在自线图片| 波野结衣二区三区在线| 午夜免费鲁丝| 国产成人午夜福利电影在线观看| 久久久国产一区二区| 欧美精品av麻豆av| 国产成人精品一,二区| 亚洲国产精品专区欧美| 午夜精品国产一区二区电影| 香蕉国产在线看| 亚洲精品第二区| 一级,二级,三级黄色视频| 美女国产高潮福利片在线看| 免费女性裸体啪啪无遮挡网站| 亚洲美女黄色视频免费看| 黄网站色视频无遮挡免费观看| 哪个播放器可以免费观看大片| 日日啪夜夜爽| 一级a做视频免费观看| 深夜精品福利| 美女内射精品一级片tv| 欧美精品人与动牲交sv欧美| 欧美成人精品欧美一级黄| 桃花免费在线播放| 另类亚洲欧美激情| 久久97久久精品| 国产精品久久久久久久电影| 狂野欧美激情性bbbbbb| 欧美人与性动交α欧美精品济南到 | xxx大片免费视频| 97超碰精品成人国产| 啦啦啦视频在线资源免费观看| 99国产综合亚洲精品| 欧美精品亚洲一区二区| 国产 一区精品| 91成人精品电影| 欧美国产精品一级二级三级| 看免费成人av毛片| www日本在线高清视频| 国产黄频视频在线观看| 久久精品国产亚洲av天美| 99久久中文字幕三级久久日本| 欧美变态另类bdsm刘玥| 美女福利国产在线| 母亲3免费完整高清在线观看 | 国产一区有黄有色的免费视频| 亚洲精品国产色婷婷电影| 少妇高潮的动态图| 曰老女人黄片| 日韩中字成人| 国产成人精品一,二区| 亚洲av电影在线观看一区二区三区| 亚洲成av片中文字幕在线观看 | 一本色道久久久久久精品综合| 国产亚洲av片在线观看秒播厂| 亚洲av男天堂| 十八禁高潮呻吟视频| 观看av在线不卡| 亚洲精品日本国产第一区| 国产精品一区二区在线观看99| av国产久精品久网站免费入址| 久久久久久久久久久免费av| 在线观看免费高清a一片| 黄色配什么色好看| 国精品久久久久久国模美| 少妇高潮的动态图| 91精品伊人久久大香线蕉| 久久久久久久久久久免费av| 亚洲成人手机| 国产男女超爽视频在线观看| 午夜老司机福利剧场| 一级毛片黄色毛片免费观看视频| 婷婷色综合www| 寂寞人妻少妇视频99o| 五月天丁香电影| 精品人妻偷拍中文字幕| 少妇 在线观看| 在线观看免费日韩欧美大片| 免费看av在线观看网站| 欧美人与性动交α欧美精品济南到 | 欧美日韩视频高清一区二区三区二| 精品国产一区二区三区久久久樱花| 亚洲五月色婷婷综合| 精品久久久久久电影网| 日韩在线高清观看一区二区三区| 最后的刺客免费高清国语| 午夜日本视频在线| 亚洲少妇的诱惑av| 国产日韩一区二区三区精品不卡| 成人亚洲欧美一区二区av| 麻豆乱淫一区二区| 日韩大片免费观看网站| 在线精品无人区一区二区三| 国产精品嫩草影院av在线观看| 一区在线观看完整版| 精品一区二区三卡| 黄色毛片三级朝国网站| 久久免费观看电影| 日日啪夜夜爽| 黑人欧美特级aaaaaa片| 免费观看在线日韩| 欧美成人精品欧美一级黄| 高清av免费在线| 五月开心婷婷网| 一二三四在线观看免费中文在 | 丝袜在线中文字幕| 在线观看三级黄色| 女人久久www免费人成看片| 亚洲欧美中文字幕日韩二区| 少妇人妻 视频| 国产精品人妻久久久久久| 国语对白做爰xxxⅹ性视频网站| 大香蕉久久成人网| 巨乳人妻的诱惑在线观看| 亚洲,一卡二卡三卡| 精品人妻偷拍中文字幕| 十八禁高潮呻吟视频| 黑人欧美特级aaaaaa片| 在线看a的网站| 欧美日韩亚洲高清精品| 97在线视频观看| 亚洲av电影在线进入| 人妻系列 视频| 日韩中文字幕视频在线看片| 国产探花极品一区二区| 国产精品一区www在线观看| 国产一区二区在线观看av| a级片在线免费高清观看视频| 美国免费a级毛片| 成人手机av| 青春草视频在线免费观看| 毛片一级片免费看久久久久| 久久99精品国语久久久| 亚洲第一av免费看| 91精品伊人久久大香线蕉| 国产极品天堂在线| 天天影视国产精品| 国产1区2区3区精品| freevideosex欧美| 欧美精品高潮呻吟av久久| 日日爽夜夜爽网站| 这个男人来自地球电影免费观看 | 三上悠亚av全集在线观看| 男女国产视频网站| 免费女性裸体啪啪无遮挡网站| 女的被弄到高潮叫床怎么办| 伦理电影免费视频| 精品熟女少妇av免费看| 成年av动漫网址| 丰满迷人的少妇在线观看| 在线看a的网站| 乱人伦中国视频| 国产亚洲最大av| 国产在线免费精品| 久久国内精品自在自线图片| 亚洲精品视频女| 国产欧美日韩一区二区三区在线| 美女国产视频在线观看| 欧美97在线视频| 9色porny在线观看| 成人亚洲欧美一区二区av| 欧美精品亚洲一区二区| 免费观看a级毛片全部| 国产一区有黄有色的免费视频| 一级毛片电影观看| 久久鲁丝午夜福利片| 国产精品久久久久成人av|