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

    西南巖溶山區(qū)某石化類場(chǎng)地地下水污染風(fēng)險(xiǎn)評(píng)價(jià)

    2022-03-29 11:54:04陳植華
    地下水 2022年1期
    關(guān)鍵詞:關(guān)嶺巖溶裂隙

    王 川,陳植華

    (1.四川省地質(zhì)工程勘察院集團(tuán)有限公司,四川 成都 610072;2.中國(guó)地質(zhì)大學(xué)(武漢)環(huán)境學(xué)院,湖北 武漢 430074)

    西南巖溶山區(qū),碳酸鹽巖分布廣泛且溶蝕作用強(qiáng)烈,水文地質(zhì)條件復(fù)雜,相比于非碳酸鹽巖地區(qū),巖溶地下水系統(tǒng)受到污染的風(fēng)險(xiǎn)更大[1],一旦受到污染,其恢復(fù)及治理的難度十分巨大。在巖溶山區(qū)進(jìn)行石化類等高污染風(fēng)險(xiǎn)項(xiàng)目的選址工作時(shí),不僅要考慮經(jīng)濟(jì)、交通以及電力等常規(guī)因素,還應(yīng)關(guān)注擬選場(chǎng)址周邊的地下水環(huán)境問(wèn)題,評(píng)價(jià)項(xiàng)目建設(shè)對(duì)地下水環(huán)境可能造成的污染風(fēng)險(xiǎn)[2-6],為項(xiàng)目科學(xué)合理的選址提供參考。

    由于巖溶地下水系統(tǒng)的特殊性,在進(jìn)行巖溶場(chǎng)區(qū)地下水污染風(fēng)險(xiǎn)評(píng)價(jià)時(shí),不僅在評(píng)價(jià)方法的選擇上會(huì)遇到難題,更重要的是如何查清場(chǎng)區(qū)巖溶水文地質(zhì)條件以及確定其主要的風(fēng)險(xiǎn)評(píng)價(jià)對(duì)象。本文以《污染場(chǎng)地風(fēng)險(xiǎn)評(píng)價(jià)導(dǎo)則(HJ25.3-2014)》[7]為基礎(chǔ),結(jié)合國(guó)內(nèi)外學(xué)者在地下水污染風(fēng)險(xiǎn)評(píng)價(jià)研究中的探索,總結(jié)分析了常規(guī)方法的適用條件及優(yōu)缺點(diǎn),在場(chǎng)地地下水污染風(fēng)險(xiǎn)識(shí)別及源項(xiàng)分析基礎(chǔ)上,建立了基于數(shù)值模擬結(jié)果的污染場(chǎng)地地下水污染風(fēng)險(xiǎn)評(píng)價(jià)體系[8]。

    1 地質(zhì)環(huán)境概況

    1.1 自然地理

    該石化類場(chǎng)地位于貴州省中部,區(qū)域地貌屬于云貴高原中部滇東高原烏蒙山系的延續(xù)部分,處于黔西山原向黔中丘原盆地的過(guò)渡帶,地形起伏較大,地勢(shì)西高東低,屬高中山山地地貌,河流強(qiáng)烈的侵蝕切割作用導(dǎo)致地形支離破碎,高山深谷、陡壁懸崖以及河澗河流等復(fù)雜地貌在研究區(qū)隨處可見(jiàn)。

    1.2 地質(zhì)條件

    場(chǎng)地出露地層主要有三疊系中統(tǒng)關(guān)嶺組(T2g)、法郎組(T2f)以及三疊系上統(tǒng)須家河組(T3x),地層巖性以碳酸鹽巖為主,主要可見(jiàn)中薄層至厚層灰?guī)r、白云巖、生物碎屑灰?guī)r及角礫白云巖,部分夾泥質(zhì)灰?guī)r及白云質(zhì)灰?guī)r,最頂部為須家河組石英砂巖。

    場(chǎng)地位于新場(chǎng)—張家寨向斜核部偏東南一翼,并與田家壩小背斜過(guò)渡。地層基本上以北傾為主,傾角從南部的30°~50°向北漸緩,但受場(chǎng)內(nèi)斷層的影響,局部產(chǎn)狀多變。發(fā)育大小斷層共7條。根據(jù)斷層分布特征及其產(chǎn)狀,判斷場(chǎng)地受NW-SE向和近EW向擠壓應(yīng)力作用,在發(fā)育一系列斷層的同時(shí)發(fā)育有多個(gè)小型褶皺。

    2 水文地質(zhì)條件

    2.1 地下水類型及其富水性

    場(chǎng)地地下水類型有第四系松散巖類孔隙水、碎屑巖裂隙水和碳酸鹽巖巖溶水三種(圖1)。

    圖1 場(chǎng)地水文地質(zhì)簡(jiǎn)圖

    第四系松散巖類孔隙水零星分布于沖溝及部分坡麓地區(qū),含水巖組以沖洪積、殘坡積松散沉積物為主,厚度在1.5~12 m之間,地下水分布不均,水位埋深淺,溝谷地區(qū)小于1.0 m,坡麓地區(qū)一般在10 m左右,大多存在于雨季,富水性差。

    碎屑巖裂隙水僅賦存于場(chǎng)區(qū)西北部須家河組(T3x)黃褐色砂巖的風(fēng)化裂隙之中,裂隙發(fā)育深度一般在30~50 m之間,裂隙發(fā)育程度隨深度增加而減弱,下部裂隙不發(fā)育或微發(fā)育具相對(duì)隔水層性質(zhì)。含水巖組整體滲透性較弱,富水性差,泉水流量一般0.5~22 m3/d。

    碳酸鹽巖巖溶水可細(xì)分為碳酸鹽巖溶蝕裂隙水和碳酸鹽巖溶洞水兩個(gè)亞類。碳酸鹽巖溶蝕裂隙水主要賦存于三疊系法郎組(T2f)、關(guān)嶺組一段(T2g1)和關(guān)嶺組三段(T2g3)的白云巖、泥質(zhì)白云巖中。受斷裂構(gòu)造影響,場(chǎng)地白云巖整體比較破碎,表現(xiàn)為網(wǎng)絡(luò)狀的節(jié)理及裂隙,局部地方發(fā)育的溶蝕裂隙、溶孔及小溶洞構(gòu)成了良好的地下水賦存空間,但深部裂隙多不發(fā)育,因此地下水多賦存于淺表層且分布不均;碳酸鹽巖溶洞水主要賦存于三疊系關(guān)嶺組二段(T2g2)的灰?guī)r、白云質(zhì)灰?guī)r中,地下水的賦存及徑流以深部管道、廊道為主,局部發(fā)育較大的溶洞,形成地下暗河,故富水性極不均勻,地下水水位埋深一般較大。

    2.2 地下水補(bǔ)徑排特征

    場(chǎng)地地下水主要接受大氣降水入滲及河流的滲漏補(bǔ)給。地下水的徑流在場(chǎng)地北側(cè)主要受到地形影響,自高向低,西北側(cè)偏西南方向徑流,東側(cè)偏西北方向徑流,南側(cè)主要是通過(guò)巖溶管道向北徑流。場(chǎng)地關(guān)嶺組三段與關(guān)嶺組二段分屬不同的水文地質(zhì)單元,地下水不連續(xù),徑流特征各有不同,關(guān)嶺組三段以淺層溶蝕裂隙水為主,關(guān)嶺組二段以深部管道巖溶地下水為主。地下水主要通過(guò)蒸發(fā)蒸騰、在地勢(shì)低洼處以泉的形式以及沿沖溝向下游的水庫(kù)排泄。

    3 地下水?dāng)?shù)值模型

    場(chǎng)地出露地層主要為關(guān)嶺組三段的白云巖,地下水類型以碳酸鹽巖溶蝕裂隙水為主,地下水污染風(fēng)險(xiǎn)評(píng)價(jià)重點(diǎn)關(guān)注淺層地下水,故本次數(shù)值建模僅考慮關(guān)嶺組三段、法郎組的淺層溶蝕裂隙水、須家河組的風(fēng)化裂隙水以及分布不連續(xù)的第四系孔隙水。將含水層概化成非均質(zhì)、水平方向各向同性垂向存在差異的空間三維穩(wěn)定地下水流系統(tǒng)。

    圖2 模型二維平面剖分結(jié)果

    本次數(shù)值模型面積6.0 km2,基于FEFLOW數(shù)值模擬軟件,將模型范圍離散為不規(guī)則的三角剖分網(wǎng)格,同時(shí)對(duì)四周邊界及內(nèi)部水系進(jìn)行加密處理,最終得到研究區(qū)平面二維剖分結(jié)果(圖2),結(jié)點(diǎn)數(shù)71 971個(gè),有限單元數(shù)133 587個(gè)。地質(zhì)模型概化為三層,第一層為第四系孔隙含水巖組,第二層為強(qiáng)溶蝕裂隙及風(fēng)化裂隙含水巖組,第三層為中溶蝕裂隙及風(fēng)化裂隙含水巖組(圖3)。

    圖3 模型三維地質(zhì)模型

    模型北部、西部及南部邊界設(shè)置為零通量邊界,東部關(guān)嶺組三段與關(guān)嶺組二段分界線概化為流量邊界,上邊界設(shè)定為為降水補(bǔ)給及蒸發(fā)邊界,下邊界以裂隙—溶蝕裂隙不發(fā)育為模型下邊界,設(shè)為隔水邊界(圖4)。

    水文地質(zhì)參數(shù)是表征含水介質(zhì)儲(chǔ)水能力、釋水能力和地下水運(yùn)動(dòng)能力的指標(biāo)[49],F(xiàn)EFLOW地下水水流模型中表征含水介質(zhì)水文地質(zhì)特征的參數(shù)主要有滲透系數(shù)(K)和給水度(μ)。研究區(qū)巖性不統(tǒng)一,各層的水文地質(zhì)參數(shù)存在差異,因此需要給出各層的水文地質(zhì)參數(shù)分區(qū),并在FEFLOW中進(jìn)行分區(qū)參數(shù)的賦值。

    圖4 模型邊界條件示意圖

    通過(guò)反復(fù)調(diào)整參數(shù),識(shí)別了模型結(jié)構(gòu)、各參數(shù)和均衡要素,并使模擬計(jì)算后的模型水量達(dá)到平衡,反算誤差達(dá)到FEFLOW模型的要求(0.001),最終識(shí)別得到了場(chǎng)地地下水流場(chǎng)(圖5)。從模擬的地下水流場(chǎng)圖中可以看出,研究區(qū)局部淺層地下水徑流方向受控于內(nèi)部溪溝和溝谷,地下水整體由西南向東北方向徑流。由此可知,研究區(qū)地下水流數(shù)值模型從空間上把握了模擬區(qū)的地下水補(bǔ)徑排特征,可以基于該水流模型進(jìn)行后續(xù)的風(fēng)險(xiǎn)評(píng)價(jià)。

    圖5 識(shí)別后的地下水流場(chǎng)圖

    4 地下水污染風(fēng)險(xiǎn)評(píng)價(jià)體系建立

    4.1 風(fēng)險(xiǎn)識(shí)別

    通過(guò)收集擬建石化類項(xiàng)目設(shè)計(jì)文件并對(duì)比國(guó)內(nèi)外同行業(yè)、同類型裝置事故統(tǒng)計(jì)資料,識(shí)別項(xiàng)目在建設(shè)、生產(chǎn)及服務(wù)期滿后各個(gè)階段可能對(duì)地下水造成風(fēng)險(xiǎn)的主要污染源。

    擬建石化類項(xiàng)目產(chǎn)污環(huán)節(jié)較多,但大部分污廢水以場(chǎng)內(nèi)設(shè)備的循環(huán)使用為主,最終排污去向?yàn)槲鬯幚韽S。污水處理廠處理達(dá)標(biāo)后,采用壓力流方式送至下游水庫(kù)出水口。通過(guò)進(jìn)一步識(shí)別,該石化類場(chǎng)地建成后主要的污染源裝置為甲醇裝置及污水處理池,按照最大風(fēng)險(xiǎn)評(píng)價(jià)原則,本次研究主要考慮在天然條件下甲醇裝置事故狀態(tài)下發(fā)生泄漏以及污水處理池正常跑冒滴漏下污染物持續(xù)泄漏對(duì)地下水造成的污染風(fēng)險(xiǎn)大小。

    4.2 源項(xiàng)分析

    正常狀況下,污水處理池跑冒滴漏下的污染物進(jìn)入地下水,污水處理池面積0.036 km2,污染物類型為COD、氨氮等,按最大風(fēng)險(xiǎn)原則,正常工況且無(wú)人工防滲措施,年存水時(shí)間按365天計(jì)算,COD濃度為1 500 mg/L、氨氮濃度為50 mg/L。

    事故狀態(tài)下,考慮爆炸破壞地表防滲結(jié)構(gòu)面積為0.35 km2,收集處理時(shí)間10 h,結(jié)合場(chǎng)區(qū)內(nèi)最大滲透系數(shù),計(jì)算得到最大下滲量為0.712 5 m3。即甲醇滲入地下體積為0.712 5 m3,泄漏面積50 m2,持續(xù)時(shí)間10 h,甲醇濃度為800 kg/m3。

    4.3 風(fēng)險(xiǎn)計(jì)算

    地下水污染風(fēng)險(xiǎn)評(píng)價(jià)包括地下水污染荷載、地下水本質(zhì)脆弱性和地下水功能價(jià)值三個(gè)部分。風(fēng)險(xiǎn)計(jì)算是在風(fēng)險(xiǎn)識(shí)別及源項(xiàng)分析的基礎(chǔ)上,綜合考慮以上三個(gè)部分的內(nèi)容對(duì)研究區(qū)地下水污染風(fēng)險(xiǎn)進(jìn)行定量化評(píng)價(jià)計(jì)算。鑒于論文研究的巖溶場(chǎng)區(qū)尺度范圍較小且該場(chǎng)地已做過(guò)多次勘察工作,巖溶水文地質(zhì)條件已基本清楚,基礎(chǔ)資料相當(dāng)豐富,因此本次評(píng)價(jià)過(guò)程選用數(shù)值模擬方法進(jìn)行地下水污染風(fēng)險(xiǎn)評(píng)價(jià)計(jì)算。

    論文建立的基于數(shù)值模擬的地下水污染風(fēng)險(xiǎn)評(píng)價(jià)不需要把各影響因素按指標(biāo)進(jìn)行加權(quán)求和,只需要根據(jù)溶質(zhì)運(yùn)移模擬得到的污染物濃度遷移時(shí)空表現(xiàn),結(jié)合地下水質(zhì)量標(biāo)準(zhǔn)按照一定的劃分原則即可得到該巖溶石化類場(chǎng)地地下水污染風(fēng)險(xiǎn)的大小。

    5 地下水污染風(fēng)險(xiǎn)評(píng)價(jià)

    5.1 地下水溶質(zhì)運(yùn)移模擬

    在前文建立的地下水水流數(shù)值模型基礎(chǔ)上,對(duì)不同情形不同污染物的滲漏進(jìn)行模擬預(yù)測(cè),模型運(yùn)移時(shí)間為該石化類場(chǎng)地的服務(wù)年限50年,通過(guò)模型計(jì)算得到了不同時(shí)刻污染物在場(chǎng)地淺層地下水中的空間分布。

    5.1.1 正常工況且無(wú)人工防滲措施下污水處理池的滲漏

    從模型預(yù)測(cè)的結(jié)果來(lái)看,正常工況、污水處理池持續(xù)滲漏條件下COD及氨氮在前5年的運(yùn)移主要以包氣帶垂向遷移為主,溶蝕裂隙含水層中污染暈的水平擴(kuò)散范圍不大,最大遷移距離僅30 m,包氣帶介質(zhì)對(duì)污染物的遷移起到了一定的阻擋作用;5 a以后COD及氨氮的運(yùn)移以溶蝕裂隙含水層中的水平遷移為主,擴(kuò)散速率變大,污染暈主要受對(duì)流作用影響,在地下水流控制下污染暈以北部為優(yōu)勢(shì)遷移方向,逐步向北擴(kuò)散,污染暈范圍越來(lái)越大,溶蝕裂隙含水層中污染物的最高濃度越來(lái)越接近污染源的濃度。表1針對(duì)四個(gè)典型時(shí)間段,統(tǒng)計(jì)了溶蝕裂隙含水層中COD及氨氮污染暈的最高濃度、遷移距離、超標(biāo)面積等預(yù)測(cè)結(jié)果。

    5.1.2 事故條件且無(wú)人工防滲措施下甲醇裝置系統(tǒng)事故

    事故條件下污染物初始濃度很大,迅速經(jīng)包氣帶進(jìn)入含水層,從模型預(yù)測(cè)的結(jié)果來(lái)看,事故發(fā)生初期,污染暈僅集中在甲醇裝置周邊一定范圍內(nèi),污染暈中心濃度很高,但由于沒(méi)有污染物的持續(xù)補(bǔ)充,因此隨著時(shí)間的推移,地下水流對(duì)污染物的稀釋作用明顯,污染暈中心濃度不斷降低,污染暈擴(kuò)散范圍不斷增大,擴(kuò)散方向受地下水流控制逐漸向東部遷移。表2針對(duì)九個(gè)典型時(shí)間段,統(tǒng)計(jì)了溶蝕裂隙含水層中甲醇污染暈的最高濃度、遷移距離、超標(biāo)面積等預(yù)測(cè)結(jié)果。

    表1 污水處理池滲漏運(yùn)移預(yù)測(cè)結(jié)果表

    表2 事故狀態(tài)下甲醇運(yùn)移預(yù)測(cè)結(jié)果表

    5.2 地下水污染風(fēng)險(xiǎn)分級(jí)原則

    按照本文地下水污染風(fēng)險(xiǎn)評(píng)價(jià)體系,地下水污染荷載及地下水本質(zhì)脆弱性已體現(xiàn)在建立的地下水水流及溶質(zhì)運(yùn)移模型中,地下水功能價(jià)值改變程度主要體現(xiàn)在污染物在地下水中的濃度分布,地下水污染風(fēng)險(xiǎn)大小可根據(jù)預(yù)測(cè)的污染暈分布情況來(lái)劃分風(fēng)險(xiǎn)級(jí)別。論文研究所選污染風(fēng)險(xiǎn)因子COD、氨氮及甲醇的評(píng)價(jià)標(biāo)準(zhǔn)分別為3.0 mg/L、0.2 mg/L和200 mg/L,如果某污染物濃度超過(guò)了標(biāo)準(zhǔn)濃度,則認(rèn)為地下水受到了污染[9]。

    圖6 污染暈擴(kuò)散范圍與風(fēng)險(xiǎn)等級(jí)關(guān)系圖

    論文結(jié)合地下水溶質(zhì)運(yùn)移模擬結(jié)果,以某風(fēng)險(xiǎn)因子標(biāo)準(zhǔn)濃度對(duì)應(yīng)的濃度等值線所圈閉的污染暈作為地下水污染風(fēng)險(xiǎn)大小的度量依據(jù),結(jié)合場(chǎng)地水文地質(zhì)條件,論文將某污染荷載對(duì)研究場(chǎng)地地下水造成的污染風(fēng)險(xiǎn)劃分為五級(jí)[10-11]。當(dāng)污染物以一定源強(qiáng)進(jìn)入含水層后,如果污染暈基本上沒(méi)有擴(kuò)散(10 m以內(nèi)),則定義污染風(fēng)險(xiǎn)級(jí)別低;污染暈擴(kuò)散范圍僅分布于污染源附近(50 m以內(nèi)),則定義污染風(fēng)險(xiǎn)級(jí)別較低;如果形成了明顯的污染暈,但沒(méi)有擴(kuò)散至場(chǎng)界范圍外,則定義污染風(fēng)險(xiǎn)級(jí)別為中;如果污染暈擴(kuò)散范圍較大,已擴(kuò)散至場(chǎng)界外但沒(méi)有到達(dá)模擬邊界,則定義污染風(fēng)險(xiǎn)級(jí)別為較高;如果污染暈擴(kuò)散范圍已經(jīng)到達(dá)了東部的模擬邊界,污染物已經(jīng)危害到了關(guān)嶺組二段的巖溶管道含水層,則定義污染風(fēng)險(xiǎn)等級(jí)為高(圖6)。

    5.3 地下水污染分級(jí)及評(píng)價(jià)

    5.3.1 污水處理池正常滲漏污染風(fēng)險(xiǎn)分級(jí)

    污水處理池發(fā)生滲漏后,COD和氨氮的污染暈擴(kuò)散面積逐漸增大,在該石化類場(chǎng)地服務(wù)年限50 a后,COD及氨氮污染暈最遠(yuǎn)遷移距離分別為396 m和379 m,污染暈最大擴(kuò)散面積分別為0.279 km2和0.263 km2,污染暈擴(kuò)散長(zhǎng)軸方向未達(dá)到擬建石化類場(chǎng)地邊界,根據(jù)前文劃定的地下水污染風(fēng)險(xiǎn)等級(jí),污水處理池COD及氨氮持續(xù)滲漏造成的地下水污染風(fēng)險(xiǎn)等級(jí)均為中。

    5.3.2 甲醇裝置系統(tǒng)事故滲漏污染風(fēng)險(xiǎn)分級(jí)

    事故狀態(tài)下甲醇裝置發(fā)生泄漏,含水層中初始甲醇濃度較高,隨著時(shí)間推移,在地下水流稀釋作用下,污染暈中心濃度逐漸降低,污染暈隨水流向東部遷移,擴(kuò)散面積逐漸增大,在該石化類場(chǎng)地服務(wù)年限50 a后,污染暈最遠(yuǎn)遷移距離365 m,污染暈最大擴(kuò)散面積為0.398 km2,污染暈擴(kuò)散長(zhǎng)軸方向未達(dá)到擬建石化類場(chǎng)區(qū)地邊界,根據(jù)前文劃定的地下水污染風(fēng)險(xiǎn)等級(jí),甲醇裝置事故狀態(tài)甲醇泄漏造成的地下水污染風(fēng)險(xiǎn)等級(jí)為中。

    6 結(jié)語(yǔ)

    (1)本文構(gòu)建了基于數(shù)值模擬結(jié)果的地下水污染風(fēng)險(xiǎn)評(píng)價(jià)體系,并以西南巖溶山區(qū)某石化類場(chǎng)地為研究對(duì)象,基于該石化類場(chǎng)地水文地質(zhì)條件特征,構(gòu)建了研究場(chǎng)地的水文地質(zhì)概念模型,最終基于地下水溶質(zhì)運(yùn)移模擬結(jié)果對(duì)場(chǎng)地地下水污染風(fēng)險(xiǎn)等級(jí)進(jìn)行了劃分。

    (2)根據(jù)擬建石化類場(chǎng)地的項(xiàng)目設(shè)計(jì)文件,識(shí)別了研究場(chǎng)地主要污染源來(lái)自于污水處理池及甲醇裝置,污水處理池持續(xù)滲漏的主要污染物為COD及氨氮,濃度分別為1 500 mg/L和50 mg/L;甲醇裝置事故狀態(tài)下的污染濃度為800 kg/m3。

    (3)本文建立的基于數(shù)值模擬結(jié)果的場(chǎng)地地下水污染風(fēng)險(xiǎn)評(píng)價(jià)體系不需要將地下水污染風(fēng)險(xiǎn)的影響因素按照常規(guī)方法進(jìn)行指標(biāo)的加權(quán)求和,只需根據(jù)溶質(zhì)運(yùn)移模擬得到的污染物濃度遷移時(shí)空表現(xiàn),結(jié)合評(píng)價(jià)標(biāo)準(zhǔn)按照一定的劃分原則即可得到評(píng)價(jià)場(chǎng)地的地下水污染風(fēng)險(xiǎn)大小。

    (4)論文根據(jù)建立的風(fēng)險(xiǎn)評(píng)價(jià)標(biāo)準(zhǔn)并結(jié)合地下水溶質(zhì)運(yùn)移的模擬結(jié)果,得出該石化類場(chǎng)地污水處理池持續(xù)滲漏以及甲醇裝置事故滲漏條件下對(duì)場(chǎng)地地下水造成的污染風(fēng)險(xiǎn)等級(jí)均為中等,即污染暈的擴(kuò)散范圍沒(méi)有遷移出廠界,更沒(méi)有到達(dá)模型東部關(guān)嶺組二段與關(guān)嶺組三段的分界線,各風(fēng)險(xiǎn)因子不會(huì)對(duì)關(guān)嶺組二段深部巖溶管道地下水造成污染。

    猜你喜歡
    關(guān)嶺巖溶裂隙
    某石灰?guī)r礦區(qū)巖溶涌水治理處理方法
    裂隙腦室綜合征的診斷治療新進(jìn)展
    古生物發(fā)現(xiàn)之旅——關(guān)嶺化石群國(guó)家地質(zhì)公園
    探討巖溶區(qū)高速公路勘察技術(shù)方法
    裂隙燈檢查的個(gè)性化應(yīng)用(下)
    冰雪小鎮(zhèn) 關(guān)嶺永寧滑雪場(chǎng)
    開(kāi)磷為關(guān)嶺縣送肥保春耕
    高密度電法在巖溶區(qū)隧道勘察中的應(yīng)用
    《老炮兒》:在時(shí)代裂隙中揚(yáng)棄焦慮
    電影新作(2016年1期)2016-02-27 09:16:48
    非完全充填裂隙的應(yīng)力應(yīng)變特性研究
    河南科技(2014年11期)2014-02-27 14:09:42
    国产一区二区在线观看日韩| 一级黄片播放器| 久久久久久久国产电影| 欧美变态另类bdsm刘玥| 国产成人精品久久久久久| 国产精品无大码| 欧美一区二区国产精品久久精品| 国产麻豆成人av免费视频| 99久久成人亚洲精品观看| 国产伦理片在线播放av一区| 亚洲精品亚洲一区二区| 久久久久久国产a免费观看| 久久国内精品自在自线图片| 亚洲不卡免费看| 中文资源天堂在线| 中文欧美无线码| 国产精品一区二区三区四区久久| 18禁裸乳无遮挡免费网站照片| 99热这里只有精品一区| 看免费成人av毛片| 久久久久精品久久久久真实原创| 69av精品久久久久久| 国产一级毛片在线| 女人久久www免费人成看片 | 日本免费a在线| 欧美色视频一区免费| 高清在线视频一区二区三区 | 九九热线精品视视频播放| 中文乱码字字幕精品一区二区三区 | 亚洲av熟女| 韩国高清视频一区二区三区| 麻豆成人av视频| 蜜臀久久99精品久久宅男| 九九热线精品视视频播放| 国产精品美女特级片免费视频播放器| 日韩欧美在线乱码| 村上凉子中文字幕在线| 久久精品久久精品一区二区三区| 男人狂女人下面高潮的视频| 白带黄色成豆腐渣| 免费av毛片视频| 国产成人a∨麻豆精品| 晚上一个人看的免费电影| 麻豆成人av视频| 亚洲欧美精品专区久久| av在线播放精品| 婷婷色av中文字幕| 亚洲色图av天堂| 黄片wwwwww| 一夜夜www| 久久久国产成人精品二区| 91狼人影院| 欧美日本亚洲视频在线播放| 久久6这里有精品| 久久草成人影院| 亚洲在久久综合| 亚洲精品,欧美精品| 中文字幕熟女人妻在线| 国产精品国产三级国产专区5o | 丝袜美腿在线中文| 老女人水多毛片| 三级男女做爰猛烈吃奶摸视频| 搡女人真爽免费视频火全软件| 色播亚洲综合网| 一卡2卡三卡四卡精品乱码亚洲| 熟妇人妻久久中文字幕3abv| 91aial.com中文字幕在线观看| 青春草亚洲视频在线观看| 人妻系列 视频| 人人妻人人看人人澡| 男女啪啪激烈高潮av片| 久久久久久久午夜电影| 成年免费大片在线观看| 大香蕉久久网| 日韩欧美三级三区| 秋霞在线观看毛片| 看免费成人av毛片| 成人国产麻豆网| 久久韩国三级中文字幕| 成人美女网站在线观看视频| 两个人的视频大全免费| 99九九线精品视频在线观看视频| a级毛色黄片| 日日摸夜夜添夜夜添av毛片| 男女下面进入的视频免费午夜| 国产精品久久视频播放| 天堂影院成人在线观看| 我的老师免费观看完整版| 久久精品影院6| 久久久久网色| 欧美日韩国产亚洲二区| 草草在线视频免费看| 亚洲高清免费不卡视频| 老师上课跳d突然被开到最大视频| 插逼视频在线观看| 麻豆国产97在线/欧美| 亚洲欧美一区二区三区国产| 国产免费男女视频| 久久99热这里只有精品18| av国产久精品久网站免费入址| 女的被弄到高潮叫床怎么办| 欧美激情久久久久久爽电影| 亚洲伊人久久精品综合 | 欧美日韩国产亚洲二区| 精品一区二区免费观看| av女优亚洲男人天堂| 国产伦精品一区二区三区四那| 2022亚洲国产成人精品| 天天躁日日操中文字幕| 欧美97在线视频| 国产精品综合久久久久久久免费| 青青草视频在线视频观看| 一级av片app| 成年版毛片免费区| 最近最新中文字幕免费大全7| 日本一二三区视频观看| 亚洲人成网站在线播| 亚洲av中文字字幕乱码综合| 亚洲精品乱码久久久久久按摩| 在线观看一区二区三区| 亚洲欧美清纯卡通| 麻豆成人av视频| 成人国产麻豆网| 日本一二三区视频观看| 久久久久免费精品人妻一区二区| 美女国产视频在线观看| 午夜激情欧美在线| 亚洲av二区三区四区| 国产一区亚洲一区在线观看| 桃色一区二区三区在线观看| 亚洲色图av天堂| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲精品乱久久久久久| 国产亚洲一区二区精品| 九九在线视频观看精品| 国产精品爽爽va在线观看网站| 女人被狂操c到高潮| 伦理电影大哥的女人| 在线观看一区二区三区| 能在线免费观看的黄片| 少妇高潮的动态图| 国产成人免费观看mmmm| 久久精品夜色国产| 黄色日韩在线| 国产视频首页在线观看| 99久久人妻综合| 精品熟女少妇av免费看| 深爱激情五月婷婷| av黄色大香蕉| 亚洲真实伦在线观看| 岛国在线免费视频观看| 久久精品熟女亚洲av麻豆精品 | 久久婷婷人人爽人人干人人爱| 男人的好看免费观看在线视频| 秋霞伦理黄片| 久久这里只有精品中国| 国产黄色视频一区二区在线观看 | 18禁裸乳无遮挡免费网站照片| 亚洲丝袜综合中文字幕| 免费播放大片免费观看视频在线观看 | 国产一级毛片在线| 1000部很黄的大片| 丝袜美腿在线中文| 午夜视频国产福利| 日日撸夜夜添| 日本黄色片子视频| 国产av码专区亚洲av| 国产国拍精品亚洲av在线观看| 成人性生交大片免费视频hd| h日本视频在线播放| 亚洲18禁久久av| 日韩一区二区三区影片| 免费观看在线日韩| 日本爱情动作片www.在线观看| 乱码一卡2卡4卡精品| 午夜福利高清视频| 国产精品久久久久久av不卡| 久久久久久久久中文| 亚洲成人中文字幕在线播放| 国产淫片久久久久久久久| 亚洲,欧美,日韩| 久久久久久伊人网av| 51国产日韩欧美| 黄色配什么色好看| 男女边吃奶边做爰视频| 男人狂女人下面高潮的视频| 国产精品综合久久久久久久免费| 好男人视频免费观看在线| 国产色婷婷99| 男女下面进入的视频免费午夜| 国产精品嫩草影院av在线观看| 日韩成人伦理影院| 国产成人免费观看mmmm| 男人狂女人下面高潮的视频| 欧美一区二区精品小视频在线| 观看美女的网站| 国产又黄又爽又无遮挡在线| 精品久久久久久久久亚洲| 国产精品嫩草影院av在线观看| 麻豆成人av视频| 国产精品国产三级国产av玫瑰| 建设人人有责人人尽责人人享有的 | 91av网一区二区| 大香蕉97超碰在线| 少妇高潮的动态图| 一级二级三级毛片免费看| 国产在线一区二区三区精 | 亚洲va在线va天堂va国产| 日韩欧美在线乱码| 国产男人的电影天堂91| 亚洲美女搞黄在线观看| 天天躁夜夜躁狠狠久久av| 狠狠狠狠99中文字幕| 天堂网av新在线| 哪个播放器可以免费观看大片| 在现免费观看毛片| 国语自产精品视频在线第100页| 国产精品日韩av在线免费观看| 中文字幕制服av| 国产极品精品免费视频能看的| 国产精品精品国产色婷婷| 国产av一区在线观看免费| 日日啪夜夜撸| 欧美色视频一区免费| 日韩,欧美,国产一区二区三区 | 国产午夜精品论理片| 免费电影在线观看免费观看| 欧美成人免费av一区二区三区| 国产高清三级在线| 91aial.com中文字幕在线观看| 亚洲最大成人手机在线| 国产精品.久久久| 麻豆国产97在线/欧美| 蜜桃久久精品国产亚洲av| 亚洲精品自拍成人| 99久国产av精品国产电影| 特级一级黄色大片| 成人特级av手机在线观看| 色尼玛亚洲综合影院| 日本五十路高清| 一本久久精品| 亚洲成色77777| or卡值多少钱| 久久这里有精品视频免费| 婷婷六月久久综合丁香| 麻豆精品久久久久久蜜桃| 国产av码专区亚洲av| 亚洲av一区综合| 国产精品麻豆人妻色哟哟久久 | 嫩草影院新地址| 成人欧美大片| 久久久成人免费电影| 简卡轻食公司| 欧美日韩综合久久久久久| 中国国产av一级| 久久久a久久爽久久v久久| 哪个播放器可以免费观看大片| 国产真实乱freesex| av福利片在线观看| 国产免费视频播放在线视频 | 国产三级中文精品| 好男人在线观看高清免费视频| 国产精品国产高清国产av| 国产亚洲最大av| 亚洲av成人精品一区久久| 97人妻精品一区二区三区麻豆| 国产人妻一区二区三区在| 欧美日韩综合久久久久久| 一个人免费在线观看电影| av专区在线播放| 国产成人免费观看mmmm| 如何舔出高潮| 精品熟女少妇av免费看| 激情 狠狠 欧美| 美女高潮的动态| 岛国毛片在线播放| 亚洲综合色惰| 水蜜桃什么品种好| 欧美人与善性xxx| 久久99热这里只有精品18| 好男人在线观看高清免费视频| 日韩欧美国产在线观看| 成人毛片60女人毛片免费| 男女视频在线观看网站免费| 人妻制服诱惑在线中文字幕| 久久久久久久久大av| 卡戴珊不雅视频在线播放| 精品人妻一区二区三区麻豆| 26uuu在线亚洲综合色| 秋霞在线观看毛片| 1000部很黄的大片| 男人和女人高潮做爰伦理| 色综合色国产| 国产黄片视频在线免费观看| 国产伦一二天堂av在线观看| 国产成人福利小说| 亚洲在线观看片| 日韩制服骚丝袜av| 欧美极品一区二区三区四区| 五月伊人婷婷丁香| 极品教师在线视频| 亚洲国产精品sss在线观看| 在线免费观看不下载黄p国产| 三级国产精品欧美在线观看| 性色avwww在线观看| 国产精品一区二区三区四区免费观看| 久久综合国产亚洲精品| 97超碰精品成人国产| 长腿黑丝高跟| 亚洲三级黄色毛片| 高清av免费在线| 丝袜喷水一区| 国产精品一及| 日韩欧美国产在线观看| 我的女老师完整版在线观看| 国产大屁股一区二区在线视频| 精品久久久久久久久亚洲| 麻豆精品久久久久久蜜桃| 日本免费一区二区三区高清不卡| 国产片特级美女逼逼视频| 欧美区成人在线视频| 亚洲av熟女| 国产免费男女视频| 我要搜黄色片| 一个人看视频在线观看www免费| 久久久久久大精品| 午夜免费男女啪啪视频观看| 一级毛片电影观看 | 日韩精品有码人妻一区| 免费黄网站久久成人精品| 五月伊人婷婷丁香| 欧美另类亚洲清纯唯美| 色噜噜av男人的天堂激情| 亚洲一区高清亚洲精品| 汤姆久久久久久久影院中文字幕 | 亚洲av电影不卡..在线观看| 欧美日本亚洲视频在线播放| 国产69精品久久久久777片| 日韩一本色道免费dvd| 亚洲欧美日韩高清专用| 亚洲国产高清在线一区二区三| 成人无遮挡网站| 精品酒店卫生间| 日日干狠狠操夜夜爽| 麻豆成人av视频| 欧美性感艳星| 啦啦啦韩国在线观看视频| 最新中文字幕久久久久| 精品免费久久久久久久清纯| 亚洲人成网站高清观看| 国产精品不卡视频一区二区| 久久久久久久久久成人| av免费观看日本| 午夜爱爱视频在线播放| 久久久久久久久久成人| 美女国产视频在线观看| 男人的好看免费观看在线视频| 白带黄色成豆腐渣| 高清日韩中文字幕在线| 日本黄色片子视频| 边亲边吃奶的免费视频| 日韩精品有码人妻一区| 国产美女午夜福利| 亚洲欧美日韩东京热| 亚洲欧洲日产国产| 国产成人福利小说| 成人二区视频| 精品人妻偷拍中文字幕| 大香蕉97超碰在线| 欧美日韩综合久久久久久| 床上黄色一级片| 亚洲美女搞黄在线观看| 中文天堂在线官网| 国产精品久久视频播放| 免费看光身美女| 亚洲av电影在线观看一区二区三区 | 最近最新中文字幕大全电影3| 人人妻人人澡欧美一区二区| 男人和女人高潮做爰伦理| 久久精品国产亚洲av涩爱| 99在线人妻在线中文字幕| 老女人水多毛片| 精华霜和精华液先用哪个| 精品人妻熟女av久视频| 国产精品.久久久| 97在线视频观看| 国产精品久久久久久久久免| 亚洲av成人精品一二三区| 亚洲精品乱码久久久v下载方式| 久久婷婷人人爽人人干人人爱| 欧美成人午夜免费资源| 国产v大片淫在线免费观看| 国产三级在线视频| 欧美成人精品欧美一级黄| 看非洲黑人一级黄片| 啦啦啦啦在线视频资源| 久久久久精品久久久久真实原创| 老师上课跳d突然被开到最大视频| 丰满人妻一区二区三区视频av| 老女人水多毛片| 国产精品一区二区三区四区久久| 麻豆一二三区av精品| 噜噜噜噜噜久久久久久91| 免费一级毛片在线播放高清视频| 最近视频中文字幕2019在线8| 亚洲成av人片在线播放无| 亚洲精品国产成人久久av| 尾随美女入室| 女人十人毛片免费观看3o分钟| 女人久久www免费人成看片 | 亚洲av中文av极速乱| 建设人人有责人人尽责人人享有的 | 精品一区二区免费观看| 亚洲最大成人手机在线| 国产毛片a区久久久久| 亚洲国产欧美人成| 丝袜喷水一区| av在线观看视频网站免费| 男人和女人高潮做爰伦理| 在线观看66精品国产| 在线观看美女被高潮喷水网站| 精品久久久噜噜| 免费av毛片视频| 免费黄网站久久成人精品| 少妇人妻一区二区三区视频| 国产一级毛片在线| 成人三级黄色视频| 久热久热在线精品观看| 亚洲精品久久久久久婷婷小说 | av线在线观看网站| 2021天堂中文幕一二区在线观| 久久国内精品自在自线图片| av在线亚洲专区| 麻豆国产97在线/欧美| 欧美成人午夜免费资源| 在线观看av片永久免费下载| 亚洲综合精品二区| 久久鲁丝午夜福利片| 麻豆国产97在线/欧美| 日韩欧美国产在线观看| 26uuu在线亚洲综合色| 免费在线观看成人毛片| 欧美日韩一区二区视频在线观看视频在线 | 亚洲aⅴ乱码一区二区在线播放| 亚洲国产日韩欧美精品在线观看| 99在线人妻在线中文字幕| 国产色爽女视频免费观看| 亚洲欧美成人综合另类久久久 | 国产一区二区在线观看日韩| 国产视频内射| 永久免费av网站大全| 最近手机中文字幕大全| 国语自产精品视频在线第100页| 亚洲av中文av极速乱| av卡一久久| 看免费成人av毛片| 午夜亚洲福利在线播放| 亚洲人成网站在线观看播放| 免费不卡的大黄色大毛片视频在线观看 | 亚洲最大成人手机在线| 亚洲成色77777| 69av精品久久久久久| 欧美人与善性xxx| 99九九线精品视频在线观看视频| 麻豆一二三区av精品| 亚洲av成人av| 成人欧美大片| 亚洲精品国产成人久久av| 尾随美女入室| 嫩草影院新地址| 国产精品.久久久| 亚洲国产最新在线播放| 美女脱内裤让男人舔精品视频| 免费观看的影片在线观看| 亚洲第一区二区三区不卡| 99在线视频只有这里精品首页| 国产乱人视频| 秋霞在线观看毛片| 久久久午夜欧美精品| 免费看av在线观看网站| 婷婷色综合大香蕉| 国产91av在线免费观看| 色噜噜av男人的天堂激情| 欧美3d第一页| 2021天堂中文幕一二区在线观| 欧美丝袜亚洲另类| 亚洲精品日韩在线中文字幕| 青春草国产在线视频| 精华霜和精华液先用哪个| 中文字幕制服av| 精品一区二区免费观看| 乱码一卡2卡4卡精品| 久久人人爽人人爽人人片va| 少妇的逼水好多| 视频中文字幕在线观看| 欧美激情久久久久久爽电影| 嫩草影院入口| 国产熟女欧美一区二区| 亚洲va在线va天堂va国产| 亚洲电影在线观看av| 亚洲av免费高清在线观看| 小说图片视频综合网站| 精品一区二区三区视频在线| 亚洲成人精品中文字幕电影| 国产av一区在线观看免费| 男女边吃奶边做爰视频| 少妇高潮的动态图| 国产精品久久久久久精品电影小说 | 免费观看的影片在线观看| 亚洲精品久久久久久婷婷小说 | 哪个播放器可以免费观看大片| 十八禁国产超污无遮挡网站| 久久久精品欧美日韩精品| 国产精品一区二区在线观看99 | 嘟嘟电影网在线观看| 日韩欧美三级三区| 国产伦在线观看视频一区| 午夜老司机福利剧场| 欧美又色又爽又黄视频| 高清在线视频一区二区三区 | 国产色婷婷99| 日韩成人av中文字幕在线观看| av.在线天堂| 啦啦啦啦在线视频资源| 亚洲av免费在线观看| 国产伦一二天堂av在线观看| 99久久精品国产国产毛片| av播播在线观看一区| 国产综合懂色| 免费观看人在逋| 精品不卡国产一区二区三区| 成人性生交大片免费视频hd| 亚洲精品成人久久久久久| 国产色婷婷99| 亚洲欧美日韩卡通动漫| 久久热精品热| 亚洲精品一区蜜桃| 青春草亚洲视频在线观看| 精品国产露脸久久av麻豆 | 亚洲不卡免费看| 久久久久久久午夜电影| 国语自产精品视频在线第100页| 成人性生交大片免费视频hd| 26uuu在线亚洲综合色| 国产淫语在线视频| 欧美激情在线99| 中文天堂在线官网| 亚洲久久久久久中文字幕| 麻豆国产97在线/欧美| 黄色配什么色好看| 国产美女午夜福利| 久久精品久久久久久久性| 亚洲伊人久久精品综合 | 国产高清视频在线观看网站| 久久鲁丝午夜福利片| 午夜爱爱视频在线播放| 国产伦精品一区二区三区视频9| 亚洲丝袜综合中文字幕| 日本免费a在线| 身体一侧抽搐| 午夜福利成人在线免费观看| 亚洲内射少妇av| 国语自产精品视频在线第100页| 非洲黑人性xxxx精品又粗又长| 国产伦精品一区二区三区视频9| 日韩欧美在线乱码| 日本一二三区视频观看| 亚洲精品成人久久久久久| av专区在线播放| 亚洲婷婷狠狠爱综合网| 自拍偷自拍亚洲精品老妇| 国产精品麻豆人妻色哟哟久久 | 婷婷色av中文字幕| 亚洲国产欧洲综合997久久,| 少妇高潮的动态图| 一级黄片播放器| 国产成人免费观看mmmm| 一个人看视频在线观看www免费| av女优亚洲男人天堂| 亚洲怡红院男人天堂| 亚洲无线观看免费| 男人和女人高潮做爰伦理| 日韩精品有码人妻一区| 97超碰精品成人国产| 毛片女人毛片| 久久久久精品久久久久真实原创| kizo精华| 国产精品熟女久久久久浪| 五月伊人婷婷丁香| 亚洲人成网站在线观看播放| 男插女下体视频免费在线播放| 日本午夜av视频| 亚洲av中文字字幕乱码综合| 国语对白做爰xxxⅹ性视频网站| 国产精品无大码| 一级毛片aaaaaa免费看小| 美女内射精品一级片tv| 国产精品无大码| 国产精品野战在线观看| 性色avwww在线观看| 亚洲精品日韩av片在线观看| 国产淫语在线视频| 国模一区二区三区四区视频| 久久国产乱子免费精品| 波多野结衣巨乳人妻| 18禁裸乳无遮挡免费网站照片| 精华霜和精华液先用哪个| 久久久欧美国产精品| 午夜精品在线福利| 午夜a级毛片| 国产成人一区二区在线| 国产老妇伦熟女老妇高清| 精品一区二区免费观看| 日韩国内少妇激情av| 在线免费观看不下载黄p国产|