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

    基于含水層-鉆孔水量交換的疏水鉆孔水量計(jì)算

    2021-10-26 02:57:22靳德武
    煤炭學(xué)報(bào) 2021年9期
    關(guān)鍵詞:鉆場(chǎng)仰角涌水量

    劉 基,靳德武,王 皓

    (1.中煤科工集團(tuán) 西安研究院有限公司,陜西 西安 710054; 2.陜西省煤礦水害防治技術(shù)重點(diǎn)實(shí)驗(yàn)室,陜西 西安 710054)

    我國(guó)西部地區(qū)侏羅紀(jì)煤炭資源豐富[1],其侏羅系延安組煤層開(kāi)采普遍受到頂板水害的威脅[2-3]。目前,采用鉆孔進(jìn)行大降深、大流量超前預(yù)疏放是煤礦頂板水害主要防治工程措施[4]。在一些特定的礦井地質(zhì)、水文地質(zhì)條件下,超前疏排煤層頂板直接或間接充水含水層水是惟一的防治手段[5]。頂板疏放水鉆孔的設(shè)計(jì)(長(zhǎng)度、仰角、數(shù)量等)直接關(guān)乎鉆探工程量、疏放時(shí)間及疏放水效果[6],同時(shí)為工作面防排水系統(tǒng)設(shè)計(jì)以及制定可控疏放方案提供直接參考。

    井下疏放水是利用鉆孔孔口標(biāo)高低于充水含水層水位標(biāo)高的天然條件,使含水層地下水從鉆孔中自然涌出,從而有效降低充水含水層水位(頭)的方法。由于承壓含水層水頭較大,大降深、大流量的井下疏放水,與地面抽水有很大區(qū)別,常規(guī)井流理論不適用于井下疏放水過(guò)程。目前,有關(guān)煤層頂板疏放水量計(jì)算多采用點(diǎn)源線(xiàn)匯理論[6-8],根據(jù)多目標(biāo)管理模型對(duì)鉆孔參數(shù)進(jìn)行優(yōu)化[9]。該方法一方面適用性未得到充分驗(yàn)證[10],另一方面未考慮疏放水鉆孔內(nèi)復(fù)雜的流態(tài)特征[11],同時(shí),不能解決鉆孔穿透多個(gè)含水層的水量分配問(wèn)題[12],導(dǎo)致計(jì)算精度不高,不能滿(mǎn)足疏放水鉆孔工程優(yōu)化設(shè)計(jì)的要求,因此目前對(duì)疏放水鉆孔參數(shù)的優(yōu)化工作開(kāi)展較少。

    事實(shí)上,傾斜鉆孔疏放水時(shí)鉆孔內(nèi)水流流態(tài)具有管道流特征[11],整個(gè)含水層-鉆孔系統(tǒng)可分解為含水層和鉆孔2個(gè)子系統(tǒng)。含水層子系統(tǒng)可以用滲流方程來(lái)刻畫(huà),而鉆孔子系統(tǒng)可以近似看作為一個(gè)小型管道系統(tǒng),通過(guò)滲流-管流耦合模型來(lái)描述。

    “滲流-管流耦合模型”最早應(yīng)用在供水水文地質(zhì)水資源量評(píng)價(jià)領(lǐng)域[13-14],該模型充分考慮了鉆孔內(nèi)多流態(tài)并存的水流特征[11,13],已廣泛用于水平井[15]、滲流井[16-17]、輻射井[18-19]等取水建(構(gòu))筑物水量計(jì)算中,效果良好。筆者借鑒“滲流-管流耦合模型”理論及能量守恒定律,以含水層-鉆孔間水量交換量為耦合點(diǎn),構(gòu)建煤層頂板復(fù)合充水含水層傾斜疏放水鉆孔定降深放水的含水層-鉆孔系統(tǒng)耦合模型。以母杜柴登礦井井下疏放水工程為例,研究鉆場(chǎng)內(nèi)不同參數(shù)(角度、數(shù)量等)下鉆孔涌水量的變化特征,確定鉆場(chǎng)內(nèi)鉆孔最佳長(zhǎng)度、仰角和數(shù)量等參數(shù),為疏放水工程優(yōu)化設(shè)計(jì)提供技術(shù)支撐。

    1 基于含水層-鉆孔水量交換的疏放鉆孔水量計(jì)算模型

    1.1 含水層子系統(tǒng)滲流模型

    疏放水鉆孔放水時(shí),含水層中地下水流呈現(xiàn)明顯的承壓三維流,其滲流模型可表示為

    (1)

    式中,Kxx,Kyy和Kzz為沿x,y,z軸方向的主滲透系數(shù),這里假設(shè)滲透系數(shù)的主方向和坐標(biāo)軸的方向一致;h為水頭高度;Ss為多孔介質(zhì)的彈性釋水率;t為時(shí)間。

    1.2 鉆孔子系統(tǒng)管流模型

    疏放水鉆孔近似看作一小型管道。當(dāng)管道完全飽和時(shí),其流動(dòng)可以假定為沿管道軸線(xiàn)的一維流。19世紀(jì)50年代末,通過(guò)對(duì)直筒形管道中水流的實(shí)驗(yàn)表明,沿管道的水頭損失與速度水頭和管道長(zhǎng)度直接相關(guān),與管道直徑成反比[20]。Darcy,Weisbach和其他水力學(xué)研究者提出以下適用于紊流和層流的一般方程:

    (2)

    在管道中v=nu,其中n=1,v=u。此時(shí)方程(2)可變?yōu)?/p>

    (3)

    式中,ΔH為沿程水頭損失;f為摩擦因數(shù);l為沿程長(zhǎng)度;d為管道直徑;u為流速;n為孔隙度;v為平均流速;g為重力加速度。

    方程(3)中的摩擦因數(shù)是流量和管道相對(duì)粗糙度的函數(shù)。方程式(3)被稱(chēng)為Darcy-Weisbach方程。方程(3)變形為

    (4)

    管道中的流量為

    (5)

    其中,A為管道的過(guò)水?dāng)嗝娣e,即管道中的流量與含水層滲流中的流量計(jì)算公式形式上是一樣;K′為管道中的“等效滲透系數(shù)”;J為管道中水力梯度。摩擦因數(shù)與雷諾數(shù)Re相關(guān),不同的雷諾數(shù)Re分別對(duì)應(yīng)不同的摩擦因數(shù)[19]。

    (1)當(dāng)Re≤3 000,水流呈層流時(shí),摩擦因數(shù)可用以下公式計(jì)算[20],即

    (6)

    式中,υ為黏滯系數(shù),那么

    (7)

    將式(7)代入式(5)得

    (8)

    (2)當(dāng)Re>3 000,水流呈紊流時(shí),摩擦因數(shù)可用以下公式計(jì)算[21],即

    (9)

    將式(9)代入式(4)[22]得

    (10)

    綜上所述,鉆孔中的流量計(jì)算公式為

    式中,Kc為鉆孔內(nèi)壁管壁的粗糙高度;μ為動(dòng)力黏滯系數(shù);τ為鉆孔內(nèi)壁粗糙度;ρ為水的密度。

    1.3 含水層-鉆孔系統(tǒng)耦合模型

    根據(jù)能量守恒定律,含水層地下水流入鉆孔的總流量等于鉆孔內(nèi)涌水量,即

    ∑Qex+Qp=0

    (11)

    其中,Qex為含水層與管道之間的交換水量,負(fù)號(hào)代表含水層流入管道,即含水層地下水進(jìn)入鉆孔的水量。

    多孔介質(zhì)地下水滲流區(qū)(含水層)與管道之間的水量交換,可根據(jù)兩區(qū)域之間的水頭差進(jìn)行線(xiàn)性計(jì)算,即

    Qex=α(hp-hi)

    (12)

    式中,α為含水層與管道交換區(qū)的水力傳導(dǎo)系數(shù);hp為交換區(qū)管道中的水頭;hi為交換區(qū)含水層中的水頭。

    那么由式(11)可推出

    (13)

    這樣,就將多孔介質(zhì)地下水的滲流與鉆孔中的管流耦合起來(lái),方程(13)即為構(gòu)建的含水層-鉆孔系統(tǒng)復(fù)雜流態(tài)水流計(jì)算耦合模型。

    方程(13)可通過(guò)改進(jìn)管道流計(jì)算程序Conduit Flow Process(CFP)進(jìn)行數(shù)值計(jì)算。CFP為在MODFLOW程序包上開(kāi)發(fā)出來(lái)的,地下水滲流區(qū)根據(jù)MODFLOW進(jìn)行離散化,利用CFP程序包將管道離散為一系列的節(jié)點(diǎn)和管子,兩相鄰節(jié)點(diǎn)用一根管子連接,地下水在相鄰節(jié)點(diǎn)間通過(guò)管子傳輸,地下水滲流區(qū)和管道之間通過(guò)水量交換量進(jìn)行耦合[22]。近年來(lái),CFP程序已經(jīng)開(kāi)始用于層流和紊流的計(jì)算[23],能夠很好地刻畫(huà)巖溶等管道流[24]以及滲流井傾斜輻射管-含水層間復(fù)雜水流[25],模擬計(jì)算效果良好。

    2 模型可靠性驗(yàn)證

    為了驗(yàn)證該耦合模型,筆者采用母杜柴登礦井臨時(shí)煤倉(cāng)放水試驗(yàn)數(shù)據(jù)進(jìn)行模型的可靠性驗(yàn)證。

    2.1 研究區(qū)概況

    母杜柴登礦井位于鄂爾多斯盆地呼吉爾特礦區(qū),礦井地層自下而上由老至新發(fā)育有:三疊系上統(tǒng)延長(zhǎng)組(T3y)、侏羅系中統(tǒng)延安組(J2y)、侏羅系中統(tǒng)直羅組(J2z)、侏羅系中統(tǒng)安定組(J2a)、白堊系下統(tǒng)志丹群(K1zh)和第四系全新統(tǒng)(Q4)(圖1)。主要的含水層為第四系松散層(Q)潛水含水層、白堊系下統(tǒng)志丹群(K1zh)孔隙裂隙承壓含水層、侏羅系直羅組(J2z)裂隙承壓含水層以及侏羅系中延安組(J2y)碎屑巖類(lèi)承壓含水層,具體簡(jiǎn)述如下:

    圖1 地層綜合柱狀Fig.1 Comprehensive stratigraphic column

    (1)第四系松散層潛水含水層。含水層結(jié)構(gòu)疏松,孔隙發(fā)育,透水性良好。含水層水位標(biāo)高1 266.65~1 315.24 m,單位涌水量0.495~12.176 L/(s·m),含水層富水性中等~極強(qiáng),滲透系數(shù)0.398 ~30.625 m/d。

    (2)白堊系下統(tǒng)志丹群孔隙裂隙承壓含水層。該含水層巖芯較為破碎,膠結(jié)程度相對(duì)較差,裂隙較發(fā)育,開(kāi)啟程度較好。含水層水位標(biāo)高1 269.22~1 317.63 m,單位涌水量0.114~0.624 L/(s·m),含水層富水性中等,滲透系數(shù)0.034 ~1.424 m/d。

    (3)侏羅系直羅組孔隙裂隙承壓含水層。該含水層分布較為連續(xù)、穩(wěn)定,同樣是由一套砂巖與泥巖含、隔水層相互疊置的組合結(jié)構(gòu)。含水層水位標(biāo)高1 261.95~1 296.8 m,單位涌水量0.014~0.203 L/(s·m),含水層富水性弱至中等,滲透系數(shù)0.031 ~0.288 m/d。

    (4)侏羅系延安組碎屑巖類(lèi)承壓含水層:延安組為含煤地層,煤層之間巖性以粉細(xì)砂巖、砂質(zhì)泥巖和泥巖為主。砂巖分布不連續(xù),砂質(zhì)泥巖、泥巖厚度大,分布較穩(wěn)定。由于埋藏深,砂巖含水層分布不連續(xù),含水層富水性弱,滲透性差。

    對(duì)礦井充水因素分析認(rèn)為,煤層開(kāi)采的直接充水含水層為侏羅系延安組三段、直羅組一段和二段的復(fù)合含水層,其中直羅組含水層富水性較強(qiáng),水壓較高,這3層含水層都是采前預(yù)疏放的對(duì)象。

    2.2 模型驗(yàn)證

    母杜柴登礦井臨時(shí)煤倉(cāng)放水試驗(yàn)是以TF1放水孔,以TG2,TG3作為井下觀(guān)測(cè)孔。3個(gè)鉆孔均從3-1煤層頂板開(kāi)口,各鉆孔參數(shù)及布設(shè)位置見(jiàn)表1。

    表1 母杜柴登礦井臨時(shí)煤倉(cāng)放水試驗(yàn)鉆孔參數(shù)Table 1 Drilling parameters of drainage test in temporary coal bunker in Muduchaideng Coal Mine

    根據(jù)放水試驗(yàn)數(shù)據(jù),其影響半徑小于5 km,因此,本次建立10 km×10 km的理想模型;由于鉆孔開(kāi)口處為3-1煤層頂板,選擇3-1煤層底板作為模型底板,模型垂向包括3-1煤層、3-1煤頂板延安組,2-2中煤層、2-2中煤頂板直羅組地層。由于本次目的主要為探討傾斜鉆孔的計(jì)算方法,因此建立的數(shù)值模型為均質(zhì)各向異性含水層理想模型。模型邊界概化為四周為變流量邊界,模型中采用GHB邊界處理,頂?shù)装鍨楦羲吔纭?/p>

    臨時(shí)煤倉(cāng)傾斜疏放水鉆孔放水試驗(yàn)期間,地下水流服從三維達(dá)西流。采用下列數(shù)學(xué)模型[26-27]進(jìn)行描述:

    式中,n為第2類(lèi)邊界外法線(xiàn)方向;q(x,y,z,t)為四周側(cè)向補(bǔ)給量;Hc為鉆孔疏放水時(shí)定水頭;Γ1為四周第2類(lèi)流量邊界;Γ2為底板邊界;D為計(jì)算區(qū)范圍。

    本次模擬范圍為10 km×10 km,模型坐標(biāo)為x為0~10 km;y為0~10 km。首先采用100 m×100 m分別平行于x軸、y軸的正交網(wǎng)格對(duì)計(jì)算域進(jìn)行平面剖分,對(duì)鉆孔布設(shè)區(qū)域進(jìn)一步精細(xì)剖分為2 m×2 m的網(wǎng)格,最終將模擬區(qū)剖分為492行×394列;為了精細(xì)刻畫(huà)鉆孔,垂向上采用水平切片方式將地層細(xì)分為22層:1~6層對(duì)應(yīng)實(shí)際的直羅組底部中粗砂巖孔隙-裂隙承壓含水層,7~21層對(duì)應(yīng)實(shí)際的延安組砂巖孔隙-裂隙承壓含水層,22層對(duì)應(yīng)實(shí)際的3-1煤層。模擬時(shí)間為1 d,應(yīng)力期劃分與現(xiàn)場(chǎng)放水試驗(yàn)觀(guān)測(cè)時(shí)間一致,共分為50個(gè)應(yīng)力期,模擬的所有含水層初始水位為未擾動(dòng)的含水層原始水位1 267.702 m。

    以耦合數(shù)學(xué)模型為基礎(chǔ),利用MODFLOW軟件進(jìn)行數(shù)值計(jì)算,以放水試驗(yàn)過(guò)程中地下水水位和水量實(shí)測(cè)數(shù)據(jù)為依據(jù),對(duì)數(shù)值模型進(jìn)行識(shí)別驗(yàn)證,最終識(shí)別的水文地質(zhì)參數(shù)見(jiàn)表2。主放水孔TF1,觀(guān)測(cè)孔TG2和TG3在模型中計(jì)算結(jié)果與實(shí)測(cè)數(shù)據(jù)擬合結(jié)果如圖2所示。由圖2可以看出,地下水水位及水量數(shù)值計(jì)算值和實(shí)際觀(guān)測(cè)值變化趨勢(shì)基本一致,表明該數(shù)值模型合理可靠,所建立的地下水流數(shù)值模型可以用于傾斜鉆孔的涌水量計(jì)算。

    表2 模型各層水文地質(zhì)參數(shù)Table 2 Hydrogeological parameters of each layer

    圖2 臨時(shí)煤倉(cāng)放水試驗(yàn)流量和降深模擬值與實(shí)測(cè)值擬 合曲線(xiàn)Fig.2 Fitting curves of flow rate and drawdown of the simulated and measured value of temporary bunker dewatering test

    3 鉆孔參數(shù)優(yōu)化

    不同的疏放水鉆孔參數(shù)(如單孔角度,單個(gè)鉆場(chǎng)的鉆孔數(shù)量)影響著總涌水量,進(jìn)而決定疏放水持續(xù)時(shí)間和工程費(fèi)用,為此需要開(kāi)展母杜柴登礦井疏放水鉆孔參數(shù)優(yōu)化研究。

    3.1 角 度

    為了探討仰角角度對(duì)鉆孔涌水量的影響,設(shè)計(jì)一個(gè)垂向長(zhǎng)度為120 m的鉆孔,仰角分別為30°,45°,60°,75°,90°,分別計(jì)算不同仰角方案下的單孔涌水量以及單位長(zhǎng)度涌水量(圖3,4)。

    圖3 不同角度下鉆孔單孔涌水量曲線(xiàn)Fig.3 Water inflow curve of single borehole at different angles

    由圖3可知,鉆孔仰角為30°時(shí),單孔穩(wěn)定涌水量約為222 m3/h;鉆孔仰角為45°時(shí),單孔穩(wěn)定涌水量約為191 m3/h;鉆孔仰角為60°時(shí),單孔穩(wěn)定涌水量約為165 m3/h;鉆孔仰角為75°時(shí),單孔穩(wěn)定涌水量約為121 m3/h;鉆孔仰角90°時(shí)涌水量與75°基本一致。根據(jù)母杜柴登礦井現(xiàn)場(chǎng)鉆孔實(shí)際施工情況,45°鉆場(chǎng)涌水量為160~180 m3/h,井底車(chē)場(chǎng)井下施工的一個(gè)仰角60°的生產(chǎn)用水鉆孔,鉆孔終孔涌水量為150~160 m3/h。由此可見(jiàn),耦合模型計(jì)算的單孔涌水量和現(xiàn)場(chǎng)實(shí)際涌水量基本一致,計(jì)算精度有保證。

    由圖4可知,隨著角度的不斷增大,單孔涌水量不斷減少。這是因?yàn)榇咕嘁欢〞r(shí),隨著角度的增大,鉆孔斜長(zhǎng)逐漸變小,有效長(zhǎng)度減少,導(dǎo)致其涌水量變小。然而,隨著角度增大,單位長(zhǎng)度鉆孔涌水量呈現(xiàn)先增加后減小的趨勢(shì)。在仰角為60°時(shí),鉆孔單位涌水量達(dá)到最大。

    圖4 不同角度單位長(zhǎng)度鉆孔涌水量變化曲線(xiàn)Fig.4 Variation curves of borehole water inflow per unit length at different angles

    不同角度下的鉆孔放水時(shí)的第1應(yīng)力期的流場(chǎng)如圖5所示。由圖5可知,地下水流場(chǎng)為以鉆孔為中心的降落漏斗,水流由含水層流向鉆孔,鉆孔內(nèi)水流由鉆孔最頂部向下流,符合現(xiàn)場(chǎng)實(shí)際情況。

    綜上所述,疏放水鉆孔的最佳角度為60°,此時(shí),鉆孔單位長(zhǎng)度涌水量最大,且主要涌水層位為直羅組一段含水層,約占總涌水量的81%(表3)。

    3.2 數(shù) 量

    在鉆孔單位長(zhǎng)度涌水量達(dá)到最大的條件下,即垂向長(zhǎng)度為120 m,仰角為60°時(shí),設(shè)計(jì)6種不同鉆孔數(shù)量組合方案,探討鉆孔數(shù)量對(duì)總涌水量及單位長(zhǎng)度涌水量的影響(圖6~8),6種不同鉆孔組合方案如圖6所示。設(shè)計(jì)鉆場(chǎng)內(nèi)施工6個(gè)鉆孔,方位角以正北為起點(diǎn),以60°遞增均勻布置,其他方案則在此基礎(chǔ)上抽稀即可。

    圖6 組合方案示意及相對(duì)位置Fig.6 Schematic diagram of combination scheme and relative position

    由圖7可知,鉆孔數(shù)量為2個(gè)時(shí),鉆孔穩(wěn)定總涌水量約為268 m3/h;鉆孔數(shù)量為3個(gè)時(shí),鉆孔穩(wěn)定總涌水量約為348 m3/h;鉆孔數(shù)量為4個(gè)時(shí),鉆孔穩(wěn)定總涌水量約為395 m3/h;鉆孔數(shù)量為5個(gè)時(shí),鉆孔穩(wěn)定總涌水量約為443 m3/h;鉆孔數(shù)量為6個(gè)時(shí),鉆孔穩(wěn)定總涌水量約為448 m3/h。鉆孔數(shù)量增加到5個(gè)以后,鉆孔總涌水量基本不再增加。

    圖7 不同方案下鉆孔總涌水量歷時(shí)曲線(xiàn)Fig.7 Curves of total water inflow of boreholes under different schemes

    由圖8可知,隨著鉆孔數(shù)量增加,鉆孔總涌水量不斷增加,但是,單位長(zhǎng)度的鉆孔涌水量呈現(xiàn)逐步減小的趨勢(shì)。由此可見(jiàn),鉆孔之間的干擾程度隨著其數(shù)量的增加不斷增大。

    圖8 不同方案下鉆孔總涌水量對(duì)比Fig.8 Comparison of total water inflow of boreholes under different schemes

    綜合以上考慮,單個(gè)鉆場(chǎng)內(nèi)部鉆孔的數(shù)量宜為5個(gè),超過(guò)5個(gè)后總涌水量增加不明顯,增加的工程量為無(wú)效工程量。

    4 600 m工作面疏放水效果

    以鉆孔仰角60°,垂向長(zhǎng)度為120 m,單個(gè)鉆場(chǎng)鉆孔5個(gè)為組合,探討150 m鉆孔間距下600 m工作面范圍內(nèi)的疏放效果。此方案下鉆場(chǎng)共布置10個(gè)鉆場(chǎng)(5鉆場(chǎng)×2巷),50個(gè)鉆孔。假設(shè)3種定降深疏降模式,分別為780,830 m以及880 m。在工作面600 m范圍內(nèi)進(jìn)行井下定水頭疏放水,分別探討其涌水量、降深以及流場(chǎng)分布情況。

    4.1 鉆場(chǎng)涌水量

    根據(jù)模型計(jì)算結(jié)果,不同降深的鉆孔總涌水量變化如圖9所示。由圖9可知,3種水頭疏降方案下,鉆孔總涌水量隨時(shí)間逐步減小,最終趨于穩(wěn)定。水頭越小,鉆孔的初始總涌水量越大,但隨著時(shí)間的不斷延續(xù),最終3種水頭下的鉆孔總涌水量相差變小。

    圖9 不同降深的鉆孔總涌水量歷時(shí)曲線(xiàn)Fig.9 Curves of total water inflow of boreholes with different drawdowns

    60 d內(nèi),780 m定水頭疏降時(shí),鉆孔總涌水量由第1天的1 700 m3/h降至170 m3/h,降幅達(dá)90%;830和880 m定水頭疏降下,鉆孔總涌水量由第1天的1 530,1 350 m3/h分別降至153,135 m3/h,降幅均為90%。

    不同鉆場(chǎng)內(nèi)的鉆孔總涌水量分布圖如圖10所示。由圖10可知,600 m范圍內(nèi)工作面兩端2個(gè)鉆場(chǎng)(1,5和6,10)鉆場(chǎng)總涌水量較大,內(nèi)側(cè)3個(gè)鉆場(chǎng)總涌水量則相對(duì)較小,位于中間第3和第8鉆場(chǎng)鉆孔總涌水量最小。780 m定水頭疏降時(shí),兩端鉆孔總涌水量約為309 m3/h,中間鉆場(chǎng)鉆孔總涌水量約為142 m3/h。由此可見(jiàn),鉆場(chǎng)內(nèi)鉆孔總涌水量受周邊鉆場(chǎng)疏降的影響,由兩端向內(nèi)側(cè)影響逐步增大。其他2個(gè)定降深疏降條件下,鉆場(chǎng)總涌水量的變化趨勢(shì)基本一致。

    圖10 第1天不同鉆場(chǎng)內(nèi)鉆孔總涌水量Fig.10 Total water inflow of different drilling field in the first day

    4.2 降 深

    不同定水頭疏降方案下,工作面內(nèi)部水頭變化情況如圖11所示。由圖11可知,隨著疏放時(shí)間的延續(xù),不同方案下,工作面內(nèi)部水頭變化趨勢(shì)均為由初始水頭1 267.702 m降至不同的定水頭(780,830和880 m),定水頭越小,疏降的時(shí)間則相對(duì)越長(zhǎng)。150 m鉆場(chǎng)間距方案下,疏降至780,830和880 m的時(shí)間分別為155,125和100 d。

    圖11 不同降深下工作面內(nèi)部水頭變化歷時(shí)曲線(xiàn)Fig.11 Duration curves of water head change in working face under different drawdown

    不同時(shí)間不同降深的范圍有所不同如圖12所示。由圖12可知,隨著時(shí)間延續(xù),不同降深范圍逐步增大(橫向)。以200 m的降深為例,10 d之前,降深均小于200 m,10 d之后,隨著時(shí)間延續(xù),降深大于200 m的范圍逐漸增大。同一時(shí)間內(nèi),小降深的范圍大于大降深的范圍(縱向)。

    圖12 定水頭880 m方案下不同時(shí)間的降深范圍曲線(xiàn)Fig.12 Curves of drawdown range at different time under 880 m constant head scheme

    4.3 流 場(chǎng)

    隨著疏放水工程推進(jìn),煤層頂板含水層的流場(chǎng)不斷發(fā)生變化。剖面和平面流場(chǎng)圖如圖13,14所示。由圖13,14可知,隨著疏放時(shí)間延續(xù),橫剖面呈現(xiàn)以?xún)A斜鉆孔為中心的降落漏斗且范圍逐步增大,平面上則呈現(xiàn)以全部鉆場(chǎng)范圍為中心的降落漏斗且范圍逐漸增大。

    圖13 定水頭880 m方案下橫剖面流場(chǎng)(中心)Fig.13 Flow field of cross section under 880 m constant head scheme (central location)

    圖14 定水頭880 m方案下侏羅系含水層平面流場(chǎng)Fig.14 Plane flow field of Jurassic aquifer under 880 m constant head scheme

    5 結(jié) 論

    (1)針對(duì)復(fù)合充水含水層傾斜疏放水鉆孔涌水量計(jì)算方法不完善的問(wèn)題,以能量守恒定律為基礎(chǔ),以“含水層-鉆孔水量交換”為耦合點(diǎn),建立了含水層-鉆孔系統(tǒng)多流態(tài)水流計(jì)算的耦合模型。

    (2)應(yīng)用母杜柴登礦井臨時(shí)煤倉(cāng)井下放水試驗(yàn)的1個(gè)主放水孔水量和2個(gè)觀(guān)測(cè)孔水位歷時(shí)數(shù)據(jù),驗(yàn)證了含水層-鉆孔系統(tǒng)耦合模型的可靠性,可以很好解決傾斜鉆孔在復(fù)合充水含水層中的水量分配問(wèn)題。

    (3)傾斜疏放水鉆孔涌水量受到其角度和數(shù)量等參數(shù)的影響。隨著角度增大,單位長(zhǎng)度鉆孔涌水量呈現(xiàn)先增加后減小的趨勢(shì),在仰角為60°時(shí),鉆孔單位涌水量最大。此外,隨著單個(gè)鉆場(chǎng)內(nèi)鉆孔數(shù)量的增多,鉆孔總涌水量不斷增加,鉆孔數(shù)量超過(guò)5個(gè)后總涌水量增加并不明顯。即鉆孔的最優(yōu)角度為60°,單個(gè)鉆場(chǎng)內(nèi)鉆孔最佳數(shù)量為5個(gè)。

    (4)150 m間距方案下600 m工作面范圍內(nèi),不同疏降水頭鉆孔總涌水量、工作面內(nèi)部降深、流場(chǎng)變化均不同。降深越大,鉆孔的初始總涌水量越大,隨著時(shí)間推移,最終鉆孔總涌水量相差變小,直至穩(wěn)定。鉆場(chǎng)內(nèi)鉆孔總涌水量受到周邊鉆場(chǎng)疏降的影響,由兩端向內(nèi)側(cè)影響逐步增大。同時(shí),疏降定水頭越小,疏降時(shí)間則相對(duì)延長(zhǎng)。

    猜你喜歡
    鉆場(chǎng)仰角涌水量
    鉆場(chǎng)大斷面錨索加強(qiáng)支護(hù)效果分析
    煤(2020年6期)2020-07-03 06:45:22
    用銳角三角函數(shù)解決仰角、俯角問(wèn)題
    采煤工作面過(guò)高位鉆場(chǎng)期間瓦斯治理技術(shù)
    小浪底引黃工程干線(xiàn)1#隧洞涌水量預(yù)測(cè)
    分段三次Hermite插值計(jì)算GNSS系統(tǒng)衛(wèi)星仰角
    科技視界(2016年15期)2016-06-30 12:48:10
    黃石互通淺埋隧道涌水量預(yù)測(cè)分析
    模糊數(shù)學(xué)在預(yù)測(cè)大孤山鐵礦深部涌水量中的應(yīng)用
    傾斜鉆孔與高位巷抽采裂隙帶的瓦斯效果試驗(yàn)研究
    山西煤炭(2013年2期)2013-10-22 06:54:34
    三維離散元法在巷道涌水量預(yù)測(cè)中的應(yīng)用
    金屬礦山(2013年12期)2013-03-11 16:55:30
    注漿充填法過(guò)頂板鉆場(chǎng)開(kāi)采實(shí)踐
    日韩中字成人| 亚洲18禁久久av| 国产亚洲欧美在线一区二区| 给我免费播放毛片高清在线观看| 国产精品久久久久久久久免 | 国产精品98久久久久久宅男小说| 久久国产乱子免费精品| 乱人视频在线观看| 亚洲成a人片在线一区二区| 日本一本二区三区精品| 99久久九九国产精品国产免费| 自拍偷自拍亚洲精品老妇| 麻豆国产97在线/欧美| 亚洲欧美清纯卡通| 夜夜躁狠狠躁天天躁| av天堂中文字幕网| 美女黄网站色视频| 国产亚洲av嫩草精品影院| 欧美另类亚洲清纯唯美| 看十八女毛片水多多多| 少妇被粗大猛烈的视频| 91久久精品国产一区二区成人| 色视频www国产| 免费一级毛片在线播放高清视频| 大型黄色视频在线免费观看| 午夜福利免费观看在线| 最新在线观看一区二区三区| 国产伦一二天堂av在线观看| 日韩中字成人| 国产精品久久视频播放| 午夜久久久久精精品| 国产乱人伦免费视频| 亚洲精品成人久久久久久| 美女xxoo啪啪120秒动态图 | 久久精品国产清高在天天线| 国产精华一区二区三区| 国产精品一区二区免费欧美| 国产免费一级a男人的天堂| 日本五十路高清| 琪琪午夜伦伦电影理论片6080| 成年人黄色毛片网站| 欧美高清性xxxxhd video| 久久久久久久久久成人| 国内揄拍国产精品人妻在线| 一区二区三区四区激情视频 | 国产高清激情床上av| 亚洲欧美日韩卡通动漫| 在线播放国产精品三级| 999久久久精品免费观看国产| 国产91精品成人一区二区三区| 色精品久久人妻99蜜桃| 别揉我奶头 嗯啊视频| 国内精品一区二区在线观看| 久久这里只有精品中国| 久久九九热精品免费| 国产在视频线在精品| 在线观看66精品国产| 亚洲美女视频黄频| 一本久久中文字幕| 亚洲人成网站高清观看| 国内少妇人妻偷人精品xxx网站| 天堂动漫精品| 少妇人妻精品综合一区二区 | 色吧在线观看| 一级作爱视频免费观看| 久久久久久久久中文| 亚洲av五月六月丁香网| 女同久久另类99精品国产91| bbb黄色大片| 啪啪无遮挡十八禁网站| 亚洲欧美清纯卡通| 久久久国产成人精品二区| 婷婷丁香在线五月| 可以在线观看的亚洲视频| 精品人妻熟女av久视频| 一本一本综合久久| 成人午夜高清在线视频| 久久亚洲真实| 国产真实乱freesex| 久久久久国内视频| 亚洲七黄色美女视频| 琪琪午夜伦伦电影理论片6080| 日本免费a在线| 制服丝袜大香蕉在线| 中亚洲国语对白在线视频| 午夜精品在线福利| 一卡2卡三卡四卡精品乱码亚洲| 欧美最黄视频在线播放免费| 国产精品人妻久久久久久| 欧美潮喷喷水| 国产成人欧美在线观看| 久久99热6这里只有精品| 欧美色视频一区免费| 嫩草影院精品99| 在线观看美女被高潮喷水网站 | 青草久久国产| 18禁黄网站禁片午夜丰满| 99久国产av精品| 欧美最黄视频在线播放免费| 99久久九九国产精品国产免费| 欧美区成人在线视频| 伦理电影大哥的女人| 91九色精品人成在线观看| 99热精品在线国产| 亚洲三级黄色毛片| 又爽又黄a免费视频| 波多野结衣巨乳人妻| 国产乱人伦免费视频| 国产乱人视频| 日韩欧美精品v在线| 久久久久久久久大av| 嫩草影院新地址| 真人一进一出gif抽搐免费| 在线免费观看不下载黄p国产 | 免费av不卡在线播放| 亚洲无线观看免费| 国产成+人综合+亚洲专区| 国产精品乱码一区二三区的特点| 在线a可以看的网站| 此物有八面人人有两片| 亚洲激情在线av| 亚洲欧美日韩卡通动漫| 99久久精品国产亚洲精品| 久久精品国产亚洲av天美| 精品久久久久久成人av| 日本黄色片子视频| 欧美中文日本在线观看视频| 最近最新中文字幕大全电影3| 国产精品永久免费网站| 亚洲不卡免费看| 美女被艹到高潮喷水动态| 少妇高潮的动态图| 悠悠久久av| 五月伊人婷婷丁香| 九九热线精品视视频播放| 久久久久九九精品影院| 久久久色成人| 亚洲性夜色夜夜综合| 国产精品av视频在线免费观看| 亚洲av.av天堂| 黄色女人牲交| 最新中文字幕久久久久| 精品午夜福利在线看| а√天堂www在线а√下载| 国产在线精品亚洲第一网站| 亚洲无线观看免费| 国产精品一区二区性色av| 香蕉av资源在线| 日本免费a在线| 免费在线观看亚洲国产| 天堂网av新在线| 国产毛片a区久久久久| av在线蜜桃| 国产精品免费一区二区三区在线| 最近中文字幕高清免费大全6 | 国产三级中文精品| 欧美+亚洲+日韩+国产| 国产爱豆传媒在线观看| 亚洲va日本ⅴa欧美va伊人久久| 国产成+人综合+亚洲专区| 色综合婷婷激情| 人人妻人人看人人澡| 757午夜福利合集在线观看| 亚洲成人免费电影在线观看| 国产三级中文精品| 久久婷婷人人爽人人干人人爱| 男女之事视频高清在线观看| 日韩有码中文字幕| 中文字幕久久专区| 国产精品亚洲一级av第二区| 午夜福利欧美成人| 亚洲精品成人久久久久久| 午夜福利在线在线| 午夜a级毛片| 高清日韩中文字幕在线| 欧美黄色片欧美黄色片| 欧美一区二区国产精品久久精品| 欧美成人一区二区免费高清观看| 人妻久久中文字幕网| 99国产精品一区二区三区| 色在线成人网| 国产一区二区三区在线臀色熟女| 在线免费观看不下载黄p国产 | 免费人成在线观看视频色| 国产一级毛片七仙女欲春2| 十八禁网站免费在线| 亚洲欧美日韩高清专用| 午夜日韩欧美国产| 日本五十路高清| 两个人视频免费观看高清| 成人亚洲精品av一区二区| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av美国av| 欧美色视频一区免费| 午夜免费男女啪啪视频观看 | 免费人成在线观看视频色| 此物有八面人人有两片| 日韩欧美免费精品| 桃色一区二区三区在线观看| 国产熟女xx| 亚洲国产精品999在线| 深爱激情五月婷婷| 国产成人aa在线观看| 国产色婷婷99| 12—13女人毛片做爰片一| 日本精品一区二区三区蜜桃| 69人妻影院| 夜夜夜夜夜久久久久| 亚洲天堂国产精品一区在线| 色综合婷婷激情| 亚洲精品在线观看二区| 床上黄色一级片| 嫩草影院精品99| 天天躁日日操中文字幕| 国产真实乱freesex| 免费黄网站久久成人精品 | 精品人妻1区二区| 亚洲激情在线av| 少妇的逼好多水| 欧美成人一区二区免费高清观看| 成人亚洲精品av一区二区| 欧美不卡视频在线免费观看| 中文字幕人成人乱码亚洲影| 欧美激情在线99| 国产精品久久电影中文字幕| 伦理电影大哥的女人| 少妇被粗大猛烈的视频| 久久久成人免费电影| 一进一出好大好爽视频| 国产黄a三级三级三级人| 成人无遮挡网站| 亚洲av电影不卡..在线观看| 99热精品在线国产| 三级男女做爰猛烈吃奶摸视频| 久久精品国产自在天天线| 99久久无色码亚洲精品果冻| 亚洲精品乱码久久久v下载方式| 免费av不卡在线播放| 国产欧美日韩一区二区三| 成人无遮挡网站| 99热这里只有是精品50| 白带黄色成豆腐渣| 在线观看66精品国产| 亚洲中文字幕日韩| 亚洲精品在线观看二区| 国产一区二区在线av高清观看| АⅤ资源中文在线天堂| 久久午夜福利片| 成熟少妇高潮喷水视频| 又黄又爽又免费观看的视频| 男人舔女人下体高潮全视频| 最后的刺客免费高清国语| 国产熟女xx| 在线观看一区二区三区| 欧美xxxx黑人xx丫x性爽| 男人和女人高潮做爰伦理| 熟女人妻精品中文字幕| 国产精品亚洲av一区麻豆| 精品久久久久久久久久久久久| 在线国产一区二区在线| 精品国产亚洲在线| 一个人看的www免费观看视频| 久久热精品热| 黄色视频,在线免费观看| 久久草成人影院| 白带黄色成豆腐渣| 亚洲av美国av| 亚洲成人久久爱视频| 国产成+人综合+亚洲专区| av欧美777| 青草久久国产| 麻豆一二三区av精品| 亚洲一区二区三区色噜噜| 级片在线观看| 成人特级av手机在线观看| 亚洲精品久久国产高清桃花| 一本久久中文字幕| 婷婷丁香在线五月| 日韩精品中文字幕看吧| 久久婷婷人人爽人人干人人爱| 日本一二三区视频观看| 国产精品野战在线观看| 精品久久久久久,| 国产伦人伦偷精品视频| 成年女人毛片免费观看观看9| 男女那种视频在线观看| 在线国产一区二区在线| 亚洲 欧美 日韩 在线 免费| 午夜影院日韩av| 十八禁网站免费在线| 国产日本99.免费观看| 欧美国产日韩亚洲一区| 看十八女毛片水多多多| 一区福利在线观看| 久久九九热精品免费| 中文字幕av在线有码专区| 亚洲国产精品sss在线观看| 国产精品久久视频播放| 成人一区二区视频在线观看| 91av网一区二区| 一级黄片播放器| 国产精品一区二区性色av| 成人无遮挡网站| 国产精品国产高清国产av| 欧美性猛交╳xxx乱大交人| 国产精品电影一区二区三区| 亚洲内射少妇av| 十八禁国产超污无遮挡网站| 草草在线视频免费看| 成年女人毛片免费观看观看9| 亚洲成人久久性| 欧美日韩乱码在线| 中文字幕人妻熟人妻熟丝袜美| 白带黄色成豆腐渣| 男人舔奶头视频| 99在线人妻在线中文字幕| 欧美成人a在线观看| 欧美另类亚洲清纯唯美| 久久久久久久久中文| 成年女人看的毛片在线观看| 色综合亚洲欧美另类图片| 亚洲国产色片| 色5月婷婷丁香| 国产伦精品一区二区三区视频9| 给我免费播放毛片高清在线观看| 91狼人影院| 五月伊人婷婷丁香| 久久久久国内视频| 国产日本99.免费观看| 亚洲av电影不卡..在线观看| 91九色精品人成在线观看| 久久热精品热| 18禁黄网站禁片免费观看直播| 精品午夜福利在线看| 成人性生交大片免费视频hd| 特级一级黄色大片| 高清日韩中文字幕在线| 国产精品美女特级片免费视频播放器| 欧美精品啪啪一区二区三区| 黄色女人牲交| 国内精品久久久久久久电影| 精品不卡国产一区二区三区| 亚洲熟妇中文字幕五十中出| 亚洲av电影在线进入| 欧美成人a在线观看| 国产淫片久久久久久久久 | 人妻久久中文字幕网| www.熟女人妻精品国产| 韩国av一区二区三区四区| 国内揄拍国产精品人妻在线| 久久久色成人| 久久精品国产亚洲av涩爱 | 免费在线观看成人毛片| 久久欧美精品欧美久久欧美| 亚洲成a人片在线一区二区| 亚洲精品在线观看二区| 99热6这里只有精品| 欧美一区二区精品小视频在线| 美女cb高潮喷水在线观看| 国产免费男女视频| 国产精品一及| 日本一二三区视频观看| 成年免费大片在线观看| 啦啦啦韩国在线观看视频| 国产午夜福利久久久久久| www.999成人在线观看| 在线免费观看不下载黄p国产 | 欧美黄色淫秽网站| 好看av亚洲va欧美ⅴa在| 香蕉av资源在线| 国产亚洲欧美98| 成人特级黄色片久久久久久久| 国内毛片毛片毛片毛片毛片| 又黄又爽又刺激的免费视频.| av天堂在线播放| 在线观看午夜福利视频| 天堂影院成人在线观看| 白带黄色成豆腐渣| 午夜a级毛片| 在线观看午夜福利视频| 十八禁国产超污无遮挡网站| 国产高潮美女av| 午夜福利免费观看在线| 午夜激情福利司机影院| 久久久久国产精品人妻aⅴ院| 美女大奶头视频| 青草久久国产| 成人国产综合亚洲| 赤兔流量卡办理| 97碰自拍视频| 一区二区三区免费毛片| netflix在线观看网站| 国产高清激情床上av| 在现免费观看毛片| 国产色爽女视频免费观看| 日韩欧美在线乱码| 美女免费视频网站| 18美女黄网站色大片免费观看| 免费人成在线观看视频色| 少妇丰满av| 国产淫片久久久久久久久 | 91在线观看av| 18禁黄网站禁片免费观看直播| 嫁个100分男人电影在线观看| 赤兔流量卡办理| 亚洲国产精品sss在线观看| 能在线免费观看的黄片| 国产毛片a区久久久久| 波多野结衣高清作品| 国内精品久久久久久久电影| 精品国内亚洲2022精品成人| 91av网一区二区| 亚洲五月婷婷丁香| 国产亚洲精品av在线| 国产高清视频在线观看网站| 人妻久久中文字幕网| 国产精品久久久久久亚洲av鲁大| a在线观看视频网站| 特级一级黄色大片| 琪琪午夜伦伦电影理论片6080| 亚洲国产精品久久男人天堂| 永久网站在线| 18禁在线播放成人免费| 真实男女啪啪啪动态图| 2021天堂中文幕一二区在线观| 欧美黑人巨大hd| 直男gayav资源| 神马国产精品三级电影在线观看| 久久6这里有精品| 真人做人爱边吃奶动态| 国产男靠女视频免费网站| 久久性视频一级片| 久久久久免费精品人妻一区二区| 日本精品一区二区三区蜜桃| 国产成人啪精品午夜网站| 精品人妻1区二区| 国产淫片久久久久久久久 | 国产不卡一卡二| 麻豆av噜噜一区二区三区| 久久人人精品亚洲av| 色av中文字幕| 丰满人妻一区二区三区视频av| av在线老鸭窝| 精品一区二区免费观看| 极品教师在线视频| 久久精品91蜜桃| 一夜夜www| 天堂网av新在线| 深夜a级毛片| 99在线视频只有这里精品首页| 欧美精品啪啪一区二区三区| 亚洲熟妇中文字幕五十中出| 欧美黄色淫秽网站| 国产亚洲精品久久久久久毛片| 欧美性猛交黑人性爽| 国产高清三级在线| 精品乱码久久久久久99久播| 精品国内亚洲2022精品成人| 国内毛片毛片毛片毛片毛片| 免费无遮挡裸体视频| 亚洲av熟女| 国产伦在线观看视频一区| 国产黄a三级三级三级人| 欧洲精品卡2卡3卡4卡5卡区| 国产精品久久久久久久电影| 久久婷婷人人爽人人干人人爱| 又黄又爽又刺激的免费视频.| 久久精品综合一区二区三区| 国产精品乱码一区二三区的特点| 十八禁人妻一区二区| 国产高清视频在线播放一区| 狂野欧美白嫩少妇大欣赏| 少妇人妻一区二区三区视频| 在线播放国产精品三级| 精品久久久久久久久av| 国产伦一二天堂av在线观看| 欧美成人一区二区免费高清观看| 一进一出抽搐gif免费好疼| 免费在线观看亚洲国产| 国产主播在线观看一区二区| 高潮久久久久久久久久久不卡| a级毛片a级免费在线| 在线十欧美十亚洲十日本专区| 久久午夜福利片| 国产高清视频在线播放一区| 久久草成人影院| 午夜a级毛片| 国产成人啪精品午夜网站| av欧美777| 欧美日韩亚洲国产一区二区在线观看| 免费在线观看影片大全网站| 最后的刺客免费高清国语| 少妇高潮的动态图| 黄色视频,在线免费观看| 久久久久久久亚洲中文字幕 | www.www免费av| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品爽爽va在线观看网站| 精品一区二区三区av网在线观看| 3wmmmm亚洲av在线观看| 首页视频小说图片口味搜索| 国产爱豆传媒在线观看| 日韩亚洲欧美综合| 波多野结衣高清无吗| 蜜桃久久精品国产亚洲av| 免费高清视频大片| 欧美+日韩+精品| 十八禁人妻一区二区| 国产三级中文精品| 日本三级黄在线观看| 亚洲精品粉嫩美女一区| 亚洲国产精品合色在线| 亚洲最大成人手机在线| 亚洲精品一区av在线观看| 亚洲成a人片在线一区二区| 日本与韩国留学比较| 很黄的视频免费| 大型黄色视频在线免费观看| 精品欧美国产一区二区三| 全区人妻精品视频| av视频在线观看入口| 亚洲,欧美,日韩| 黄色视频,在线免费观看| 丁香六月欧美| 十八禁网站免费在线| 国产在视频线在精品| 午夜a级毛片| 内地一区二区视频在线| av国产免费在线观看| 免费看日本二区| www日本黄色视频网| 国产成人啪精品午夜网站| 色哟哟·www| 亚洲av五月六月丁香网| 18禁黄网站禁片午夜丰满| 精品国产三级普通话版| 少妇的逼水好多| 日本免费a在线| 天天一区二区日本电影三级| 欧美zozozo另类| 俺也久久电影网| .国产精品久久| 18美女黄网站色大片免费观看| 变态另类成人亚洲欧美熟女| 性插视频无遮挡在线免费观看| 国产精品1区2区在线观看.| 一区二区三区免费毛片| 午夜福利高清视频| 久久伊人香网站| 国产一区二区三区在线臀色熟女| 一本综合久久免费| 久久久久久久午夜电影| 午夜免费激情av| 亚洲精品一卡2卡三卡4卡5卡| 国模一区二区三区四区视频| 一区二区三区激情视频| 性色avwww在线观看| 女同久久另类99精品国产91| 成人无遮挡网站| 啪啪无遮挡十八禁网站| 欧美绝顶高潮抽搐喷水| 欧美zozozo另类| 网址你懂的国产日韩在线| 最近最新免费中文字幕在线| 在线观看66精品国产| 国产精华一区二区三区| 中文亚洲av片在线观看爽| 免费观看的影片在线观看| 国产极品精品免费视频能看的| 一本一本综合久久| 12—13女人毛片做爰片一| 18禁黄网站禁片午夜丰满| av在线观看视频网站免费| 蜜桃亚洲精品一区二区三区| 久久久久久久久久黄片| 亚洲最大成人中文| 性插视频无遮挡在线免费观看| 在线天堂最新版资源| 欧美成狂野欧美在线观看| 成人毛片a级毛片在线播放| 欧美xxxx性猛交bbbb| av在线观看视频网站免费| 老司机午夜福利在线观看视频| 亚洲精品日韩av片在线观看| 91麻豆av在线| 国产蜜桃级精品一区二区三区| 免费在线观看日本一区| 色播亚洲综合网| 嫩草影视91久久| 国产色爽女视频免费观看| av视频在线观看入口| 亚洲精品亚洲一区二区| 国产在视频线在精品| 男女之事视频高清在线观看| 国内少妇人妻偷人精品xxx网站| 99久久成人亚洲精品观看| 少妇裸体淫交视频免费看高清| 亚洲一区二区三区色噜噜| 国产视频内射| 国内揄拍国产精品人妻在线| 99热这里只有是精品50| 变态另类丝袜制服| 久久久久久久久中文| 午夜亚洲福利在线播放| 婷婷亚洲欧美| a在线观看视频网站| 欧美日韩黄片免| 亚洲成人久久爱视频| 成人国产一区最新在线观看| 首页视频小说图片口味搜索| 美女大奶头视频| 尤物成人国产欧美一区二区三区| 亚洲国产欧洲综合997久久,| 9191精品国产免费久久| 成人午夜高清在线视频|