• <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
    婷婷亚洲欧美| 88av欧美| 午夜福利在线观看吧| 夜夜夜夜夜久久久久| 免费观看精品视频网站| 永久网站在线| 五月玫瑰六月丁香| 亚洲精品久久国产高清桃花| 尤物成人国产欧美一区二区三区| 日本成人三级电影网站| 日韩欧美一区二区三区在线观看| 免费搜索国产男女视频| 女的被弄到高潮叫床怎么办 | 少妇裸体淫交视频免费看高清| 日韩 亚洲 欧美在线| 欧美一级a爱片免费观看看| 又紧又爽又黄一区二区| 国产精品,欧美在线| 夜夜爽天天搞| 一本精品99久久精品77| 黄色日韩在线| 色吧在线观看| 精品欧美国产一区二区三| 中文在线观看免费www的网站| 国产乱人伦免费视频| 精华霜和精华液先用哪个| 九色成人免费人妻av| 久久久久久久精品吃奶| 熟女人妻精品中文字幕| 欧美高清性xxxxhd video| 日韩一区二区视频免费看| 身体一侧抽搐| 91久久精品国产一区二区三区| 在线免费十八禁| 亚洲人成网站高清观看| 老熟妇仑乱视频hdxx| 免费高清视频大片| 午夜激情欧美在线| 内射极品少妇av片p| 日本与韩国留学比较| 最近视频中文字幕2019在线8| 别揉我奶头~嗯~啊~动态视频| 亚洲图色成人| 国产亚洲欧美98| 少妇的逼好多水| 色综合亚洲欧美另类图片| 日本精品一区二区三区蜜桃| 美女高潮喷水抽搐中文字幕| 婷婷六月久久综合丁香| 搡女人真爽免费视频火全软件 | 老熟妇乱子伦视频在线观看| 色综合亚洲欧美另类图片| 精品久久久久久久末码| 国产大屁股一区二区在线视频| 免费av不卡在线播放| 久久亚洲精品不卡| 久久九九热精品免费| 九色国产91popny在线| 国产熟女欧美一区二区| 亚洲熟妇熟女久久| 97热精品久久久久久| 乱系列少妇在线播放| 国产高清三级在线| 亚洲一区二区三区色噜噜| www日本黄色视频网| 国产成年人精品一区二区| 日韩av在线大香蕉| 毛片女人毛片| 国产激情偷乱视频一区二区| 国产三级在线视频| 亚洲黑人精品在线| 蜜桃久久精品国产亚洲av| 国产乱人伦免费视频| 99久久精品热视频| 好男人在线观看高清免费视频| 黄色欧美视频在线观看| 91久久精品电影网| 久久精品综合一区二区三区| 又爽又黄无遮挡网站| 人人妻人人看人人澡| 亚洲中文日韩欧美视频| 春色校园在线视频观看| 天美传媒精品一区二区| 热99在线观看视频| 国产精品人妻久久久影院| 国产男人的电影天堂91| 我要看日韩黄色一级片| 国产精品久久视频播放| 人妻少妇偷人精品九色| 成人亚洲精品av一区二区| 欧美精品啪啪一区二区三区| 亚洲av.av天堂| 亚洲四区av| 免费人成在线观看视频色| ponron亚洲| 九九热线精品视视频播放| x7x7x7水蜜桃| 成人国产综合亚洲| 波多野结衣高清作品| 99在线视频只有这里精品首页| 国产av麻豆久久久久久久| 久久欧美精品欧美久久欧美| 免费大片18禁| 精品久久久久久久久久久久久| 国产亚洲精品久久久久久毛片| 在线免费观看的www视频| 久久久色成人| 成人欧美大片| 国产真实伦视频高清在线观看 | 中国美白少妇内射xxxbb| 淫妇啪啪啪对白视频| 亚洲人与动物交配视频| 亚洲欧美日韩高清在线视频| 欧美+亚洲+日韩+国产| 听说在线观看完整版免费高清| 麻豆精品久久久久久蜜桃| 搞女人的毛片| 久久99热这里只有精品18| 99riav亚洲国产免费| 精品一区二区三区av网在线观看| 亚洲真实伦在线观看| 观看美女的网站| 中文字幕免费在线视频6| 亚洲精品一卡2卡三卡4卡5卡| 国产精品久久久久久精品电影| 亚洲欧美清纯卡通| 联通29元200g的流量卡| 波多野结衣高清无吗| 日韩欧美免费精品| 亚洲天堂国产精品一区在线| 亚洲真实伦在线观看| 91在线精品国自产拍蜜月| 欧美高清性xxxxhd video| 中文亚洲av片在线观看爽| 我的女老师完整版在线观看| 欧美成人a在线观看| 国产久久久一区二区三区| 成人永久免费在线观看视频| 全区人妻精品视频| 日韩欧美精品v在线| 麻豆国产97在线/欧美| 男人和女人高潮做爰伦理| av在线老鸭窝| 日韩欧美精品v在线| 亚洲精华国产精华精| 99热这里只有是精品50| 少妇人妻一区二区三区视频| 国产亚洲av嫩草精品影院| 可以在线观看毛片的网站| 精品人妻1区二区| 欧美成人性av电影在线观看| 大又大粗又爽又黄少妇毛片口| 久久久午夜欧美精品| 亚洲七黄色美女视频| 又爽又黄a免费视频| 91久久精品电影网| 亚洲av一区综合| 亚洲精品乱码久久久v下载方式| 亚洲国产高清在线一区二区三| 国产一区二区激情短视频| 精品一区二区三区视频在线| 18+在线观看网站| 婷婷色综合大香蕉| 99久久精品国产国产毛片| 赤兔流量卡办理| 亚洲精品日韩av片在线观看| 性色avwww在线观看| 精品国内亚洲2022精品成人| 国产高潮美女av| 日韩欧美国产一区二区入口| 色尼玛亚洲综合影院| 少妇的逼水好多| 午夜精品在线福利| 中文字幕免费在线视频6| 成人一区二区视频在线观看| 亚洲色图av天堂| 此物有八面人人有两片| 国产成年人精品一区二区| 男女视频在线观看网站免费| 可以在线观看的亚洲视频| 亚洲无线观看免费| 久久人人爽人人爽人人片va| 亚洲avbb在线观看| 亚洲一区二区三区色噜噜| 免费看美女性在线毛片视频| 黄色欧美视频在线观看| 少妇人妻一区二区三区视频| 精品人妻熟女av久视频| 九九在线视频观看精品| 狠狠狠狠99中文字幕| 国产高清三级在线| 成人三级黄色视频| 1024手机看黄色片| 大又大粗又爽又黄少妇毛片口| 日韩欧美 国产精品| 国产亚洲91精品色在线| 嫁个100分男人电影在线观看| 亚洲av不卡在线观看| 午夜免费男女啪啪视频观看 | 免费无遮挡裸体视频| 男女下面进入的视频免费午夜| 久久天躁狠狠躁夜夜2o2o| 真人做人爱边吃奶动态| 热99在线观看视频| xxxwww97欧美| 老熟妇仑乱视频hdxx| 久久6这里有精品| 午夜久久久久精精品| 精品一区二区三区视频在线观看免费| 久久精品91蜜桃| 99久久精品一区二区三区| 精品一区二区三区视频在线观看免费| 丰满的人妻完整版| 免费观看人在逋| 又紧又爽又黄一区二区| 99久久中文字幕三级久久日本| 18禁裸乳无遮挡免费网站照片| 亚洲人成网站在线播| 国产精品亚洲一级av第二区| 国产亚洲91精品色在线| x7x7x7水蜜桃| 日日摸夜夜添夜夜添小说| 久久久久久九九精品二区国产| 天堂√8在线中文| 亚洲国产色片| 蜜桃久久精品国产亚洲av| 亚洲成人久久性| 精品免费久久久久久久清纯| 人妻久久中文字幕网| 国产白丝娇喘喷水9色精品| 亚洲美女视频黄频| 啪啪无遮挡十八禁网站| 99久久成人亚洲精品观看| 午夜福利在线观看免费完整高清在 | 亚洲av美国av| 久9热在线精品视频| 波多野结衣巨乳人妻| 人人妻人人澡欧美一区二区| 嫩草影院精品99| 欧美3d第一页| 久久久久久久久久成人| 午夜福利18| 日韩欧美免费精品| 中国美女看黄片| 精品日产1卡2卡| 国产精品98久久久久久宅男小说| 亚洲综合色惰| 欧美色视频一区免费| 又黄又爽又刺激的免费视频.| 精品不卡国产一区二区三区| 欧美日韩精品成人综合77777| 欧美日韩乱码在线| 成人鲁丝片一二三区免费| 国产精品亚洲美女久久久| 少妇人妻一区二区三区视频| 一个人免费在线观看电影| 欧美极品一区二区三区四区| 久久欧美精品欧美久久欧美| 给我免费播放毛片高清在线观看| 精品不卡国产一区二区三区| 嫩草影院入口| 成年版毛片免费区| 最近中文字幕高清免费大全6 | 成人国产麻豆网| 欧美xxxx性猛交bbbb| 亚洲中文字幕日韩| 亚洲精品久久国产高清桃花| 国产黄片美女视频| 两个人视频免费观看高清| 人妻丰满熟妇av一区二区三区| 99国产极品粉嫩在线观看| 亚洲内射少妇av| 久久精品人妻少妇| 中文字幕熟女人妻在线| 人妻少妇偷人精品九色| 搡老熟女国产l中国老女人| 大型黄色视频在线免费观看| 琪琪午夜伦伦电影理论片6080| 99在线视频只有这里精品首页| 精品久久久噜噜| 午夜福利高清视频| 舔av片在线| 99久国产av精品| 在线观看66精品国产| 久久精品影院6| 91麻豆av在线| 国产真实伦视频高清在线观看 | 两个人的视频大全免费| 在现免费观看毛片| 国产一级毛片七仙女欲春2| 嫩草影院入口| 亚洲经典国产精华液单| 69av精品久久久久久| 搡老熟女国产l中国老女人| 国产一区二区在线av高清观看| 真实男女啪啪啪动态图| 国产 一区精品| 狂野欧美白嫩少妇大欣赏| 欧美日韩黄片免| 最近最新免费中文字幕在线| 网址你懂的国产日韩在线| 少妇裸体淫交视频免费看高清| 精品久久久久久久人妻蜜臀av| 亚洲欧美日韩无卡精品| 国产伦人伦偷精品视频| 国内揄拍国产精品人妻在线| av在线天堂中文字幕| 日本黄色片子视频| 国产精品国产高清国产av| 天堂动漫精品| 观看免费一级毛片| 国产精品一及| 九色国产91popny在线| 国产伦精品一区二区三区视频9| 久久久国产成人精品二区| 嫩草影视91久久| 精品一区二区三区视频在线| 一区二区三区高清视频在线| 最近在线观看免费完整版| 熟女电影av网| 男女视频在线观看网站免费| 国产精品1区2区在线观看.| 亚洲av熟女| 国产伦一二天堂av在线观看| 在线观看美女被高潮喷水网站| 日韩欧美在线乱码| 亚洲性夜色夜夜综合| 真人做人爱边吃奶动态| 露出奶头的视频| a级一级毛片免费在线观看| 亚洲18禁久久av| 最近最新中文字幕大全电影3| 久久精品国产亚洲av香蕉五月| 欧美极品一区二区三区四区| 久久精品久久久久久噜噜老黄 | 日日啪夜夜撸| aaaaa片日本免费| 又粗又爽又猛毛片免费看| 午夜爱爱视频在线播放| 久久精品人妻少妇| 国产日本99.免费观看| 男人舔奶头视频| 能在线免费观看的黄片| 日本免费一区二区三区高清不卡| 久99久视频精品免费| 亚洲成人久久性| 国产欧美日韩一区二区精品| 91久久精品电影网| 午夜激情欧美在线| 自拍偷自拍亚洲精品老妇| 亚洲人成伊人成综合网2020| av在线观看视频网站免费| 国内精品久久久久久久电影| 午夜福利18| 啦啦啦啦在线视频资源| 亚洲成人精品中文字幕电影| 国产欧美日韩一区二区精品| 日日摸夜夜添夜夜添av毛片 | 在线观看66精品国产| 久久精品夜夜夜夜夜久久蜜豆| 男人的好看免费观看在线视频| 天堂av国产一区二区熟女人妻| 一区福利在线观看| 午夜精品久久久久久毛片777| 久久久久国产精品人妻aⅴ院| 国产精品一区二区三区四区久久| 亚洲精品456在线播放app | 乱系列少妇在线播放| 成人av一区二区三区在线看| 亚洲人与动物交配视频| 国产男靠女视频免费网站| 在线免费观看不下载黄p国产 | 国内久久婷婷六月综合欲色啪| 国产伦精品一区二区三区视频9| 日本三级黄在线观看| 国产人妻一区二区三区在| 中文亚洲av片在线观看爽| 欧美xxxx黑人xx丫x性爽| 国产高清不卡午夜福利| 久久九九热精品免费| 不卡一级毛片| 国产成人av教育| 午夜福利在线在线| 欧美中文日本在线观看视频| 男女下面进入的视频免费午夜| 亚洲男人的天堂狠狠| 草草在线视频免费看| 久久久久精品国产欧美久久久| 亚洲av.av天堂| 日本一本二区三区精品| 国产一区二区激情短视频| 久久久久久久久中文| 最近视频中文字幕2019在线8| netflix在线观看网站| 国产精品自产拍在线观看55亚洲| 熟妇人妻久久中文字幕3abv| 一个人观看的视频www高清免费观看| 一个人看的www免费观看视频| 午夜福利在线观看吧| 在线天堂最新版资源| 日本欧美国产在线视频| 嫩草影视91久久| 久久99热6这里只有精品| 精品人妻视频免费看| 午夜福利成人在线免费观看| 内射极品少妇av片p| 久久6这里有精品| 亚洲精品一卡2卡三卡4卡5卡| 性色avwww在线观看| 午夜福利在线观看免费完整高清在 | 蜜桃亚洲精品一区二区三区| 少妇人妻一区二区三区视频| 国内久久婷婷六月综合欲色啪| 美女黄网站色视频| 天堂av国产一区二区熟女人妻| 舔av片在线| 久久精品国产亚洲av涩爱 | 舔av片在线| 搞女人的毛片| 一区二区三区免费毛片| 两个人的视频大全免费| 两人在一起打扑克的视频| 久久国产乱子免费精品| 黄色丝袜av网址大全| 久久久久久久久久黄片| 波多野结衣巨乳人妻| 国内毛片毛片毛片毛片毛片| 亚洲经典国产精华液单| 精品欧美国产一区二区三| 男女啪啪激烈高潮av片| 久久精品国产99精品国产亚洲性色| 国产一区二区在线观看日韩| 久久人妻av系列| 亚洲熟妇熟女久久| 人人妻人人看人人澡| 亚洲性夜色夜夜综合| 欧美日韩瑟瑟在线播放| 高清毛片免费观看视频网站| 午夜福利高清视频| 蜜桃亚洲精品一区二区三区| 亚洲av.av天堂| 免费人成在线观看视频色| 我要看日韩黄色一级片| 国产白丝娇喘喷水9色精品| 桃红色精品国产亚洲av| 嫩草影视91久久| 亚洲精品456在线播放app | 真人做人爱边吃奶动态| 嫩草影院新地址| av福利片在线观看| 真人一进一出gif抽搐免费| 一个人看的www免费观看视频| 亚洲美女视频黄频| 国产精品久久久久久久久免| 亚洲专区国产一区二区| 亚洲av第一区精品v没综合| netflix在线观看网站| 国产黄a三级三级三级人| 内地一区二区视频在线| 美女免费视频网站| 免费看av在线观看网站| 日韩欧美 国产精品| 精品欧美国产一区二区三| 亚洲国产高清在线一区二区三| 岛国在线免费视频观看| 久久久久免费精品人妻一区二区| 有码 亚洲区| 国产探花在线观看一区二区| 波多野结衣高清作品| 又爽又黄a免费视频| 国内揄拍国产精品人妻在线| 99久久成人亚洲精品观看| 变态另类丝袜制服| 欧美绝顶高潮抽搐喷水| 嫁个100分男人电影在线观看| 中文字幕av成人在线电影| 久久草成人影院| 国产91精品成人一区二区三区| 成人欧美大片| 亚洲av熟女| 精品一区二区三区视频在线观看免费| 99热这里只有是精品在线观看| 午夜a级毛片| 啦啦啦观看免费观看视频高清| 99九九线精品视频在线观看视频| 国产高清不卡午夜福利| 少妇人妻精品综合一区二区 | 看免费成人av毛片| 在现免费观看毛片| 少妇人妻一区二区三区视频| av在线老鸭窝| 免费大片18禁| 夜夜看夜夜爽夜夜摸| x7x7x7水蜜桃| www.色视频.com| 天堂av国产一区二区熟女人妻| 狠狠狠狠99中文字幕| 亚洲电影在线观看av| 麻豆精品久久久久久蜜桃| 婷婷六月久久综合丁香| a级一级毛片免费在线观看| 日日干狠狠操夜夜爽| 午夜福利在线观看免费完整高清在 | 成人特级av手机在线观看| 成人高潮视频无遮挡免费网站| 国产亚洲精品久久久com| 2021天堂中文幕一二区在线观| 国产av一区在线观看免费| 色综合站精品国产| 一个人看的www免费观看视频| 成人av在线播放网站| 搡老岳熟女国产| 成人美女网站在线观看视频| 观看美女的网站| 久久精品91蜜桃| 日本黄色片子视频| 3wmmmm亚洲av在线观看| 欧美另类亚洲清纯唯美| 欧美日本亚洲视频在线播放| 女人被狂操c到高潮| 久久国产乱子免费精品| 如何舔出高潮| 久久这里只有精品中国| 欧美性猛交╳xxx乱大交人| 天堂网av新在线| 久久精品综合一区二区三区| 亚洲性久久影院| 国产成人福利小说| 日日夜夜操网爽| 变态另类丝袜制服| 亚洲中文日韩欧美视频| 97超视频在线观看视频| 免费一级毛片在线播放高清视频| 女生性感内裤真人,穿戴方法视频| 国产激情偷乱视频一区二区| 综合色av麻豆| 又黄又爽又刺激的免费视频.| 国产亚洲精品久久久久久毛片| 1024手机看黄色片| 日本五十路高清| 精品午夜福利视频在线观看一区| 日日撸夜夜添| 日日摸夜夜添夜夜添小说| 免费大片18禁| 三级男女做爰猛烈吃奶摸视频| 久久久国产成人免费| 村上凉子中文字幕在线| 精品一区二区三区视频在线| 中文字幕av成人在线电影| 女人被狂操c到高潮| 男女做爰动态图高潮gif福利片| 天堂影院成人在线观看| 国产成人福利小说| 中文资源天堂在线| 国产精品一及| 中文字幕av在线有码专区| 哪里可以看免费的av片| 熟妇人妻久久中文字幕3abv| 色精品久久人妻99蜜桃| 国产免费一级a男人的天堂| 干丝袜人妻中文字幕| 久久婷婷人人爽人人干人人爱| 国产综合懂色| 久久精品影院6| 内射极品少妇av片p| 免费搜索国产男女视频| 搞女人的毛片| 春色校园在线视频观看| 欧美+亚洲+日韩+国产| 99久国产av精品| 国产精品,欧美在线| 观看免费一级毛片| 我要看日韩黄色一级片| 国产在线精品亚洲第一网站| 美女黄网站色视频| 免费av观看视频| 亚洲黑人精品在线| 精品一区二区三区人妻视频| 欧美xxxx黑人xx丫x性爽| 人妻夜夜爽99麻豆av| 91久久精品国产一区二区成人| av.在线天堂| 在线免费观看不下载黄p国产 | 久久6这里有精品| 麻豆成人av在线观看| 午夜免费激情av| 日韩精品青青久久久久久| 蜜桃亚洲精品一区二区三区| 丰满人妻一区二区三区视频av| 午夜视频国产福利| 色在线成人网| 欧美成人免费av一区二区三区| 国产蜜桃级精品一区二区三区| 精品久久久噜噜| 舔av片在线| 在线a可以看的网站| 国产真实乱freesex| 22中文网久久字幕| 国产单亲对白刺激| 99久久中文字幕三级久久日本| 麻豆一二三区av精品| 国产伦人伦偷精品视频| av.在线天堂| 久久99热6这里只有精品| 国产主播在线观看一区二区| 国产色婷婷99| 在线观看美女被高潮喷水网站| 99视频精品全部免费 在线| 国产淫片久久久久久久久| 午夜福利18| 精品日产1卡2卡| 美女高潮喷水抽搐中文字幕| 在线国产一区二区在线| 两个人视频免费观看高清| 欧美绝顶高潮抽搐喷水|