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

    天橋泉域某礦奧灰含水層水文地質(zhì)特征研究

    2012-09-16 07:21:08李國(guó)紅
    焦作大學(xué)學(xué)報(bào) 2012年2期
    關(guān)鍵詞:水化學(xué)井田灰?guī)r

    李國(guó)紅

    (河南省煤田地質(zhì)局二隊(duì),河南洛陽(yáng)471023)

    天橋泉域某礦奧灰含水層水文地質(zhì)特征研究

    李國(guó)紅

    (河南省煤田地質(zhì)局二隊(duì),河南洛陽(yáng)471023)

    奧陶紀(jì)灰?guī)r(奧灰)巖溶水是華北型煤田下組煤安全開(kāi)采的主要威脅,查明奧灰水文地質(zhì)條件是下組煤開(kāi)采的關(guān)鍵和前提。文章通過(guò)水化學(xué)分析,獲取了研究區(qū)奧灰?guī)r溶水水化學(xué)參數(shù),探討了井田奧灰?guī)r溶水水動(dòng)力條件及其演化規(guī)律。通過(guò)群孔抽水試驗(yàn),研究了奧灰含水層水力聯(lián)系特征,并驗(yàn)證了水化學(xué)分析結(jié)論。同時(shí)運(yùn)用GMS數(shù)值模擬軟件進(jìn)行了群孔抽水試驗(yàn)含水層流場(chǎng)特征模擬,反演出滲透系數(shù)及貯水系數(shù)等水文地質(zhì)參數(shù)。研究成果為該礦下組煤開(kāi)采奧灰?guī)r溶水防治提供了可靠依據(jù)。

    奧灰含水層;水文地質(zhì)特征

    1. 引言

    華北型煤田主要水害類(lèi)型為底板突水,充水含水層是奧灰?guī)r溶含水層及煤系內(nèi)薄層灰?guī)r含水層。目前華北型礦井平均采深已超過(guò)650m,并以12~20m/a的速度向深部拓展[1],逐步開(kāi)發(fā)下組煤。與淺部開(kāi)采相比,下組煤深受奧灰含水層威脅,且之間隔水層厚度小,隔水層結(jié)構(gòu)受到多種構(gòu)造因素的破壞而阻水能力嚴(yán)重不足[2],礦井水害隱患大大增加。只有充分掌握和查明奧灰含水層水文地質(zhì)特征,才能進(jìn)行正確的下組煤開(kāi)采安全性評(píng)價(jià),保證煤炭資源的有效回收和采掘工人的人身安全。因此,奧灰含水層水文地質(zhì)特征研究是解放下組煤煤炭資源的關(guān)鍵所在[3]。

    山西省某礦為設(shè)計(jì)生產(chǎn)能力1000余萬(wàn)噸的特大型礦井。井田地處河?xùn)|煤田北部,構(gòu)造簡(jiǎn)單,呈平緩的單斜構(gòu)造形態(tài),稀疏小斷層發(fā)育。目前主采煤層為上組煤8煤,接續(xù)煤層為下組煤13煤,直接突水威脅來(lái)自底板奧灰含水層。為保障 13煤安全開(kāi)采,本文對(duì)井田奧灰含水層水文地質(zhì)特征進(jìn)行了研究,為下組煤開(kāi)采提供技術(shù)資料。

    2. 井田奧灰含水層概況

    2.1 奧灰含水層層段劃分

    井田內(nèi)奧陶系自上而下劃分為峰峰組和馬家溝組:

    (1)峰峰組灰?guī)r巖溶含水層:根據(jù)鉆孔資料,該含水層厚度為89.50~125.75m,平均104.12m;巖性以泥灰?guī)r、白云質(zhì)泥灰?guī)r及角礫狀泥灰?guī)r為主,上部夾有幾層薄層泥巖,下部賦存石膏層。

    (2)馬家溝組灰?guī)r巖溶含水層:鉆孔僅揭露了上馬家溝組,揭露厚度33.87~108.62m。主要巖性為泥灰?guī)r、白云質(zhì)泥灰?guī)r、角礫狀泥灰?guī)r和豹皮狀灰?guī)r。本段含水層巖溶發(fā)育,以溶隙、溶孔和小型溶洞為主,小溶洞多被充填。

    2.2 垂向富水性變化規(guī)律

    3. 井田奧灰含水層水化學(xué)特征分析

    3.1 區(qū)域巖溶水補(bǔ)逕排條件

    研究區(qū)處于天橋泉域西側(cè),泉域巖溶水系統(tǒng)為一補(bǔ)給、徑流、排泄完整、獨(dú)立的全排型巖溶水水文地質(zhì)單元[4],由天橋泉、龍口泉及老牛灣泉三個(gè)巖溶水子系統(tǒng)組成,子系統(tǒng)間以地下水分水嶺相間。泉域接受東部裸露巖溶區(qū)大氣降水入滲補(bǔ)給、黃河及其支流的滲漏補(bǔ)給[5],總體向西部中段匯流,最終在黃河沿岸以溢流泉形式排泄。

    3.2 奧灰含水層水化學(xué)三線圖分析

    選取井田及外圍12個(gè)鉆孔(供水井)獲取12個(gè)巖溶水水樣,并進(jìn)行了水化學(xué)成分分析,得到了該區(qū)奧灰?guī)r溶水主要化學(xué)特征值,并依據(jù)舒卡列夫分類(lèi)方法確定了各組水樣的水質(zhì)類(lèi)型,成果如表1所示。

    表1 奧灰?guī)r溶水水化學(xué)特征值

    將奧灰含水層巖溶水所有水樣的陰陽(yáng)離子毫克當(dāng)量百分?jǐn)?shù)投影到 Piper三線圖上,得到如圖2所示的水化學(xué)三線圖。Piper三線圖可以填入大量的分析結(jié)果,從而查明所分析水中陰陽(yáng)離子的對(duì)比關(guān)系,這種表示簡(jiǎn)捷有效[6]。

    圖2 奧灰含水層水化學(xué)三線圖

    根據(jù)Piper三線圖圖解分區(qū),菱形區(qū)共劃分為 9個(gè)分區(qū),落在不同分區(qū)的水樣具有不同的化學(xué)特征。分析圖2可知,奧灰水樣分布分散,主要分布區(qū)為 5區(qū)、6區(qū)、7區(qū)及9區(qū)。5區(qū)表示奧灰水以 Ca2++Mg2+離子和為主,酸性較弱,硬度較小;6區(qū)表示奧灰水中陰離子以為主,Na+(K+)逐漸增多,硬度逐漸增大;7區(qū)表示奧灰水中陰離子以 Cl-為主,陽(yáng)離子以 Na+(K+)為主,硬度大;9區(qū)表示區(qū)內(nèi)水樣為6區(qū)與5區(qū)及7區(qū)的過(guò)渡區(qū),奧灰水中含有的陰陽(yáng)離子多樣,可視為不同演化階段巖溶水的過(guò)渡地下水。

    3.3 井田奧灰含水層巖溶水演化特征分析

    根據(jù)3.2分析,所取的12個(gè)水樣水化學(xué)成分差別較大,說(shuō)明水樣選取點(diǎn)奧灰?guī)r溶水徑流條件及演化階段的不同。研究表明,奧灰?guī)r溶水隨徑流條件變化,其化學(xué)成分的演化主要受以下三項(xiàng)化學(xué)溶解反應(yīng)影響:

    受以上三個(gè)平衡化學(xué)式影響,奧灰?guī)r溶水補(bǔ)給區(qū)CO2供給充足可以溶解方解石、石膏和白云石,沿巖溶水徑流路徑,Mg2+、Ca2+、SO2-4及 HCO-3離子含量均呈增加趨勢(shì),礦化度也隨之增加[7];隨著巖溶水徑流條件變差,CO2逐漸消耗,則方解石趨于飽和、不再溶解,奧灰水中HCO-3離子含量逐漸穩(wěn)定,而石膏和白云石仍在溶解,SO42-離子含量仍在增加;巖溶水徑流條件進(jìn)一步變差,白云石趨于飽和、停止溶解,Mg2+離子含量逐漸趨于穩(wěn)定,當(dāng)石膏也趨于飽和、停止溶解后,SO2-4離子含量也逐漸穩(wěn)定,此時(shí)巖溶水處于滯流狀態(tài)。隨著巖溶水徑流路徑的延長(zhǎng),濃縮作用不斷增強(qiáng),加之煤系中淋濾作用,巖溶水中 Na+(K+)、Cl-離子含量隨徑流條件的改變不斷增加。

    (1)靜脈炎。與導(dǎo)管相關(guān)性靜脈炎性癥狀或體征:持續(xù)疼痛(靜脈輸液結(jié)束后持續(xù)疼痛大于2小時(shí));發(fā)紅;腫脹;硬化(條索狀);膿性分泌物。分析采用美國(guó)靜脈輸液護(hù)理學(xué)會(huì)(Infusion Nursing Society,INS)分級(jí)系統(tǒng) [6]和 Tager’s[7]評(píng)分綜合考慮,分為三級(jí)。

    圖3 井田奧灰?guī)r溶水徑流路徑及演化過(guò)程

    根據(jù)以上分析,可將取樣按照水化學(xué)特征劃分為以下5組(如表1所示):

    (1)Ⅰ組:包括2#、4#及5#供水井,該組奧灰水水質(zhì)類(lèi)型為Na(Na·Ca)型,反映了徑流區(qū)水質(zhì)特征,處于井田以東區(qū)域;

    (2)Ⅱ組:包括1#及6#供水井,該組奧灰水水質(zhì)類(lèi)型為HCO3·Cl-Na型,離子比例減少,Cl-離子比例明顯增加,反映了徑流區(qū)向弱徑流區(qū)轉(zhuǎn)化的水質(zhì)特征,處于井田北部;

    (3)Ⅲ組:包括CK6孔及3#供水井,該組奧灰水水質(zhì)類(lèi)型為HCO3·Cl-Na·Ca(Ca·Mg)型,巖溶水徑流條件減弱,反映了弱徑流區(qū)水質(zhì)特征,處于井田東部邊界;

    (4)Ⅳ組:包括CK3孔,奧灰水水質(zhì)為Cl·SO4-Ca型離子為主要陰離子,巖溶水徑流條件進(jìn)一步減弱,反映了弱徑流區(qū)向滯流區(qū)轉(zhuǎn)化的水質(zhì)特征,處于井田東南部;

    (5)Ⅴ組:包括CK2、CK4、CK5及CK1孔,該組奧灰水水質(zhì)為Cl-Na,反映出滯流奧灰?guī)r溶水特點(diǎn),處于井田中部。

    各組水化學(xué)特征變化反映了井田內(nèi)不同區(qū)段奧灰?guī)r溶水水動(dòng)力條件及演化路徑,即井田內(nèi)巖溶水總體上由東北向西南徑流,依次可劃分為弱徑流區(qū)、弱徑流 -滯流過(guò)渡區(qū)及滯流區(qū)。隨著巖溶水演化過(guò)程的延續(xù),地下水徑流路徑不斷延長(zhǎng),礦化度逐漸升高,水質(zhì)類(lèi)型總體呈HCO3-Na(Na·Ca)→HCO3·Cl-Na·Ca(Ca· Mg)→Cl-Na的演化順序。此規(guī)律與天橋泉域內(nèi)水化學(xué)特征基本一致[8、9]。

    4. 奧灰含水層群孔抽水試驗(yàn)

    群孔抽水試驗(yàn)是研究含水層地下水連通性、流場(chǎng)特征及水動(dòng)力條件的有效手段[10]。為檢驗(yàn)水化學(xué)分析結(jié)論,并進(jìn)一步查明奧灰含水層流場(chǎng)特征,獲取水文地質(zhì)參數(shù),本礦首采區(qū)進(jìn)行群孔抽水試驗(yàn)。

    4.1 抽水孔及觀測(cè)孔布置

    受井田地面及井下條件所限,不能施工新奧灰孔,只能以現(xiàn)有鉆孔為基礎(chǔ)進(jìn)行群孔抽水試驗(yàn)。抽水孔及觀測(cè)孔位置如圖 5所示。在首采區(qū)中部選取抽水孔 3個(gè),分別為CK1、CK2及CK4,觀測(cè)孔5個(gè),位于抽水孔外圍,分別為CK8、CK21、CK10、CK22和CK3。

    4.2 抽水試驗(yàn)結(jié)果及分析

    抽水試驗(yàn)從2009年 12月 4日 12時(shí)開(kāi)始,同時(shí)打開(kāi)CK 1、CK 2及CK 4抽水孔,水量控制閥開(kāi)至最大,同步進(jìn)行觀測(cè)孔水位觀測(cè)。13日8時(shí)關(guān)閉抽水孔閥門(mén),并開(kāi)始觀測(cè)各觀測(cè)孔水位恢復(fù),直至16日11時(shí)各孔奧灰水位恢復(fù)至抽水前水平。抽水試驗(yàn)結(jié)束,歷時(shí)17220min,其中抽水歷時(shí)12720min,恢復(fù)水位觀測(cè)歷時(shí)4500min。整個(gè)抽水試驗(yàn)過(guò)程中抽水孔水量變化不大,在抽水初期稍有波動(dòng),后期基本穩(wěn)定,抽水量及其變化見(jiàn)圖4。試驗(yàn)過(guò)程中,各觀測(cè)孔水位變化都不大(如圖6所示),說(shuō)明群孔抽水試驗(yàn)過(guò)程中,地下水位降落漏斗擴(kuò)散范圍普遍較小,橫向上奧灰含水層水力聯(lián)系弱。

    圖4 抽水孔水量歷時(shí)變化曲線圖

    5. 奧灰含水層地下水流場(chǎng)數(shù)值模擬

    GMS軟件是國(guó)際上最受歡迎的地下水模擬軟件[11],在地下水流場(chǎng)模擬及水文地質(zhì)參數(shù)反演方面等領(lǐng)域具有較大的應(yīng)用前景。本文基于群孔抽水試驗(yàn)結(jié)果,利用該軟件進(jìn)行水文地質(zhì)參數(shù)反演,研究奧灰含水層抽水試驗(yàn)過(guò)程中流場(chǎng)變化特征。

    5.1 水文地質(zhì)概念模型

    研究區(qū)域?yàn)榫锛安糠滞鈹U(kuò)區(qū)域,由于區(qū)域內(nèi)無(wú)較大的裂隙構(gòu)造,邊界均為人為邊界,即由井田邊界外擴(kuò)1800m。根據(jù)地質(zhì)資料和抽水試驗(yàn)資料,含水層內(nèi)部結(jié)構(gòu)概化為非均質(zhì)各向異性的非穩(wěn)定承壓二維流;邊界定為一類(lèi)補(bǔ)給邊界;補(bǔ)給項(xiàng)為大氣降水,排泄項(xiàng)為抽水孔排泄抽水。

    5.2 空間及時(shí)間離散

    全區(qū)離散采用矩形網(wǎng)格剖分,劃分網(wǎng)格時(shí)在抽水孔及觀測(cè)孔結(jié)點(diǎn)處適當(dāng)細(xì)化,共計(jì)剖分單元 5512個(gè),有效單元5082個(gè)(圖5)。模擬時(shí)段采用整個(gè)抽水試驗(yàn)過(guò)程,本次模擬時(shí)間離散時(shí)段劃分較細(xì),共69個(gè)應(yīng)力期,離散結(jié)果列于表2中。

    圖5 研究區(qū)網(wǎng)格剖分圖

    表2 模擬時(shí)間離散結(jié)果

    5.3 水文地質(zhì)模型識(shí)別與檢驗(yàn)

    由于研究模型邊界均為人為劃定,模型中存在較多不確定量,給參數(shù)識(shí)別帶來(lái)較大困難。為克服這個(gè)難題,本文采用試估法進(jìn)行模型識(shí)別,識(shí)別內(nèi)容包括滲透系數(shù)(Kx、Ky)和貯水系數(shù)μ*。首先進(jìn)行水文地質(zhì)參數(shù)預(yù)估算,即采用配線法利用 CK4孔的抽水量及 CK21孔和CK10孔的水位觀測(cè)資料對(duì)含水層導(dǎo)水系數(shù)、滲透系數(shù)及貯水系數(shù)等進(jìn)行初步估算,并作為數(shù)值模擬求參的初始值。其次將觀測(cè)孔的水位動(dòng)態(tài)實(shí)測(cè)資料作為擬合對(duì)象,以此確定含水層的水文地質(zhì)參數(shù),觀測(cè)孔水頭的擬合程度基本上可以反映所求參數(shù)的精度。設(shè)目標(biāo)函數(shù):

    式中,p1,p2,pn——待求的水文地質(zhì)參數(shù);

    m——用于水頭對(duì)比的觀測(cè)孔個(gè)數(shù);

    n——水頭對(duì)比時(shí)段數(shù);

    H(i,j)——模型計(jì)算的 i點(diǎn)觀測(cè)孔 j時(shí)段的水位(m);

    Hg(i,j)——i點(diǎn)觀測(cè)孔j時(shí)段的實(shí)測(cè)水位(m)。

    根據(jù)各參數(shù)的約束條件,求取目標(biāo)函數(shù) E最小。若此時(shí)的水位擬合差未達(dá)到足夠小,則進(jìn)一步分析其原因,必要時(shí)對(duì)約束條件及模型定解條件進(jìn)行必要的調(diào)整,再次求取目標(biāo)函數(shù)最小。如此反復(fù)計(jì)算、分析,直到求得滿意的結(jié)果為止。

    識(shí)別選取抽水試驗(yàn)中所有抽水孔及觀測(cè)孔資料,水位擬合針對(duì)所有觀測(cè)孔,擬合數(shù)據(jù)時(shí)間段為整個(gè)抽水試驗(yàn)過(guò)程,其中包括觀測(cè)孔水位恢復(fù)階段。經(jīng)過(guò)模型參數(shù)識(shí)別,最終確定奧灰含水層水文地質(zhì)參數(shù)為:Kx=0.00152m/ min,Ky=0.00182m/min,μ* =0.45×10-5m-1。

    在整個(gè)模型識(shí)別、驗(yàn)證過(guò)程中觀測(cè)孔水位計(jì)算值與實(shí)測(cè)值的誤差統(tǒng)計(jì)見(jiàn)表3,觀測(cè)孔水位擬合曲線見(jiàn)圖6。由結(jié)果可以看出奧灰含水層的實(shí)測(cè)水位與計(jì)算水位相對(duì)誤差集中在1%內(nèi),符合《中華人民共和國(guó)國(guó)家標(biāo)準(zhǔn)GB/T14497-93地下水資源管理模型工作要求》的規(guī)定,模擬效果較理想。

    表3 奧灰含水層地下水位計(jì)算值與觀測(cè)值擬合統(tǒng)計(jì)

    圖6 觀測(cè)孔水位擬合曲線

    5.4 模擬結(jié)果分析

    圖7為奧灰含水層數(shù)值模擬初始及結(jié)束時(shí)地下水流場(chǎng)分布圖。對(duì)比(a)、(b)可知,抽水試驗(yàn)過(guò)程中,地下水流場(chǎng)分布及流向變化總體變化不大。同時(shí),除CK8及CK22兩孔降深相對(duì)較大外,其他觀測(cè)孔水位沒(méi)有明顯波動(dòng),說(shuō)明奧灰含水層橫向上溶(裂)隙連通性差、水力聯(lián)系差。究其原因,研究區(qū)奧灰含水層處于巖溶水弱徑流-滯流區(qū),水動(dòng)力條件較差,溶隙、溶孔等含導(dǎo)水構(gòu)造水力聯(lián)系差,致使地下水運(yùn)移緩慢,水力交替周期長(zhǎng),橫向滲透性及富水性差別較大。

    圖7 模擬初始及結(jié)束時(shí)奧灰含水層地下水流場(chǎng)

    6. 結(jié)論

    奧灰含水層水文地質(zhì)條件研究是下組煤開(kāi)采安全評(píng)價(jià)的前提和關(guān)鍵。本文在研究區(qū)選取了12個(gè)奧灰?guī)r溶水取樣點(diǎn),通過(guò)對(duì)水樣進(jìn)行水化學(xué)分析,得到了奧灰?guī)r溶水水化學(xué)參數(shù),并采用 Piper三線圖方法分析了地下水的水質(zhì)特征。將水化學(xué)資料整理分組,發(fā)現(xiàn)各組水化學(xué)特征反映了井田不同區(qū)段巖溶水演化階段,據(jù)此將井田奧灰流場(chǎng)劃分為弱徑流區(qū)、弱徑流 -滯流過(guò)渡區(qū)及滯流區(qū)。井田奧灰群孔抽水試驗(yàn)驗(yàn)證了水化學(xué)特征分析結(jié)論,即井田奧灰含水層橫向水力聯(lián)系弱,水動(dòng)力條件差。同時(shí),運(yùn)用地下水模擬軟件 GMS對(duì)井田奧灰含水層地下水流場(chǎng)變化規(guī)律進(jìn)行了模擬和水文地質(zhì)參數(shù)反演。結(jié)果表明群孔抽水影響范圍小;且?jiàn)W灰含水層滲透系數(shù)及貯水系數(shù)量級(jí)較小,分別為 Kx=0.00152m/ min,Ky=0.00182m/min,μ* =0.45×10-5m-1。本文初步掌握了該礦奧灰含水層水文地質(zhì)特征,為下組煤開(kāi)采奧灰?guī)r溶水防治提供了可靠依據(jù)。

    [1] 虎維岳,田干.我國(guó)煤礦水害類(lèi)型及其防治對(duì)策[J].煤炭科學(xué)技術(shù),2010,38(1):92-96.

    [2] 虎維岳,何滿潮.深部煤炭資源及開(kāi)發(fā)地質(zhì)條件[M].北京:煤炭工業(yè)出版社,2008.

    [3] 王希良,鄭世書(shū),郭勵(lì)生.邢臺(tái)煤礦下組煤開(kāi)采水文地質(zhì)條件評(píng)價(jià)及突水危險(xiǎn)性預(yù)測(cè)[J].煤田地質(zhì)與勘探,2000,28 (3):43-45.

    [4] 韓行瑞,魯榮安,李慶松,等.巖溶水系統(tǒng) -山西巖溶大泉研究[M].北京:地質(zhì)出版社,1993.306-313.

    [5] 牛威.天橋泉域巖溶水水資源評(píng)價(jià)[J].山西水利科技,2007,37(2):58-60.

    [6] 徐慧珍,段秀銘,高贊東.濟(jì)南泉域排泄區(qū)巖溶地下水水化學(xué)特征[J].水文地質(zhì)工程地質(zhì),2007,51(3):15-19.

    [7] 裴捍華,楊親民,郭振中,等.山西巖溶水強(qiáng)徑流帶的成因類(lèi)型及其水文地質(zhì)特[J].中國(guó)巖溶,2003,22(3):219-224.

    [8] [9]曹金亮,韓穎,袁新華,等.天橋泉域巖溶水系統(tǒng)水動(dòng)力場(chǎng)、水化學(xué)場(chǎng)特征分析[J].中國(guó)巖溶,2005,24(4):312-317.

    [9] 徐連利,牛志剛,高建中.大口徑群孔抽水試驗(yàn)成果在評(píng)價(jià)趙固一井水文地質(zhì)條件中的應(yīng)用[J].中國(guó)煤田地質(zhì),2005,17(6):17-19.

    [10] 祝曉彬.地下水模擬系統(tǒng)(GMS)軟件[J].水文地質(zhì)工程地質(zhì),2003,30(5):53-55.

    (責(zé)任編輯 陳永康)

    book=65,ebook=65

    P641

    A

    1008-7257(2012)02-0075-05

    2012-03-22

    李國(guó)紅(1974-),男,河南洛陽(yáng)人,河南省煤田地質(zhì)局二隊(duì)助理工程師。

    猜你喜歡
    水化學(xué)井田灰?guī)r
    幸福湖泉域地下水水化學(xué)特征及其演化機(jī)制
    人民黃河(2025年1期)2025-01-16 00:00:00
    灰?guī)r礦礦石分級(jí)利用的應(yīng)用實(shí)踐
    山東冶金(2022年2期)2022-08-08 01:51:38
    秦皇島北部地表水和地下水水化學(xué)特征及水質(zhì)評(píng)價(jià)
    煤礦地面水害區(qū)域治理工程鉆進(jìn)層位控制
    關(guān)于大平井田煤層缺失帶成因的探討
    晉煤成莊井田煤儲(chǔ)層特征及其對(duì)煤層氣成藏的控制
    衡水市新近系館陶組地?zé)崴瘜W(xué)特征與分布規(guī)律研究
    筠連煤層氣井產(chǎn)出水化學(xué)特征及意義
    金湖凹陷西斜坡阜寧組灰?guī)r段混積儲(chǔ)層特征
    不同圍壓下茅口灰?guī)r漸進(jìn)性破壞的試驗(yàn)研究
    五月玫瑰六月丁香| 国产精品精品国产色婷婷| 黑人高潮一二区| 久久人人爽av亚洲精品天堂 | 久久午夜福利片| 午夜爱爱视频在线播放| 亚洲在线观看片| 免费看光身美女| 久久热精品热| 18禁裸乳无遮挡动漫免费视频 | 看黄色毛片网站| 欧美+日韩+精品| 国国产精品蜜臀av免费| 日韩人妻高清精品专区| 黄色欧美视频在线观看| 久久久a久久爽久久v久久| 色综合色国产| 视频区图区小说| 嫩草影院精品99| 亚洲av免费高清在线观看| 日本一二三区视频观看| 在线观看国产h片| 九九久久精品国产亚洲av麻豆| 国产熟女欧美一区二区| 99re6热这里在线精品视频| 久久这里有精品视频免费| 国产精品不卡视频一区二区| 亚洲精品视频女| 亚州av有码| 国产av国产精品国产| 欧美日韩综合久久久久久| 精品99又大又爽又粗少妇毛片| 国产中年淑女户外野战色| 韩国高清视频一区二区三区| 亚洲精品成人av观看孕妇| 亚洲av在线观看美女高潮| 尤物成人国产欧美一区二区三区| av国产久精品久网站免费入址| 亚洲欧洲日产国产| 亚洲精品国产av成人精品| 狠狠精品人妻久久久久久综合| 国产在线男女| av在线观看视频网站免费| 有码 亚洲区| 插阴视频在线观看视频| 欧美一级a爱片免费观看看| 白带黄色成豆腐渣| 天天一区二区日本电影三级| 国产有黄有色有爽视频| 欧美精品国产亚洲| 少妇 在线观看| 欧美97在线视频| 全区人妻精品视频| 五月伊人婷婷丁香| 一二三四中文在线观看免费高清| 久久99热6这里只有精品| 看黄色毛片网站| 在线观看av片永久免费下载| av又黄又爽大尺度在线免费看| 免费观看无遮挡的男女| 亚洲精品乱码久久久v下载方式| 少妇的逼水好多| 九草在线视频观看| 久久久久精品久久久久真实原创| 在线观看国产h片| 亚洲精品第二区| av专区在线播放| 国产免费福利视频在线观看| 国产亚洲午夜精品一区二区久久 | 女人久久www免费人成看片| 少妇人妻一区二区三区视频| 搡女人真爽免费视频火全软件| 黄片无遮挡物在线观看| 国产有黄有色有爽视频| 黄色配什么色好看| 久久久久久久亚洲中文字幕| 日本一本二区三区精品| 男女那种视频在线观看| av线在线观看网站| 99热这里只有精品一区| 国产综合懂色| 三级国产精品片| 欧美 日韩 精品 国产| 91精品一卡2卡3卡4卡| 亚洲欧美精品自产自拍| 国产成人精品婷婷| 成年免费大片在线观看| 国产 一区 欧美 日韩| 不卡视频在线观看欧美| 亚洲精品亚洲一区二区| 深爱激情五月婷婷| 晚上一个人看的免费电影| 午夜免费男女啪啪视频观看| 91aial.com中文字幕在线观看| 我的女老师完整版在线观看| 搡老乐熟女国产| 久久久久国产精品人妻一区二区| 国产日韩欧美亚洲二区| 亚洲av日韩在线播放| 26uuu在线亚洲综合色| 久久精品久久精品一区二区三区| 女人被狂操c到高潮| 特级一级黄色大片| eeuss影院久久| 性插视频无遮挡在线免费观看| 久久久久网色| 国产精品嫩草影院av在线观看| 免费av观看视频| 精品久久久久久久人妻蜜臀av| 成人一区二区视频在线观看| 免费看光身美女| 亚洲第一区二区三区不卡| 久久久午夜欧美精品| 婷婷色综合大香蕉| 久久久精品欧美日韩精品| 嫩草影院精品99| 午夜福利视频精品| 国产视频内射| 国产视频内射| 中文字幕亚洲精品专区| 99热这里只有是精品在线观看| 五月天丁香电影| 成人欧美大片| 日韩av不卡免费在线播放| 成人特级av手机在线观看| 最近的中文字幕免费完整| 免费av观看视频| 哪个播放器可以免费观看大片| 一级毛片电影观看| 99久久精品国产国产毛片| 青春草国产在线视频| 国产成人精品一,二区| 日韩国内少妇激情av| 久久女婷五月综合色啪小说 | 网址你懂的国产日韩在线| 国产黄片美女视频| 夜夜看夜夜爽夜夜摸| 亚洲av欧美aⅴ国产| 亚洲真实伦在线观看| 大话2 男鬼变身卡| 国内揄拍国产精品人妻在线| 一本色道久久久久久精品综合| 国产亚洲91精品色在线| 日本一二三区视频观看| 国产亚洲一区二区精品| 在线观看一区二区三区激情| 亚洲在久久综合| 国产精品一区二区性色av| 亚洲成人精品中文字幕电影| 国产精品av视频在线免费观看| 嫩草影院精品99| 亚洲欧美清纯卡通| 最近最新中文字幕大全电影3| 成人亚洲欧美一区二区av| 亚洲av福利一区| 日韩亚洲欧美综合| 新久久久久国产一级毛片| 一级毛片 在线播放| 哪个播放器可以免费观看大片| 97人妻精品一区二区三区麻豆| 又粗又硬又长又爽又黄的视频| 大陆偷拍与自拍| 亚洲精品aⅴ在线观看| 国产伦精品一区二区三区视频9| 精品国产三级普通话版| 丰满少妇做爰视频| 国产成人午夜福利电影在线观看| 成人黄色视频免费在线看| 亚洲欧美日韩卡通动漫| 欧美极品一区二区三区四区| 精品久久久久久久人妻蜜臀av| 街头女战士在线观看网站| 久久精品久久精品一区二区三区| 午夜精品国产一区二区电影 | 天堂中文最新版在线下载 | 2021少妇久久久久久久久久久| 国产色爽女视频免费观看| 亚洲av欧美aⅴ国产| 日韩成人av中文字幕在线观看| 久久精品久久久久久久性| 搞女人的毛片| 草草在线视频免费看| 男插女下体视频免费在线播放| 欧美一级a爱片免费观看看| 我的老师免费观看完整版| 一级毛片电影观看| 国产综合懂色| 欧美亚洲 丝袜 人妻 在线| 免费观看性生交大片5| 极品少妇高潮喷水抽搐| 成人亚洲欧美一区二区av| 免费不卡的大黄色大毛片视频在线观看| 亚洲综合精品二区| 91精品一卡2卡3卡4卡| av国产免费在线观看| 少妇的逼好多水| a级毛色黄片| 成人亚洲欧美一区二区av| 亚洲自偷自拍三级| 亚洲欧美一区二区三区国产| 男人添女人高潮全过程视频| 大香蕉久久网| 亚洲av日韩在线播放| 日韩欧美精品v在线| 18禁裸乳无遮挡动漫免费视频 | 免费人成在线观看视频色| kizo精华| 国产午夜精品一二区理论片| 亚洲无线观看免费| 麻豆成人av视频| 亚洲国产日韩一区二区| 亚洲精品自拍成人| 少妇人妻 视频| av在线播放精品| 日本欧美国产在线视频| 九九爱精品视频在线观看| 白带黄色成豆腐渣| 国产精品国产三级专区第一集| 国产成人福利小说| 午夜福利视频精品| 国产一级毛片在线| 欧美一级a爱片免费观看看| 欧美老熟妇乱子伦牲交| 日本与韩国留学比较| 97超视频在线观看视频| 国国产精品蜜臀av免费| 国产精品嫩草影院av在线观看| 国产爽快片一区二区三区| 国产免费又黄又爽又色| 久久久精品免费免费高清| 午夜视频国产福利| 国产熟女欧美一区二区| 国产伦精品一区二区三区视频9| 亚洲高清免费不卡视频| 日韩国内少妇激情av| 亚洲精品成人久久久久久| 亚洲欧美中文字幕日韩二区| 国产高清不卡午夜福利| 在线观看美女被高潮喷水网站| 高清毛片免费看| 国产乱人视频| av在线播放精品| 人妻 亚洲 视频| 亚洲av成人精品一二三区| 午夜激情福利司机影院| 久久久久久久久久成人| 亚洲av免费在线观看| 日韩 亚洲 欧美在线| 网址你懂的国产日韩在线| xxx大片免费视频| 18禁动态无遮挡网站| 22中文网久久字幕| 丰满少妇做爰视频| 白带黄色成豆腐渣| 在线观看一区二区三区| 日本一二三区视频观看| 国产成人freesex在线| 哪个播放器可以免费观看大片| 免费观看性生交大片5| 大又大粗又爽又黄少妇毛片口| 自拍偷自拍亚洲精品老妇| 婷婷色综合www| 免费看光身美女| 亚洲av.av天堂| 亚洲成人一二三区av| 边亲边吃奶的免费视频| 一本一本综合久久| 欧美+日韩+精品| 亚洲内射少妇av| 久久精品国产亚洲av天美| 亚洲精品视频女| 美女cb高潮喷水在线观看| 精品少妇黑人巨大在线播放| 搡老乐熟女国产| 99热这里只有是精品50| 成人黄色视频免费在线看| 亚洲成人久久爱视频| 在线播放无遮挡| 国产亚洲91精品色在线| 男女那种视频在线观看| 亚洲欧美成人综合另类久久久| 国产精品熟女久久久久浪| 精品酒店卫生间| 中国国产av一级| 青春草视频在线免费观看| 亚洲怡红院男人天堂| 又爽又黄无遮挡网站| 免费观看av网站的网址| 六月丁香七月| 久久久久久久久久久免费av| 国产高清有码在线观看视频| 久久久久久久亚洲中文字幕| 人妻制服诱惑在线中文字幕| 欧美日韩视频精品一区| 69av精品久久久久久| av黄色大香蕉| 成人亚洲欧美一区二区av| 婷婷色av中文字幕| 国产午夜精品一二区理论片| 国产成人精品久久久久久| 大片电影免费在线观看免费| 少妇 在线观看| 亚洲精华国产精华液的使用体验| 国产精品久久久久久久电影| 中国三级夫妇交换| 99久久精品国产国产毛片| 可以在线观看毛片的网站| 神马国产精品三级电影在线观看| av在线天堂中文字幕| 亚洲人成网站高清观看| 免费观看a级毛片全部| 欧美xxxx性猛交bbbb| 免费观看在线日韩| 免费观看的影片在线观看| 亚洲人成网站在线观看播放| 国产精品三级大全| 在线观看一区二区三区激情| av免费在线看不卡| 亚洲美女搞黄在线观看| 亚洲欧美成人综合另类久久久| 丝袜脚勾引网站| 夫妻性生交免费视频一级片| 国产成人免费无遮挡视频| 亚洲,欧美,日韩| 国产亚洲最大av| 成人一区二区视频在线观看| 国产精品久久久久久精品电影| 性色av一级| 在线观看av片永久免费下载| 亚洲欧美一区二区三区国产| 久久久亚洲精品成人影院| 亚洲av免费高清在线观看| 欧美最新免费一区二区三区| 国精品久久久久久国模美| 免费av观看视频| 亚洲国产欧美人成| 日本一本二区三区精品| 国产成人免费观看mmmm| 尤物成人国产欧美一区二区三区| 3wmmmm亚洲av在线观看| 香蕉精品网在线| freevideosex欧美| 日本欧美国产在线视频| 人妻一区二区av| 自拍欧美九色日韩亚洲蝌蚪91 | 日本色播在线视频| 久久99热6这里只有精品| 婷婷色麻豆天堂久久| 人妻少妇偷人精品九色| 日韩一本色道免费dvd| 国产精品久久久久久精品电影小说 | 亚洲在线观看片| 欧美bdsm另类| 中文字幕亚洲精品专区| 少妇被粗大猛烈的视频| 欧美极品一区二区三区四区| 欧美亚洲 丝袜 人妻 在线| 51国产日韩欧美| 亚洲欧美一区二区三区黑人 | 2018国产大陆天天弄谢| 成人综合一区亚洲| 最近最新中文字幕免费大全7| 好男人在线观看高清免费视频| 久热这里只有精品99| 日韩欧美精品v在线| 成人漫画全彩无遮挡| 精品人妻一区二区三区麻豆| 天堂中文最新版在线下载 | 夫妻午夜视频| 成人鲁丝片一二三区免费| 精品少妇黑人巨大在线播放| 成人黄色视频免费在线看| 性插视频无遮挡在线免费观看| 亚洲图色成人| 日本免费在线观看一区| 亚洲欧洲日产国产| 超碰97精品在线观看| 高清午夜精品一区二区三区| 久久影院123| 免费黄网站久久成人精品| 夫妻性生交免费视频一级片| 99久久中文字幕三级久久日本| 久久久色成人| 成人综合一区亚洲| 久久热精品热| 久久综合国产亚洲精品| 免费黄频网站在线观看国产| 天天躁夜夜躁狠狠久久av| 亚洲精品国产色婷婷电影| 女的被弄到高潮叫床怎么办| 国产毛片在线视频| 国产精品av视频在线免费观看| 好男人在线观看高清免费视频| 精品久久久久久久人妻蜜臀av| www.色视频.com| 久久精品国产a三级三级三级| 亚洲成人中文字幕在线播放| 麻豆乱淫一区二区| 哪个播放器可以免费观看大片| 午夜激情福利司机影院| 老司机影院毛片| 色播亚洲综合网| 不卡视频在线观看欧美| 天天躁日日操中文字幕| 亚洲欧美成人精品一区二区| 久久99热这里只频精品6学生| 免费看不卡的av| 欧美少妇被猛烈插入视频| 人妻 亚洲 视频| 搡老乐熟女国产| 91精品一卡2卡3卡4卡| 免费大片18禁| 成人午夜精彩视频在线观看| 国产伦精品一区二区三区视频9| 女人被狂操c到高潮| 深夜a级毛片| 国产v大片淫在线免费观看| 男女边摸边吃奶| 狂野欧美白嫩少妇大欣赏| 国产女主播在线喷水免费视频网站| 色播亚洲综合网| 国产欧美亚洲国产| 王馨瑶露胸无遮挡在线观看| 国产成人精品久久久久久| 亚洲成人久久爱视频| 国产有黄有色有爽视频| 下体分泌物呈黄色| 亚洲精品久久久久久婷婷小说| 美女内射精品一级片tv| 亚洲欧美日韩东京热| 激情 狠狠 欧美| 亚洲精品一二三| 免费大片18禁| 高清视频免费观看一区二区| 久久人人爽av亚洲精品天堂 | 国产av码专区亚洲av| 国产精品熟女久久久久浪| 久久热精品热| 国产一级毛片在线| 久久久久久伊人网av| 亚洲av成人精品一区久久| 美女国产视频在线观看| 午夜精品一区二区三区免费看| 亚洲综合色惰| 成人欧美大片| 嫩草影院精品99| 国产欧美亚洲国产| 能在线免费看毛片的网站| 亚洲av日韩在线播放| 少妇 在线观看| 亚洲熟女精品中文字幕| 精品人妻一区二区三区麻豆| 一级a做视频免费观看| 久久久欧美国产精品| 午夜亚洲福利在线播放| 九九在线视频观看精品| 最后的刺客免费高清国语| 两个人的视频大全免费| 最近手机中文字幕大全| 日韩欧美 国产精品| 久久精品国产a三级三级三级| 国产精品福利在线免费观看| 女人被狂操c到高潮| 免费看a级黄色片| 亚洲成色77777| 97人妻精品一区二区三区麻豆| 观看美女的网站| 亚洲在线观看片| 成年女人看的毛片在线观看| 精品午夜福利在线看| 亚洲图色成人| av播播在线观看一区| 久久精品夜色国产| 亚洲精品日韩av片在线观看| 国产在线一区二区三区精| 国产精品人妻久久久影院| 久久这里有精品视频免费| 听说在线观看完整版免费高清| 日韩在线高清观看一区二区三区| av在线天堂中文字幕| 日韩欧美一区视频在线观看 | 22中文网久久字幕| 国产精品久久久久久精品古装| 2021天堂中文幕一二区在线观| 香蕉精品网在线| 美女国产视频在线观看| 真实男女啪啪啪动态图| 日韩大片免费观看网站| 成人欧美大片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品人妻久久久久久| 欧美另类一区| 国产黄片美女视频| 一二三四中文在线观看免费高清| 亚洲成色77777| 国产片特级美女逼逼视频| 久久久午夜欧美精品| 亚洲美女搞黄在线观看| 国产亚洲午夜精品一区二区久久 | 性插视频无遮挡在线免费观看| 欧美三级亚洲精品| 国产精品一区www在线观看| 少妇人妻一区二区三区视频| 尤物成人国产欧美一区二区三区| 亚洲内射少妇av| 日韩,欧美,国产一区二区三区| 亚洲国产最新在线播放| 亚洲激情五月婷婷啪啪| 别揉我奶头 嗯啊视频| 精品久久国产蜜桃| 国产欧美亚洲国产| 麻豆成人午夜福利视频| 中文欧美无线码| 国产综合懂色| 综合色av麻豆| 成年版毛片免费区| 国产美女午夜福利| 亚洲精品日韩av片在线观看| 免费少妇av软件| 一级二级三级毛片免费看| 国产老妇女一区| 99热全是精品| 我的女老师完整版在线观看| 国产伦理片在线播放av一区| 久久久久九九精品影院| 亚洲欧美精品专区久久| 免费黄频网站在线观看国产| 国产精品女同一区二区软件| 一级毛片黄色毛片免费观看视频| 亚洲精品第二区| 日韩制服骚丝袜av| 日韩在线高清观看一区二区三区| 欧美成人a在线观看| 日韩av不卡免费在线播放| 免费黄网站久久成人精品| 国产精品三级大全| 人人妻人人澡人人爽人人夜夜| 国产片特级美女逼逼视频| 亚洲精品成人久久久久久| 大又大粗又爽又黄少妇毛片口| 听说在线观看完整版免费高清| 国产黄片美女视频| 亚洲成人久久爱视频| 日本午夜av视频| 联通29元200g的流量卡| 日日摸夜夜添夜夜爱| 永久免费av网站大全| 极品教师在线视频| 色哟哟·www| 色综合色国产| 51国产日韩欧美| 国产精品久久久久久精品电影| a级一级毛片免费在线观看| av卡一久久| 亚洲精品亚洲一区二区| 尤物成人国产欧美一区二区三区| 大码成人一级视频| 一级毛片aaaaaa免费看小| 男女国产视频网站| 国产 一区 欧美 日韩| 日韩人妻高清精品专区| 日韩av免费高清视频| av线在线观看网站| 春色校园在线视频观看| 欧美高清性xxxxhd video| 亚洲精品乱码久久久久久按摩| av又黄又爽大尺度在线免费看| 中文字幕制服av| 3wmmmm亚洲av在线观看| 听说在线观看完整版免费高清| 人妻少妇偷人精品九色| 日韩av在线免费看完整版不卡| 欧美成人一区二区免费高清观看| 看十八女毛片水多多多| 欧美精品国产亚洲| 精品一区二区三区视频在线| 2021少妇久久久久久久久久久| 国产成人免费观看mmmm| 爱豆传媒免费全集在线观看| 韩国高清视频一区二区三区| 国产有黄有色有爽视频| 欧美性感艳星| 如何舔出高潮| 国产爱豆传媒在线观看| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品影视一区二区三区av| 午夜精品国产一区二区电影 | 精品久久久精品久久久| 一级爰片在线观看| 又大又黄又爽视频免费| 免费黄网站久久成人精品| 久久久久国产精品人妻一区二区| 精品一区二区三卡| 久久热精品热| 性色avwww在线观看| .国产精品久久| 夫妻午夜视频| 国产亚洲91精品色在线| 九九爱精品视频在线观看| 我的老师免费观看完整版| 国产成人免费观看mmmm| 中文字幕人妻熟人妻熟丝袜美| 一区二区三区精品91| 直男gayav资源| 亚洲欧美日韩东京热| 亚洲av.av天堂| 成人特级av手机在线观看| 亚洲av不卡在线观看| 一个人看的www免费观看视频| 亚州av有码| 国产成人精品婷婷| 亚洲欧美精品专区久久| 内地一区二区视频在线| 一区二区三区四区激情视频| 免费在线观看成人毛片|