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

    酸法地浸采鈾多井系統(tǒng)中滲透系數(shù)時(shí)空演化模擬

    2022-07-07 03:10:42陳經(jīng)明周澤超陳茜茜李尋羅躍
    有色金屬科學(xué)與工程 2022年3期
    關(guān)鍵詞:注液浸液滲透系數(shù)

    陳經(jīng)明, 周澤超, 陳茜茜, 李尋, 羅躍

    (1.廣東華南環(huán)保產(chǎn)業(yè)技術(shù)研究院有限公司,廣州 511466;2.核工業(yè)二九〇研究所,廣東 韶關(guān) 512029;3.江西省勘察設(shè)計(jì)研究院,南昌 330095;4.東華理工大學(xué),南昌,330013)

    開(kāi)發(fā)鈾資源,需要有高效采鈾工藝技術(shù)。經(jīng)過(guò)多年實(shí)踐證明,原地浸出采鈾(簡(jiǎn)稱地浸采鈾)技術(shù)具有顯著優(yōu)勢(shì),該技術(shù)不需井巷工程、不需要建立礦石、廢石場(chǎng),生產(chǎn)成本低、對(duì)地表生態(tài)環(huán)境破壞小、能適應(yīng)于低品位鈾礦開(kāi)采,是國(guó)內(nèi)外鈾礦采冶的主流技術(shù)之一[1-5]。由于地浸采鈾技術(shù)優(yōu)勢(shì)明顯,不適合地浸的南方硬巖型鈾礦,已相繼政策性停產(chǎn)。新疆、內(nèi)蒙古等地相繼發(fā)現(xiàn)多個(gè)適合地浸的砂巖型鈾礦后[6],我國(guó)鈾資源采冶已基本采用地浸方式采鈾。因此地浸采鈾已成為我國(guó)提高鈾產(chǎn)量的重要技術(shù),甚至是唯一技術(shù)。

    地浸采鈾技術(shù)基本原理是:首先通過(guò)注液井向含礦含水層中注入溶浸液,讓其與礦石充分接觸,使鈾礦溶解于地下水中,然后通過(guò)抽液井將富含溶解態(tài)鈾的地下水抽出,最后將抽出的地下水(稱為浸出液)輸送至浸出液處理廠提取鈾,尾液將作為溶浸液又注入含水層中,形成地下水抽注循環(huán)(圖1和圖2)。按照溶浸液使用的試劑不同將地浸分為酸法、堿法和中性,其中酸法使用硫酸作為溶浸液(圖2),使用較普遍。

    圖1 地浸采鈾原理示意Fig.1 Principle of in-situ leaching uranium

    圖2 酸法地浸采鈾過(guò)程中地下水抽注循環(huán)示意Fig.2 Schematic diagram of groundwater extraction cycle during acid in-situ leaching uranium

    溶浸范圍內(nèi)抽注井網(wǎng)布設(shè)與鉆井生產(chǎn)能力、礦層連續(xù)性及滲透性等因素有關(guān),目前各個(gè)鈾礦抽注孔數(shù)量主要根據(jù)單孔抽注液量及水文地質(zhì)試驗(yàn)結(jié)果來(lái)確定[7-8]。地浸井場(chǎng)鉆孔布設(shè)類型主要有網(wǎng)格式和行列式2種,網(wǎng)格式主要有五點(diǎn)型(四注一抽)和七點(diǎn)型(六注一抽)2種,行列式抽注井各自成行平行排列[9]。

    目前國(guó)內(nèi)外關(guān)于地浸采鈾模擬的研究較少,大部分學(xué)者進(jìn)行的是現(xiàn)場(chǎng)或者室內(nèi)試驗(yàn),在實(shí)驗(yàn)的基礎(chǔ)上建立模型研究。趙春虎等較早提出了將數(shù)值模擬應(yīng)用于地浸采鈾中,利用數(shù)值模擬軟件Visual MODFLOW模擬研究了我國(guó)某地浸礦山中水動(dòng)力場(chǎng)與溶質(zhì)運(yùn)移過(guò)程[10-11]。周義朋等通過(guò)結(jié)合現(xiàn)場(chǎng)試驗(yàn)利用數(shù)值模擬研究了某砂巖型鈾礦床地浸過(guò)程中NO3-的運(yùn)移情況[12]。田亮等模擬研究了某礦山地浸試驗(yàn)流場(chǎng)SO42-的運(yùn)移規(guī)律[13]。何智等通過(guò)成熟的地下水?dāng)?shù)值模擬軟件GMS模擬研究了某鈾礦床中地下水溶質(zhì)運(yùn)移規(guī)律,同時(shí)預(yù)測(cè)了未來(lái)一段時(shí)間內(nèi)該采區(qū)浸出鈾濃度變化[14]。焦友軍等對(duì)比模擬研究了地下水含水層中六價(jià)鈾吸附反應(yīng)模型[15]。陳茜茜等通過(guò)PHT3D對(duì)酸法地浸過(guò)程中在低含氧情況下鈾水解與遷移的時(shí)空演化規(guī)律進(jìn)行了模擬研究[16]。王兵等利用PHT3D模擬分析了酸法地浸采鈾中不同酸度對(duì)鈾礦浸出效率影響,及雜質(zhì)礦物對(duì)溶浸液的消耗降低鈾礦浸出率[17]。邱文杰等利用TOUGHREACT模擬研究了地浸采鈾過(guò)程中溶浸液與礦層水-巖相互作用及反應(yīng)動(dòng)力學(xué)過(guò)程[18]。

    大部分學(xué)者對(duì)于地下水的模擬研究更多的著重點(diǎn)是地下水滲透系數(shù)的模擬。桂勇等模擬研究了單孔原地浸礦的滲流量計(jì)算模型,結(jié)果表明計(jì)算模型穩(wěn)定有效[19]。張將偉等利用不確定分析構(gòu)建了地下水與地表水耦合模型,確定了滲透系數(shù)等參數(shù)在模型中的靈敏度較高[20]。趙勇勝等模擬研究了低滲透系數(shù)非均質(zhì)含水層中污染物遷移規(guī)律[21]。金磊等通過(guò)開(kāi)發(fā)多孔介質(zhì)滲流耦合模擬平臺(tái),揭示了介質(zhì)滲透性降低的本質(zhì)是滲流與介質(zhì)顆粒骨架耦合作用結(jié)果[22]。馬鵬飛等構(gòu)建的近場(chǎng)多孔介質(zhì)滲透模型,可應(yīng)用于多孔介質(zhì)滲流過(guò)程模擬[23]。

    巴彥烏拉鈾礦是內(nèi)蒙古開(kāi)發(fā)的第一個(gè)酸法地浸砂巖型鈾礦山,自2015年8月正式投產(chǎn)運(yùn)營(yíng)以來(lái),礦層滲透性顯著降低,導(dǎo)致鈾浸出量下降,浸鈾成本增高,資源回收率降低,嚴(yán)重阻礙了該礦床的鈾資源開(kāi)發(fā)。本文以巴彥烏拉鈾礦床礦層堵塞為背景,采用數(shù)值模擬法,計(jì)算地浸采鈾過(guò)程中礦層滲透系數(shù)的時(shí)空演化,預(yù)測(cè)堵塞發(fā)生的時(shí)間和位置,為礦區(qū)開(kāi)采過(guò)程中解決堵塞等問(wèn)題提供科學(xué)依據(jù)。

    1 研究方法

    1.1 模擬軟件

    TOUGH(Trransport of Unsaturated Groundwater and Heat)是由美國(guó)勞倫斯-伯克利國(guó)家實(shí)驗(yàn)室基于MULKOM模擬程序開(kāi)發(fā)的數(shù)值模擬程序;是一個(gè)可以模擬裂隙或孔隙介質(zhì)中一維到三維的多組分、多相流及非等溫的水流和熱運(yùn)移的數(shù)值模擬程序,不同程序及所適用的領(lǐng)域如表1所列。本文所用到的數(shù)值模擬軟件TOUGHREACT是由許天福等基于TOUGH基礎(chǔ)上開(kāi)發(fā)的可應(yīng)用于裂隙或孔隙介質(zhì)中一維到三維的多組分、多相流及非等溫的水流和熱反應(yīng)運(yùn)移的數(shù)值模擬程序[25-26],不同程序及所適用的領(lǐng)域如表2所列。TOUGH系列軟件自開(kāi)發(fā)以來(lái)被應(yīng)用于飽和-非飽和帶地下水科學(xué)、環(huán)境影響評(píng)價(jià)和修復(fù)、地?zé)衢_(kāi)發(fā)工程、CO2地質(zhì)封存、核廢物處置、白云石化、頁(yè)巖氣及石油開(kāi)采等領(lǐng)域[27-34]。

    表1 TOUGH家族代碼的發(fā)展歷史(自施小清(2009)修[24])Table 1 Development history of TOUGH family codes(repaired from Shi(2009)[24])

    表2 TOUGHREACT模塊功能(自XU(2004)修[35])Tab le 2 TOUGHREACT module function(repaired from XU(2004)[35])

    1.2 控制方程

    本文相應(yīng)控制方程均引TOUGHREACT模擬軟件使用說(shuō)明[25]。在TOUGHREACT中的質(zhì)量守恒示意圖如圖3。

    圖3 質(zhì)量(能量)守恒示意Fig.3 Conservation of mass(energy)

    質(zhì)量守恒方程可以寫(xiě)成以下形式:

    式(1)中:Vn為 研究區(qū)內(nèi)的任何積分子區(qū)域;Гn為積分區(qū)域表面,m2;M為單位體積的能量或質(zhì)量;k表示不同的組分,如:水、氣、溶質(zhì)等;F為質(zhì)量(能量)對(duì)流量,m3/s;q表示源匯相;n表示通過(guò)表面積為dΓn的法向量。

    式(2)中:k表示組分;Φ表示孔隙度;Sβ表示β相的飽和度;ρβ表示β相的密度,kg/m3;Xkβ表示在相β中組分k的質(zhì)量分?jǐn)?shù);β表示液相。

    式(5)中:μβ表示β相滲流速度,m/s;k表示絕對(duì)滲透率,m/s;krβ表示相對(duì)滲透率,m/s;μβ表示液體黏滯系數(shù),Pa;g表示重力加速度常量,g/s2。

    式(6)中:Pβ表示液相壓強(qiáng),Pa;P為參考系壓強(qiáng),Pa;Pcβ為毛管壓強(qiáng),Pa。

    所采用的反應(yīng)速率表達(dá)式來(lái)自Lasaga等的研究成果[36]。

    式(7)中:n表示礦物;rn表示反應(yīng)速率,(mol/(L·s-1)),正表示溶解,負(fù)表示沉淀;kn表示隨溫度變化的速率常數(shù)(mol/(m2·s-1));An表示每千克水的礦物比表面積(cm2/g);Ωn表示礦物飽和比率;θ與η由實(shí)驗(yàn)得出,通常取1。

    式(8)中:k25表示25℃的速率常數(shù),mol/(m2·s);Ea為活化能,kJ/mol;R為氣體常數(shù),kJ/(mol·K),通常取8.314 51;T表示絕對(duì)溫度,K。

    滲透率k(Permeability)表征巖層對(duì)不同流體的固有滲透性能,滲透率僅取決于巖石的空隙性質(zhì),與滲流的液體性質(zhì)無(wú)關(guān)。滲透率與滲透系數(shù)K的關(guān)系為[37]:

    式(9)中:ρ為液體密度,kg/m3;g為重力加速度,m/s2;μ為液體動(dòng)力黏滯系數(shù),N·s/m2;k為滲透率,常用單位為μm2(即達(dá)西)或cm2。

    在進(jìn)行多場(chǎng)耦合建模時(shí)需要給定不同物理參數(shù)之間的關(guān)系,如孔隙度-滲透率之間的關(guān)系[38]。孔隙度-滲透率之間復(fù)雜的關(guān)系取決于孔徑分布、孔隙形狀及孔隙間的連通性等因素的相互作用,無(wú)法僅僅通過(guò)簡(jiǎn)單的立方根定理和Konzeny-Carman孔隙度-滲透率方程反映[39]。有室內(nèi)實(shí)驗(yàn)表明,產(chǎn)生適量的沉淀物會(huì)導(dǎo)致孔隙度減少且滲透系數(shù)會(huì)大幅降低[40]。現(xiàn)場(chǎng)實(shí)驗(yàn)表明介質(zhì)滲透率對(duì)于孔喉變化較為敏感[41],這與在天然情況下滲流通道的收斂-發(fā)散性一致,由于孔喉被沉淀物堵塞,導(dǎo)致介質(zhì)中存在不連續(xù)的孔隙空間[39]。滲透率降低不僅取決于總孔隙度的減少,同時(shí)也取決于孔隙空間幾何形狀和沉淀物的空間分布有關(guān)[42]。模型中采用VERNA等改進(jìn)的孔隙度-滲透率關(guān)系式[39]。

    Verman—Pruess模型:

    式(10)中:ki和φi分別表示初始條件的滲透率和孔隙度;φc表示臨界孔隙度,當(dāng)孔隙度達(dá)到這個(gè)值時(shí)滲透率為0;n表示冪指數(shù);φc和n的取值由介質(zhì)的性質(zhì)決定。

    根據(jù)式(9)、式(10)可得滲透系數(shù)與孔隙度之間的關(guān)系,如式(11):

    本文中孔隙度初始值取0.2;前人的研究認(rèn)為當(dāng)多孔介質(zhì)中孔隙度約20%時(shí),滲透系數(shù)會(huì)降低幾個(gè)數(shù)量級(jí)[43],因此φc臨界孔隙度取初始值的80%即0.16;滲透率為8.03×10-12m2;由此可確定模型中滲透系數(shù)(K,m/d)與孔隙度的關(guān)系式為:

    滲透系數(shù)與孔隙度關(guān)系曲線見(jiàn)圖4。

    圖4 滲透系數(shù)—孔隙度關(guān)系Fig.4 Relationship between permeability coefficient and porosity

    2 數(shù)值模擬

    2.1 模型概化

    基于內(nèi)蒙古某砂巖型鈾礦的多個(gè)抽注系統(tǒng),模擬酸法地浸采鈾過(guò)程中出現(xiàn)滲透系數(shù)減小的現(xiàn)象。該場(chǎng)地的地下水動(dòng)力場(chǎng)以及地浸采鈾的堵塞物質(zhì)已有詳細(xì)研究[44-45]。模型共有6個(gè)注液孔(Z1-Z6),2個(gè)抽液孔(C1、C2)(圖5),構(gòu)成兩組五點(diǎn)型二維地浸抽注系統(tǒng);在確定模擬區(qū)域時(shí)分別模擬了邊長(zhǎng)為2 000、1 000、500、200 m正方形區(qū)域,模擬結(jié)果顯示溶浸液的滲流范圍較?。▓D6模擬結(jié)果流場(chǎng)圖),為了減少模擬時(shí)所需的計(jì)算量最終確定將邊界大小為200 m較為合適。

    圖5 二維模型概念示意Fig.5 Schematic diagram of the 2-D model

    圖6 二維流場(chǎng)矢量示意Fig.6 Vector diagram of two-bit flow field

    2.2 模型條件

    本研究的模擬條件來(lái)自鄭和秋野等對(duì)巴彥烏拉鈾礦現(xiàn)場(chǎng)地浸的條件試驗(yàn)和地浸場(chǎng)水動(dòng)力模擬。模型邊界大小確定為200 m,含礦含水層厚度為58.3 m,每個(gè)注液孔與抽液孔的距離為30 m,每個(gè)注液孔的注入流量為4 m3/h,抽液孔抽出流量為12 m3/h;含礦層厚度58.3 m,含礦層的初始孔隙度假設(shè)為0.2,初始滲透系數(shù)為7.6 m/d(滲透率為8.03×10-12m2),總模擬時(shí)間為690 d。在地浸采鈾中將尾液中的有用元素(鈾)沉淀之后,再次作為溶浸液注入注液孔,這種工藝也稱為循環(huán)抽注系統(tǒng)。這種工藝容易造成礦層滲透系數(shù)減?。ㄉ踔炼氯?,原因是回灌到礦層中的地浸尾液礦化度較高。通過(guò)對(duì)區(qū)域地浸尾液室內(nèi)試驗(yàn)數(shù)據(jù)的分析得出:研究區(qū)域地浸開(kāi)采100 d左右后,浸出液中的主要離子濃度趨于穩(wěn)定(尤其是Ca2+、Mg2+),因此模型中的邊界條件各類離子濃度采用100 d后的平均值。初始條件與邊界條件的各離子濃度列于表3。根據(jù)已有的資料整理得到模擬區(qū)域雜質(zhì)礦物含量列于表4,用于計(jì)算礦物動(dòng)力學(xué)速率常數(shù)的參數(shù)列于表5。

    表3 模型初始/邊界條件Table 3 Model initial/boundary conditions

    表4 模擬區(qū)域雜質(zhì)礦物含量Tab l e 4 Impurity mineral content in the simulated area

    表5 用于計(jì)算礦物動(dòng)力學(xué)速率常數(shù)的參數(shù)Tab le 5 Parameters for calculating the kinetic rate contents of minerals

    3 模擬結(jié)果與討論

    基于模型概化,選取模擬結(jié)果中各因子變化較明顯時(shí)刻,進(jìn)行討論分析;對(duì)模擬結(jié)果中滲透系數(shù)(K)的時(shí)空演化規(guī)律進(jìn)行討論分析。模擬結(jié)果中各因子變化較明顯的時(shí)間分別為:6.9、41.4、200、283、469、690 d;同時(shí)選取剖面A、B對(duì)滲透系數(shù)的空間演化結(jié)果進(jìn)行討論分析。

    3.1 不同時(shí)刻滲透系數(shù)演化結(jié)果

    在t=6.9 d時(shí),模擬區(qū)域內(nèi)的滲透系數(shù)均大于初始值(7.6 m/d),滲透系數(shù)介于7.6~7.69 m/d之間,變化量較小,相對(duì)初始值最大僅增大1.2%。滲透系數(shù)最大值均位于注液孔附近,與注液孔距離約4 m。在抽液孔附近的滲透系數(shù)幾乎不發(fā)生變化。從圖7中可以看出滲透系數(shù)發(fā)生變化的區(qū)域也主要集中在注液孔附近,最大變化半徑約為10 m,且變化方向朝著注液孔一側(cè)。

    圖7 不同時(shí)刻滲透系數(shù)演化結(jié)果Fig.7 Evolution results of permeability coefficient at different times

    在溶浸液(H2SO4)反應(yīng)前鋒作用下與注液孔距離越大,滲透系數(shù)越接近初始值,不同注液孔周圍的滲透系數(shù)演化情況相似,沒(méi)有井群干擾現(xiàn)象。

    與上一個(gè)時(shí)刻(6.9 d)對(duì)比,在t=41.1 d時(shí)模擬區(qū)域內(nèi)滲透系數(shù)發(fā)生變化的面積有所增大,且該時(shí)刻開(kāi)始出現(xiàn)滲透系數(shù)小于初始值的網(wǎng)格點(diǎn),以注液孔為圓心滲透系數(shù)發(fā)生變化區(qū)域最大半徑為22.63 m。

    根據(jù)圖7可知,t=41.1 d時(shí)模擬區(qū)域滲透系數(shù)介于6.65~7.75 m/d之間;滲透系數(shù)最大值和最小值較初始值分別增大2.0%和-12.5%。

    說(shuō)明隨著時(shí)間的推移模擬區(qū)域內(nèi)的滲透系數(shù)區(qū)間會(huì)不斷擴(kuò)大,并且主要的變化是滲透系數(shù)最小值會(huì)不斷減小,從圖7中可以看出,在滲透系數(shù)改變的區(qū)域內(nèi),主要表現(xiàn)為滲透系數(shù)不斷減小,僅有少量緊鄰注液孔的位置滲透系數(shù)增大。

    根據(jù)圖7給出的t=200 d模擬區(qū)域滲透系數(shù)演化結(jié)果可知,在該時(shí)刻模擬區(qū)域內(nèi)滲透系數(shù)介于5.99~9.08 m/d之間,相對(duì)初始值最小值減少21.14%,最大值增大19.54%。

    從圖7中可以看出,模擬區(qū)域內(nèi)滲透系數(shù)增大的范圍僅出現(xiàn)在靠近注液孔周圍的極小區(qū)域內(nèi)(半徑約為4 m),沿溶浸液滲流范圍內(nèi)的其他區(qū)域均出現(xiàn)了滲透系數(shù)減小現(xiàn)象。通過(guò)與之前的幾個(gè)時(shí)刻(6.9 d、41.4 d)對(duì)比,說(shuō)明在模擬區(qū)域內(nèi)靠近注液孔處的礦層隨著時(shí)間的延長(zhǎng)不斷與溶浸液發(fā)生溶解反應(yīng),且在滲流作用和反應(yīng)動(dòng)力學(xué)共同作用下,在滲流路徑中不斷會(huì)有沉淀物生成,同時(shí)也有溶解反應(yīng)發(fā)生。

    根據(jù)模擬結(jié)果顯示,該時(shí)刻滲透系數(shù)小于初始值的大部分礦層內(nèi)滲透系數(shù)均介于5.99~6.77 m/d之間,少部分區(qū)域滲透系數(shù)介于7.15~7.45 m/d之間。

    從模擬結(jié)果可知在抽液孔附近(半徑小于2 m)滲透系數(shù)已經(jīng)開(kāi)始發(fā)生變化。說(shuō)明模擬到200 d時(shí),溶浸液與雜質(zhì)礦物發(fā)生的沉淀-溶解反應(yīng)已經(jīng)開(kāi)始影響抽液周圍的滲透系數(shù)。且與注液孔距離越遠(yuǎn)滲透系數(shù)減少值則越小,說(shuō)明在滲流通道上沉淀-溶解量不斷發(fā)生變化;與注液孔距離近的區(qū)域沉淀物的量較多,溶解量較少;靠近抽液孔處沉淀量相對(duì)減少,溶解量也相對(duì)增加。

    圖7給出了t=283 d模擬區(qū)域滲透系數(shù)演化結(jié)果,t=283d時(shí)模擬區(qū)域內(nèi)滲透系數(shù)區(qū)間為5.88~9.44 m/d。最大值相對(duì)初始值增大了30.34%,最小值相對(duì)初始值減小了22.61%。對(duì)比前文中討論的幾個(gè)時(shí)刻,該時(shí)刻最大的區(qū)別是模擬區(qū)域內(nèi)滲透系數(shù)的最大值變化量比最小值變化量大。

    滲透系數(shù)最大值(9.94 m/d)均出現(xiàn)在注液孔附近,且滲透系數(shù)發(fā)生變化的范圍幾乎沒(méi)有變化;從模擬結(jié)果圖可以看出該區(qū)域內(nèi)滲透系數(shù)均大于7.9 m/d。說(shuō)明在注液孔附近滲透系數(shù)增大的區(qū)域內(nèi),不斷發(fā)生溶解反應(yīng),產(chǎn)生的離子隨滲流通斷不斷運(yùn)移;在該區(qū)域內(nèi)礦物溶解量遠(yuǎn)大于沉淀量。

    從圖7可知該時(shí)刻滲透系數(shù)發(fā)生改變的區(qū)域內(nèi)只有少部分區(qū)域滲透系數(shù)區(qū)間為[7.4 m/d,7.9 m/d],說(shuō)明在模擬達(dá)到283 d時(shí),礦層內(nèi)發(fā)生沉淀和溶解的區(qū)域分界線較明顯。

    在滲透系數(shù)發(fā)生改變的區(qū)域中滲透系數(shù)介于區(qū)間為[5.88 m/d,6.9 m/d]的區(qū)域占比較大。說(shuō)明在這些區(qū)域內(nèi),反應(yīng)生成的沉淀物量大于溶解量;且該區(qū)域內(nèi)大部分位置滲透系數(shù)均小于6.39 m/d,只有少量靠近抽液孔處區(qū)域的滲透系數(shù)為6.9 m/d。

    隨著模擬時(shí)間不斷增加,滲流場(chǎng)-化學(xué)場(chǎng)耦合作用對(duì)模擬區(qū)域內(nèi)的滲透系數(shù)以及其他因子影響越明顯。

    圖7給出496 d時(shí)的滲透系數(shù)演化結(jié)果,該時(shí)刻礦層內(nèi)滲透系數(shù)最大值為12.01 m/d,相對(duì)初始值增大了57.99%;滲透系數(shù)最小值為5.65 m/d,相對(duì)初始值減小了25.71%。

    在前文中討論到283 d時(shí)模擬區(qū)域滲透系數(shù)在兩組五點(diǎn)型抽注系統(tǒng)的影響范圍有公共區(qū)域,從模擬結(jié)果看出,該時(shí)刻(496 d)滲透系數(shù)在模擬區(qū)域內(nèi)分布較均一。在滲透系數(shù)減小的區(qū)域內(nèi),不同抽注液孔之間幾乎沒(méi)有差異;兩組抽注系統(tǒng)的滲透系數(shù)空間分布關(guān)于直線X=100對(duì)稱,這是之前幾個(gè)時(shí)刻沒(méi)有出現(xiàn)的情況;說(shuō)明經(jīng)過(guò)496 d的反應(yīng)運(yùn)移,滲流場(chǎng)-化學(xué)場(chǎng)相互影響已經(jīng)開(kāi)始趨于穩(wěn)定狀態(tài),至少模擬區(qū)域內(nèi)的沉淀物空間分布差異性開(kāi)始減小。

    對(duì)于區(qū)域內(nèi)滲透系數(shù)增大和減小的區(qū)域,與前文的分析相似。區(qū)別在于該時(shí)刻礦層內(nèi)滲透系數(shù)減小的范圍有所增大,這是由于不斷注入的溶浸液與伴生礦物反應(yīng)增加了沉淀物的量,同時(shí)也由于反應(yīng)動(dòng)力學(xué)的作用,沉淀物會(huì)在滲流路徑中生成,同時(shí)可以解釋,礦層內(nèi)滲透系數(shù)小于初始值的區(qū)域內(nèi)滲透系數(shù)差異性較小。

    根據(jù)圖7給出的模擬結(jié)束時(shí)(690 d)的滲透系數(shù)演化結(jié)果,該時(shí)刻礦層內(nèi)滲透系數(shù)最大值為14.68 m/d,相對(duì)初始值增大了93.18%;滲透系數(shù)最小值為5.37 m/d,相對(duì)初始值減小了29.31%。

    隨著時(shí)間的推移,礦層內(nèi)滲透系數(shù)改變的區(qū)域不斷增大。從模擬結(jié)果可以看出,模擬區(qū)域內(nèi)滲透系數(shù)發(fā)生改變的區(qū)域主要集中在抽注液孔之間,在注液孔外側(cè)也有部分區(qū)域滲透系數(shù)減小。從圖7可知,注液孔Z1、Z3、Z4、Z6演化情況一致,Z2、Z5演化情況一致。

    滲透系數(shù)區(qū)間值不斷增大,模擬結(jié)束時(shí)(690 d)滲透系數(shù)最大值(14.68 m/d)是最小值(5.37 m/d)的2.73倍,且滲透系數(shù)增大的區(qū)域除了集中在注液孔附近,在外側(cè)也有少量位置滲透系數(shù)所有增大;滲透系數(shù)減小的區(qū)域不斷增大。由于不斷注入的溶浸液與雜質(zhì)礦物發(fā)生反應(yīng),在注液孔處大量的雜質(zhì)礦物發(fā)生溶解反應(yīng),但是在反應(yīng)動(dòng)力學(xué)的作用下沉淀反應(yīng)不會(huì)馬上發(fā)生,因此沉淀物會(huì)分布在滲流通道上。對(duì)于滲透系數(shù)減小的區(qū)域內(nèi),也會(huì)有溶解反應(yīng)生成,但是模擬結(jié)果顯示生成的沉淀物量大于溶解量。

    對(duì)比前文中討論的幾個(gè)時(shí)刻,模擬結(jié)束時(shí)最大的特點(diǎn)是在滲透系數(shù)減少的區(qū)域外側(cè),有少量區(qū)域的滲透系數(shù)相對(duì)初始值有所增大(8.86 m/d)。分析可能的原因是這些區(qū)域也是滲流場(chǎng)影響的范圍內(nèi),運(yùn)移到這些網(wǎng)格的溶浸液內(nèi)還有少量的H+未反應(yīng),與礦層發(fā)生反應(yīng)導(dǎo)致滲透系數(shù)增大。

    圖7給出了不同時(shí)刻(6.9、41.4、200、690 d)二維模型滲透系數(shù)演化結(jié)果。從演化結(jié)果來(lái)看,滲透系數(shù)發(fā)生變化的區(qū)域隨時(shí)間不斷增大,模擬區(qū)域內(nèi)滲透系數(shù)范圍隨時(shí)間不斷增大。

    在滲流作用和化學(xué)反應(yīng)前鋒不斷向抽液孔推進(jìn),井群干擾現(xiàn)象越明顯,6個(gè)注液孔周圍的滲透系數(shù)演化情況略有差異。Z1、Z3、Z4、Z6 4個(gè)注液孔演化情況相似,從圖中可以看出滲透系數(shù)變化主要方向不斷向C1或C2抽液孔推進(jìn);Z2、Z5 2個(gè)注液孔演化情況相似,從圖7中可以看出滲透系數(shù)變化主要方向不斷同時(shí)向C1和C2抽液孔推進(jìn)。出現(xiàn)這種情況的原因是C1、C2共同影響Z2、Z5注液孔周圍的流場(chǎng)分布,其他4個(gè)注液孔周圍的流場(chǎng)僅受C1或者C2單獨(dú)影響。

    3.2 不同空間滲透系數(shù)演化結(jié)果

    圖8給出了兩個(gè)剖面(A、B)滲透系數(shù)在t=6.9 d、t=41.4 d、t=200 d、t=283 d、t=469 d及t=690 d時(shí)的演化結(jié)果。0 m表示注液孔位置,30 m表示抽液孔位置。

    從圖8中可以看出在A、B兩個(gè)剖面中,從注液孔到抽液孔之間的滲透系數(shù)演化情況相同之處是:滲透系數(shù)隨距離的總體趨勢(shì)均是不斷減?。辉赿=17.23 m到注液孔處兩組抽注孔之間的滲透系數(shù)演化結(jié)果相同。

    圖8 不同空間位置滲透系數(shù)演化結(jié)果Fig.8 Evolution results of permeability coefficients at different spatial locations

    不同之處:A和B兩個(gè)剖面上滲透系數(shù)在該時(shí)刻演化結(jié)果有所差異。A剖面上滲透系數(shù)區(qū)間為[7.60 m/d,7.68 m/d],B剖面上滲透系數(shù)區(qū)間為[7.60 m/d,7.69 m/d]。

    在d=17.23 m之前兩組抽注孔之間的滲透系數(shù)演化雖然是沿滲流路徑不斷減小,但是相同距離下的滲透系數(shù)均不相同;在d=17.23 m之前的滲流路徑內(nèi)大部分(占抽注液孔間距的81.7%)位置中Z2-C1之間的滲透系數(shù)大于Z4-C1之間的滲透系數(shù)。

    根據(jù)圖8給出的兩個(gè)剖面(A、B)滲透系數(shù)在t=41.4 d時(shí)演化結(jié)果。對(duì)比2個(gè)剖面滲透系數(shù)演化結(jié)果,相同點(diǎn)是:兩個(gè)剖面上滲透系數(shù)沿滲流路徑的演化趨勢(shì)相同,均是先減小后增大;兩個(gè)剖面上該時(shí)刻滲透系數(shù)最小值出現(xiàn)的位置相同(均是于注液孔相距5.58 m處);兩個(gè)剖面在14.82 m之后的空間位置上滲透系數(shù)差異較小。

    在A、B兩個(gè)剖面上從注液孔到滲透系數(shù)最小值處(5.63 m處)之間區(qū)域內(nèi)也有所差異,A剖面上滲透系數(shù)沿滲流路徑不斷減小,B剖面從注液孔(7.50 m/d)到2.83 m(7.53 m/d)之間有所增大且均小于初始值。

    兩個(gè)剖面上滲透系數(shù)演化的不同點(diǎn):兩個(gè)剖面上滲透系數(shù)最小值有所不同,其中A剖面上滲透系數(shù)最小值為7.02 m/d,B剖面上滲透系數(shù)最小值為6.99 m/d。

    A剖面在靠近注液孔處的滲透系數(shù)大于初始值(7.75 m/d),在B剖面上靠近注液孔處的滲透系數(shù)小于初始值(7.50 m/d)。

    兩個(gè)剖面從注液孔到5.88 m處兩個(gè)剖面上滲透系數(shù)差異性較大,主要體現(xiàn)在B剖面不同位置滲透系數(shù)演化較復(fù)雜,從注液孔處出發(fā)有少量位置的滲透系數(shù)相對(duì)注液孔出有所增大;從5.58 m至抽液孔處滲透系數(shù)演化也有所差異,B剖面大部分空間位置的滲透系數(shù)均小于A剖面相同位置。

    根據(jù)圖8給出的兩個(gè)剖面(A、B)滲透系數(shù)在t=200 d時(shí)演化結(jié)果。A、B兩個(gè)剖面滲透系數(shù)最小值均出現(xiàn)在5.66 m處,該時(shí)刻兩個(gè)剖面上滲透系數(shù)最小值差異較小(A剖面滲透系數(shù)最小值為5.99 m/d,B剖面滲透系數(shù)最小值為6.05 m/d,相差1%);兩個(gè)剖面上滲透系數(shù)演化結(jié)果均出現(xiàn)了大于初始值和小于初始值的空間位置。

    A和B兩個(gè)剖面上滲透系數(shù)在該時(shí)刻演化結(jié)果有所差異。A剖面上滲透系數(shù)區(qū)間為[5.99 m/d,9.08 m/d],B剖面上滲透系數(shù)區(qū)間為[6.05 m/d,8.41 m/d],B剖面上出現(xiàn)了滲透系數(shù)明顯大于初始值的位置,A剖面上滲透系數(shù)區(qū)間明顯大于B剖面上滲透系數(shù)區(qū)間,說(shuō)明B剖面上的滲透系數(shù)較A剖面滲透系數(shù)演化較慢。

    模擬結(jié)果看出在該時(shí)刻A、B剖面大部分位置處,B剖面上的滲透系數(shù)率大于A剖面上的滲透系數(shù)(注液孔除外);14.15~30 m之間的區(qū)域內(nèi),B剖面上滲透系數(shù)值曲線一直在A剖面滲透系數(shù)值曲線上方;在5.63~14.15 m之間存在少部分區(qū)域(約2.81 m)B剖面上滲透系數(shù)值曲線一直在A剖面滲透系數(shù)值曲線下方。

    圖8給出的兩個(gè)剖面(A、B)滲透系數(shù)在t=283 d時(shí)演化結(jié)果。A、B兩個(gè)剖面滲透系數(shù)最小值均出現(xiàn)在5.66 m處,該時(shí)刻兩個(gè)剖面上滲透系數(shù)最小值差異較?。ˋ剖面滲透系數(shù)最小值為5.88 m/d,B剖面滲透系數(shù)最小值為6.94 m/d,相差1%)。

    A和B兩個(gè)剖面上滲透系數(shù)在該時(shí)刻演化結(jié)果有所差異。A剖面上滲透系數(shù)區(qū)間為[5.88 m/d,9.94 m/d],B剖面上滲透系數(shù)區(qū)間為[5.94 m/d,9.17 m/d],B剖面上出現(xiàn)了滲透系數(shù)明顯大于初始值的位置,A剖面上滲透系數(shù)區(qū)間明顯大于B剖面上滲透系數(shù)區(qū)間,說(shuō)明B剖面上的滲透系數(shù)較A剖面滲透系數(shù)演化較慢;在該時(shí)刻,兩個(gè)剖面在滲透系數(shù)較小的區(qū)域內(nèi)(5.66~16 m)之間的演化情況一致。

    兩個(gè)剖面滲透系數(shù)演化的差異性體現(xiàn)在靠近注液孔處,以及16~30 m這兩段區(qū)域內(nèi);在16~30 m之間,B剖面滲透系數(shù)均大于A剖面上相同位置處的滲透系數(shù),注液孔到5.66 m之間兩個(gè)剖面上滲透系數(shù)的差異與前文中討論一致,就不贅述。

    圖8給出的兩個(gè)剖面(A、B)滲透系數(shù)在t=469 d時(shí)演化結(jié)果。A、B兩個(gè)剖面滲透系數(shù)最小值均出現(xiàn)在5.66 m處,該時(shí)刻兩個(gè)剖面上滲透系數(shù)最小值差異較小(A剖面滲透系數(shù)最小值為5.88 m/d,B剖面滲透系數(shù)最小值為6.94 m/d,相差0.8%);在該時(shí)刻,兩個(gè)剖面從滲透系數(shù)最小值到抽液孔處(5.66~30 m)之間的演化情況一致。

    A和B兩個(gè)剖面上滲透系數(shù)在該時(shí)刻演化結(jié)果有所差異。A剖面上滲透系數(shù)區(qū)間為[5.65 m/d,12 m/d],B剖面上滲透系數(shù)區(qū)間為[5.70 m/d,10.94 m/d],B剖面上出現(xiàn)了滲透系數(shù)明顯大于初始值的位置,A剖面上滲透系數(shù)區(qū)間明顯大于B剖面上滲透系數(shù)區(qū)間,說(shuō)明B剖面上的滲透系數(shù)較A剖面滲透系數(shù)演化較慢。

    通過(guò)兩個(gè)具有代表性的剖面滲透系數(shù)演化結(jié)果可知,在該時(shí)刻,模擬區(qū)域內(nèi)大部分空間位置的滲透系數(shù)均小于初始值。

    根據(jù)圖8給出的兩個(gè)剖面(A、B)滲透系數(shù)在t=690 d(模擬結(jié)束)時(shí)的演化結(jié)果。A、B兩個(gè)剖面滲透系數(shù)最小值均出現(xiàn)在5.66 m處,該時(shí)刻兩個(gè)剖面上滲透系數(shù)最小值差異較?。ˋ剖面滲透系數(shù)最小值為5.37 m/d,B剖面滲透系數(shù)最小值為5.42 m/d,相差1.11%);在該時(shí)刻,兩個(gè)剖面在滲透系數(shù)較小的從滲透系數(shù)最小值到抽液孔處(5.66~30 m)之間的演化情況一致。

    A和B兩個(gè)剖面上滲透系數(shù)在該時(shí)刻演化結(jié)果有所差異。A剖面上滲透系數(shù)區(qū)間為[5.37 m/d,14.68 m/d],B剖面上滲透系數(shù)區(qū)間為[5.42 m/d,13.20 m/d],B剖面上出現(xiàn)了滲透系數(shù)明顯大于初始值的位置,A剖面上滲透系數(shù)區(qū)間明顯大于B剖面上滲透系數(shù)區(qū)間,說(shuō)明B剖面上的滲透系數(shù)較A剖面滲透系數(shù)演化較慢。

    在該時(shí)刻,模擬區(qū)域內(nèi)僅有注液孔出的滲透系數(shù)大于初始值,且從模擬結(jié)果看出,在溶液滲流通道內(nèi)大部分區(qū)域的滲透系數(shù)均在6 m/d左右。

    通過(guò)對(duì)比模擬區(qū)域內(nèi)兩個(gè)剖面內(nèi)滲透系數(shù)演化的差異可知,在多組抽注系統(tǒng)共同作用下,不同抽注井之間滲透系數(shù)演化差異性隨時(shí)間推移不斷減小。在文中給的抽注系統(tǒng)中,滲透系數(shù)最小值出現(xiàn)在與注液孔相距大約5.66 m處;受單個(gè)抽液孔影響的注液孔周圍滲透系數(shù)在滲流路徑上演化較為單一,受多個(gè)(文中僅有兩個(gè))抽液孔影響的注液孔附近滲透系數(shù)在模擬初期演化較為多樣性;

    隨著模擬時(shí)間的推移,溶浸液滲流通內(nèi)滲透系數(shù)會(huì)不斷較小,且區(qū)域內(nèi)滲透系數(shù)會(huì)區(qū)域一致。

    從選取的兩個(gè)特征剖面來(lái)看,不同抽注系統(tǒng)內(nèi)滲透系數(shù)存在一定的相互影響的現(xiàn)象,隨著時(shí)間的推移影響越明顯,主要體現(xiàn)在模擬末期兩個(gè)剖面上的滲透系數(shù)演化趨于一致,差異性較小。

    4 結(jié)論與展望

    本文基于前人的研究基礎(chǔ),運(yùn)用數(shù)值模擬的方法,對(duì)某砂巖型鈾礦在原位地浸采鈾過(guò)程中滲透系數(shù)的時(shí)空分布規(guī)律進(jìn)行了模擬研究。通過(guò)討論分析反應(yīng)動(dòng)力學(xué)參數(shù)、地下水滲流場(chǎng)對(duì)滲透系數(shù)及溶質(zhì)運(yùn)移的影響,得到了以下主要結(jié)論:

    模擬區(qū)域內(nèi)滲透系數(shù)區(qū)間隨時(shí)間推移增大,模擬結(jié)束時(shí)滲透系數(shù)最大值為14.7 m/d,最小值為5.37 m/d。從滲透系數(shù)空間演化規(guī)律可知滲透系數(shù)最大值出現(xiàn)在注液孔處,滲透系數(shù)最小值出現(xiàn)在與注液孔中心相距約5.66 m處的圓周上;模擬結(jié)果表明,在模擬初期不同抽注系統(tǒng)之間溶質(zhì)運(yùn)移受井群干擾作用較明顯。

    本文基于參考文獻(xiàn)[44-45]的研究基礎(chǔ),建立多場(chǎng)耦合模型對(duì)內(nèi)蒙古某砂巖型鈾礦開(kāi)采中滲透系數(shù)進(jìn)行模擬研究,對(duì)提升鈾礦開(kāi)采效率具有一定的科學(xué)意義。

    由于水文地質(zhì)單元概化及模擬程序限制等因素的影響,本文在酸法地浸采鈾過(guò)程的模擬研究中未對(duì)介質(zhì)的非均一性、微生物過(guò)程、溶浸液中不同離子相互作用及溶浸液循環(huán)抽注等因素對(duì)結(jié)果影響。

    猜你喜歡
    注液浸液滲透系數(shù)
    扣式電池自動(dòng)注液機(jī)方案設(shè)計(jì)
    新型蠶蛾研磨批量自動(dòng)加注器的研發(fā)應(yīng)用
    “疊壓供水”技術(shù)在地浸采鈾設(shè)計(jì)中的應(yīng)用研究
    鈾礦冶(2022年4期)2022-11-11 08:03:46
    基于Origin的滲透系數(shù)衰減方程在地?zé)崴毓嘀械膽?yīng)用
    鉑鈀精礦酸浸與還原工藝的試驗(yàn)探索
    石煤釩礦酸浸液中萃取提釩的研究進(jìn)展與前景
    多孔材料水滲透系數(shù)預(yù)測(cè)的隨機(jī)行走法
    輸水渠防滲墻及基巖滲透系數(shù)敏感性分析
    河北平原新近系熱儲(chǔ)層滲透系數(shù)規(guī)律性分析
    從低濃度含鈀硝浸液中富集鈀工藝研究
    国产久久久一区二区三区| 99国产极品粉嫩在线观看| 91av网一区二区| 美女大奶头视频| 亚洲成人精品中文字幕电影| 精品一区二区三区人妻视频| 亚洲经典国产精华液单 | 天美传媒精品一区二区| 嫩草影院新地址| 噜噜噜噜噜久久久久久91| 观看美女的网站| 亚洲精华国产精华精| 国产淫片久久久久久久久 | 欧美最黄视频在线播放免费| 脱女人内裤的视频| 99热精品在线国产| 免费观看精品视频网站| 99精品在免费线老司机午夜| 丁香六月欧美| 乱码一卡2卡4卡精品| 婷婷丁香在线五月| 极品教师在线视频| 久久午夜福利片| 亚洲aⅴ乱码一区二区在线播放| 精品国产三级普通话版| 亚洲综合色惰| 天堂√8在线中文| 亚洲精品色激情综合| 亚洲熟妇熟女久久| 搡女人真爽免费视频火全软件 | 又黄又爽又刺激的免费视频.| 亚洲专区国产一区二区| 久久精品久久久久久噜噜老黄 | 美女 人体艺术 gogo| 国产亚洲精品久久久com| 免费在线观看亚洲国产| 国产亚洲av嫩草精品影院| 熟女人妻精品中文字幕| 51国产日韩欧美| 日本黄色片子视频| 欧美高清性xxxxhd video| 乱码一卡2卡4卡精品| 中文字幕久久专区| a级一级毛片免费在线观看| 99久久99久久久精品蜜桃| 国产色爽女视频免费观看| 成年女人毛片免费观看观看9| 给我免费播放毛片高清在线观看| 99久久无色码亚洲精品果冻| 色精品久久人妻99蜜桃| 亚洲精品粉嫩美女一区| 久久中文看片网| 在线观看一区二区三区| 午夜福利在线在线| 国产一区二区三区视频了| 亚洲色图av天堂| 在线国产一区二区在线| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日韩有码中文字幕| 中国美女看黄片| 免费大片18禁| 亚洲无线在线观看| 免费高清视频大片| 国产欧美日韩一区二区三| 成人三级黄色视频| 久久精品人妻少妇| 精品一区二区三区人妻视频| 日本成人三级电影网站| 伦理电影大哥的女人| 中文资源天堂在线| 国产伦精品一区二区三区四那| 精品久久久久久久久久久久久| 精品日产1卡2卡| bbb黄色大片| 亚洲精品久久国产高清桃花| 两个人视频免费观看高清| av在线天堂中文字幕| 国产黄片美女视频| 国产精品久久电影中文字幕| 两人在一起打扑克的视频| 日韩有码中文字幕| 99精品久久久久人妻精品| www.www免费av| 丝袜美腿在线中文| 美女xxoo啪啪120秒动态图 | 欧美成人免费av一区二区三区| 亚洲成人久久爱视频| 真人一进一出gif抽搐免费| 嫩草影院入口| 老鸭窝网址在线观看| 69人妻影院| 国产精品爽爽va在线观看网站| www.999成人在线观看| 色综合亚洲欧美另类图片| 久久精品国产自在天天线| 深夜a级毛片| 亚洲av成人不卡在线观看播放网| ponron亚洲| 日韩欧美精品v在线| 丰满人妻一区二区三区视频av| 最近视频中文字幕2019在线8| 99国产精品一区二区蜜桃av| 午夜老司机福利剧场| 真人一进一出gif抽搐免费| 久久久久性生活片| 51国产日韩欧美| 国产真实乱freesex| 国产淫片久久久久久久久 | 亚洲午夜理论影院| 国产精品一区二区免费欧美| 亚洲美女搞黄在线观看 | 亚洲精品亚洲一区二区| 人妻夜夜爽99麻豆av| 亚洲电影在线观看av| 久久香蕉精品热| 日本成人三级电影网站| .国产精品久久| 亚洲经典国产精华液单 | 国产精品女同一区二区软件 | 女生性感内裤真人,穿戴方法视频| 深夜精品福利| 亚洲精品色激情综合| 亚洲av美国av| 免费电影在线观看免费观看| 在线免费观看的www视频| 99热只有精品国产| 一进一出好大好爽视频| 精品一区二区免费观看| 香蕉av资源在线| 黄片小视频在线播放| 亚洲av美国av| 国产国拍精品亚洲av在线观看| 国产一区二区亚洲精品在线观看| 国产一区二区亚洲精品在线观看| 国产一区二区激情短视频| 国产精品三级大全| 90打野战视频偷拍视频| 最近视频中文字幕2019在线8| 免费看美女性在线毛片视频| 欧美日韩乱码在线| 女人十人毛片免费观看3o分钟| 两个人视频免费观看高清| 亚洲黑人精品在线| 久久午夜福利片| 亚洲av免费高清在线观看| 国产精品嫩草影院av在线观看 | 亚洲欧美日韩无卡精品| 色噜噜av男人的天堂激情| 久久婷婷人人爽人人干人人爱| 国产伦精品一区二区三区视频9| 色噜噜av男人的天堂激情| 久久久成人免费电影| 免费av毛片视频| 日本精品一区二区三区蜜桃| 亚州av有码| 九色成人免费人妻av| 色综合站精品国产| 免费看美女性在线毛片视频| 青草久久国产| 色av中文字幕| 色5月婷婷丁香| 国产高清视频在线播放一区| 亚洲片人在线观看| 色综合站精品国产| 男插女下体视频免费在线播放| 蜜桃亚洲精品一区二区三区| 91麻豆av在线| 丰满人妻熟妇乱又伦精品不卡| 国产色爽女视频免费观看| 中文字幕久久专区| 少妇熟女aⅴ在线视频| 亚洲av五月六月丁香网| 国内精品一区二区在线观看| 亚洲最大成人中文| 国产精品乱码一区二三区的特点| 韩国av一区二区三区四区| 日日摸夜夜添夜夜添av毛片 | 亚洲无线观看免费| 成人高潮视频无遮挡免费网站| 亚洲精品一卡2卡三卡4卡5卡| 色哟哟哟哟哟哟| 97碰自拍视频| 亚洲精品日韩av片在线观看| 久久婷婷人人爽人人干人人爱| 99久久精品热视频| 亚洲熟妇中文字幕五十中出| 琪琪午夜伦伦电影理论片6080| 偷拍熟女少妇极品色| 一个人观看的视频www高清免费观看| 男插女下体视频免费在线播放| 欧美黄色淫秽网站| 看免费av毛片| 黄色日韩在线| 国产野战对白在线观看| 亚洲精品色激情综合| 亚洲国产精品成人综合色| 国产亚洲精品久久久com| 99热这里只有是精品在线观看 | av中文乱码字幕在线| 久久这里只有精品中国| 久久人人爽人人爽人人片va | 国产成人av教育| 亚洲美女搞黄在线观看 | 97超视频在线观看视频| 18禁在线播放成人免费| 国产成人影院久久av| 三级毛片av免费| 欧美成人免费av一区二区三区| 国产高清三级在线| 亚洲国产日韩欧美精品在线观看| 亚洲色图av天堂| 人妻夜夜爽99麻豆av| 久久国产乱子免费精品| 久久久久久久久大av| 色综合站精品国产| 69av精品久久久久久| 成年人黄色毛片网站| 美女高潮的动态| 日本在线视频免费播放| av专区在线播放| 搞女人的毛片| 欧美三级亚洲精品| 日本五十路高清| 久久99热这里只有精品18| 亚洲国产欧美人成| 欧美zozozo另类| 久久午夜亚洲精品久久| 夜夜夜夜夜久久久久| 日本一本二区三区精品| 亚洲av五月六月丁香网| 一级a爱片免费观看的视频| 亚洲av.av天堂| 亚洲狠狠婷婷综合久久图片| 久久天躁狠狠躁夜夜2o2o| 精品无人区乱码1区二区| 直男gayav资源| 亚洲中文日韩欧美视频| 黄色日韩在线| 午夜免费成人在线视频| 亚洲av.av天堂| 99国产综合亚洲精品| 一本久久中文字幕| 色视频www国产| av在线老鸭窝| 国产亚洲精品av在线| 亚洲人成网站在线播放欧美日韩| 好男人电影高清在线观看| 少妇被粗大猛烈的视频| 精品久久国产蜜桃| 最近最新免费中文字幕在线| 一进一出抽搐动态| 欧洲精品卡2卡3卡4卡5卡区| 人人妻人人看人人澡| 免费无遮挡裸体视频| 日韩精品青青久久久久久| 久久久久久久久大av| 午夜亚洲福利在线播放| 国产淫片久久久久久久久 | 搡老熟女国产l中国老女人| 亚洲人成网站在线播放欧美日韩| 午夜影院日韩av| 天堂√8在线中文| 久久亚洲真实| 在线免费观看不下载黄p国产 | 国产单亲对白刺激| 免费人成在线观看视频色| 欧美最新免费一区二区三区 | 国产欧美日韩精品亚洲av| 欧美zozozo另类| 成年女人毛片免费观看观看9| 99国产精品一区二区三区| 搡老岳熟女国产| 性插视频无遮挡在线免费观看| 日本三级黄在线观看| 一个人看的www免费观看视频| 免费av毛片视频| 一本一本综合久久| 国产在线男女| 久久国产乱子伦精品免费另类| 天美传媒精品一区二区| 亚洲久久久久久中文字幕| 全区人妻精品视频| 黄色女人牲交| 久久午夜福利片| 九色成人免费人妻av| 国模一区二区三区四区视频| 婷婷丁香在线五月| 黄色配什么色好看| 午夜a级毛片| 午夜福利在线观看吧| 中文字幕免费在线视频6| 亚洲va日本ⅴa欧美va伊人久久| 12—13女人毛片做爰片一| 少妇被粗大猛烈的视频| 国产一区二区三区在线臀色熟女| 午夜福利视频1000在线观看| 校园春色视频在线观看| 国产亚洲精品综合一区在线观看| 欧美日韩黄片免| 国产蜜桃级精品一区二区三区| 午夜亚洲福利在线播放| 亚洲人与动物交配视频| 精品久久久久久久久久免费视频| 国产乱人伦免费视频| 国产真实伦视频高清在线观看 | 久久精品国产清高在天天线| 亚洲七黄色美女视频| 97超级碰碰碰精品色视频在线观看| 精品国内亚洲2022精品成人| 国产老妇女一区| 日韩中字成人| 观看美女的网站| 免费在线观看成人毛片| 国产欧美日韩精品一区二区| 亚洲国产高清在线一区二区三| 搡女人真爽免费视频火全软件 | 欧美日韩福利视频一区二区| 亚洲av成人精品一区久久| 免费av不卡在线播放| 久久精品久久久久久噜噜老黄 | 好看av亚洲va欧美ⅴa在| 国内少妇人妻偷人精品xxx网站| 精品午夜福利在线看| 欧美乱妇无乱码| 男女那种视频在线观看| 国产精品98久久久久久宅男小说| 一级a爱片免费观看的视频| 美女xxoo啪啪120秒动态图 | 亚洲精品在线美女| 高清日韩中文字幕在线| 午夜福利在线观看吧| 别揉我奶头 嗯啊视频| 国内精品一区二区在线观看| 一级a爱片免费观看的视频| av在线天堂中文字幕| 亚洲欧美激情综合另类| 人人妻人人看人人澡| 嫩草影视91久久| 国产精品爽爽va在线观看网站| 简卡轻食公司| 淫妇啪啪啪对白视频| 在线观看免费视频日本深夜| 国产伦精品一区二区三区视频9| 狠狠狠狠99中文字幕| 精品午夜福利视频在线观看一区| 在线十欧美十亚洲十日本专区| 精品人妻1区二区| 欧美性猛交╳xxx乱大交人| 精品人妻熟女av久视频| 天堂av国产一区二区熟女人妻| 成年女人看的毛片在线观看| 久久精品国产亚洲av香蕉五月| 成人午夜高清在线视频| 欧美丝袜亚洲另类 | 国产精品久久久久久久电影| 九九热线精品视视频播放| 国产成人啪精品午夜网站| 91久久精品电影网| 亚洲欧美日韩卡通动漫| 淫妇啪啪啪对白视频| 97人妻精品一区二区三区麻豆| 美女cb高潮喷水在线观看| 日日摸夜夜添夜夜添av毛片 | 国产国拍精品亚洲av在线观看| 亚洲精品456在线播放app | 国产视频一区二区在线看| www.999成人在线观看| 国产视频内射| 亚洲av二区三区四区| 亚洲av成人av| 国内精品久久久久久久电影| 人妻制服诱惑在线中文字幕| 国产精品,欧美在线| 给我免费播放毛片高清在线观看| 久久国产精品人妻蜜桃| 亚洲精品在线观看二区| 毛片女人毛片| 特级一级黄色大片| 2021天堂中文幕一二区在线观| 国产精品女同一区二区软件 | 久99久视频精品免费| 丁香六月欧美| 又粗又爽又猛毛片免费看| a级毛片a级免费在线| 精品人妻熟女av久视频| 精品久久久久久久人妻蜜臀av| 亚洲精品亚洲一区二区| 我要看日韩黄色一级片| 国产精品亚洲av一区麻豆| 老司机深夜福利视频在线观看| 精华霜和精华液先用哪个| 丝袜美腿在线中文| 欧美国产日韩亚洲一区| 在线观看66精品国产| 国产欧美日韩一区二区三| 国产亚洲av嫩草精品影院| 啦啦啦观看免费观看视频高清| 日韩欧美 国产精品| 精品久久国产蜜桃| 亚洲欧美激情综合另类| 亚洲,欧美,日韩| 欧美在线黄色| 老鸭窝网址在线观看| 国产又黄又爽又无遮挡在线| 亚洲电影在线观看av| 床上黄色一级片| 两个人视频免费观看高清| 欧美高清性xxxxhd video| 午夜两性在线视频| 天堂√8在线中文| 中文字幕精品亚洲无线码一区| netflix在线观看网站| 欧美在线黄色| 桃红色精品国产亚洲av| 美女被艹到高潮喷水动态| 亚洲国产精品sss在线观看| 亚洲精品久久国产高清桃花| 脱女人内裤的视频| 免费看a级黄色片| 亚洲av免费高清在线观看| 在线a可以看的网站| 男女之事视频高清在线观看| 欧美最黄视频在线播放免费| 一a级毛片在线观看| 国产视频内射| ponron亚洲| 午夜福利视频1000在线观看| 白带黄色成豆腐渣| 综合色av麻豆| 欧美一级a爱片免费观看看| 香蕉av资源在线| 色av中文字幕| 久久国产乱子伦精品免费另类| 国产美女午夜福利| 精品国产亚洲在线| 欧美xxxx黑人xx丫x性爽| 国产免费男女视频| 乱人视频在线观看| 午夜激情福利司机影院| 午夜福利在线观看免费完整高清在 | 国产精品久久久久久亚洲av鲁大| 国产欧美日韩一区二区精品| 国产精品久久久久久久电影| 国产精品久久视频播放| 欧美不卡视频在线免费观看| 婷婷精品国产亚洲av| 51午夜福利影视在线观看| 国产精品不卡视频一区二区 | 看片在线看免费视频| 国产国拍精品亚洲av在线观看| 九色国产91popny在线| 99久久精品热视频| 内地一区二区视频在线| 五月伊人婷婷丁香| 桃色一区二区三区在线观看| av天堂在线播放| 一本久久中文字幕| 日韩欧美在线乱码| 97超视频在线观看视频| 欧美丝袜亚洲另类 | 精品国内亚洲2022精品成人| 中文字幕人成人乱码亚洲影| 亚洲成人久久爱视频| 亚洲中文字幕日韩| 欧美zozozo另类| 在线国产一区二区在线| 亚州av有码| 欧洲精品卡2卡3卡4卡5卡区| 久久久久亚洲av毛片大全| 高清毛片免费观看视频网站| 欧美激情国产日韩精品一区| www.999成人在线观看| 中文字幕免费在线视频6| 亚洲 国产 在线| 日韩高清综合在线| 欧美激情国产日韩精品一区| 长腿黑丝高跟| 亚洲美女搞黄在线观看 | 精品一区二区三区视频在线| 少妇高潮的动态图| 国产一区二区亚洲精品在线观看| 国产精品永久免费网站| 婷婷精品国产亚洲av在线| 亚洲五月婷婷丁香| 亚洲国产精品合色在线| 精品久久久久久久久av| 99热这里只有是精品50| 精品国产三级普通话版| 国产精品亚洲av一区麻豆| 精品人妻偷拍中文字幕| 色精品久久人妻99蜜桃| 中文字幕av在线有码专区| 免费av毛片视频| 亚洲一区高清亚洲精品| 亚洲电影在线观看av| 亚洲av第一区精品v没综合| 亚洲三级黄色毛片| 最新中文字幕久久久久| 午夜激情福利司机影院| 午夜两性在线视频| 少妇裸体淫交视频免费看高清| 人妻丰满熟妇av一区二区三区| 久久精品久久久久久噜噜老黄 | 99久久久亚洲精品蜜臀av| 久久久久久久精品吃奶| 中文字幕高清在线视频| 久久中文看片网| 亚洲第一欧美日韩一区二区三区| 午夜福利在线观看免费完整高清在 | 国产精品综合久久久久久久免费| 一个人观看的视频www高清免费观看| 久久久久久久久中文| 久久久久久久久久黄片| 欧美日韩亚洲国产一区二区在线观看| 此物有八面人人有两片| 给我免费播放毛片高清在线观看| 亚洲国产精品合色在线| 国产精品爽爽va在线观看网站| av在线老鸭窝| avwww免费| 成人鲁丝片一二三区免费| 久久久久久九九精品二区国产| 久久久久国内视频| 亚洲欧美日韩东京热| a级毛片免费高清观看在线播放| 中文在线观看免费www的网站| 婷婷色综合大香蕉| 在线看三级毛片| 国产高清视频在线播放一区| 1024手机看黄色片| 麻豆久久精品国产亚洲av| 美女cb高潮喷水在线观看| 国产精品久久久久久久电影| 亚洲成av人片免费观看| 精品一区二区三区视频在线观看免费| 久久久色成人| 91字幕亚洲| 99热这里只有是精品在线观看 | 在线看三级毛片| 美女大奶头视频| 蜜桃久久精品国产亚洲av| 色综合亚洲欧美另类图片| 国产aⅴ精品一区二区三区波| 欧美黑人欧美精品刺激| av天堂中文字幕网| 我要搜黄色片| 又黄又爽又免费观看的视频| 啪啪无遮挡十八禁网站| 欧美区成人在线视频| 欧美黑人欧美精品刺激| 啦啦啦韩国在线观看视频| a级毛片a级免费在线| 老司机深夜福利视频在线观看| 免费人成视频x8x8入口观看| 国产亚洲精品久久久久久毛片| 麻豆国产av国片精品| 亚洲精品在线观看二区| 婷婷丁香在线五月| 欧美午夜高清在线| 免费在线观看日本一区| a级毛片a级免费在线| 国产免费男女视频| av视频在线观看入口| 亚洲色图av天堂| 色精品久久人妻99蜜桃| 国产老妇女一区| 国产人妻一区二区三区在| 欧美又色又爽又黄视频| 成年女人毛片免费观看观看9| avwww免费| 中文字幕人妻熟人妻熟丝袜美| 国产欧美日韩精品亚洲av| 国产精品人妻久久久久久| а√天堂www在线а√下载| .国产精品久久| 国产精品嫩草影院av在线观看 | 国产午夜精品论理片| 国产高清有码在线观看视频| 亚洲色图av天堂| 欧美乱妇无乱码| 午夜福利在线在线| 免费黄网站久久成人精品 | 精品久久久久久久久久久久久| 丰满人妻一区二区三区视频av| 怎么达到女性高潮| 一个人免费在线观看的高清视频| 在线国产一区二区在线| 制服丝袜大香蕉在线| 国产白丝娇喘喷水9色精品| 男人的好看免费观看在线视频| 成人av一区二区三区在线看| 在线看三级毛片| 97人妻精品一区二区三区麻豆| 亚洲精品乱码久久久v下载方式| 久久这里只有精品中国| 国产又黄又爽又无遮挡在线| 一二三四社区在线视频社区8| 国产激情偷乱视频一区二区| 欧美性感艳星| 国产精品1区2区在线观看.| 少妇人妻精品综合一区二区 | 中文字幕精品亚洲无线码一区| 久久中文看片网| 欧美精品啪啪一区二区三区| 狂野欧美白嫩少妇大欣赏| 美女黄网站色视频| 欧美激情久久久久久爽电影| 亚洲人成伊人成综合网2020| 深夜精品福利| 国产乱人伦免费视频| 色视频www国产| 床上黄色一级片| 国产精品亚洲一级av第二区| 婷婷丁香在线五月| 噜噜噜噜噜久久久久久91|