• <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)

    為更準(zhǔn)確地模擬核素的非菲克遷移行為并尋求一種更合理的建模方法,根據(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ì)處置是目前公認(rèn)的較安全可靠的高放射性廢物處置方法,但在較大的時間尺度上將面臨處置庫退化導(dǎo)致放射性核素釋放進入地質(zhì)環(huán)境的潛在風(fēng)險;核素若隨地下水遷移進入生物圈,將對生態(tài)環(huán)境和人類健康產(chǎn)生潛在危害,因此,需要研究核素在地下水中的遷移行為這一重要問題.利用數(shù)值模擬手段對核素遷移行為進行定量研究是一種重要的方法.然而由于地質(zhì)介質(zhì)的非均質(zhì)性且受到水解反應(yīng),絡(luò)合反應(yīng),氧化還原反應(yīng)以及膠體作用的共同影響,核素在地下水中的遷移通常呈現(xiàn)非菲克性.基于對流-彌散方程(ADE)的溶質(zhì)運移模型無法準(zhǔn)確描述核素遷移過程中普遍存在的非菲克效應(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)時間到達(dá)距離處另一位置的遷移行為.示蹤劑在實驗室尺度土柱中遷移的經(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在地下水中的實際遷移速度和距離遠(yuǎn)高于計算值.因此本文選擇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)門。當(dāng)年,她家的牛肉面辣得登峰造極——我獨喜歡門口小爐子里燜煮的鹵干。

    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)境,沒有遭受破壞的當(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é)"等,不斷挖掘瓊中什寒旅游新熱點、新賣點,成功打造什寒鄉(xiāng)村休閑旅游品牌。

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

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

    式中: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)象,而穿透時間對的敏感度較小,不同的值下濃度到達(dá)峰值的時間相近.如圖1d, BTC對截斷時間2的敏感度低,當(dāng)2>103后,2呈指數(shù)級增長時BTC卻基本無變化.如圖1e,穿透曲線對表征參與吸附和解吸作用的參數(shù)平均吸附率較敏感,值較大時遷移初始速度較低,峰值濃度小,而達(dá)到穿透時間以后,隨著的增大,拖尾的漸近坡度幾乎不再變化.如圖1f,參數(shù)表征冪律形式下的拖尾效應(yīng),值越大,遷移速度更快,拖尾效應(yīng)更明顯.

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

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

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

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

    針對當(dāng)前目的地感官營銷實踐面臨的問題,結(jié)合上述目的地感官營銷的3種運用方式分析,本文認(rèn)為現(xià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進行評價:

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

    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次??蛇x高鈣液肥有黃腐酸生態(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ù),可以簡化建模過程.

    “王九,我說,你瞧,大爺大姑娘不來,先生可來了。好,咱們動手,先生不會走的。你小心別讓趙四小子扔倒。先生幫咱們繃個場面,看你摔趙四這小子,先生準(zhǔn)不走?!?/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] 夏 源,吳吉春.張 勇.改進時間分?jǐn)?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)

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

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

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機模型
    提煉模型 突破難點
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達(dá)及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
    3D打印中的模型分割與打包
    人妻久久中文字幕网| 久久热在线av| 亚洲伊人色综图| 国产淫语在线视频| 51午夜福利影视在线观看| 亚洲人成电影免费在线| 91精品国产国语对白视频| 一级a爱视频在线免费观看| 国产深夜福利视频在线观看| 成人亚洲精品一区在线观看| 超色免费av| 亚洲欧洲日产国产| 欧美日韩精品网址| 热99久久久久精品小说推荐| 久久久久久久精品吃奶| 午夜久久久在线观看| 高清欧美精品videossex| 国产深夜福利视频在线观看| 久久中文字幕一级| 日本wwww免费看| 免费在线观看日本一区| 三上悠亚av全集在线观看| 深夜精品福利| 欧美成狂野欧美在线观看| 大片电影免费在线观看免费| 高清欧美精品videossex| 午夜视频精品福利| 日本精品一区二区三区蜜桃| 欧美人与性动交α欧美精品济南到| 久久久久国产一级毛片高清牌| 久久精品国产综合久久久| 狠狠精品人妻久久久久久综合| 欧美黑人欧美精品刺激| 99香蕉大伊视频| 亚洲欧美日韩高清在线视频 | 亚洲五月婷婷丁香| 天堂8中文在线网| 一边摸一边做爽爽视频免费| 男女床上黄色一级片免费看| 欧美精品一区二区免费开放| 18禁黄网站禁片午夜丰满| 汤姆久久久久久久影院中文字幕| 老司机影院毛片| 国产免费现黄频在线看| 丝袜美腿诱惑在线| 国产成人精品无人区| 人成视频在线观看免费观看| 久久亚洲真实| 日本av免费视频播放| av片东京热男人的天堂| 国产精品久久电影中文字幕 | 欧美在线一区亚洲| 老司机午夜福利在线观看视频 | 嫩草影视91久久| 亚洲精品粉嫩美女一区| 久久久国产成人免费| 午夜免费鲁丝| 狠狠精品人妻久久久久久综合| 国产一区二区三区综合在线观看| 中文字幕人妻熟女乱码| 高潮久久久久久久久久久不卡| 黄色成人免费大全| 欧美午夜高清在线| 久久精品91无色码中文字幕| 免费av中文字幕在线| 手机成人av网站| 久久精品熟女亚洲av麻豆精品| 久久99一区二区三区| av又黄又爽大尺度在线免费看| 欧美精品啪啪一区二区三区| 国产真人三级小视频在线观看| cao死你这个sao货| 成人18禁在线播放| 另类精品久久| 91成年电影在线观看| 热99re8久久精品国产| 国精品久久久久久国模美| 欧美 日韩 精品 国产| 国产精品.久久久| 国产成人av教育| 精品少妇内射三级| 啦啦啦免费观看视频1| 国产福利在线免费观看视频| 91大片在线观看| 大陆偷拍与自拍| 亚洲中文av在线| 狠狠精品人妻久久久久久综合| 国产在线观看jvid| 精品乱码久久久久久99久播| 亚洲中文av在线| 高清黄色对白视频在线免费看| 男女边摸边吃奶| 免费在线观看黄色视频的| 欧美日韩一级在线毛片| 国产成人欧美在线观看 | 一边摸一边抽搐一进一小说 | 久久久久国产一级毛片高清牌| 亚洲欧美日韩高清在线视频 | 交换朋友夫妻互换小说| 人人妻人人澡人人看| 色尼玛亚洲综合影院| 大片免费播放器 马上看| 久久国产精品影院| 亚洲伊人色综图| 一二三四社区在线视频社区8| 欧美精品人与动牲交sv欧美| 亚洲人成伊人成综合网2020| 国产aⅴ精品一区二区三区波| 亚洲中文字幕日韩| 一个人免费在线观看的高清视频| 99热网站在线观看| 欧美激情久久久久久爽电影 | 少妇猛男粗大的猛烈进出视频| 黑人巨大精品欧美一区二区蜜桃| 亚洲九九香蕉| 国产高清videossex| 中文亚洲av片在线观看爽 | 亚洲中文av在线| 波多野结衣一区麻豆| 怎么达到女性高潮| 久久国产精品男人的天堂亚洲| 啪啪无遮挡十八禁网站| 黄色a级毛片大全视频| 欧美日韩亚洲高清精品| 亚洲欧美一区二区三区黑人| 制服人妻中文乱码| 老熟妇仑乱视频hdxx| 亚洲人成77777在线视频| 成人国语在线视频| 午夜两性在线视频| 国产精品自产拍在线观看55亚洲 | 五月天丁香电影| 午夜福利一区二区在线看| 高清毛片免费观看视频网站 | 免费观看av网站的网址| 亚洲国产av新网站| 高清欧美精品videossex| h视频一区二区三区| 成人手机av| 中文字幕av电影在线播放| 国产精品影院久久| 狠狠婷婷综合久久久久久88av| av网站在线播放免费| 亚洲专区国产一区二区| 国产男女超爽视频在线观看| 国产男女内射视频| 国产欧美日韩一区二区三区在线| 久久久久久久精品吃奶| 国产99久久九九免费精品| 日韩制服丝袜自拍偷拍| 国产精品 国内视频| 9色porny在线观看| av天堂在线播放| 在线看a的网站| 精品第一国产精品| 999久久久国产精品视频| 国产淫语在线视频| 精品视频人人做人人爽| 亚洲国产欧美网| 色播在线永久视频| av天堂在线播放| 久久久久久久久久久久大奶| 国产免费视频播放在线视频| 精品人妻在线不人妻| 制服诱惑二区| 欧美日韩国产mv在线观看视频| 精品一区二区三卡| 最黄视频免费看| 热99国产精品久久久久久7| 在线永久观看黄色视频| av福利片在线| 热99国产精品久久久久久7| 黄片小视频在线播放| 菩萨蛮人人尽说江南好唐韦庄| 国产精品久久久久久精品电影小说| 亚洲第一青青草原| 久久久久视频综合| 欧美成人午夜精品| 欧美激情高清一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 国产又爽黄色视频| 在线观看免费午夜福利视频| 欧美日韩av久久| 日韩三级视频一区二区三区| 法律面前人人平等表现在哪些方面| 不卡一级毛片| 纵有疾风起免费观看全集完整版| 少妇猛男粗大的猛烈进出视频| 男女无遮挡免费网站观看| 国产色视频综合| 黑人操中国人逼视频| 国产精品香港三级国产av潘金莲| 视频区图区小说| 欧美日韩av久久| 岛国毛片在线播放| 欧美另类亚洲清纯唯美| 嫩草影视91久久| 纵有疾风起免费观看全集完整版| 国精品久久久久久国模美| 在线观看一区二区三区激情| 国产av又大| av电影中文网址| 亚洲国产看品久久| 女警被强在线播放| 成年版毛片免费区| 大陆偷拍与自拍| 久久国产精品人妻蜜桃| 久久国产精品男人的天堂亚洲| 亚洲熟女精品中文字幕| 国产精品98久久久久久宅男小说| 狠狠婷婷综合久久久久久88av| 黑人欧美特级aaaaaa片| 天天影视国产精品| 午夜精品久久久久久毛片777| 中文字幕精品免费在线观看视频| 国产在线视频一区二区| 制服人妻中文乱码| 久久亚洲精品不卡| 日韩中文字幕视频在线看片| 国产黄频视频在线观看| 女人爽到高潮嗷嗷叫在线视频| av视频免费观看在线观看| 18禁黄网站禁片午夜丰满| 一二三四在线观看免费中文在| 成人免费观看视频高清| 精品国产乱子伦一区二区三区| 国产一区有黄有色的免费视频| 久久精品亚洲av国产电影网| 9热在线视频观看99| 国产免费现黄频在线看| 国产精品一区二区在线不卡| 亚洲欧美色中文字幕在线| 香蕉国产在线看| 亚洲va日本ⅴa欧美va伊人久久| 久久久水蜜桃国产精品网| 自拍欧美九色日韩亚洲蝌蚪91| 999精品在线视频| 精品福利永久在线观看| 超碰成人久久| 久久中文字幕人妻熟女| 少妇猛男粗大的猛烈进出视频| 欧美变态另类bdsm刘玥| 老司机福利观看| 亚洲国产成人一精品久久久| 亚洲av国产av综合av卡| 欧美性长视频在线观看| 国产精品.久久久| √禁漫天堂资源中文www| 亚洲 国产 在线| 老鸭窝网址在线观看| 黄色毛片三级朝国网站| 国产91精品成人一区二区三区 | 久久久久久久大尺度免费视频| cao死你这个sao货| 不卡av一区二区三区| 99在线人妻在线中文字幕 | 国产精品1区2区在线观看. | 国产色视频综合| 免费日韩欧美在线观看| 亚洲性夜色夜夜综合| 国产一区二区 视频在线| 日韩 欧美 亚洲 中文字幕| 国产精品美女特级片免费视频播放器 | 99国产综合亚洲精品| 最近最新中文字幕大全电影3 | 免费在线观看完整版高清| 欧美成人免费av一区二区三区 | 国产又色又爽无遮挡免费看| 一本色道久久久久久精品综合| 999精品在线视频| 日韩人妻精品一区2区三区| 亚洲av片天天在线观看| 国产精品九九99| 久久久久精品国产欧美久久久| 亚洲成人免费av在线播放| 这个男人来自地球电影免费观看| 一边摸一边抽搐一进一出视频| 老汉色av国产亚洲站长工具| 黑人猛操日本美女一级片| 人人妻,人人澡人人爽秒播| 亚洲人成电影观看| 精品福利观看| 丝袜美腿诱惑在线| 久久青草综合色| 久久国产精品男人的天堂亚洲| 国产精品秋霞免费鲁丝片| 国产日韩欧美视频二区| 这个男人来自地球电影免费观看| 99精品在免费线老司机午夜| 黄色视频在线播放观看不卡| 91麻豆精品激情在线观看国产 | 久久久国产欧美日韩av| 午夜久久久在线观看| 午夜成年电影在线免费观看| 丰满人妻熟妇乱又伦精品不卡| 日日夜夜操网爽| 757午夜福利合集在线观看| 亚洲精品在线美女| 免费黄频网站在线观看国产| 另类亚洲欧美激情| 天天添夜夜摸| 日韩大码丰满熟妇| 一进一出抽搐动态| 法律面前人人平等表现在哪些方面| 欧美久久黑人一区二区| 99国产精品一区二区三区| 久久久久久久大尺度免费视频| 精品视频人人做人人爽| 男人操女人黄网站| 搡老熟女国产l中国老女人| 国产亚洲精品第一综合不卡| 亚洲人成77777在线视频| 一区在线观看完整版| 久久久久视频综合| 亚洲成人手机| 亚洲五月色婷婷综合| 精品视频人人做人人爽| 国产男女超爽视频在线观看| 又大又爽又粗| 成年女人毛片免费观看观看9 | 国产av一区二区精品久久| 欧美日韩亚洲综合一区二区三区_| 国产精品成人在线| 大片电影免费在线观看免费| 精品少妇久久久久久888优播| 久久精品人人爽人人爽视色| a级片在线免费高清观看视频| 国产午夜精品久久久久久| 欧美日韩视频精品一区| 制服人妻中文乱码| 中文亚洲av片在线观看爽 | 动漫黄色视频在线观看| 国产精品国产av在线观看| a在线观看视频网站| 熟女少妇亚洲综合色aaa.| 黑人巨大精品欧美一区二区mp4| 成年动漫av网址| 电影成人av| 欧美一级毛片孕妇| 在线 av 中文字幕| 亚洲精品国产区一区二| 亚洲国产欧美网| 久久久国产欧美日韩av| 男女边摸边吃奶| 日韩一卡2卡3卡4卡2021年| 99久久99久久久精品蜜桃| videos熟女内射| 亚洲性夜色夜夜综合| 一区二区av电影网| 99精品在免费线老司机午夜| 精品国产超薄肉色丝袜足j| 久久精品熟女亚洲av麻豆精品| 视频区欧美日本亚洲| 在线观看免费日韩欧美大片| 国产精品一区二区免费欧美| 精品一区二区三区视频在线观看免费 | 正在播放国产对白刺激| 中文字幕精品免费在线观看视频| videos熟女内射| 欧美黄色片欧美黄色片| 欧美黄色片欧美黄色片| 性色av乱码一区二区三区2| 久久久国产一区二区| 亚洲午夜精品一区,二区,三区| 丁香六月天网| 亚洲人成77777在线视频| 亚洲成a人片在线一区二区| 丰满少妇做爰视频| 日本欧美视频一区| 人人妻人人澡人人爽人人夜夜| 午夜老司机福利片| 王馨瑶露胸无遮挡在线观看| 十八禁高潮呻吟视频| 大陆偷拍与自拍| 亚洲伊人久久精品综合| 黑人操中国人逼视频| 国产深夜福利视频在线观看| videosex国产| 亚洲伊人久久精品综合| 777久久人妻少妇嫩草av网站| 日韩欧美三级三区| 精品国产乱码久久久久久小说| 不卡av一区二区三区| 手机成人av网站| 国产精品久久久人人做人人爽| 成人黄色视频免费在线看| 丝袜美足系列| 国产一区有黄有色的免费视频| 亚洲国产精品一区二区三区在线| 黑人猛操日本美女一级片| 一边摸一边做爽爽视频免费| 欧美午夜高清在线| 菩萨蛮人人尽说江南好唐韦庄| 亚洲欧美色中文字幕在线| 露出奶头的视频| 亚洲综合色网址| 一级黄色大片毛片| 欧美日本中文国产一区发布| 精品午夜福利视频在线观看一区 | 悠悠久久av| 日韩熟女老妇一区二区性免费视频| 国产淫语在线视频| 五月开心婷婷网| 一本—道久久a久久精品蜜桃钙片| 精品少妇一区二区三区视频日本电影| 久久精品人人爽人人爽视色| 一级片'在线观看视频| 不卡一级毛片| 免费观看av网站的网址| 久久青草综合色| 免费少妇av软件| 亚洲人成电影观看| 午夜91福利影院| 欧美午夜高清在线| 国产欧美日韩综合在线一区二区| 99九九在线精品视频| 美女高潮喷水抽搐中文字幕| 多毛熟女@视频| 丁香六月天网| 欧美成狂野欧美在线观看| 99国产精品一区二区三区| 岛国毛片在线播放| 人人妻人人澡人人看| 变态另类成人亚洲欧美熟女 | 69av精品久久久久久 | av有码第一页| 欧美另类亚洲清纯唯美| 悠悠久久av| 色精品久久人妻99蜜桃| 午夜福利影视在线免费观看| 精品乱码久久久久久99久播| 91麻豆av在线| 欧美+亚洲+日韩+国产| 国产精品久久久久久人妻精品电影 | 99国产精品一区二区三区| 亚洲精品国产精品久久久不卡| www.熟女人妻精品国产| 国产97色在线日韩免费| 亚洲精品在线美女| 精品福利永久在线观看| 午夜91福利影院| 日本黄色视频三级网站网址 | av超薄肉色丝袜交足视频| 老熟妇乱子伦视频在线观看| 中文欧美无线码| 欧美日韩亚洲高清精品| 最黄视频免费看| 天堂俺去俺来也www色官网| 在线看a的网站| 久久精品国产综合久久久| 国产不卡av网站在线观看| 18禁裸乳无遮挡动漫免费视频| 久久久久久久国产电影| 桃红色精品国产亚洲av| 一边摸一边抽搐一进一出视频| 最近最新免费中文字幕在线| 飞空精品影院首页| 亚洲成国产人片在线观看| 日韩欧美一区视频在线观看| 一级毛片女人18水好多| 亚洲性夜色夜夜综合| 99国产精品一区二区三区| 视频在线观看一区二区三区| 18禁观看日本| 日韩中文字幕视频在线看片| 热99久久久久精品小说推荐| 亚洲 欧美一区二区三区| 国产男女超爽视频在线观看| 精品福利观看| 美女主播在线视频| 91成人精品电影| 欧美国产精品一级二级三级| 岛国毛片在线播放| 天天躁日日躁夜夜躁夜夜| 老鸭窝网址在线观看| 国产精品久久久久久精品古装| 这个男人来自地球电影免费观看| 777久久人妻少妇嫩草av网站| 高清毛片免费观看视频网站 | 久久性视频一级片| av天堂在线播放| 午夜成年电影在线免费观看| 国产欧美日韩精品亚洲av| 午夜91福利影院| 亚洲国产欧美一区二区综合| 日韩免费av在线播放| 精品欧美一区二区三区在线| 日本精品一区二区三区蜜桃| 亚洲成人手机| 最近最新免费中文字幕在线| 天天躁夜夜躁狠狠躁躁| 亚洲精品国产一区二区精华液| 一本一本久久a久久精品综合妖精| 王馨瑶露胸无遮挡在线观看| 中亚洲国语对白在线视频| 999精品在线视频| 成人国语在线视频| kizo精华| 亚洲av成人一区二区三| 午夜福利在线观看吧| 午夜视频精品福利| av网站在线播放免费| 欧美黄色片欧美黄色片| 亚洲专区字幕在线| 99riav亚洲国产免费| 亚洲av成人不卡在线观看播放网| 亚洲精品国产一区二区精华液| 精品少妇黑人巨大在线播放| 日韩免费高清中文字幕av| 啪啪无遮挡十八禁网站| 12—13女人毛片做爰片一| 精品久久蜜臀av无| 久久av网站| av电影中文网址| 一级,二级,三级黄色视频| 大型黄色视频在线免费观看| 免费久久久久久久精品成人欧美视频| 国产一区二区三区综合在线观看| 超碰成人久久| 如日韩欧美国产精品一区二区三区| 色在线成人网| 成人18禁高潮啪啪吃奶动态图| 老司机影院毛片| 黑人巨大精品欧美一区二区mp4| 精品国产一区二区久久| 亚洲精品久久成人aⅴ小说| 中文字幕av电影在线播放| 日本a在线网址| 动漫黄色视频在线观看| 亚洲伊人久久精品综合| 久久天堂一区二区三区四区| 亚洲国产欧美网| 这个男人来自地球电影免费观看| 亚洲少妇的诱惑av| 国产精品一区二区在线不卡| 丝袜美腿诱惑在线| 亚洲欧美激情在线| 首页视频小说图片口味搜索| 欧美激情 高清一区二区三区| 国产亚洲精品一区二区www | 黄色怎么调成土黄色| 国内毛片毛片毛片毛片毛片| 午夜两性在线视频| 国产成人欧美在线观看 | 国产精品偷伦视频观看了| 亚洲国产欧美日韩在线播放| 国产熟女午夜一区二区三区| 久久久久久免费高清国产稀缺| 欧美精品一区二区免费开放| 在线亚洲精品国产二区图片欧美| 中文字幕人妻熟女乱码| 伊人久久大香线蕉亚洲五| 亚洲伊人久久精品综合| 无遮挡黄片免费观看| 最新在线观看一区二区三区| 国产淫语在线视频| 久久影院123| 高清视频免费观看一区二区| 91老司机精品| 男女床上黄色一级片免费看| 窝窝影院91人妻| 亚洲精品国产精品久久久不卡| 成人精品一区二区免费| 国产精品一区二区在线观看99| 午夜日韩欧美国产| 国产伦人伦偷精品视频| 9热在线视频观看99| av网站免费在线观看视频| 国产精品久久久久久精品电影小说| 精品免费久久久久久久清纯 | 精品少妇内射三级| 大陆偷拍与自拍| 亚洲va日本ⅴa欧美va伊人久久| 人人妻人人澡人人爽人人夜夜| 女人爽到高潮嗷嗷叫在线视频| 黄色视频,在线免费观看| 又黄又粗又硬又大视频| 国产av精品麻豆| 免费在线观看视频国产中文字幕亚洲| 久久99一区二区三区| 岛国在线观看网站| 精品福利永久在线观看| 国产精品99久久99久久久不卡| 啪啪无遮挡十八禁网站| 在线观看66精品国产| 制服诱惑二区| 999精品在线视频| 午夜精品久久久久久毛片777| 久久久精品区二区三区| 成人黄色视频免费在线看| 国产亚洲一区二区精品| 亚洲熟女精品中文字幕| 久久人人97超碰香蕉20202| 国产成人欧美在线观看 | 国产淫语在线视频| 亚洲精品国产色婷婷电影| 国产不卡av网站在线观看| av天堂久久9| 国产在线视频一区二区| 窝窝影院91人妻| 人人妻人人澡人人爽人人夜夜| 9191精品国产免费久久| 成在线人永久免费视频| 咕卡用的链子| 免费观看人在逋| 色婷婷久久久亚洲欧美| 电影成人av| 在线观看www视频免费| 国产色视频综合| 女人被躁到高潮嗷嗷叫费观| 视频区图区小说| 涩涩av久久男人的天堂| 自拍欧美九色日韩亚洲蝌蚪91| 少妇被粗大的猛进出69影院| 国产精品一区二区免费欧美|