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

    細粒層存在條件下礦巖散體內的溶液流動特性

    2017-10-13 14:36:13尹升華王雷鳴潘晨陽謝芳芳
    中國有色金屬學報 2017年3期
    關鍵詞:壓力降礦巖細粒

    尹升華,王雷鳴,潘晨陽,陳 勛, 謝芳芳

    ?

    細粒層存在條件下礦巖散體內的溶液流動特性

    尹升華1, 2,王雷鳴1, 2,潘晨陽1,陳 勛1, 2, 謝芳芳1

    (1. 北京科技大學土木與資源工程學院金屬礦山高效開采與安全教育部重點實驗室,北京 100083;2. 北京金誠信礦山技術研究院有限公司,北京101500)

    柱浸過程中礦巖散體的不均勻分布影響著溶液的滲流規(guī)律。為考察細粒層的存在下溶液流動特性,利用CT技術與COMSOL Multiphysics多相耦合模擬軟件,開展礦巖散體內溶液滲流模擬研究,獲得細粒層存在條件下的流體流動軌跡、滲流速度場及壓力場分布規(guī)律。結果表明:細粒層影響著溶液的流動軌跡和優(yōu)先流的形成,具體表現(xiàn)為大部分溶液繞過細粒層形成優(yōu)先流,細小支流橫穿細粒層,在細粒層內部溶液難以到達的位置,溶液流速較低,形成溶液流動的停滯區(qū)。細粒層內溶液速度降與壓力降較小,在相鄰細粒層的孔喉位置處,溶液流速及該處壓力明顯增高。

    滲流;偏析;優(yōu)先流;數(shù)值模擬;流動行為;CT技術

    隨著經(jīng)濟的飛速發(fā)展,礦產(chǎn)資源的需求日益枯竭,對低品位、難處理的銅等貴重金屬的開采成為采礦未來發(fā)展的重要方向[1]。溶浸采礦能夠有效地回收復雜難處理的低品位礦石,并具有經(jīng)濟性、高效性與環(huán)境友好性等優(yōu)點,因而被廣泛應用于提取難采選礦石及廢石中的有用元素[2?6]。

    近年來,核磁共振技術(MRI)、計算機斷層掃描技術(CT)以及粒子圖像測速法(PIV)等技術,因其具有無損地內部孔隙結構、溶液分布及演化規(guī)律等優(yōu)點,被廣泛應用礦巖散體內溶液滲流規(guī)律研究。比如:FAGAN等[7]基于MRI與重力研究滴灌堆體內部水體及微生物運移繁殖規(guī)律,吳愛祥等[8?9]利用MRI技術和CT技術,對柱浸系統(tǒng)中礦巖散體結構及溶液分布情況開展研究。DHAWAN等[10]對CT技術在堆浸系統(tǒng)分析及應用的最新進展進行總結。尹升華等[11]利用PIV技術對礦堆溶液滲流場進行了無損探測,并考察了噴淋強度等對滲流規(guī)律的影響。堆浸堆體內部礦石與孔隙分布分布具有不均勻性和隨機性,MILLER 等[12]和LIN等[13?14]基于CT技術與結合Lattice Boltzmann(LB),對柱浸體系內流體滲濾過程中的孔隙結構內的流體進行了模擬,并分析獲得了孔隙結構與柱內有效流之間的關系。王貽明等[15]利用超聲波技術,探究了其強化浸出低品位硫化銅礦的可行性和影響機理。利用COMSOL Multiphysics等多相多場耦合模擬軟件分析浸礦孔隙內部流體流動特征[16]。此外,DHAWAN等[17]研究也證實了不同尺寸礦石顆粒的不均勻分布對堆體內溶液的擴散和有用元素的浸出具有明顯作用。葉勇軍等[18]開展柱浸試驗,考察了不同粒徑分布分維數(shù)條件下的鈾礦石浸出規(guī)律。吳愛祥等[19]基于顆粒群組構特性,對氧化銅礦柱浸滲流行為進行分析。但是,對于細粒區(qū)域存在條件下柱浸體系內溶液流動軌跡、滲流速度場及壓力場特性的細觀研究仍較為匱乏。

    本文作者針對堆浸筑堆過程中常見的礦石不均勻分布(顆粒偏析)這一現(xiàn)象,將CT技術與COMSOL Multiphysics耦合模擬軟件結合,對礦巖散體內的溶液流動軌跡、滲流速度場及出口邊界速度分布、滲流壓力場分布特性等進行研究,研究結果對于改善堆浸筑堆方式、提高浸出效果等具有很好的借鑒 意義。

    1 實驗

    1.1 CT圖像獲取及二值化處理

    本實驗采用福建某次生硫化銅礦,采用德國 Siemens AG X型線電子計算機斷層掃描儀,獲取若干柱浸系統(tǒng)的橫縱截面的CT圖像。利用Adore Illustrator圖像處理軟件對二維灰度圖像進行二值化處理,如圖1所示。其中,黑色部分表示礦石間孔隙,白色部分表示礦石顆粒。此外,材料的密度是通過材料的衰減系數(shù)[20],以X射線表示。當X射線穿過被檢物,透過光線的明亮程度依據(jù)式(1)所示:

    式中:0是X射線穿透物體前的強度;是X射線穿透物體后的強度;m是每檢測到的對象的單位質量吸收系數(shù)。由圖1可見,在礦巖散體內存在3個以細粒層為主導的細粒區(qū)域,將其分別命名為細粒區(qū)域、和??梢钥吹郊氼w粒主要位于柱體的核部,而粗顆粒主要位于柱體的周邊位置,由此導致礦石筑堆過程中不同粒徑礦石顆粒的分布不均勻,即偏析現(xiàn)象。

    圖1 柱體縱截面CT圖像獲取及圖像的二值化處理

    Fig. 1 CT images and its binarized images of column longitudinal sections: (a) CT scanning apparatus; (b) Pretreatment of processing Binarized

    1.2 模型構建及邊界條件設定

    利用COMSOL Multiphysics耦合模擬軟件,構建二維滲流模型并設定邊界條件。將二值化圖像導入Adobe Illustrator進行圖像矢量化處理,獲取DXF文件,再導入COMSOL Multiphysics耦合軟件,在二維多孔介質模型環(huán)境下,通過常微分方程和邊界條件等控制模擬過程,研究細粒區(qū)域存在條件下,礦石孔隙之間流體的流動軌跡、滲流速度場及壓力場分布特性。

    其中,多孔介質中的慢速流動由Brinkman方程控制,如式(2)及(3)所示;礦石間孔隙中的飽和流由Navier-Stokes方程控制,如式(4)所示。

    (3)

    (4)

    基于Darcy定律及水流連續(xù)性方程的一般滲流偏微分方程,如式(5)所示:

    式中:是壓力;是達西速度場;是溶液的動態(tài)黏度;P是孔隙度;為介質的磁導率;為溶液密度;為位置水頭;m為源匯項;為孔隙率。

    入口邊界的約束條件如式(6)所示:

    出口邊界的約束條件如式(7)所示:

    (7)

    式中:0為0.715 Pa。此外,具體的滲流參數(shù)如表1所列。

    此外,建設模型前,需要作出以下假設:

    1) 溶液流動僅發(fā)生在孔隙中,未滲透到礦石顆粒內部,忽略毛細作用;

    2) 礦石顆粒的粒度是不變化的;

    3) 孔隙中的溶液是不可壓縮、連續(xù)的。

    最終,通過設定的邊界條件與初始值等,導入COMSOL Mutliphysics多場多相耦合模擬軟件,獲得滲流物理模型、網(wǎng)格模型以及計算網(wǎng)格的局部放大情況,如圖2所示。

    表1 滲流模擬的關鍵參數(shù)

    圖2 滲流物理模型及網(wǎng)格劃分

    2 結果與討論

    2.1 基于流線的溶液流動軌跡模擬

    滲流速度場內流體流線模擬結果如圖3所示。設置流線模擬的起點為柱體上表面,藍色曲線即為流線,用于表征礦巖散體內的溶液流動軌跡及優(yōu)先流通道等滲流特性。其中,流線為粗細均一的曲線,流線的數(shù)量與流體的流量呈正比,流線間距與該處流體流量呈正比,其積分容差為0.001,穩(wěn)態(tài)點停止容差與循環(huán)容差均為0.01,流線數(shù)量為20。

    由圖3(a)可見,在礦體顆粒非均勻分布的作用下,柱體上部的流線分布較分散,間距較大;下部流線分布較為集中,并形成了以兩股優(yōu)先流為主導的流動軌跡。將流線模擬結果圖進行局部放大,得到圖3(b)。由圖3(b)可見,在流線不斷向下延伸傳遞的過程中,位于柱體核部的細粒區(qū)域對流線分布具有明顯的分流作用。依據(jù)細粒區(qū)域內顆粒的尺寸與外形的差異,礦體顆粒對流線分布的影響程度略有差別,流線間不斷的匯聚與分開表征著溶液的不斷匯合與分流。此外,在毛細吸力的作用下,尺寸較細的顆粒傾向于被粘結在一起。在這種橋接效應下,顆粒處于擺動狀態(tài),該現(xiàn)象也在前人的研究中得以驗證[21]。

    在柱浸體系下部區(qū)域,礦巖粒徑加大且分布隨機,相比細粒區(qū)域,下部的優(yōu)先流更加明顯。由此表明,由于細粒區(qū)域內若干微細裂隙的存在,在大部分溶液繞過細粒區(qū)域的同時,部分溶液被分流形成若干局部優(yōu)先流,以不斷繞過細顆粒的方式向下滲流,縱穿細粒區(qū)域,并在溶液難以到達的部位形成了數(shù)量較多的停滯區(qū)。這一現(xiàn)象也通過后續(xù)滲流速度場、壓力場的研究得以驗證。

    2.2 礦巖孔隙內部及出口邊界的速度場分布

    滲流速度場及壓力場分布情況的模擬結果如圖4和5所示。其中,亮藍色表示流體的流速較高,暗藍色表示速度較低。如圖4(a)所示,在礦石孔隙中的不同位置有若干矢量箭頭,其方向標明了溶液流動方向,矢量箭頭的長度與該處流體流速成正比。

    由圖4(a)可知,溶液集中分布于柱體中間部位,在柱體近邊壁位置的溶液量較小,流速較低,形成明顯的浸出盲區(qū)。礦巖散體孔隙內的流體逐步形成優(yōu)先流,將圖4(a)局部放大后獲得圖4(b),可見在細粒區(qū)域內,溶液矢量箭頭較短、分布均勻且為暗藍色,表明該處溶液流速較小,溶液流動緩慢,即細粒區(qū)域導致了溶液流動停滯區(qū)的生成,如圖4(b)紅色虛線圈出部分所示。在礦巖散體的狹小通道區(qū)域內,流速矢量箭頭較長且呈亮藍色,在礦石顆粒的分流作用下,箭頭方向不同,連通的狹小管狀通道內溶液增速明顯,隨滲流孔道的直徑減小,滲流速度明顯增大。因此,細粒區(qū)域導致溶液流動停滯區(qū)的形成,并以此間接地控制著溶液流動軌跡及滲流場特性。

    圖3 滲流場內溶液流動軌跡模擬

    圖4 細粒區(qū)域存在條件下礦巖散體內滲流速度場分布

    流體流速分布的高度表達式是以流速值為因變量,凸起的高度為因變量的關系曲面,如圖5所示。其中,流體流速與凸起的高度成正比。研究表明:在柱體的核部形成了優(yōu)勢凸起,這表明在浸礦過程中,礦巖散體中的優(yōu)先流逐漸形成并主要集中于柱體核部區(qū)域,并且該優(yōu)先流在柱體下部得到了明顯體現(xiàn)。優(yōu)先流繞過大塊礦巖顆粒向下流動,在柱體左上部一狹小通道內達局部流體速度的最高值,為1.98×10?5 m/s。由圖5(b)可知,柱體前段凸起較少,而末端出口位置凸起較多,表明流體進入柱體后,在細粒區(qū)域的作用下流體流速較為穩(wěn)定且流速值較低;在流體通過細粒區(qū)域后,流速呈現(xiàn)普遍增加的趨勢且變化程度增加。這表明細粒層的存在是導致形成溶液流動停滯區(qū)的重要因素。在停滯區(qū)內,溶液量稀少且溶液流動緩慢,處于相對停滯的狀態(tài),流體交換主要依靠液體的橫向毛細作用和滲透壓作用[22]。因此,停滯區(qū)中的浸出作用緩慢且效果較差,有用元素回收率較低。

    圖5 基于高度表達式的滲流速度場分布

    圖6 出口邊界速度場分布

    圖6所示為模型下部出口邊界的速度分布情況,出口邊界的位置已于圖4中標明。其中,橫坐標表示與出口邊界左端的距離,藍線表示流體流速的分布規(guī)律,紅線表示礦石顆粒,該處流體流速為0 m/s。

    由圖6可知,在溶液的出口邊界處,速度分布曲線呈現(xiàn)若干不同程度的拋物線狀凸起的形態(tài),曲線凸起的高度與該處溶液的流速成正比,凸起曲線斜率與孔道直徑減小的幅度成正比。對于由多尺寸礦石構成的滲流通道,在黏滯力的作用下,礦石表面溶液流速較低;隨取值點與孔道中心的距離減小,溶液流速逐漸增加。即在孔喉的中心處流速最大,為1.882×10?6 m/s。隨著流體與礦石顆粒邊界之間距離的縮短,流體流速逐漸減小,并于固體?水界面處降至0 m/s??缀碇行某隽魉倏臁⑦叡谔幜魉俾牧魉俜植家?guī)律,證實了礦石表面不動水層的存在及流動孔道邊壁的黏滯作用[23]。此外,孔隙的有效連通率與筑堆礦石的尺寸成正比[24],該現(xiàn)象也被本次模擬所證實。

    2.3 滲流壓力場的分布特性

    滲流壓力場分布情況如圖7所示。圖7存在以紅、橙、黃及藍4色區(qū)分的4個壓力分布區(qū),壓力值呈現(xiàn)非均勻遞減的規(guī)律。依據(jù)柱體內壓力降變化規(guī)律,可將柱體內的壓力分布自上而下分為4個區(qū)域,分別命名為1、2、3與4,其中,1~3中均存在細粒區(qū)域。

    由圖7(a)可知,柱體上部的3個壓力分布區(qū)1~3內的等壓線間距較大,分布較為稀疏;最底部的壓力分布區(qū)4內的等壓線間距較小,分布較為稠密。這表明1~3區(qū)域內壓力值較為均一,壓力降較小,而下部區(qū)域壓力降較大。這表明細粒區(qū)域的存在對階梯狀的壓力降區(qū)域的形成起到了決定性作用。具體而言,在細粒區(qū)域內,壓力分布較為均勻且壓力降較小,細粒層內的溶液擴散主要依靠毛細作用,處于非飽和狀態(tài)。由圖7可見,狹小孔道集中于相鄰細粒區(qū)域結合部,當粗顆粒礦石所占比例增加時,該區(qū)域的孔隙率增加,導致壓力下降的速率增大。由圖7(b)可見,1、2及3三者的壓力降均較小,4壓力降較大,1~3壓力降總和基本等同于4壓力降。這表明細粒區(qū)域內的流體壓力下降緩慢,即細粒區(qū)域對流體壓力的下降具有延遲作用。當流體遇到礦石顆粒時會形成一個局部的壓力升高區(qū),等壓線繞過礦石繼續(xù)向下傳遞。流體通過最下部的細粒區(qū)域后,出現(xiàn)一個壓力降紊亂區(qū),該區(qū)域壓力降不再受細粒區(qū)域干擾,主要是相鄰礦巖顆粒耦合作用的結果,進一步證明了柱體內部礦石顆粒分布的不均勻性與隨機性對壓力降起到的重要作用。細粒區(qū)域的存在控制壓力降,進而控制孔隙內溶液的流速。

    圖7 細粒區(qū)域存在條件下礦巖散體內滲流壓力場分布

    3 結論

    1) 細粒層對溶液的流動軌跡和優(yōu)先流的形成具有明顯影響。由于細粒層內有效連通的孔道數(shù)量較少,大部分溶液繞過細粒層,在細粒層外表面形成優(yōu)先流,極少量溶液流經(jīng)細粒層內的微細孔隙,自上而下緩慢地穿過細粒層,該部分溶液流過細粒層后再次匯聚形成優(yōu)先流。

    2) 細粒層內的流速及壓力降均明顯小于粗顆粒礦石區(qū)域,單位時間內流速及壓力下降速率穩(wěn)定且較小。細粒層內溶液難以到達的區(qū)域形成了溶液流動的停滯區(qū),該區(qū)域內溶液擴散主要依靠毛細作用,溶液流速較低且溶液的壓力降幾乎為0,處于非飽和狀態(tài)。

    3) 溶液滲流場分布整體上具有明顯的非均勻性,在細粒層內溶液流動緩慢且局部流向相對均一,細粒層的存在延緩了溶液滲流速率及壓力下降速率。在黏滯力作用下,礦石顆粒表面的流體流速近似為0;隨著孔喉中心與礦石表面的距離增大,孔道內流體流速壓力逐漸增大并于孔喉中心處取最大值。

    [1] WATLING H R. The bioleaching of sulphide minerals with emphasis on copper sulphides—A review[J]. Hydrometallurgy, 2006, 84(1): 81?108.

    [2] PANDA S, BISWAL A, MISHRA S, PANDA P K, PRADHAN N, MOHAPATRA U, SUKLA L B, B K MISHRA, A AKCIL. Reductive dissolution by waste newspaper for enhanced meso-acidophilic bioleaching of copper from low grade chalcopyrite: A new concept of biohydrometallurgy[J]. Hydrometallurgy, 2015, 153(3): 98?105.

    [3] 賴紹師, 覃文慶, 楊聰仁, 王 軍, 張雁生, 張 博, 常自勇, 匡浩華. 低品位硫化銅礦的細菌浸出[J]. 中國有色金屬學報, 2011, 21(6): 1473?1479. LAI Shao-shi, QIN Wen-qing, YANG Cong-ren, WANG Jun, ZHANG Yan-sheng, ZHANG Bo, CHANG Zi-yong, KUANG Hao-hua. Bioleaching of low grade copper sulfide ore[J]. The Chinese Journal of Nonferrous Metals, 2011, 21(6): 1473?1479.

    [4] MISHRA S, PANDA P P, PRADHAN N, SATAPATHY D, SUBUDHI U, BISWAI S K, MISHRA B K. Effect of native bacteria Sinomonas flava 1C and Acidithiobacillus ferrooxidans on desulphurization of Meghalaya coal and its combustion properties[J]. Fuel, 2014, 117(1): 415?421.

    [5] WU A X, YIN S H, WANG H J, QIN W Q, QIU G Z. Technological assessment of a mining-waste dump at the Dexing copper mine, China, for possible conversion to an in situ bioleaching operation[J]. Bioresource Technology, 2009, 100(6): 1931?1936.

    [6] 王少勇, 吳愛祥, 王洪江, 尹升華, 顧曉春. 高含泥氧化銅礦水洗?分級堆浸工藝[J]. 中國有色金屬學報, 2013, 23(1): 229?237. WANG Shao-yong, WU Ai-xiang, WANG Hong-jiang, YIN Sheng-hua, GU Xiao-chun. Craft of wash-classification of heap leaching used in high-clay copper oxide ore[J]. The Chinese Journal of Nonferrous Metals, 2013, 23(1): 229?237.

    [7] FAGAN M A, NGOMA I E, CHIUME R A, MINNAAR S, SEDERMAN A J, JOHNS M L, HARRION S T L. MRI and gravimetric studies of hydrology in drip irrigated heaps and its effect on the propagation of bioleaching micro-organisms[J]. Hydrometallurgy, 2014, 150(12): 210?221.

    [8] 吳愛祥, 薛振林, 尹升華, 王 恒. 基于核磁共振技術的礦巖散體結構及溶液分布研究[J]. 中國礦業(yè)大學學報, 2014, 43(4): 582?587. WU Ai-xiang, XUE Zhen-lin, YIN Sheng-hua, WANG Heng. Study of the granular medium structure of mine rock and solution distribution based on MRI technology[J]. Journal of China University of Mining & Technology, 2014, 43(4): 582?587.

    [9] 薛振林, 吳愛祥, 尹升華, 劉 超. 氧化銅礦堆浸雙重介質分形滲流演化規(guī)律[J]. 中國有色金屬學報, 2014, 24(9): 2373?2379. XUE Zhen-lin, WU Ai-xiang, YIN Sheng-hua, LIU Chao. Evolution rules of dual medium fractal seepage in copper oxide ore heap leaching[J]. The Chinese Journal of Nonferrous Metals, 2014, 24(9): 2373?2379.

    [10] DHAWAN N, SAFARZADEH M S, MILLER J D, MOATS M S, RAJAMANI R K, LIN C L. Recent advances in the application of X-ray computed tomography in the analysis of heap leaching systems[J]. Minerals Engineering, 2012, 35(6): 75?86.

    [11] 尹升華, 陳 勛, 吳愛祥, 王雷鳴, 劉 超. 基于PIV技術的細觀礦堆溶液滲流場無擾動測試[J]. 中南大學學報(自然科學版), 2015, 46(7): 2597?2604. YIN Sheng-hua, CHEN Xun, WU Ai-xiang, WANG Lei-ming, LIU Chao. Undisturbed test on mesoscopic solution seepage field in ore heap based on PIV[J].Journal of Central South University (Science and Technology), 2015, 46(7): 2597?2604.

    [12] MILLER J D, LIN C L, GARCIA C, ARIAS H. Ultimate recovery in heap leaching operations as established from mineral exposure analysis by X-ray microtomography[J]. International Journal of Mineral Processing, 2003, 72(3): 331?340.

    [13] LIN C L, MILLER J D, GARCIA C. Saturated flow characteristics in column leaching as described by LB simulation[J]. Minerals Engineering, 2005, 18(10): 1045?1051.

    [14] LIN C L, MILLER J D. Pore structure analysis of particle beds for fluid transport simulation during filtration[J]. International Journal of Mineral Processing, 2004, 73(2): 281?294.

    [15] 王貽明, 吳愛祥, 艾純明. 低品位硫化銅礦超聲強化浸出實驗與機理分析[J]. 中國有色金屬學報, 2013, 23(7): 2019?2025. WANG Yi-ming, WU Ai-xiang, AI Chun-ming. Experiment and mechanism analysis on leaching process of low grade copper sulfide intensified by ultrasonic wave[J]. The Chinese Journal of Nonferrous Metals, 2013, 23(7): 2019?2025.

    [16] FOURIE W, SAID R, YOUNG P, BARNES D L. The simulation of pore scale fluid flow with real world geometries obtained from X-ray computed tomography[C]//Proceedings of the COMSOL Conference. Boston: COMSOL, 2007.

    [17] DHAWAN N, SAFARZADEH M S, MILLER J D, MOATS M S, RAJAMANI R K. Crushed ore agglomeration and its control for heap leach operations[J]. Minerals Engineering, 2013, 41(2): 53?70.

    [18] 葉勇軍, 丁德馨, 李廣悅, 扶海鷹, 宋鍵斌, 胡 南. 不同粒徑分布分維數(shù)鈾礦石的浸出規(guī)律[J]. 中國有色金屬學報, 2013, 23(10): 2921?2927. YE Yong-jun, DING De-xin, LI Guang-yue, FU Hai-ying, SONG Jian-bin, HU Nan. Leaching behavior of uranium ore with different fractal dimensions of particle size distribution[J]. The Chinese Journal of Nonferrous Metals, 2013, 23(10): 2921?2927.

    [19] 吳愛祥, 姚高輝, 薛振林, 王貽明, 胡凱建. 基于顆粒群組構特性的氧化銅礦柱浸滲流行為分析[J]. 中南大學學報(自然科學版), 2014, 45(5): 1605?1611. WU Ai-xiang, YAO Gao-hui, XUE Zhen-lin, WANG Yi-ming, HU Kai-jian. Seepage behavior of column leaching of copper oxide ore based on particles fabric characteristics[J]. Journal of Central South University (Science and Technology), 2014, 45(5): 1605?1611.

    [20] YANG B H, WU A X, MIAO X X, LIU J Z. 3D characterization and analysis of pore structure of packed ore particle beds based on computed tomography images[J]. Transactions of Nonferrous Metals Society of China, 2014, 24(3): 833?838.

    [21] KRISTENSEN H G, SCHAEFER T. Granulation: A review on pharmaceutical wet-granulation[J]. Drug Development and Industrial Pharmacy, 1987, 13(4/5): 803?872.

    [22] WU A X, YIN S H, QIN W Q, LIU J Z, QIU G Z. The effect of preferential flow on extraction and surface morphology of copper sulphides during heap leaching[J]. Hydrometallurgy, 2009, 95(1/2): 76?81.

    [23] WU A X, YANG B H, ZHOU X. Fractal analysis of granular ore media based on computed tomography image processing[J]. Transactions of Nonferrous Metals Society of China, 2008, 18(6): 1523?1528.

    [24] YIN S H, WU A, HU K J, WANG Y M, XUE Z L. Visualization of flow behavior during bioleaching of waste rock dumps under saturated and unsaturated conditions[J]. Hydrometallurgy, 2013, 133(2): 1?6.

    (編輯 王 超)

    Fluid flowing characteristics in ore granular with fine interlayers existed

    YIN Sheng-hua1, 2, WANG Lei-ming1, 2, PAN Chen-yang1, CHEN Xun1, 2, XIE Fang-fang1

    (1. Key Laboratory of Ministry of Education for High-Efficient Mining and Safety of Metal,School of Civil and Resource Engineering, University of Science and Technology Beijing, Beijing 100083, China;2. Beijing JCHX Mine Technology Research Institute Co., LTD, Beijing 101500, China)

    The uneven distribution of ore granular media affects the solution seepage law in heap dumping process. To investigate the solution seepage law with fine interlayers existed, the seepage simulation of column with fines interlayers existed was carried out by using the computed tomography (CT) technology and COMSOL Multiphysics modeling software. The fluid flow trajectory, distribution law of seepage velocity and pressure field with fine interlayers existed were ascertained. The results show that fine interlayers influence flow trajectory and formation of preferential flow. In detail, the preferential flow is emerged after the most solution bypassed the fine interlayers, and the plenty of small tributaries flow across fine interlayers. Stagnant regions of fluid flow are emerged where the solution is hard to reach in the internal of fine interlayers. The drops of speed and pressure are small in the fine interlayers, and the speed and pressure increase significantly when the fluid passes through the pore throats connecting the adjacent fine interlayers.

    seepage; segregation; preferential flow; numerical simulation; flow behavior; computed tomography technology

    Project (51374035) supported by the National Natural Science Foundation of China; Project (NCET-13-0669) supported by the Program for New Century Excellent Talents in University, China; Project (201351) supported by the Foundation for the Author of National Excellent Doctoral Dissertation of China

    2015-10-20; Accepted date:2016-09-23

    YIN Sheng-hua; Tel: +86-10-62332750; E-mail: csuysh@126.com

    10.19476/j.ysxb.1004.0609.2017.03.016

    1004-0609(2017)-03-0574-08

    TD862

    A

    國家自然科學基金資助項目(51374035);新世紀優(yōu)秀人才支持計劃資助(NCET-13-0669);全國優(yōu)秀博士學位論文作者專項基金資助項目(201351)

    2015-10-20;

    2016-09-23

    尹升華,教授,博士;電話:010-62332750;E-mail:csuysh@126.com

    猜你喜歡
    壓力降礦巖細粒
    陽極壓力降在PEMFC故障診斷中的應用
    電源技術(2023年11期)2023-11-29 03:00:04
    垂直溜井貯礦段放礦中礦巖運動速度變化特征
    金屬礦山(2021年9期)2021-10-20 09:54:56
    MG-9#捕收劑在極細粒煤泥浮選中的應用
    選煤技術(2021年3期)2021-10-13 07:33:36
    細粒級尾砂高濃度膠結充填試驗研究與工業(yè)應用
    傾斜溜井中的礦巖運動特征及其對井壁的損傷與破壞
    金屬礦山(2020年9期)2020-10-26 13:51:54
    礦巖初始運動對其沖擊溜井井壁規(guī)律的影響
    金屬礦山(2020年9期)2020-10-26 13:51:52
    機載火焰抑制器流通性能研究
    濟陽陸相斷陷湖盆泥頁巖細粒沉積層序初探
    微細粒磁鐵礦提鐵降硫試驗研究
    氣體熱載體干餾爐內壓力降的研究
    亚洲人成电影观看| 亚洲aⅴ乱码一区二区在线播放 | 亚洲 欧美一区二区三区| 一级毛片女人18水好多| 色综合婷婷激情| 国产成+人综合+亚洲专区| 看黄色毛片网站| 欧美绝顶高潮抽搐喷水| 日韩精品中文字幕看吧| 免费在线观看影片大全网站| 一级黄色大片毛片| 老熟妇仑乱视频hdxx| 999精品在线视频| 99国产极品粉嫩在线观看| 精品国产一区二区三区四区第35| 巨乳人妻的诱惑在线观看| 亚洲激情在线av| 免费在线观看完整版高清| 成在线人永久免费视频| 黄色视频,在线免费观看| 夜夜爽天天搞| 波多野结衣一区麻豆| 欧美丝袜亚洲另类 | 99久久国产精品久久久| 精品欧美国产一区二区三| 操美女的视频在线观看| 91大片在线观看| 看免费av毛片| 欧美乱码精品一区二区三区| 久久精品亚洲熟妇少妇任你| 麻豆成人av在线观看| 中文亚洲av片在线观看爽| 亚洲色图av天堂| 每晚都被弄得嗷嗷叫到高潮| 乱人伦中国视频| 嫩草影视91久久| 亚洲自拍偷在线| 午夜免费观看网址| 久久精品影院6| 香蕉丝袜av| 最近最新免费中文字幕在线| 一级作爱视频免费观看| 国产97色在线日韩免费| 午夜久久久在线观看| 一本综合久久免费| 国产精品99久久99久久久不卡| 亚洲五月色婷婷综合| 女同久久另类99精品国产91| 悠悠久久av| 亚洲一区二区三区色噜噜| 伊人久久大香线蕉亚洲五| 国产成人系列免费观看| 亚洲片人在线观看| 黄色视频不卡| 国产精品爽爽va在线观看网站 | 91精品三级在线观看| 久久久久国内视频| 可以在线观看毛片的网站| 老汉色∧v一级毛片| 亚洲男人天堂网一区| 十八禁人妻一区二区| 国产精品一区二区在线不卡| 亚洲三区欧美一区| 久久国产乱子伦精品免费另类| 黄片小视频在线播放| 视频区欧美日本亚洲| 欧美激情极品国产一区二区三区| 久久人妻福利社区极品人妻图片| 欧美在线一区亚洲| 国产精华一区二区三区| 99久久精品国产亚洲精品| 美女扒开内裤让男人捅视频| 久久人人爽av亚洲精品天堂| 国产精品乱码一区二三区的特点 | 视频在线观看一区二区三区| 国产精品美女特级片免费视频播放器 | 最近最新免费中文字幕在线| 法律面前人人平等表现在哪些方面| 成年版毛片免费区| 国产91精品成人一区二区三区| 在线观看午夜福利视频| 国产激情欧美一区二区| 国产野战对白在线观看| 夜夜爽天天搞| 久久久久久久久免费视频了| 天天躁夜夜躁狠狠躁躁| 满18在线观看网站| 亚洲国产精品久久男人天堂| 久久香蕉精品热| 黄片小视频在线播放| 真人一进一出gif抽搐免费| 啦啦啦韩国在线观看视频| 一区二区三区高清视频在线| av在线天堂中文字幕| 在线av久久热| 老熟妇乱子伦视频在线观看| 久久中文看片网| 精品国产亚洲在线| 午夜福利,免费看| 如日韩欧美国产精品一区二区三区| 我的亚洲天堂| 亚洲自拍偷在线| 国产成人精品久久二区二区91| svipshipincom国产片| 精品少妇一区二区三区视频日本电影| 美女高潮喷水抽搐中文字幕| 国产精品美女特级片免费视频播放器 | 91麻豆精品激情在线观看国产| 十八禁人妻一区二区| 午夜免费鲁丝| 天天躁狠狠躁夜夜躁狠狠躁| 免费搜索国产男女视频| 亚洲欧洲精品一区二区精品久久久| 亚洲熟女毛片儿| 亚洲熟妇熟女久久| 一进一出抽搐动态| 搞女人的毛片| 久久中文看片网| 别揉我奶头~嗯~啊~动态视频| 99国产精品99久久久久| 巨乳人妻的诱惑在线观看| 亚洲精华国产精华精| 国产aⅴ精品一区二区三区波| 亚洲一区二区三区色噜噜| av有码第一页| 成熟少妇高潮喷水视频| 黄网站色视频无遮挡免费观看| 变态另类丝袜制服| 自拍欧美九色日韩亚洲蝌蚪91| 黄色成人免费大全| 国产三级黄色录像| 精品久久久久久成人av| 悠悠久久av| 大陆偷拍与自拍| 午夜福利成人在线免费观看| 夜夜看夜夜爽夜夜摸| 国产黄a三级三级三级人| 国产成人av教育| 午夜福利影视在线免费观看| av视频免费观看在线观看| 在线观看免费日韩欧美大片| 热re99久久国产66热| 亚洲欧美日韩无卡精品| 亚洲色图 男人天堂 中文字幕| 精品高清国产在线一区| 日韩精品免费视频一区二区三区| 国产精华一区二区三区| 欧美 亚洲 国产 日韩一| 一区二区日韩欧美中文字幕| 一区在线观看完整版| 9191精品国产免费久久| 亚洲成人国产一区在线观看| 精品国产一区二区久久| 麻豆成人av在线观看| 久久久水蜜桃国产精品网| 日本vs欧美在线观看视频| 久久久久国产精品人妻aⅴ院| av视频免费观看在线观看| 99在线视频只有这里精品首页| 搡老熟女国产l中国老女人| 一a级毛片在线观看| 又大又爽又粗| 一本大道久久a久久精品| 国产高清激情床上av| 91国产中文字幕| 精品卡一卡二卡四卡免费| 在线观看一区二区三区| 久久香蕉激情| 国产精品1区2区在线观看.| 国产精品一区二区三区四区久久 | 多毛熟女@视频| 欧美亚洲日本最大视频资源| 欧美日韩乱码在线| 亚洲 欧美一区二区三区| 亚洲国产毛片av蜜桃av| 亚洲人成电影免费在线| 九色亚洲精品在线播放| 国内毛片毛片毛片毛片毛片| 最近最新中文字幕大全电影3 | 国产精品香港三级国产av潘金莲| 在线观看日韩欧美| 国产精品1区2区在线观看.| 欧美激情极品国产一区二区三区| 亚洲男人的天堂狠狠| 丝袜在线中文字幕| 一夜夜www| 精品第一国产精品| 免费少妇av软件| 亚洲五月婷婷丁香| 又黄又爽又免费观看的视频| av有码第一页| 久久人妻av系列| 久久欧美精品欧美久久欧美| 亚洲五月天丁香| 精品久久久久久久人妻蜜臀av | 非洲黑人性xxxx精品又粗又长| 国产精品国产高清国产av| 天堂动漫精品| www.精华液| 亚洲专区中文字幕在线| 亚洲精品久久国产高清桃花| 美国免费a级毛片| 国产欧美日韩一区二区三| tocl精华| 亚洲精品一卡2卡三卡4卡5卡| 国产亚洲欧美在线一区二区| 日韩欧美在线二视频| 女人被躁到高潮嗷嗷叫费观| 电影成人av| 真人做人爱边吃奶动态| 在线免费观看的www视频| 日本vs欧美在线观看视频| 一边摸一边做爽爽视频免费| 日本免费a在线| 黑人欧美特级aaaaaa片| 久久精品aⅴ一区二区三区四区| 麻豆国产av国片精品| 757午夜福利合集在线观看| 久久精品影院6| 午夜福利成人在线免费观看| 亚洲av成人不卡在线观看播放网| 99riav亚洲国产免费| 亚洲国产欧美网| 一进一出好大好爽视频| 成人三级做爰电影| 亚洲精品一区av在线观看| 视频在线观看一区二区三区| 欧美乱码精品一区二区三区| 亚洲男人天堂网一区| 精品久久蜜臀av无| 老汉色∧v一级毛片| 亚洲av熟女| 亚洲色图综合在线观看| 欧美激情 高清一区二区三区| 女人被狂操c到高潮| 激情视频va一区二区三区| 制服丝袜大香蕉在线| 黄色a级毛片大全视频| 免费少妇av软件| 可以在线观看的亚洲视频| 亚洲午夜精品一区,二区,三区| 国产午夜精品久久久久久| 欧美国产日韩亚洲一区| 国产xxxxx性猛交| 99riav亚洲国产免费| 精品久久久久久成人av| 高潮久久久久久久久久久不卡| 亚洲人成网站在线播放欧美日韩| 国产精品香港三级国产av潘金莲| 淫秽高清视频在线观看| 别揉我奶头~嗯~啊~动态视频| 欧美一级a爱片免费观看看 | 午夜精品在线福利| 日本三级黄在线观看| 欧美精品亚洲一区二区| 18禁国产床啪视频网站| 日韩大码丰满熟妇| 国产aⅴ精品一区二区三区波| 国产亚洲av高清不卡| 欧美黄色片欧美黄色片| 十八禁人妻一区二区| 淫妇啪啪啪对白视频| 国产99白浆流出| 19禁男女啪啪无遮挡网站| 国产国语露脸激情在线看| 中文字幕精品免费在线观看视频| 美女大奶头视频| 久久久国产成人免费| 久久国产精品影院| 99国产精品一区二区蜜桃av| 咕卡用的链子| 久久久精品欧美日韩精品| xxx96com| 久久中文看片网| 校园春色视频在线观看| 久久久久亚洲av毛片大全| 一区二区三区精品91| 亚洲va日本ⅴa欧美va伊人久久| 国产亚洲欧美在线一区二区| 免费观看精品视频网站| av视频免费观看在线观看| 少妇的丰满在线观看| 亚洲全国av大片| 日韩大码丰满熟妇| 岛国在线观看网站| 免费看美女性在线毛片视频| 精品卡一卡二卡四卡免费| 国产欧美日韩精品亚洲av| 午夜成年电影在线免费观看| 搡老岳熟女国产| 日韩一卡2卡3卡4卡2021年| 麻豆久久精品国产亚洲av| 国产精品 国内视频| 一本大道久久a久久精品| 国产精品日韩av在线免费观看 | 国产一区二区三区在线臀色熟女| 日日摸夜夜添夜夜添小说| 国产私拍福利视频在线观看| 国产av一区在线观看免费| 亚洲国产精品999在线| 国产国语露脸激情在线看| 亚洲精品在线美女| 一本久久中文字幕| 精品不卡国产一区二区三区| 满18在线观看网站| 国产精品国产高清国产av| 国内毛片毛片毛片毛片毛片| 国产亚洲精品一区二区www| 女生性感内裤真人,穿戴方法视频| 色尼玛亚洲综合影院| 精品久久久久久,| 亚洲五月天丁香| 日日夜夜操网爽| 美女高潮喷水抽搐中文字幕| 国产成人欧美在线观看| 激情视频va一区二区三区| 最近最新免费中文字幕在线| 一二三四在线观看免费中文在| 两个人视频免费观看高清| 国产亚洲精品一区二区www| 亚洲无线在线观看| 男男h啪啪无遮挡| 天堂√8在线中文| 黑人欧美特级aaaaaa片| 人人妻人人爽人人添夜夜欢视频| 国产aⅴ精品一区二区三区波| 婷婷六月久久综合丁香| 久久久久精品国产欧美久久久| 国产xxxxx性猛交| 欧美中文日本在线观看视频| 亚洲人成网站在线播放欧美日韩| 天天添夜夜摸| 俄罗斯特黄特色一大片| 国产精品,欧美在线| 久久天躁狠狠躁夜夜2o2o| 成年女人毛片免费观看观看9| 日韩精品中文字幕看吧| 久久久久久久午夜电影| 欧美国产日韩亚洲一区| 亚洲专区中文字幕在线| 亚洲七黄色美女视频| 91大片在线观看| 亚洲片人在线观看| 91av网站免费观看| 日日爽夜夜爽网站| 在线观看免费视频网站a站| 麻豆国产av国片精品| 亚洲va日本ⅴa欧美va伊人久久| 亚洲av美国av| 成人国语在线视频| 亚洲中文字幕一区二区三区有码在线看 | 一边摸一边做爽爽视频免费| 久久香蕉国产精品| 首页视频小说图片口味搜索| 亚洲三区欧美一区| 色综合站精品国产| 多毛熟女@视频| 久久精品91蜜桃| 亚洲无线在线观看| 久久午夜综合久久蜜桃| 日韩免费av在线播放| 久久人妻福利社区极品人妻图片| www.熟女人妻精品国产| 成人国产一区最新在线观看| 免费搜索国产男女视频| 国产日韩一区二区三区精品不卡| 午夜免费成人在线视频| 一边摸一边抽搐一进一小说| 高潮久久久久久久久久久不卡| 美女大奶头视频| 国产av又大| 激情在线观看视频在线高清| 91精品三级在线观看| 久久草成人影院| 欧美日韩黄片免| 男女之事视频高清在线观看| 国产在线精品亚洲第一网站| av在线播放免费不卡| 日本撒尿小便嘘嘘汇集6| 黄网站色视频无遮挡免费观看| 日韩高清综合在线| 欧美成人午夜精品| 国产av又大| 少妇粗大呻吟视频| 国产麻豆成人av免费视频| 久久人人97超碰香蕉20202| 三级毛片av免费| 嫩草影视91久久| 老司机午夜十八禁免费视频| 搡老岳熟女国产| 色av中文字幕| 波多野结衣高清无吗| 如日韩欧美国产精品一区二区三区| 久久久久久久久久久久大奶| 丰满人妻熟妇乱又伦精品不卡| 午夜福利18| 亚洲成av人片免费观看| 国产成+人综合+亚洲专区| 免费av毛片视频| 黄色片一级片一级黄色片| 日本 欧美在线| 亚洲中文日韩欧美视频| 免费av毛片视频| 日本撒尿小便嘘嘘汇集6| 亚洲三区欧美一区| 国产亚洲精品久久久久5区| 久久久久精品国产欧美久久久| av中文乱码字幕在线| 午夜福利视频1000在线观看 | 亚洲视频免费观看视频| 日本vs欧美在线观看视频| 欧美日韩黄片免| 桃红色精品国产亚洲av| 精品电影一区二区在线| 女人精品久久久久毛片| 午夜a级毛片| 级片在线观看| 亚洲熟妇中文字幕五十中出| 亚洲国产欧美网| 日本免费a在线| 激情在线观看视频在线高清| 亚洲五月婷婷丁香| 91字幕亚洲| 日本vs欧美在线观看视频| 国产成人系列免费观看| 久久久久精品国产欧美久久久| 88av欧美| 亚洲国产欧美网| 51午夜福利影视在线观看| 国产精品99久久99久久久不卡| 99riav亚洲国产免费| 欧美亚洲日本最大视频资源| aaaaa片日本免费| 国产在线观看jvid| 午夜成年电影在线免费观看| 国产精品一区二区三区四区久久 | 女性被躁到高潮视频| 久久精品91无色码中文字幕| 亚洲国产欧美日韩在线播放| 很黄的视频免费| 美女免费视频网站| av天堂久久9| 天堂√8在线中文| 亚洲av成人一区二区三| 69av精品久久久久久| 中文亚洲av片在线观看爽| av福利片在线| 国产成人系列免费观看| 中文字幕最新亚洲高清| 国产黄a三级三级三级人| 韩国av一区二区三区四区| 亚洲aⅴ乱码一区二区在线播放 | 久久国产精品影院| 中文字幕人妻熟女乱码| 久久国产乱子伦精品免费另类| 侵犯人妻中文字幕一二三四区| a在线观看视频网站| 每晚都被弄得嗷嗷叫到高潮| 夜夜夜夜夜久久久久| 亚洲精品久久成人aⅴ小说| 亚洲情色 制服丝袜| 久久久国产精品麻豆| 欧美性长视频在线观看| 这个男人来自地球电影免费观看| 国产欧美日韩一区二区三区在线| 成年女人毛片免费观看观看9| 国产免费av片在线观看野外av| 日韩大尺度精品在线看网址 | 亚洲欧美日韩另类电影网站| 精品熟女少妇八av免费久了| 搡老岳熟女国产| 中文字幕最新亚洲高清| 欧美日韩中文字幕国产精品一区二区三区 | 成人av一区二区三区在线看| 精品一区二区三区四区五区乱码| 99re在线观看精品视频| 9色porny在线观看| 亚洲av美国av| 精品一区二区三区av网在线观看| 琪琪午夜伦伦电影理论片6080| 超碰成人久久| 51午夜福利影视在线观看| 又黄又粗又硬又大视频| 人成视频在线观看免费观看| 国产成人欧美| 欧美黄色淫秽网站| 欧美激情 高清一区二区三区| 久久久久久久久中文| 久久久久久久精品吃奶| 亚洲国产欧美日韩在线播放| 老司机午夜福利在线观看视频| 国产精品,欧美在线| 看黄色毛片网站| 色播亚洲综合网| 免费在线观看影片大全网站| 久久久久国产一级毛片高清牌| 宅男免费午夜| 女人爽到高潮嗷嗷叫在线视频| 1024视频免费在线观看| 每晚都被弄得嗷嗷叫到高潮| 欧美日韩福利视频一区二区| 12—13女人毛片做爰片一| 波多野结衣一区麻豆| 女同久久另类99精品国产91| 母亲3免费完整高清在线观看| 又黄又粗又硬又大视频| 久久人妻熟女aⅴ| 久久精品亚洲精品国产色婷小说| 成人亚洲精品一区在线观看| 国产精品电影一区二区三区| 久久中文看片网| 最近最新中文字幕大全电影3 | 国产一级毛片七仙女欲春2 | 国产成人av教育| 亚洲国产精品合色在线| 精品免费久久久久久久清纯| 女人爽到高潮嗷嗷叫在线视频| 精品久久久久久久毛片微露脸| 精品久久久久久久人妻蜜臀av | 夜夜夜夜夜久久久久| 大型黄色视频在线免费观看| 免费在线观看黄色视频的| 亚洲电影在线观看av| 黄片大片在线免费观看| 999精品在线视频| 久久精品国产99精品国产亚洲性色 | 夜夜爽天天搞| 大型黄色视频在线免费观看| 久久伊人香网站| 久久精品人人爽人人爽视色| 国产乱人伦免费视频| 一级作爱视频免费观看| 日韩精品免费视频一区二区三区| 一级毛片精品| a级毛片在线看网站| 好男人电影高清在线观看| √禁漫天堂资源中文www| av有码第一页| 一级a爱视频在线免费观看| bbb黄色大片| 亚洲片人在线观看| 久久久久久久久中文| 中国美女看黄片| 亚洲伊人色综图| 久久午夜综合久久蜜桃| 91av网站免费观看| 色哟哟哟哟哟哟| 精品日产1卡2卡| 久久久水蜜桃国产精品网| 一边摸一边抽搐一进一出视频| 日韩有码中文字幕| 少妇被粗大的猛进出69影院| 丝袜美足系列| 午夜免费激情av| 国产野战对白在线观看| 国产成人av激情在线播放| 国产高清有码在线观看视频 | 中文字幕高清在线视频| 一边摸一边抽搐一进一小说| 欧美日本亚洲视频在线播放| 夜夜躁狠狠躁天天躁| av视频在线观看入口| 精品电影一区二区在线| 黄色丝袜av网址大全| 99国产综合亚洲精品| 久久久久久免费高清国产稀缺| 欧美一级毛片孕妇| 老汉色av国产亚洲站长工具| 欧美日韩精品网址| 免费不卡黄色视频| 天堂动漫精品| 亚洲av五月六月丁香网| 很黄的视频免费| 亚洲男人天堂网一区| www日本在线高清视频| 免费在线观看黄色视频的| av视频免费观看在线观看| 搡老熟女国产l中国老女人| 亚洲成av人片免费观看| av在线播放免费不卡| 色综合婷婷激情| 禁无遮挡网站| av在线播放免费不卡| 国产成人精品久久二区二区免费| 久久久久九九精品影院| 搡老岳熟女国产| 久久国产亚洲av麻豆专区| 又大又爽又粗| 巨乳人妻的诱惑在线观看| 麻豆一二三区av精品| 国产成人影院久久av| 亚洲欧美日韩无卡精品| 亚洲精品av麻豆狂野| 久久精品成人免费网站| 国产在线观看jvid| 波多野结衣巨乳人妻| 亚洲五月婷婷丁香| 一a级毛片在线观看| 9191精品国产免费久久| 精品欧美国产一区二区三| 在线观看免费视频网站a站| 在线观看舔阴道视频| 欧美色视频一区免费| 久久精品国产亚洲av香蕉五月| 搡老熟女国产l中国老女人| 麻豆国产av国片精品| 丰满的人妻完整版| 女人爽到高潮嗷嗷叫在线视频| 久久午夜亚洲精品久久| 久久精品aⅴ一区二区三区四区| 国产熟女xx| 99香蕉大伊视频| 久久天堂一区二区三区四区| 国产免费男女视频| 18禁黄网站禁片午夜丰满|