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

    裂隙介質(zhì)中Pu遷移的連續(xù)時間隨機游走模擬

    2024-01-06 02:04:44黃詩棋胡立堂劉東旭林劍鋒
    中國環(huán)境科學(xué) 2023年12期
    關(guān)鍵詞:模型

    黃詩棋,胡立堂*,劉東旭,林劍鋒,田 蕾,王 道

    裂隙介質(zhì)中Pu遷移的連續(xù)時間隨機游走模擬

    黃詩棋1,胡立堂1*,劉東旭2,林劍鋒2,田 蕾1,王 道1

    (1.北京師范大學(xué)水科學(xué)研究院,地下水污染控制與修復(fù)教育部工程研究中心,北京 100875;2.西北核技術(shù)研究所,陜西 西安 710024)

    為更準確地模擬核素的非菲克遷移行為并尋求一種更合理的建模方法,根據(jù)一維裂隙巖柱中Pu(钚)遷移試驗資料,利用基于不同聯(lián)合概率密度的連續(xù)時間隨機游走(CTRW)模型和傳統(tǒng)對流-彌散方程(ADE)對Pu遷移的穿透曲線(BTC)進行了過程模擬,對比分析了ADE模型、CTRW-TPL模型、CTRW-TPLA模型和CTRW-ETA模型的擬合程度及適用性,進而選擇出最優(yōu)的模型對不同膠體濃度下Pu在裂隙介質(zhì)中的遷移進行模擬研究.結(jié)果表明,相較于傳統(tǒng)ADE模型,參數(shù)擬合后考慮吸附作用的CTRW-TPL模型能有效地描述非菲克遷移行為特征,CTRW-ETA擬合過程更易實現(xiàn),但能為遷移行為的刻畫提供的信息較少.研究還發(fā)現(xiàn)膠體的存在顯著提高了Pu在裂隙介質(zhì)中的遷移性,但膠體濃度超過某個臨界值后,隨著濃度升高Pu的遷移性反而減弱,在CTRW理論框架中,該過程表現(xiàn)為表征廣義運移速度和彌散系數(shù)的參數(shù)值先增大后減小,而介質(zhì)對溶質(zhì)顆粒的平均吸附率先減小后增大.

    核素;非菲克遷移;連續(xù)時間隨機游走(CTRW);對流-彌散方程(ADE);Pu;膠體促進遷移

    深地質(zhì)處置是目前公認的較安全可靠的高放射性廢物處置方法,但在較大的時間尺度上將面臨處置庫退化導(dǎo)致放射性核素釋放進入地質(zhì)環(huán)境的潛在風險;核素若隨地下水遷移進入生物圈,將對生態(tài)環(huán)境和人類健康產(chǎn)生潛在危害,因此,需要研究核素在地下水中的遷移行為這一重要問題.利用數(shù)值模擬手段對核素遷移行為進行定量研究是一種重要的方法.然而由于地質(zhì)介質(zhì)的非均質(zhì)性且受到水解反應(yīng),絡(luò)合反應(yīng),氧化還原反應(yīng)以及膠體作用的共同影響,核素在地下水中的遷移通常呈現(xiàn)非菲克性.基于對流-彌散方程(ADE)的溶質(zhì)運移模型無法準確描述核素遷移過程中普遍存在的非菲克效應(yīng),盡管目前已有很多基于ADE的改進模型[1],如動水-不動水兩區(qū)模型和多重速率傳輸模型,能模擬物理-化學(xué)作用下的異常遷移,但作為確定性模型對復(fù)雜非均質(zhì)介質(zhì)中溶質(zhì)遷移問題的處理仍存在一定的局限性.隨機理論方法為解決地下水中的溶質(zhì)運移問題提供了另一種思路,連續(xù)時間隨機游走(CTRW)是一種概率模型方法[2],基于歐拉-拉格朗日方法,將溶質(zhì)在非均質(zhì)介質(zhì)中的遷移建模為無數(shù)粒子在空間和時間中的隨機游走,采用適宜的時間空間聯(lián)合概率密度函數(shù)(,)量化這種隨機過程,便可以描述粒子在對流-彌散以及與膠體?固相介質(zhì)發(fā)生復(fù)雜吸附/解吸等作用下從某空間位置經(jīng)時間到達距離處另一位置的遷移行為.示蹤劑在實驗室尺度土柱中遷移的經(jīng)典實驗及CTRW模擬結(jié)果[3-4]表明穿透曲線(BTCs)的不規(guī)律早到和拖尾現(xiàn)象在CTRW理論中能得到很好的解釋.近年來CTRW方法不斷被改進以應(yīng)用于更多場景:Margolin等[5]在CTRW框架中加入了平衡和動力吸附過程,可以處理運移過程中的吸附行為,使得CTRW模型能用于更多復(fù)雜機制下的遷移行為描述;在求解方法上,CTRW在偏微分方程解和粒子跟蹤(PT)兩方面得到了發(fā)展和應(yīng)用,兩種方法對非反應(yīng)性化學(xué)物質(zhì)的結(jié)果基本相同[6];此外,Ben-Zvi[7]給出相對更通用的CTRW微積分方程有限元解法(CTRW-FEM),適用于多維問題,以及化學(xué)反應(yīng)項引起的非線性方程求解;目前CTRW-PT方法通過與MODFLOW結(jié)合,已能被應(yīng)用于場地尺度存在的非菲克遷移的捕捉[8].但是基于這些改進的CTRW方法的應(yīng)用研究仍然很少,CTRW理論主要應(yīng)用于保守示蹤劑的遷移模擬中[9-12],也有少數(shù)應(yīng)用于重金屬離子[13]、有機污染物[14]和放射性核素[15]等水體污染物遷移的模擬研究.總體而言,CTRW在放射性核素遷移研究中的應(yīng)用較少.

    Pu(钚)是高放廢物處置研究的關(guān)注元素之一[16],屬于錒系元素,原子序數(shù)為94,其同位素的原子量從232~246不等.在劑量計算中,比較關(guān)注的同位素是238Pu,239Pu和242Pu,半衰期分別為87.8a,2.41′104a和3.73′105a.膠體對放射性核素遷移行為的影響相對重要[17],地下水中Pu的膠體存在形式主要包括真膠體和假膠體,真膠體主要由Pu(Ⅳ)發(fā)生水解或絡(luò)合反應(yīng)產(chǎn)生的,由于真膠體自身穩(wěn)定性相關(guān)的擴散作用及其動力吸附作用,其穿透曲線具有明顯的拖尾效應(yīng);假膠體是Pu吸附于地下水中各種天然膠體上形成的,Pu對膠體的吸附過程迅速,而吸附于膠體的Pu具有很低的解吸動力學(xué)常數(shù),這種緩慢的解吸過程使吸附在可遷移膠體上的Pu遷移能力增強.傳統(tǒng)的溶質(zhì)運移模型中Pu被視為遷移緩慢的難溶強吸附性元素,但實際的場地監(jiān)測結(jié)果表明,裂隙介質(zhì)能為Pu-膠體物種提供優(yōu)先通道,使得Pu在地下水中的實際遷移速度和距離遠高于計算值.因此本文選擇Pu作為代表性元素,以室內(nèi)一維裂隙柱試驗結(jié)果為例,分別建立ADE模型、CTRW-TPL(概率密度函數(shù)為截斷冪律模型的CTRW方法)模型、含吸附作用項的CTRW-TPL模型(在本文中簡記為CTRW-TPLA)以及CTRW-ETA (概率密度函數(shù)為修正指數(shù)模型的CTRW方法)模型,描述裂隙巖柱實驗中的Pu遷移行為,對實測的穿透曲線進行擬合,識別相應(yīng)的模型參數(shù),對比不同模型的擬合效果,然后在此基礎(chǔ)上選擇最優(yōu)的模型對不同膠體濃度下裂隙介質(zhì)中Pu遷移行為進行模擬,分析膠體對Pu遷移行為的影響以及相應(yīng)的參數(shù)變化規(guī)律.以期為放射性核素在裂隙介質(zhì)的遷移行為研究提供參考.

    雙桐巷旁邊的和平大戲院不見了,唯有青弋江畔的新百大廈尚在。這眼前的建筑也不知翻修了多少遍,唯有赤豆酒釀的滋味永恒不滅。我們?nèi)サ猛砹?,老奶奶牛肉面館已歇業(yè)關(guān)門。當年,她家的牛肉面辣得登峰造極——我獨喜歡門口小爐子里燜煮的鹵干。

    1 研究方法

    1.1 裂隙巖柱實驗

    本文進行數(shù)值模擬研究使用的數(shù)據(jù)來自Tran[18]和林劍鋒[19-20]等發(fā)表的裂隙巖柱實驗資料,分別用于分析不同模型對膠體作用下Pu在裂隙巖柱中遷移的穿透曲線(BTC)的擬合優(yōu)度,以及不同濃度的膠體對Pu遷移行為的影響及CTRW模型對應(yīng)參數(shù)的變化規(guī)律.

    1.2 對流-彌散模型(ADE)

    在均質(zhì)介質(zhì)和穩(wěn)定飽和流場的條件下,考慮平衡吸附和放射性衰變的一維對流-彌散方程為[21-22]:

    推進全域旅游是海南旅游發(fā)展的必然趨勢,必須突破產(chǎn)業(yè)界限,深入加強產(chǎn)業(yè)融合,促進旅游產(chǎn)業(yè)與各個產(chǎn)業(yè)的的融合發(fā)展,全域都是景區(qū),全域協(xié)同發(fā)展,保護完好的自然生態(tài)環(huán)境,沒有遭受破壞的當?shù)刂脖缓妥匀痪坝^,濃郁的熱帶民族風情和原汁原味的鄉(xiāng)土文化,熱情友善的當?shù)鼐用?,打造誠信可靠、便捷周到的服務(wù)體系。什寒村的營銷推廣必須充分利用瓊中縣整體的旅游營銷活動,通過節(jié)慶活動、旅游聯(lián)盟、媒體聯(lián)盟等多渠道宣傳,通過舉辦各類主題活動,比如"奔格內(nèi)"瓊中綠橙旅游季、好美瓊中·"奔格內(nèi)"鄉(xiāng)村旅游季、"三月三"黎苗文化旅游節(jié)、"鄉(xiāng)村旅游音樂節(jié)"等,不斷挖掘瓊中什寒旅游新熱點、新賣點,成功打造什寒鄉(xiāng)村休閑旅游品牌。

    CTRW模型是一種脫離ADE基本假設(shè)而更加具有一般性的隨機方法,在CTRW框架下多孔介質(zhì)和裂隙介質(zhì)中溶質(zhì)質(zhì)點在任意時間和空間位置的概率分布函數(shù),即溶質(zhì)濃度的時間與空間分布的拉普拉斯表達式為:

    1.自覺學(xué)習的習慣。在日常教學(xué)中,制定詳細可行的規(guī)則,要求學(xué)生在規(guī)定時間內(nèi)完成規(guī)定的任務(wù),并且及時進行反饋和適當激勵,引導(dǎo)學(xué)生逐步克服惰性,養(yǎng)成良好的學(xué)習習慣。

    式中:0為源層中核素的初始質(zhì)量,g;為體積流量,mL/h;L為浸出率系數(shù),h-1;為有效孔隙度,V為源層體積,cm3.

    1.3 連續(xù)時間隨機游走模型(CTRW)

    作為我國本土學(xué)者對跨文化能力模式的研究成果,對于我國國內(nèi)學(xué)者的研究具有特別的指導(dǎo)意義,也為筆者的研究提供了重要的理論參考。

    修正指數(shù)模型(ETA):

    漸近模型(ASY):

    截斷冪律模型(TPL):

    為了描述核素隨時間發(fā)生放射性衰變以及在水流作用下從源層不斷浸出而造成的源層核素釋放濃度減少,Baes[23]對淋濾過程中復(fù)雜的物理化學(xué)過程進行簡化,源層釋放濃度如下:

    “你跟那什么嬌,奧,唐嬌分手了嗎?”何守四看何西走了就過來問何北,還沒等何北回話,他又接著說“你可別給我鬧出這事兒來!”看三伯也走了過來,何北趕緊以上廁所為名溜了。

    2 應(yīng)用研究

    2.1 CTRW模型參數(shù)靈敏度分析

    根據(jù)表1所示的條件進行數(shù)值模擬實驗,參數(shù)、、2、和的靈敏度結(jié)果示于圖1.如圖1a所示,隨著增大?初始遷移速度增大,?穿透曲線的峰值濃度也相應(yīng)增大,?而隨著減小,BTC的不對稱性更強且拖尾現(xiàn)象更加明顯.?如圖1b,隨著值增大,BTC早到和拖尾現(xiàn)象更明顯,峰值濃度減小,而較小的導(dǎo)致穿透時間延遲但峰值濃度較高.如圖1c,冪律指數(shù)對BTC峰值濃度的影響較大,較大的值導(dǎo)致更大的峰值濃度,較小的值(<1)導(dǎo)致更明顯的拖尾現(xiàn)象,而穿透時間對的敏感度較小,不同的值下濃度到達峰值的時間相近.如圖1d, BTC對截斷時間2的敏感度低,當2>103后,2呈指數(shù)級增長時BTC卻基本無變化.如圖1e,穿透曲線對表征參與吸附和解吸作用的參數(shù)平均吸附率較敏感,值較大時遷移初始速度較低,峰值濃度小,而達到穿透時間以后,隨著的增大,拖尾的漸近坡度幾乎不再變化.如圖1f,參數(shù)表征冪律形式下的拖尾效應(yīng),值越大,遷移速度更快,拖尾效應(yīng)更明顯.

    表1 靈敏度分析數(shù)值實驗條件

    注:CTRW Toolbox計算中為無量綱單位,輸入時需分別除以土柱長度和2.

    圖1 CTRW模型參數(shù)靈敏度分析

    2.2 Pu運移穿透曲線的參數(shù)反演方法

    針對當前目的地感官營銷實踐面臨的問題,結(jié)合上述目的地感官營銷的3種運用方式分析,本文認為現(xiàn)階段目的地感官營銷研究所要重點解決的問題包括:

    ②敷貼藥物組方及配制方法:包括炒白芥子、炙甘遂、細辛、延胡索、肉桂、干姜、生麻黃、沉香、冰片、麝香,按4∶4∶2∶2∶2∶2∶2∶2∶1∶1配成,上述藥物研末并篩出細粉后用新鮮生姜汁、蜂蜜等溶劑調(diào)成膏狀,制成1 cm×1 cm×1 cm大小藥球,置于5 cm×5 cm的防過敏貼中央,為保證最大療效,敷貼皆為現(xiàn)場調(diào)制。

    ADE模型和CTRW模型模擬精度采用可決系數(shù)2和均方根誤差RMSE進行評價:

    李自成撤離北京后,永王與其兄逃出嘉定府失散。永王逃到河南鄉(xiāng)村耕田自食其力。清軍聞訊追捕到河南,永王又逃往安徽鳳陽地祖宗朱元璋的故鄉(xiāng),被王姓鄉(xiāng)紳收留,易名王士元。王鄉(xiāng)紳去世后王士元被其后人猜忌,無奈流浪到湖北大別山落發(fā)為僧,云游到野村谷搭茅庵修行,懸壺濟世。孰料當朝皇帝又追殺過來,再逃到浙江余姚,在一廢棄古廟棲身……

    3 結(jié)果與討論

    3.1 ADE和CTRW模型擬合效果對比分析

    圖2 基于ADE和CTRW-TPL模型擬合的Pu穿透曲線

    實測數(shù)據(jù)來自Tran等[18]

    實驗數(shù)據(jù)表明Pu遷移行為呈現(xiàn)強烈的非菲克性,實測BTC表現(xiàn)出早到和拖尾效應(yīng).早到現(xiàn)象的主要原因是含有優(yōu)勢裂隙的介質(zhì)為Pu膠體遷移提供了優(yōu)先流通道,大部分溶質(zhì)粒子在對流作用下迅速遷移到出口,導(dǎo)致溶質(zhì)穿透時間提前,而與Pu真膠體自身穩(wěn)定性相關(guān)的擴散作用和沉降過程導(dǎo)致了明顯的拖尾效應(yīng).ADE模型對穿透曲線的擬合結(jié)果如圖2所示,BTC的早到和拖尾現(xiàn)象難以被很好地刻畫,擬合的可決系數(shù)2小于0.8,均方根誤差RMSE大于0.3,這是由于傳統(tǒng)ADE模型的參數(shù)平均流速和彌散系數(shù)及阻滯系數(shù)R都具有尺度效應(yīng),是與空間有關(guān)的函數(shù)而非固定值.R本質(zhì)上是基于等溫吸附模型的分配系數(shù)得出的計算值,僅能反應(yīng)簡化的吸附過程,對實際問題中復(fù)雜的吸附-解吸?膠體遷移機制的描述有很大的局限.相比之下,CTRW-TPL模型擬合的可決系數(shù)2大于0.9,均方根誤差RMSE小于0.2,擬合精度更高,能更好地擬合不對稱的BTC,同時TPL模型還能通過參數(shù)可以描述菲克遷移向非菲克遷移轉(zhuǎn)變的時間.

    圖3 基于CTRW-TPL、CTRW-TPLA、CTRW-ETA模型擬合的Pu穿透曲線

    實測數(shù)據(jù)來自Tran等[18]

    圖4 ADE和CTRW模型對穿透曲線的擬合效果

    CTRW-TPLA模型對穿透曲線的擬合精度如圖4所示,2大于0.9,RMSE小于0.2,與CTRW-TPL一樣都能較好地擬合一維柱實驗穿透曲線,兩種模型反演識別的參數(shù)均小于1,揭示了柱實驗中Pu在膠體作用下遷移的高度非菲克性.CTRW-TPLA模型的參數(shù)平均吸附率擬合值較小,說明介質(zhì)對Pu的吸附作用較弱,這是因為相較于游離態(tài)Pu,裂隙介質(zhì)對Pu真膠體的吸附性有很大程度的降低.僅從擬合曲線和擬合精度來看,兩種模型都能得到較好的模擬結(jié)果,這可歸結(jié)于CTRW各個參數(shù)之間的內(nèi)部相關(guān)性,導(dǎo)致擬合參數(shù)組合不唯一,即使不考慮吸附作用最小化算法也能收斂到局部最小值,但在吸附作用較強烈的情況下可能使參數(shù)識別結(jié)果產(chǎn)生較大的誤差.

    CTRW-ETA模型對穿透曲線的擬合結(jié)果如圖3所示,2大于0.9,RMSE小于0.2,從擬合BTC和實測數(shù)據(jù)吻合程度上來說,ETA模型模擬效果較好,且由于模型參數(shù)很少,采用試錯法進行參數(shù)擬合比其他兩種模型更易實現(xiàn),但是單一的參數(shù)能為遷移行為的分析提供的信息有限.此外,ETA模型的參數(shù)處于0.06~0.1的范圍內(nèi)時均能得出多組相似的擬合結(jié)果,難以排除在局部收斂的參數(shù)組合.

    表2 ADE和CTRW模型擬合的參數(shù)值

    注:CXTFIT求解ADE模型參數(shù)在95%的置信區(qū)間內(nèi).

    3.2 膠體濃度對Pu遷移模擬的影響分析

    表3 CTRW-TPLA模型擬合的參數(shù)值

    圖5 基于CTRW-TPLA模型擬合的不同膠體濃度下Pu穿透曲線

    實測數(shù)據(jù)來自文獻[19-20]

    4)落花后10~15天,全樹噴高鈣液肥,每隔10~15天噴1次,共噴3~5次。可選高鈣液肥有黃腐酸生態(tài)鈣、甾醇鈣、全能鈣、硝酸鈣、氯化鈣等。

    4 結(jié)論

    4.1 傳統(tǒng)ADE模型對膠體作用下Pu遷移的穿透曲線模擬效果較差,而三種基于不同聯(lián)合概率密度函數(shù)(,)的CTRW模型都能實現(xiàn)較好的擬合,其中含吸附項的CTRW-TPLA的參數(shù)能較全面地描述非菲克遷移特征.CTRW建模的優(yōu)勢在于,在CTRW理論的框架下溶質(zhì)的遷移行為僅與粒子在介質(zhì)中的隨機游走有關(guān),通過建立合適的(,)能較好描述復(fù)雜物理-化學(xué)作用下溶質(zhì)運移的宏觀規(guī)律,且在實際應(yīng)用層面,通過(,)有重要特征的較少參數(shù)就可定量分析觀測數(shù)據(jù),可以簡化建模過程.

    “王九,我說,你瞧,大爺大姑娘不來,先生可來了。好,咱們動手,先生不會走的。你小心別讓趙四小子扔倒。先生幫咱們繃個場面,看你摔趙四這小子,先生準不走?!?/p>

    4.2 CTRW參數(shù)反演可能出現(xiàn)最小化算法在局部收斂的情況,導(dǎo)致擬合參數(shù)不唯一,因此本文先使用參數(shù)較少的ADE擬合吸附性很弱的保守示蹤劑氚的穿透曲線,得到與水流特征有關(guān)的參數(shù)和,盡管其物理意義不同于CTRW模型參數(shù),但對溶質(zhì)遷移行為的影響具有相似的規(guī)律,參考該結(jié)果和實驗條件來確定部分CTRW初始參數(shù)值,可以縮短試錯時間,并在一定程度上減小局部收斂造成的誤差.

    4.3 膠體的存在顯著提高了Pu在裂隙介質(zhì)中的遷移性,但膠體濃度超過某一閾值后,隨著濃度升高,遷移性反而降低,這可歸結(jié)于膠體濃度變大導(dǎo)致膠體碰撞概率增加從而發(fā)生沉積的可能性變大,這個過程中在CTRW參數(shù)中表現(xiàn)為值先增大后減小,介質(zhì)對溶質(zhì)顆粒的平均吸附率先減小后增大.目前采用的CTRW模型主要考慮了膠體的吸附作用,膠體對溶質(zhì)的作用機制仍待進一步深入研究.

    [1] 夏 源,吳吉春.張 勇.改進時間分數(shù)階模型模擬非Fick溶質(zhì)運移 [J]. 水科學(xué)進展, 2013,24(3):349-357.

    Xia Y, Wu J C, Zhang Y. Tempered time-fractional advection- dispersion equation for modeling non-Fickian transport [J]. Advances in Water Science, 2013,24(3):349-357.

    [2] Berkowitz B, Scher H. The role of probabilistic approaches to transport theory in heterogeneous media [J]. Transport in Porous Media, 2001,42(1):241-263.

    [3] Cortis A, Berkowitz B. Anomalous transport in "classical" soil and sand columns [J]. Soil Science Society of America Journal, 2004, 68(5):1539-1548.

    [4] Raveh-Rubin S, Edery Y, Dror I, et al. Nickel migration and retention dynamics in natural soil columns [J]. Water Resources Research, 2015,51(9):7702-7722.

    [5] Margolin G, Dentz M, Berkowitz B. Continuous time random walk and multirate mass transfer modeling of sorption [J]. Chemical Physics, 2003,295(1):71-80.

    [6] Dentz M, Cortis A, Scher H, et al. Time behavior of solute transport in heterogeneous media: transition from anomalous to normal transport [J]. Advances in Water Resources, 2004,27(2):155-173.

    [7] Ben-zvi R, Jiang S D, Scher H, et al. Finite-Element Method Solution of Non-Fickian Transport in Porous Media: The CTRW-FEM Package [J]. Groundwater, 2019,57(3):479-484.

    [8] Hansen S K, Berkowitz B. Aurora: A non-Fickian (and Fickian) particle tracking package for modeling groundwater contaminant transport with MODFLOW [J]. Environmental Modelling & Software, 2020,134:104871.

    [9] Muljadi B P, Bijeljic B, Blunt M J, et al. Modelling and upscaling of transport in carbonates during dissolution: Validation and calibration with NMR experiments [J]. Journal of Contaminant Hydrology, 2018, 212:85-95.

    [10] 孫立偉.一維河道穿透曲線的拖尾現(xiàn)象和多峰現(xiàn)象的模擬研究 [D]. 廣州:暨南大學(xué), 2019.

    Sun L W, Simulation study on tailing and multi-peaking phenomena of breakthrough curves in 1D channel [D]. Guangzhou: Jinan University, 2019.

    [11] Kurotori T, Zahasky C, Benson S M, et al.Description of chemical transport in laboratory rock cores using the continuous random walk formalism [J]. Water Resources Research, 2020,56(9):1-19.

    [12] Goeppert N, Goldscheider N, Berkowitz B. Experimental and modeling evidence of kilometer-scale anomalous tracer transport in an alpine karst aquifer [J]. Water Research, 2020,178:115755.

    [13] Li N, Ren L, Shang J. Application of continuous time random walk theory to simulate cadmium transport in the soil-water environment [J]. Journal of Hydraulic Engineering, 2012,43(8):906-913,925.

    [14] Li Y, Bian J, Wang Q, et al. Experiment and simulation of non- reactive solute transport in porous media [J]. Groundwater, 2022, 60(3):330-343.

    [15] 劉東旭,司高華,李 哲,等.非均質(zhì)條件下鍶遷移的反向隨機模擬 [J]. 原子能科學(xué)技術(shù), 2017,51(4):609-616.

    Liu D X, Si G H, Li Z, et al. Random inverse modelling of strontium migration in heterogeneous media [J]. Atomic Energy Science and Technology, 2017,51(4):609-616.

    [16] 司高華,于 靜,王青海,等.钚在粘土中的遷移實驗 [J]. 核化學(xué)與放射化學(xué), 2013,35(1):29-33.

    Si G H, Yu J, Wang Q H, et al. Migration experiment of Plutonium in clay [J]. Journal of Nuclear and Radiochemistry, 2013,35(1):29-33.

    [17] 劉 晨,方 升,許強偉.Am(Ⅲ)膠體穩(wěn)定性及其在膨潤土膠體上的吸附行為 [J]. 核化學(xué)與放射化學(xué), 2022,44(6):610-618.

    Liu C, Fang S, Xu W Q, Stability of Am(III) Colloid and its adsorption behavior on bentonite colloids [J]. Journal of Nuclear and Radiochemistry, 2022,44(6):610-618.

    [18] Tran E, Zavrin M, Kersting A B, et al. Colloid-facilitated transport of238Pu,233U and137Cs through fractured chalk: Laboratory experiments, modelling, and implications for nuclear waste disposal [J]. Science of the Total Environment, 2020,757:143818.

    [19] 林劍鋒.花崗巖裂隙介質(zhì)中α-FeOOH對钚遷移行為的影響研究[D]. 西安:西北核技術(shù)研究所, 2014.

    Lin J F, Effects of α-FeOOH on Plutonium transport in fractured granite media [D]. Xi’an:Northwest Institute of Nuclear Technology, 2014.

    [20] Lin J F, Dang H J, Xie J C, et al. Goethite colloid enhanced Pu transport through a single saturated fracture in granite [J]. Journal of Contaminant Hydrology, 2014,164:251-258.

    [21] 劉東旭.基于非菲克行為、膠體作用和地球化學(xué)演化的核素遷移模擬與反應(yīng)機制研究 [D]. 北京:北京師范大學(xué), 2020.

    Liu D X, Radionuclide transport modeling and reactive mechanisms based on non-Fickian behavior, Colloid Effects and Geochemical Evolution [D]. Beijing: Beijing Normal University, 2020.

    [22] Liu D X, Jivkov A P, Wang L C, et al. Non-Fickian dispersive transport of strontium in laboratory-scale columns: Modelling and evaluation [J]. Journal of Hydrology, 2017,549:1-11.

    [23] Baes C F, Sharp R D. A proposal for estimation of soil leaching and leaching constants for use in assessment models [J]. Journal of Environmental Quality, 1983,12(1):17-28.

    [24] Berkowitz B, Cortis A, Dentz M, et al. Modeling non-Fickian transport in geological formation as a continuous time random walk[J]. Reviews of Geophysics, 2006,44(2):RG2003(1-49).

    [25] Hansen S K, Berkowitz B. Interpretation and nonuniqueness of CTRW transition distributions: Insights from an alternative solute transport formulation [J]. Advances in Water Resources, 2014,74:54-63.

    [26] Toride N F, Leij F J, Genuchten M. The CXTFIT code for estimating transport parameters from laboratory or field tracer experiments [R]. California: U. S. Salinity Laboratory, 1995.

    [27] Cortis A, Berkowitz B. Computing anomalous contaminant transport in porous media: the CTRW MATLAB Toolbox [J]. Ground Water, 2010,43(6):947-950.

    [28] Simunek J, Genuchten M T V, Sejna M, et al. The STANMOD computer software for evaluating solute transport in porous media using analytical solutions of convection-dispersion equation [R]. Colorado: Colorado School of Mines International Ground Water Modeling Center, 1999.

    [29] 武貴斌,錢天偉,杜曉麗,等.非飽和帶核素遷移過程中滯留因子Rd意義的探討 [J]. 輻射防護通訊, 2003,(6):30-32.

    Wu G B, Qian T W, Du X L, et al. Discussion on Signification of Retardation Factor Rd in the Course of Nuclide Migration in Unsaturated Zone [J]. Radiation Protection Bulletin, 2003,(6):30-32.

    Transport simulation of Pu through fractured media based on continuous time random walk model.

    HUANG Shi-qi1, HU Li-tang1*, LIU Dong-xu2, LIN Jian-feng2, TIAN Lei1, WANG Dao1

    (1.College of Water Sciences, Beijing Normal University, Engineering Research Center of Groundwater Pollution Control and Remediation of Ministry of Education, Beijing 100875, China;2.Northwest Institute of Nuclear Technology, Xi'an 710024, China)., 2023,43(12):6321~6328

    Non-Fickian dispersion occurs in the migration of nuclides in groundwater, which is influenced by geological media heterogeneity and physicochemical processes. To improve the accuracy of the simulation of this phenomenon, an alternative model is required. To analyze Pu (plutonium) migration in one-dimensional fractured rock columns, we employed the continuous time random walk (CTRW) model together with the traditional advection-dispersion equation (ADE) to effectively represent Pu breakthrough curves (BTCs). The goodness of fit and feasibility of the ADE, CTRW-TPL, CTRW-TPLA, and CTRW-ETA models were later compared to identify the superior model for simulating Pu transport in fractured media with varying colloidal concentrations. The results indicated that the CTRW-TPL model, which incorporated adsorption, outperformed the traditional ADE in describing non-Fickian migration characteristics. Although the CTRW-ETA model lacked comprehensive migration characterization, its inversion process was more feasible. Furthermore, the migration capacity of Pu in fractured media was significantly enhanced due to colloidal effects. However, beyond a certain concentration threshold, the migration capacity of Pu decreased with increasing colloidal concentrations. Regarding the theoretical framework of the CTRW model, our analysis showed that the generalized transport velocity and dispersion coefficient first increased and then decreased, while the average adsorption rate of solute particles first decreased and then increased. These findings serve as a valuable reference for future research on the parameters of the CTRW model in fractured media.

    nuclides;non-Fickian;continuous time random walk (CTRW);advection-dispersion equation(ADE);Pu;colloid- facilitated transport

    X591

    A

    1000-6923(2023)12-6321-08

    黃詩棋,胡立堂,劉東旭,等.裂隙介質(zhì)中Pu遷移的連續(xù)時間隨機游走模擬 [J]. 中國環(huán)境科學(xué), 2023,43(12):6321-6328.

    Huang S Q, Hu L T, Liu D X, et al. Transport simulation of Pu through fractured media based on continuous time random walk model [J]. China Environmental Science, 2023,43(12):6321-6328.

    2023-03-29

    國家自然科學(xué)基金核技術(shù)創(chuàng)新聯(lián)合基金資助項目(U2167211)

    * 責任作者, 教授, litanghu@bnu.edu.cn

    黃詩棋(1999-),女,四川眉山人,北京師范大學(xué)碩士研究生,從事地下水數(shù)值模擬研究.發(fā)表論文1篇.sq_huang177@126.com.

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機模型
    提煉模型 突破難點
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    3D打印中的模型分割與打包
    亚洲中文字幕日韩| 国产精品 欧美亚洲| 大陆偷拍与自拍| videos熟女内射| 日韩中文字幕欧美一区二区 | 波多野结衣av一区二区av| 99国产精品一区二区三区| 国产精品三级大全| 国产一区二区在线观看av| 久久免费观看电影| 欧美在线一区亚洲| 欧美人与善性xxx| 一个人免费看片子| 人妻人人澡人人爽人人| 色婷婷av一区二区三区视频| 狠狠婷婷综合久久久久久88av| 久久毛片免费看一区二区三区| 国产在线一区二区三区精| 亚洲欧美中文字幕日韩二区| √禁漫天堂资源中文www| 欧美日韩视频精品一区| 你懂的网址亚洲精品在线观看| 捣出白浆h1v1| 搡老岳熟女国产| 亚洲精品日本国产第一区| 青春草亚洲视频在线观看| 国产97色在线日韩免费| 欧美精品一区二区免费开放| 熟女av电影| 又粗又硬又长又爽又黄的视频| av网站免费在线观看视频| 国产精品国产av在线观看| 国产黄色视频一区二区在线观看| 啦啦啦中文免费视频观看日本| 桃花免费在线播放| 欧美亚洲 丝袜 人妻 在线| 丝袜美腿诱惑在线| 久久亚洲国产成人精品v| 国产日韩欧美亚洲二区| 秋霞在线观看毛片| 亚洲精品日本国产第一区| 亚洲成人手机| 丁香六月欧美| 男女边吃奶边做爰视频| 老司机午夜十八禁免费视频| 一级片免费观看大全| 国产成人一区二区在线| 国产在视频线精品| 夫妻午夜视频| 一级毛片 在线播放| 亚洲美女黄色视频免费看| 欧美变态另类bdsm刘玥| 亚洲专区中文字幕在线| www.熟女人妻精品国产| 女人精品久久久久毛片| 亚洲美女黄色视频免费看| 叶爱在线成人免费视频播放| 18禁国产床啪视频网站| 亚洲精品在线美女| 国产一区二区 视频在线| 色婷婷av一区二区三区视频| 亚洲欧美精品自产自拍| 国产精品 国内视频| 1024视频免费在线观看| 少妇裸体淫交视频免费看高清 | 热re99久久精品国产66热6| 菩萨蛮人人尽说江南好唐韦庄| 日韩精品免费视频一区二区三区| 亚洲伊人色综图| 精品亚洲乱码少妇综合久久| 搡老岳熟女国产| 一个人免费看片子| 久热这里只有精品99| 久久久欧美国产精品| 精品一区二区三卡| 巨乳人妻的诱惑在线观看| 国产免费又黄又爽又色| 免费黄频网站在线观看国产| 男女午夜视频在线观看| 亚洲图色成人| 在线亚洲精品国产二区图片欧美| 欧美成人精品欧美一级黄| 国产成人欧美| 男人舔女人的私密视频| 婷婷色av中文字幕| 晚上一个人看的免费电影| 午夜精品国产一区二区电影| 国产精品一二三区在线看| 热re99久久精品国产66热6| 亚洲av在线观看美女高潮| 999精品在线视频| 国产片特级美女逼逼视频| 久久人人爽人人片av| 亚洲三区欧美一区| 亚洲国产精品成人久久小说| 日韩中文字幕欧美一区二区 | 狠狠婷婷综合久久久久久88av| 色精品久久人妻99蜜桃| 午夜av观看不卡| e午夜精品久久久久久久| 国产激情久久老熟女| 精品一区在线观看国产| 久久久久久亚洲精品国产蜜桃av| 老司机午夜十八禁免费视频| 精品一区二区三区四区五区乱码 | 欧美精品av麻豆av| 亚洲av电影在线观看一区二区三区| 国产日韩一区二区三区精品不卡| 久久热在线av| 日日摸夜夜添夜夜爱| 人妻人人澡人人爽人人| 一区二区三区乱码不卡18| 色婷婷久久久亚洲欧美| 免费看av在线观看网站| 国产成人欧美| 日韩一区二区三区影片| 久久久精品国产亚洲av高清涩受| 亚洲人成电影免费在线| 18禁裸乳无遮挡动漫免费视频| 777米奇影视久久| www.熟女人妻精品国产| 国精品久久久久久国模美| 欧美人与性动交α欧美软件| 国产精品一区二区精品视频观看| 国产亚洲一区二区精品| 国产深夜福利视频在线观看| 国产一区有黄有色的免费视频| 国产一卡二卡三卡精品| 亚洲av日韩精品久久久久久密 | 国产黄色免费在线视频| 欧美精品一区二区免费开放| 欧美黄色淫秽网站| av网站在线播放免费| 久久久久久久久免费视频了| 女人爽到高潮嗷嗷叫在线视频| 亚洲欧美精品综合一区二区三区| 国产亚洲一区二区精品| 嫩草影视91久久| av欧美777| 国产欧美亚洲国产| 国产女主播在线喷水免费视频网站| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩一卡2卡3卡4卡2021年| av天堂久久9| av电影中文网址| 一二三四在线观看免费中文在| 亚洲成人免费电影在线观看 | 亚洲美女黄色视频免费看| 一本色道久久久久久精品综合| 中文字幕高清在线视频| 侵犯人妻中文字幕一二三四区| 欧美成人精品欧美一级黄| 国产日韩欧美在线精品| 看免费成人av毛片| 国产女主播在线喷水免费视频网站| 亚洲一卡2卡3卡4卡5卡精品中文| 成人国语在线视频| 国产成人免费无遮挡视频| av不卡在线播放| 亚洲精品第二区| 97在线人人人人妻| 国产精品国产av在线观看| 90打野战视频偷拍视频| 久久天躁狠狠躁夜夜2o2o | 在线亚洲精品国产二区图片欧美| 欧美97在线视频| 中文字幕人妻丝袜一区二区| 一级毛片黄色毛片免费观看视频| 欧美成狂野欧美在线观看| bbb黄色大片| 国产片内射在线| 婷婷色综合www| 美女视频免费永久观看网站| av在线app专区| 久久亚洲精品不卡| 欧美老熟妇乱子伦牲交| 成在线人永久免费视频| 纵有疾风起免费观看全集完整版| 国产成人影院久久av| 成人黄色视频免费在线看| 精品一区在线观看国产| 99国产精品一区二区三区| 欧美日韩黄片免| 国产成人系列免费观看| 亚洲,欧美,日韩| 国产免费福利视频在线观看| 免费高清在线观看视频在线观看| 爱豆传媒免费全集在线观看| 日本av手机在线免费观看| 国产精品一区二区精品视频观看| 欧美av亚洲av综合av国产av| 性高湖久久久久久久久免费观看| 这个男人来自地球电影免费观看| 99香蕉大伊视频| 大码成人一级视频| 啦啦啦视频在线资源免费观看| 狂野欧美激情性bbbbbb| 老熟女久久久| 深夜精品福利| 久久精品国产亚洲av高清一级| 在线观看免费日韩欧美大片| av又黄又爽大尺度在线免费看| 亚洲欧美日韩高清在线视频 | 中文字幕人妻丝袜一区二区| 免费观看av网站的网址| 久久性视频一级片| 国产亚洲欧美精品永久| 日日夜夜操网爽| 亚洲精品成人av观看孕妇| 国产亚洲av片在线观看秒播厂| 黄片小视频在线播放| 波野结衣二区三区在线| 色精品久久人妻99蜜桃| 可以免费在线观看a视频的电影网站| 国产一区二区 视频在线| 热99国产精品久久久久久7| 亚洲精品在线美女| av欧美777| 亚洲av日韩在线播放| 国产成人av教育| 久久99热这里只频精品6学生| 久热爱精品视频在线9| 免费黄频网站在线观看国产| 国产欧美亚洲国产| 久久精品成人免费网站| 国产成人91sexporn| 黄片播放在线免费| 国产男女内射视频| 一级毛片女人18水好多 | 国产欧美日韩精品亚洲av| 亚洲av国产av综合av卡| av网站在线播放免费| 婷婷色综合www| 一区二区三区激情视频| 一区福利在线观看| 免费高清在线观看日韩| 最黄视频免费看| 国产人伦9x9x在线观看| 精品福利观看| 搡老岳熟女国产| 99国产综合亚洲精品| 99精品久久久久人妻精品| 久久国产精品人妻蜜桃| 亚洲国产精品成人久久小说| 国产亚洲精品第一综合不卡| 视频区欧美日本亚洲| 麻豆国产av国片精品| svipshipincom国产片| 久久精品久久精品一区二区三区| 丝袜美足系列| 女性被躁到高潮视频| 女性被躁到高潮视频| 在线观看人妻少妇| 亚洲精品美女久久久久99蜜臀 | 99久久综合免费| 18禁观看日本| 女性生殖器流出的白浆| 久久精品亚洲av国产电影网| 亚洲精品国产色婷婷电影| a级毛片黄视频| 99久久人妻综合| 久久精品国产亚洲av涩爱| 这个男人来自地球电影免费观看| 可以免费在线观看a视频的电影网站| 精品高清国产在线一区| 欧美亚洲日本最大视频资源| 深夜精品福利| 制服诱惑二区| av一本久久久久| 曰老女人黄片| avwww免费| 久久精品国产亚洲av高清一级| 免费不卡黄色视频| 啦啦啦中文免费视频观看日本| 日韩中文字幕欧美一区二区 | 在线av久久热| 一区二区三区乱码不卡18| av天堂在线播放| 欧美97在线视频| 国产av国产精品国产| 热re99久久国产66热| 两个人看的免费小视频| 伊人久久大香线蕉亚洲五| 国产精品99久久99久久久不卡| 午夜视频精品福利| 日韩熟女老妇一区二区性免费视频| 波多野结衣一区麻豆| 脱女人内裤的视频| 纯流量卡能插随身wifi吗| 国产三级黄色录像| 精品熟女少妇八av免费久了| 一级毛片女人18水好多 | 嫁个100分男人电影在线观看 | 另类亚洲欧美激情| 最新的欧美精品一区二区| av天堂久久9| 欧美成人午夜精品| 亚洲第一青青草原| 成人手机av| 亚洲 欧美一区二区三区| 黄色 视频免费看| 国产欧美亚洲国产| 亚洲欧美成人综合另类久久久| 精品少妇久久久久久888优播| 久久久久久久大尺度免费视频| 黄色视频不卡| 别揉我奶头~嗯~啊~动态视频 | 日本av手机在线免费观看| 婷婷成人精品国产| 国产一区二区激情短视频 | 久9热在线精品视频| 人人妻人人添人人爽欧美一区卜| 亚洲精品国产区一区二| 国产在线视频一区二区| 国产精品一国产av| 久久精品久久精品一区二区三区| 亚洲专区中文字幕在线| 成人免费观看视频高清| 精品一区二区三区四区五区乱码 | 国产在线一区二区三区精| 90打野战视频偷拍视频| 人人妻人人澡人人看| 中文字幕av电影在线播放| 亚洲图色成人| 国产成人精品久久久久久| 黑人猛操日本美女一级片| 99热全是精品| 亚洲第一青青草原| 亚洲一卡2卡3卡4卡5卡精品中文| 国精品久久久久久国模美| 亚洲av在线观看美女高潮| 久久人妻熟女aⅴ| 成人影院久久| 每晚都被弄得嗷嗷叫到高潮| 亚洲av日韩精品久久久久久密 | 精品高清国产在线一区| 精品人妻一区二区三区麻豆| 99精国产麻豆久久婷婷| 亚洲精品国产一区二区精华液| 亚洲欧美中文字幕日韩二区| 国产精品二区激情视频| 国产色视频综合| 久久久久久久久免费视频了| 亚洲免费av在线视频| 日日爽夜夜爽网站| 波多野结衣av一区二区av| 久久精品国产亚洲av高清一级| 大片电影免费在线观看免费| 国产黄色视频一区二区在线观看| 高清视频免费观看一区二区| 亚洲欧洲国产日韩| 丝袜美足系列| 欧美av亚洲av综合av国产av| 交换朋友夫妻互换小说| 9色porny在线观看| 国产无遮挡羞羞视频在线观看| 深夜精品福利| 99香蕉大伊视频| 精品卡一卡二卡四卡免费| 青春草视频在线免费观看| 十八禁人妻一区二区| 黄色毛片三级朝国网站| 国产野战对白在线观看| 人成视频在线观看免费观看| av一本久久久久| 纯流量卡能插随身wifi吗| 欧美xxⅹ黑人| 久久热在线av| 久久人人97超碰香蕉20202| 一级毛片黄色毛片免费观看视频| 亚洲国产中文字幕在线视频| 在线观看www视频免费| 大香蕉久久网| 大码成人一级视频| 一级黄片播放器| 宅男免费午夜| 91国产中文字幕| 亚洲成人手机| 成年动漫av网址| 成年人黄色毛片网站| netflix在线观看网站| 欧美97在线视频| 两性夫妻黄色片| 美女福利国产在线| 丝袜喷水一区| 免费高清在线观看视频在线观看| 9191精品国产免费久久| 亚洲成人免费av在线播放| 亚洲一码二码三码区别大吗| 啦啦啦在线观看免费高清www| 亚洲成人免费电影在线观看 | 成人18禁高潮啪啪吃奶动态图| 免费在线观看影片大全网站 | 精品国产一区二区久久| 无限看片的www在线观看| 蜜桃在线观看..| bbb黄色大片| 777久久人妻少妇嫩草av网站| 日韩精品免费视频一区二区三区| 一区在线观看完整版| 99国产精品一区二区三区| 亚洲男人天堂网一区| e午夜精品久久久久久久| 亚洲专区国产一区二区| 看十八女毛片水多多多| 色综合欧美亚洲国产小说| 欧美亚洲日本最大视频资源| 欧美xxⅹ黑人| 色视频在线一区二区三区| 亚洲欧美色中文字幕在线| 两性夫妻黄色片| 国产精品久久久久成人av| 人妻一区二区av| 亚洲免费av在线视频| 建设人人有责人人尽责人人享有的| 丝袜美腿诱惑在线| 色视频在线一区二区三区| 丁香六月天网| 欧美日韩视频精品一区| 国产片特级美女逼逼视频| 岛国毛片在线播放| 国产激情久久老熟女| 午夜免费成人在线视频| 久久天躁狠狠躁夜夜2o2o | 亚洲人成网站在线观看播放| 中文字幕人妻熟女乱码| 色婷婷久久久亚洲欧美| 亚洲久久久国产精品| 欧美国产精品一级二级三级| 乱人伦中国视频| 欧美在线一区亚洲| 真人做人爱边吃奶动态| 日韩欧美一区视频在线观看| 中文字幕人妻熟女乱码| 男人爽女人下面视频在线观看| 精品人妻熟女毛片av久久网站| 日本猛色少妇xxxxx猛交久久| 中文字幕制服av| 婷婷色麻豆天堂久久| 人妻人人澡人人爽人人| 中文字幕人妻丝袜一区二区| 国产激情久久老熟女| 老熟女久久久| 如日韩欧美国产精品一区二区三区| 亚洲欧美清纯卡通| 老汉色av国产亚洲站长工具| 男女边吃奶边做爰视频| 久久精品国产a三级三级三级| 大香蕉久久成人网| 老司机亚洲免费影院| 中文字幕制服av| 日韩欧美一区视频在线观看| 免费观看a级毛片全部| 一本大道久久a久久精品| 捣出白浆h1v1| 亚洲国产中文字幕在线视频| 久久人人97超碰香蕉20202| 1024视频免费在线观看| 亚洲成人手机| 宅男免费午夜| 日韩av免费高清视频| 男女午夜视频在线观看| 波多野结衣一区麻豆| a 毛片基地| 国产亚洲欧美在线一区二区| 一级毛片电影观看| 亚洲欧美色中文字幕在线| 别揉我奶头~嗯~啊~动态视频 | 亚洲欧美日韩另类电影网站| 国产成人精品无人区| 看免费av毛片| 精品人妻在线不人妻| 巨乳人妻的诱惑在线观看| 操美女的视频在线观看| 欧美成人精品欧美一级黄| 制服诱惑二区| 亚洲欧美日韩另类电影网站| 国产男女超爽视频在线观看| 久久影院123| 亚洲国产毛片av蜜桃av| 久久性视频一级片| 久久久久网色| 国产高清不卡午夜福利| 黄色a级毛片大全视频| 亚洲第一av免费看| 久久精品久久久久久噜噜老黄| 97精品久久久久久久久久精品| 美女脱内裤让男人舔精品视频| 欧美大码av| 一区福利在线观看| 十分钟在线观看高清视频www| 飞空精品影院首页| 亚洲五月婷婷丁香| 亚洲欧美中文字幕日韩二区| 少妇裸体淫交视频免费看高清 | 自拍欧美九色日韩亚洲蝌蚪91| 少妇粗大呻吟视频| 久久久久久久久久久久大奶| av电影中文网址| 性高湖久久久久久久久免费观看| 久久人人爽av亚洲精品天堂| 亚洲成人免费av在线播放| 老司机午夜十八禁免费视频| 亚洲国产欧美日韩在线播放| 丝袜人妻中文字幕| 我要看黄色一级片免费的| 女性生殖器流出的白浆| 50天的宝宝边吃奶边哭怎么回事| 人妻 亚洲 视频| 一级片'在线观看视频| 男女午夜视频在线观看| 麻豆国产av国片精品| 国产精品一区二区精品视频观看| 亚洲七黄色美女视频| 大片免费播放器 马上看| 国产精品二区激情视频| svipshipincom国产片| 免费少妇av软件| 少妇人妻久久综合中文| av天堂在线播放| 精品国产一区二区久久| 欧美成人精品欧美一级黄| 欧美成狂野欧美在线观看| 亚洲熟女毛片儿| 欧美精品人与动牲交sv欧美| 亚洲欧美一区二区三区久久| 久久精品人人爽人人爽视色| 黄色a级毛片大全视频| 人人妻人人爽人人添夜夜欢视频| 国产高清国产精品国产三级| 国产精品久久久人人做人人爽| 多毛熟女@视频| 人妻 亚洲 视频| 日本猛色少妇xxxxx猛交久久| 久久精品熟女亚洲av麻豆精品| av网站在线播放免费| 巨乳人妻的诱惑在线观看| 亚洲免费av在线视频| 老熟女久久久| 精品国产一区二区久久| 亚洲av成人不卡在线观看播放网 | 久久精品亚洲av国产电影网| 一区二区三区激情视频| 久久99精品国语久久久| 国产精品久久久久久精品电影小说| 精品国产一区二区三区久久久樱花| 国产成人a∨麻豆精品| 国产在线免费精品| 脱女人内裤的视频| 捣出白浆h1v1| 一区二区三区四区激情视频| 高清视频免费观看一区二区| 久久综合国产亚洲精品| 免费在线观看日本一区| 亚洲,欧美,日韩| xxxhd国产人妻xxx| 亚洲精品美女久久久久99蜜臀 | av欧美777| 久久久久国产精品人妻一区二区| 欧美日韩av久久| 久久中文字幕一级| 国产成人影院久久av| 中文字幕色久视频| 精品亚洲成a人片在线观看| 久久久精品国产亚洲av高清涩受| 中文字幕亚洲精品专区| 麻豆国产av国片精品| 自线自在国产av| 日本av免费视频播放| 成年av动漫网址| 在线观看免费午夜福利视频| 天天躁夜夜躁狠狠躁躁| 少妇人妻 视频| 亚洲精品中文字幕在线视频| 亚洲欧洲国产日韩| 亚洲国产精品国产精品| 国产成人免费观看mmmm| 黑人巨大精品欧美一区二区蜜桃| 性少妇av在线| 亚洲国产毛片av蜜桃av| 少妇人妻久久综合中文| av在线app专区| 成人国产av品久久久| 日本黄色日本黄色录像| 啦啦啦在线免费观看视频4| 精品少妇黑人巨大在线播放| 久久久久久久久免费视频了| 午夜福利,免费看| xxxhd国产人妻xxx| 日韩,欧美,国产一区二区三区| 91精品伊人久久大香线蕉| 国产一区二区激情短视频 | netflix在线观看网站| 波野结衣二区三区在线| 婷婷色综合www| 最近最新中文字幕大全免费视频 | 99精国产麻豆久久婷婷| 亚洲自偷自拍图片 自拍| 最近最新中文字幕大全免费视频 | 人体艺术视频欧美日本| 我要看黄色一级片免费的| 亚洲一码二码三码区别大吗| 精品福利观看| 中文欧美无线码| 美女视频免费永久观看网站| 视频区图区小说| 精品人妻1区二区| 啦啦啦在线观看免费高清www| 亚洲人成网站在线观看播放| 久久人人爽av亚洲精品天堂| 国产一区二区激情短视频 | 少妇猛男粗大的猛烈进出视频| 精品久久久久久久毛片微露脸 | 亚洲免费av在线视频| 精品一品国产午夜福利视频|