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

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

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

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

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

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

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

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

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

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

    Pu(钚)是高放廢物處置研究的關(guān)注元素之一[16],屬于錒系元素,原子序數(shù)為94,其同位素的原子量從232~246不等.在劑量計(jì)算中,比較關(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)的擴(kuò)散作用及其動力吸附作用,其穿透曲線具有明顯的拖尾效應(yīng);假膠體是Pu吸附于地下水中各種天然膠體上形成的,Pu對膠體的吸附過程迅速,而吸附于膠體的Pu具有很低的解吸動力學(xué)常數(shù),這種緩慢的解吸過程使吸附在可遷移膠體上的Pu遷移能力增強(qiáng).傳統(tǒng)的溶質(zhì)運(yùn)移模型中Pu被視為遷移緩慢的難溶強(qiáng)吸附性元素,但實(shí)際的場地監(jiān)測結(jié)果表明,裂隙介質(zhì)能為Pu-膠體物種提供優(yōu)先通道,使得Pu在地下水中的實(shí)際遷移速度和距離遠(yuǎn)高于計(jì)算值.因此本文選擇Pu作為代表性元素,以室內(nèi)一維裂隙柱試驗(yàn)結(jié)果為例,分別建立ADE模型、CTRW-TPL(概率密度函數(shù)為截?cái)鄡缏赡P偷腃TRW方法)模型、含吸附作用項(xiàng)的CTRW-TPL模型(在本文中簡記為CTRW-TPLA)以及CTRW-ETA (概率密度函數(shù)為修正指數(shù)模型的CTRW方法)模型,描述裂隙巖柱實(shí)驗(yàn)中的Pu遷移行為,對實(shí)測的穿透曲線進(jìn)行擬合,識別相應(yīng)的模型參數(shù),對比不同模型的擬合效果,然后在此基礎(chǔ)上選擇最優(yōu)的模型對不同膠體濃度下裂隙介質(zhì)中Pu遷移行為進(jìn)行模擬,分析膠體對Pu遷移行為的影響以及相應(yīng)的參數(shù)變化規(guī)律.以期為放射性核素在裂隙介質(zhì)的遷移行為研究提供參考.

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

    1 研究方法

    1.1 裂隙巖柱實(shí)驗(yàn)

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

    1.2 對流-彌散模型(ADE)

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

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

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

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

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

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

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

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

    漸近模型(ASY):

    截?cái)鄡缏赡P?TPL):

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

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

    2 應(yīng)用研究

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

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

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

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

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

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

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

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

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

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

    3 結(jié)果與討論

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

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

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

    實(shí)驗(yàn)數(shù)據(jù)表明Pu遷移行為呈現(xiàn)強(qiáng)烈的非菲克性,實(shí)測BTC表現(xiàn)出早到和拖尾效應(yīng).早到現(xiàn)象的主要原因是含有優(yōu)勢裂隙的介質(zhì)為Pu膠體遷移提供了優(yōu)先流通道,大部分溶質(zhì)粒子在對流作用下迅速遷移到出口,導(dǎo)致溶質(zhì)穿透時間提前,而與Pu真膠體自身穩(wěn)定性相關(guān)的擴(kuò)散作用和沉降過程導(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ù)得出的計(jì)算值,僅能反應(yīng)簡化的吸附過程,對實(shí)際問題中復(fù)雜的吸附-解吸?膠體遷移機(jī)制的描述有很大的局限.相比之下,CTRW-TPL模型擬合的可決系數(shù)2大于0.9,均方根誤差RMSE小于0.2,擬合精度更高,能更好地?cái)M合不對稱的BTC,同時TPL模型還能通過參數(shù)可以描述菲克遷移向非菲克遷移轉(zhuǎn)變的時間.

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

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

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

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

    CTRW-ETA模型對穿透曲線的擬合結(jié)果如圖3所示,2大于0.9,RMSE小于0.2,從擬合BTC和實(shí)測數(shù)據(jù)吻合程度上來說,ETA模型模擬效果較好,且由于模型參數(shù)很少,采用試錯法進(jìn)行參數(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í)測數(shù)據(jù)來自文獻(xiàn)[19-20]

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

    4 結(jié)論

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

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

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

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

    [1] 夏 源,吳吉春.張 勇.改進(jìn)時間分?jǐn)?shù)階模型模擬非Fick溶質(zhì)運(yùn)移 [J]. 水科學(xué)進(jìn)展, 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ī)模擬 [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] 司高華,于 靜,王青海,等.钚在粘土中的遷移實(shí)驗(yàn) [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] 劉 晨,方 升,許強(qiáng)偉.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)機(jī)制研究 [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]. 輻射防護(hù)通訊, 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ī)游走模擬 [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)合基金資助項(xiàng)目(U2167211)

    * 責(zé)任作者, 教授, litanghu@bnu.edu.cn

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

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機(jī)模型
    提煉模型 突破難點(diǎn)
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達(dá)及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    国产私拍福利视频在线观看| 国产亚洲精品一区二区www| 窝窝影院91人妻| 中文字幕最新亚洲高清| 欧美大码av| 国内精品久久久久精免费| 男男h啪啪无遮挡| 精品久久久久久久久久免费视频| 男女那种视频在线观看| 亚洲 欧美一区二区三区| 中文字幕av在线有码专区| 久久久久国产精品人妻aⅴ院| 精品福利观看| 在线观看www视频免费| 大型黄色视频在线免费观看| 日韩欧美国产一区二区入口| 97碰自拍视频| 国产成人系列免费观看| 中文字幕人成人乱码亚洲影| 国产精品久久久人人做人人爽| 久久久久久久久久黄片| 亚洲国产欧美人成| 91av网站免费观看| 精品久久久久久久末码| 99riav亚洲国产免费| 麻豆成人午夜福利视频| 桃色一区二区三区在线观看| 午夜福利成人在线免费观看| 亚洲国产中文字幕在线视频| 久久国产精品人妻蜜桃| 亚洲精品av麻豆狂野| 亚洲专区字幕在线| 午夜激情av网站| 欧美日韩黄片免| 欧美又色又爽又黄视频| 女生性感内裤真人,穿戴方法视频| 一进一出好大好爽视频| 亚洲18禁久久av| 婷婷六月久久综合丁香| 日韩 欧美 亚洲 中文字幕| 午夜影院日韩av| 丁香欧美五月| 少妇被粗大的猛进出69影院| 国产精品,欧美在线| 日韩欧美一区二区三区在线观看| 国产精品久久久人人做人人爽| 国产一区二区在线观看日韩 | 黄色片一级片一级黄色片| 窝窝影院91人妻| 一区二区三区高清视频在线| 黄色视频,在线免费观看| 非洲黑人性xxxx精品又粗又长| 日日爽夜夜爽网站| 男人舔奶头视频| videosex国产| 亚洲人与动物交配视频| 特级一级黄色大片| 午夜a级毛片| 99国产精品一区二区蜜桃av| 校园春色视频在线观看| 一二三四在线观看免费中文在| 亚洲美女视频黄频| 一本久久中文字幕| 久久久精品大字幕| 久久久久国产精品人妻aⅴ院| 熟女电影av网| 欧美三级亚洲精品| 久久欧美精品欧美久久欧美| av中文乱码字幕在线| 国产成人系列免费观看| 成人三级黄色视频| 97人妻精品一区二区三区麻豆| av天堂在线播放| 亚洲熟女毛片儿| 少妇粗大呻吟视频| 巨乳人妻的诱惑在线观看| 亚洲av美国av| 我要搜黄色片| 日本 欧美在线| 国产熟女午夜一区二区三区| 很黄的视频免费| 法律面前人人平等表现在哪些方面| 看片在线看免费视频| 欧美性长视频在线观看| 精品福利观看| 日韩精品青青久久久久久| 亚洲美女视频黄频| 久久精品综合一区二区三区| 久久精品91无色码中文字幕| 欧美日韩一级在线毛片| 色av中文字幕| 国产激情偷乱视频一区二区| 免费无遮挡裸体视频| 天天躁夜夜躁狠狠躁躁| 一a级毛片在线观看| 九色成人免费人妻av| 国产成人精品久久二区二区91| 好男人在线观看高清免费视频| 亚洲成a人片在线一区二区| √禁漫天堂资源中文www| 精品高清国产在线一区| 在线观看66精品国产| 黄色 视频免费看| 国产真人三级小视频在线观看| 九色成人免费人妻av| 最近最新中文字幕大全免费视频| 亚洲成av人片免费观看| 国产单亲对白刺激| 脱女人内裤的视频| 手机成人av网站| 国产野战对白在线观看| 亚洲,欧美精品.| 性色av乱码一区二区三区2| 欧美人与性动交α欧美精品济南到| 欧美一级毛片孕妇| 黄色毛片三级朝国网站| 女同久久另类99精品国产91| 亚洲成人精品中文字幕电影| 十八禁网站免费在线| 国产精品 国内视频| aaaaa片日本免费| 久久热在线av| 欧美性长视频在线观看| 亚洲全国av大片| 日韩欧美一区二区三区在线观看| av在线天堂中文字幕| 久久九九热精品免费| 一卡2卡三卡四卡精品乱码亚洲| 国产在线精品亚洲第一网站| 国产成人精品久久二区二区免费| 天天一区二区日本电影三级| 国产一区二区三区视频了| 国产一区二区在线观看日韩 | 欧美国产日韩亚洲一区| 欧美日本视频| 动漫黄色视频在线观看| 国产三级在线视频| 亚洲国产欧洲综合997久久,| 欧美在线一区亚洲| 精品一区二区三区四区五区乱码| xxx96com| 亚洲国产精品合色在线| 久久精品人妻少妇| 草草在线视频免费看| av在线天堂中文字幕| 99在线视频只有这里精品首页| 99热这里只有精品一区 | 99国产极品粉嫩在线观看| 亚洲成人国产一区在线观看| 51午夜福利影视在线观看| 99久久国产精品久久久| 日本在线视频免费播放| 欧美乱色亚洲激情| 午夜福利高清视频| 亚洲男人天堂网一区| 日韩av在线大香蕉| 欧美黄色淫秽网站| 一级毛片精品| 久久精品91蜜桃| av欧美777| 日韩欧美免费精品| 国产麻豆成人av免费视频| 久久午夜综合久久蜜桃| 天天一区二区日本电影三级| 国产三级中文精品| 1024香蕉在线观看| av视频在线观看入口| aaaaa片日本免费| 亚洲精品av麻豆狂野| 国产精品久久久久久久电影 | 国产精品野战在线观看| 国产精品免费视频内射| www日本黄色视频网| 露出奶头的视频| 国产av又大| 天堂影院成人在线观看| 欧美激情久久久久久爽电影| 2021天堂中文幕一二区在线观| 国产亚洲欧美在线一区二区| 午夜精品在线福利| 国产久久久一区二区三区| 成人永久免费在线观看视频| 又粗又爽又猛毛片免费看| 亚洲精品美女久久av网站| 叶爱在线成人免费视频播放| 中国美女看黄片| 成熟少妇高潮喷水视频| 国产午夜精品论理片| 久久精品人妻少妇| 国产av不卡久久| 毛片女人毛片| 国产片内射在线| 亚洲午夜精品一区,二区,三区| 色在线成人网| 久久中文字幕人妻熟女| 一进一出抽搐动态| 怎么达到女性高潮| 亚洲熟女毛片儿| 香蕉av资源在线| 日本 欧美在线| 日本三级黄在线观看| 久久久久国产精品人妻aⅴ院| 国产成人系列免费观看| 久久中文字幕一级| 国内精品久久久久久久电影| 国产激情偷乱视频一区二区| 色哟哟哟哟哟哟| 亚洲男人天堂网一区| 一本久久中文字幕| 午夜久久久久精精品| 黄色成人免费大全| 一区福利在线观看| 国产成人精品久久二区二区91| 老司机午夜福利在线观看视频| 欧美乱码精品一区二区三区| 亚洲 国产 在线| 国产午夜精品久久久久久| 久久久国产成人免费| 人妻丰满熟妇av一区二区三区| 熟女电影av网| 丰满人妻一区二区三区视频av | 精品久久久久久久人妻蜜臀av| 黑人巨大精品欧美一区二区mp4| 国产免费av片在线观看野外av| 大型av网站在线播放| 成熟少妇高潮喷水视频| 久久久久国产精品人妻aⅴ院| 一级a爱片免费观看的视频| 免费看美女性在线毛片视频| e午夜精品久久久久久久| 99久久99久久久精品蜜桃| 亚洲真实伦在线观看| 激情 狠狠 欧美| 色视频www国产| 欧美激情国产日韩精品一区| 亚洲欧美日韩无卡精品| 午夜亚洲福利在线播放| 一级av片app| 九草在线视频观看| 男人舔女人下体高潮全视频| 女人十人毛片免费观看3o分钟| 校园春色视频在线观看| 三级经典国产精品| 亚洲综合色惰| 美女内射精品一级片tv| 99久久成人亚洲精品观看| 女同久久另类99精品国产91| 色吧在线观看| 99久国产av精品国产电影| 啦啦啦观看免费观看视频高清| 日韩欧美 国产精品| 国产一区二区亚洲精品在线观看| 精品人妻熟女av久视频| 中文字幕av在线有码专区| 欧美又色又爽又黄视频| 色综合站精品国产| 非洲黑人性xxxx精品又粗又长| 国产亚洲精品av在线| 亚洲成人精品中文字幕电影| 日韩成人av中文字幕在线观看| 国产精品久久久久久久久免| 人妻制服诱惑在线中文字幕| 欧美bdsm另类| 欧美成人a在线观看| 尤物成人国产欧美一区二区三区| 在线观看免费视频日本深夜| 丰满人妻一区二区三区视频av| 一本一本综合久久| 久久久久久久久中文| 亚洲乱码一区二区免费版| 搡老妇女老女人老熟妇| 成人亚洲精品av一区二区| 亚洲成人中文字幕在线播放| av在线播放精品| 亚洲国产精品久久男人天堂| 插逼视频在线观看| 中文精品一卡2卡3卡4更新| 欧美色视频一区免费| 国产成人精品一,二区 | 卡戴珊不雅视频在线播放| 日本三级黄在线观看| 日本黄色视频三级网站网址| 亚洲熟妇中文字幕五十中出| 日韩 亚洲 欧美在线| 91精品一卡2卡3卡4卡| av免费观看日本| 国产黄片视频在线免费观看| 啦啦啦啦在线视频资源| 中国美白少妇内射xxxbb| 国产真实伦视频高清在线观看| 女人十人毛片免费观看3o分钟| 国产av一区在线观看免费| 不卡一级毛片| 一级毛片电影观看 | 国产精品人妻久久久影院| 乱系列少妇在线播放| 亚洲性久久影院| 亚洲欧洲国产日韩| 99久国产av精品国产电影| 赤兔流量卡办理| 欧美日韩在线观看h| 91久久精品国产一区二区成人| 午夜久久久久精精品| 日韩欧美 国产精品| 成人漫画全彩无遮挡| 亚洲国产精品国产精品| 亚洲欧美日韩东京热| 毛片一级片免费看久久久久| 中国美女看黄片| 国产亚洲av嫩草精品影院| 夜夜爽天天搞| 69av精品久久久久久| 国产淫片久久久久久久久| 欧美高清性xxxxhd video| 日韩在线高清观看一区二区三区| 特级一级黄色大片| 成人永久免费在线观看视频| www.av在线官网国产| 久久久久免费精品人妻一区二区| 欧美人与善性xxx| 国产乱人视频| av又黄又爽大尺度在线免费看 | 插逼视频在线观看| 亚洲精品亚洲一区二区| 男女边吃奶边做爰视频| 熟女人妻精品中文字幕| 91精品一卡2卡3卡4卡| 少妇高潮的动态图| 国产午夜福利久久久久久| 久久国产乱子免费精品| 国产精品爽爽va在线观看网站| 日本欧美国产在线视频| 亚洲不卡免费看| 国内精品一区二区在线观看| 亚洲真实伦在线观看| 欧美一区二区国产精品久久精品| 最后的刺客免费高清国语| 亚洲不卡免费看| 成年女人看的毛片在线观看| 国产一区二区三区av在线 | 国产成人午夜福利电影在线观看| 国产伦在线观看视频一区| 日本黄大片高清| 日日撸夜夜添| 精品免费久久久久久久清纯| 熟女人妻精品中文字幕| 丝袜喷水一区| 国产精品.久久久| 寂寞人妻少妇视频99o| eeuss影院久久| 天美传媒精品一区二区| 午夜福利成人在线免费观看| 免费人成在线观看视频色| 日本成人三级电影网站| 亚洲欧美精品自产自拍| 97在线视频观看| 亚洲国产色片| 亚洲人成网站在线观看播放| 在线观看av片永久免费下载| 国产男人的电影天堂91| 亚洲四区av| 99精品在免费线老司机午夜| 亚洲成人久久性| 久久这里只有精品中国| 免费看av在线观看网站| 熟妇人妻久久中文字幕3abv| 99国产极品粉嫩在线观看| 久久久久久大精品| 亚洲av免费高清在线观看| 久久人人爽人人片av| 日韩大尺度精品在线看网址| 国产亚洲精品久久久久久毛片| 国产伦精品一区二区三区视频9| 一区二区三区免费毛片| 中文字幕人妻熟人妻熟丝袜美| 精品不卡国产一区二区三区| 搡女人真爽免费视频火全软件| 日日干狠狠操夜夜爽| 女的被弄到高潮叫床怎么办| 国产精品一区二区三区四区免费观看| 国产精品av视频在线免费观看| 国产精品精品国产色婷婷| 秋霞在线观看毛片| 欧美+日韩+精品| 搡女人真爽免费视频火全软件| 美女cb高潮喷水在线观看| 亚洲美女搞黄在线观看| 亚洲国产精品合色在线| 2021天堂中文幕一二区在线观| 亚洲av第一区精品v没综合| 最近的中文字幕免费完整| 日韩视频在线欧美| 午夜福利高清视频| 99九九线精品视频在线观看视频| 中文字幕制服av| 在线观看免费视频日本深夜| 久久亚洲精品不卡| 日韩强制内射视频| 国产成人a∨麻豆精品| 国产久久久一区二区三区| 日本成人三级电影网站| 国产成人精品婷婷| 非洲黑人性xxxx精品又粗又长| 亚洲丝袜综合中文字幕| 天堂中文最新版在线下载 | 51国产日韩欧美| 九九爱精品视频在线观看| 日本熟妇午夜| 日韩制服骚丝袜av| 女人十人毛片免费观看3o分钟| 非洲黑人性xxxx精品又粗又长| 免费观看精品视频网站| 亚洲电影在线观看av| 亚洲av男天堂| 少妇猛男粗大的猛烈进出视频 | 国产精品嫩草影院av在线观看| 啦啦啦观看免费观看视频高清| 亚洲成人久久性| 亚州av有码| 给我免费播放毛片高清在线观看| 久久久精品大字幕| 一区二区三区高清视频在线| 一级毛片电影观看 | 男的添女的下面高潮视频| 免费电影在线观看免费观看| 禁无遮挡网站| 国产在视频线在精品| 人人妻人人澡人人爽人人夜夜 | 日韩欧美 国产精品| 男女边吃奶边做爰视频| 久久精品国产亚洲网站| 一个人看视频在线观看www免费| 久久精品久久久久久噜噜老黄 | 男女下面进入的视频免费午夜| 午夜福利在线观看吧| 久久久久久久久中文| 成人国产麻豆网| 国产淫片久久久久久久久| 国产成人a区在线观看| 精品熟女少妇av免费看| 亚洲一级一片aⅴ在线观看| 午夜免费激情av| 成人综合一区亚洲| 国产午夜精品久久久久久一区二区三区| 日韩av在线大香蕉| 国产精品电影一区二区三区| 日韩欧美 国产精品| 亚洲精品粉嫩美女一区| 午夜激情福利司机影院| 老女人水多毛片| 天堂av国产一区二区熟女人妻| 国产午夜精品久久久久久一区二区三区| 亚洲高清免费不卡视频| av专区在线播放| 在线观看一区二区三区| 麻豆国产av国片精品| 国产亚洲精品久久久com| 欧美变态另类bdsm刘玥| 国内揄拍国产精品人妻在线| 精品人妻熟女av久视频| 婷婷亚洲欧美| 少妇的逼好多水| 在线播放国产精品三级| 国产乱人视频| 成人无遮挡网站| 波多野结衣巨乳人妻| 色噜噜av男人的天堂激情| 国产精品嫩草影院av在线观看| 特级一级黄色大片| av又黄又爽大尺度在线免费看 | 国内揄拍国产精品人妻在线| 日本一二三区视频观看| 国产亚洲精品av在线| 国产精品一及| 日本撒尿小便嘘嘘汇集6| 在线观看66精品国产| 国产视频首页在线观看| 别揉我奶头 嗯啊视频| 久久久午夜欧美精品| 欧美成人一区二区免费高清观看| 免费av观看视频| 菩萨蛮人人尽说江南好唐韦庄 | 精品一区二区免费观看| a级毛色黄片| 亚洲中文字幕日韩| 国国产精品蜜臀av免费| 91精品国产九色| 国产高潮美女av| 日日摸夜夜添夜夜爱| 联通29元200g的流量卡| 亚洲成人中文字幕在线播放| 亚洲欧美清纯卡通| 国产伦一二天堂av在线观看| 高清日韩中文字幕在线| 国产日本99.免费观看| 熟女电影av网| 午夜亚洲福利在线播放| 亚洲三级黄色毛片| av免费观看日本| 中文字幕精品亚洲无线码一区| 欧美最新免费一区二区三区| 久久综合国产亚洲精品| 成年av动漫网址| 国产三级中文精品| 亚洲五月天丁香| 看黄色毛片网站| 色吧在线观看| 国产久久久一区二区三区| 日韩欧美精品v在线| 国产中年淑女户外野战色| 色哟哟哟哟哟哟| 成人无遮挡网站| 成人永久免费在线观看视频| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品国产成人久久av| 一卡2卡三卡四卡精品乱码亚洲| 91午夜精品亚洲一区二区三区| 免费搜索国产男女视频| 少妇的逼水好多| 久久精品夜夜夜夜夜久久蜜豆| 精品人妻熟女av久视频| 身体一侧抽搐| 最后的刺客免费高清国语| 亚洲色图av天堂| 青青草视频在线视频观看| 成人一区二区视频在线观看| 丰满的人妻完整版| 麻豆av噜噜一区二区三区| 亚洲一区二区三区色噜噜| 亚洲婷婷狠狠爱综合网| 国产成人一区二区在线| 久久精品影院6| 大又大粗又爽又黄少妇毛片口| 亚洲美女视频黄频| 精品日产1卡2卡| 国产真实伦视频高清在线观看| 久久久国产成人精品二区| 亚洲精品日韩在线中文字幕 | h日本视频在线播放| 亚洲av第一区精品v没综合| 欧美日本视频| 男女边吃奶边做爰视频| 成人美女网站在线观看视频| 国产精品一及| 国产久久久一区二区三区| 国产一区二区在线观看日韩| 人妻系列 视频| 成人毛片a级毛片在线播放| 老司机影院成人| 中文亚洲av片在线观看爽| 在线观看午夜福利视频| 国产黄片美女视频| 91久久精品国产一区二区成人| 国产亚洲精品av在线| 久久久久久大精品| 亚洲av成人精品一区久久| 久久婷婷人人爽人人干人人爱| 女的被弄到高潮叫床怎么办| 国产极品天堂在线| 秋霞在线观看毛片| 欧美+亚洲+日韩+国产| 可以在线观看的亚洲视频| 91久久精品电影网| 男插女下体视频免费在线播放| www.av在线官网国产| 日本爱情动作片www.在线观看| 亚洲无线观看免费| 噜噜噜噜噜久久久久久91| 麻豆一二三区av精品| 国产精品三级大全| 国产欧美日韩精品一区二区| 长腿黑丝高跟| 免费观看精品视频网站| 久久久久久久久久黄片| 国产av一区在线观看免费| 少妇被粗大猛烈的视频| 国产精品蜜桃在线观看 | 欧美精品国产亚洲| 国产成人精品一,二区 | 桃色一区二区三区在线观看| 精品人妻偷拍中文字幕| 国产精品久久久久久精品电影小说 | 三级国产精品欧美在线观看| 性色avwww在线观看| 在线观看66精品国产| 毛片一级片免费看久久久久| 色综合色国产| 91久久精品国产一区二区成人| 天天一区二区日本电影三级| 在线观看美女被高潮喷水网站| 日韩成人伦理影院| av在线天堂中文字幕| 亚洲性久久影院| 欧美性感艳星| 日韩成人伦理影院| 欧美性感艳星| 特级一级黄色大片| 日本av手机在线免费观看| 国产免费男女视频| 精品免费久久久久久久清纯| 伊人久久精品亚洲午夜| а√天堂www在线а√下载| 久久精品国产亚洲av涩爱 | 精品国内亚洲2022精品成人| АⅤ资源中文在线天堂| 国产伦精品一区二区三区四那| a级毛片免费高清观看在线播放| 午夜免费男女啪啪视频观看| 国产一区二区三区在线臀色熟女| 国产午夜精品论理片| 99九九线精品视频在线观看视频| 日日摸夜夜添夜夜添av毛片| 中文字幕熟女人妻在线| 国产片特级美女逼逼视频| 午夜福利在线观看免费完整高清在 |