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

    基于137Cs本底值計(jì)算模型和地理加權(quán)回歸克里格對(duì)137Cs本底值空間分布的預(yù)測(cè)

    2018-11-15 08:03:06文安邦
    中國(guó)水土保持科學(xué) 2018年5期
    關(guān)鍵詞:區(qū)域模型研究

    李 豪, 文安邦, 劉 濤, 李 婷

    (1.四川農(nóng)業(yè)大學(xué)資源學(xué)院, 611130,成都;2.中國(guó)科學(xué)院水利部成都山地災(zāi)害與環(huán)境研究所,610041,成都)

    土壤侵蝕造成土地資源退化,生產(chǎn)力下降,是危害人類生存與發(fā)展的世界性環(huán)境問(wèn)題之一。深入了解土壤侵蝕規(guī)律、獲取可靠的土壤侵蝕信息是開(kāi)展水土保持工作的基礎(chǔ)。核示蹤技術(shù)是土壤侵蝕研究的一種重要方法,具有成本低、高效可靠等優(yōu)點(diǎn),其中又以137Cs示蹤技術(shù)應(yīng)用最為廣泛和成功[1]。自20世紀(jì)60年代以來(lái),在國(guó)際原子能機(jī)構(gòu)(IAEA,International Atomic Energy Agency)等國(guó)際組織的大力推動(dòng)下,其技術(shù)體系逐漸完善,測(cè)定的土壤侵蝕數(shù)據(jù)的準(zhǔn)確性得到廣泛驗(yàn)證和認(rèn)可。

    在一定區(qū)域內(nèi),未受擾動(dòng)、無(wú)侵蝕亦無(wú)沉積發(fā)生的土壤條件下137Cs的面積濃度代表該區(qū)域內(nèi)每m2上137Cs大氣沉降的累積量,稱為該區(qū)域的137Cs本底值(137Cs reference inventory,Bq/m2)。137Cs示蹤法研究土壤侵蝕的原理是分析研究區(qū)某一土壤剖面的137Cs含量與137Cs本底值之間的差異,選擇合適的計(jì)算模型,定量計(jì)算得到該處的土壤侵蝕量或沉積量,因此準(zhǔn)確的137Cs本底值是運(yùn)用137Cs示蹤技術(shù)研究土壤侵蝕的關(guān)鍵因素。目前,獲取137Cs本底值數(shù)據(jù)的主要手段是野外選擇符合沒(méi)有人為擾動(dòng)、無(wú)侵蝕和沉積發(fā)生等條件的樣地采集土壤樣品,然后在實(shí)驗(yàn)室測(cè)試得到;但由于近幾十年來(lái)人類活動(dòng)強(qiáng)度大,在許多區(qū)域難以找到完全符合標(biāo)準(zhǔn)的采樣地,該方法的實(shí)現(xiàn)存在一定難度。有學(xué)者指出,由于該方法在采樣過(guò)程中普遍存在的一系列等問(wèn)題,導(dǎo)致獲取的137Cs本底值數(shù)據(jù)的可靠性存在疑問(wèn)[2-3]。

    以實(shí)測(cè)數(shù)據(jù)為基礎(chǔ)、基于GIS的空間插值技術(shù)已廣泛用于獲取土壤、環(huán)境和氣象等要素的高分辨率空間分布數(shù)據(jù)。137Cs具有與上述要素類似的空間分布特性,空間插值技術(shù)亦適用于137Cs本底值研究;但由于現(xiàn)階段我國(guó)137Cs本底值的實(shí)測(cè)樣點(diǎn)數(shù)量偏少(近30年來(lái),全國(guó)137Cs本底值研究樣點(diǎn)僅有100余個(gè)[4]),且空間分布不均勻,并不適合通過(guò)137Cs本底值實(shí)測(cè)數(shù)據(jù)直接進(jìn)行空間插值。采用GIS的空間插值技術(shù)獲取區(qū)域137Cs本底值資料,須另辟途徑。

    筆者以長(zhǎng)江上游的四川省為研究對(duì)象,嘗試采用137Cs本底值理論計(jì)算模型、主成分分析與地理加權(quán)回歸克里格插值相結(jié)合的技術(shù)方法,對(duì)137Cs本底值空間分布進(jìn)行研究,探討通過(guò)氣象站點(diǎn)相對(duì)豐富的實(shí)測(cè)數(shù)據(jù)、基于本底值理論計(jì)算模型和GIS空間插值技術(shù)獲取區(qū)域137Cs本底值數(shù)據(jù)的合適方法,以期獲取更高精度的區(qū)域137Cs本底值資料,為在該區(qū)開(kāi)展土壤侵蝕的核示蹤研究、進(jìn)行水土流失治理等工作提供有效的數(shù)據(jù)支持。

    1 研究區(qū)概況

    四川省位于我國(guó)西南,地處長(zhǎng)江上游,E 92°21′~108°12′,N 26°03′~34°19′,區(qū)域面積48.5萬(wàn)km2。本區(qū)處于第1級(jí)青藏高原和第2級(jí)長(zhǎng)江中下游平原的過(guò)渡帶,區(qū)內(nèi)海拔高差懸殊,地貌復(fù)雜多樣,區(qū)域以山地地貌為主,全省可分為東部四川盆地、川西北高原和川西南山地等3大部分。區(qū)內(nèi)季風(fēng)氣候明顯,氣候類型多樣,差異顯著:東部為亞熱帶濕潤(rùn)、半濕潤(rùn)氣候,雨熱同季,水熱條件好,年降水量1 000~1 200 mm,但多云霧、日照時(shí)間較短;西部為高山高原高寒氣候,氣候垂直變化明顯,總體上以寒溫帶氣候?yàn)橹?,冬寒夏涼,水熱不足,年降水?00~900 mm。四川省土壤類型豐富,區(qū)域分異明顯,據(jù)第2次土壤普查,全省土壤類型共有25個(gè)土類、66個(gè)亞類,以水稻土、紫色土、黃壤等為主,本區(qū)是我國(guó)紫色土分布最集中的地區(qū)。

    四川省是我國(guó)土壤侵蝕較嚴(yán)重的省區(qū)之一,土壤侵蝕分布廣,面積大,危害嚴(yán)重。本區(qū)土壤侵蝕面積和數(shù)量約占整個(gè)長(zhǎng)江上游總量的70%,長(zhǎng)江上游的泥沙主要來(lái)源于地處金沙江和嘉陵江中下游的四川盆地丘陵地區(qū)。根據(jù)第1次全國(guó)水利普查公報(bào)數(shù)據(jù),全省水蝕面積達(dá)11.4萬(wàn)km2,占幅員面積的23.54%,其中中度以上土壤侵蝕面積占比達(dá)到了57.63%[5]。

    圖1 研究區(qū)域及氣象站點(diǎn)分布圖Fig.1 Spatial distribution of meteorological stations in the studied area

    2 數(shù)據(jù)來(lái)源與研究方法

    2.1 數(shù)據(jù)來(lái)源

    筆者采用的四川省147個(gè)氣象站點(diǎn)的經(jīng)緯度、年均降水量(觀測(cè)精度為0.1 mm)等數(shù)據(jù)來(lái)自國(guó)家氣象信息中心提供的“中國(guó)基本、基準(zhǔn)和一般地面氣象觀測(cè)站1981—2010年累年值年值數(shù)據(jù)集”。各氣象站點(diǎn)的空間分布如圖1所示。雖然該數(shù)據(jù)集與本文的研究對(duì)象在時(shí)間尺度上并不完全相同;但是,該數(shù)據(jù)集具有較長(zhǎng)的時(shí)間序列(30年),均為地面實(shí)測(cè)數(shù)據(jù),可信度高,且站點(diǎn)數(shù)量多,能保證插值的精度,滿足研究要求。在參照類似研究的技術(shù)方法[4]后,我們采用該數(shù)據(jù)集作為本研究的降水?dāng)?shù)據(jù)。

    2.2 數(shù)據(jù)處理

    D.E.Walling等[6]建立的137Cs本底值全球空間分布模型(以下簡(jiǎn)稱W&H模型)是基于核爆產(chǎn)物大氣沉降數(shù)據(jù)、建立在認(rèn)識(shí)137Cs沉降過(guò)程的基礎(chǔ)上的理論模型,物理意義明確,目前應(yīng)用最為廣泛。本文各氣象站點(diǎn)的137Cs本底值模型計(jì)算值(以下簡(jiǎn)稱137Cs本底計(jì)算值)通過(guò)該模型計(jì)算得到,本文涉及的137Cs本底值數(shù)據(jù)均衰變校正至2016年1月1日。

    W&H模型[6]按照10°緯度帶劃分計(jì)算帶,對(duì)各帶分別計(jì)算137Cs本底值數(shù)據(jù),而四川省地處N 26°~34°之間,跨越了N 20°~30°和30°~40° 2個(gè)計(jì)算帶。為提高預(yù)測(cè)結(jié)果的精度,筆者對(duì)四川省進(jìn)行了細(xì)分,以N30°為界劃分為2個(gè)研究區(qū)域(圖1)。其中,低于N30°區(qū)域(以下簡(jiǎn)稱L30區(qū)域)包含63個(gè)站點(diǎn),高于N 30°區(qū)域(以下簡(jiǎn)稱H30區(qū)域)包含84個(gè)站點(diǎn),對(duì)這2個(gè)研究區(qū)域分別預(yù)測(cè)其137Cs本底計(jì)算值的空間分布。

    降水、地理位置及進(jìn)行大氣核試驗(yàn)的場(chǎng)所等是影響我國(guó)137Cs本底值分布的主要因素[7]。而我國(guó)的大氣核試驗(yàn)均在新疆羅布泊地區(qū)進(jìn)行,對(duì)本研究區(qū)域的137Cs本底值分布基本沒(méi)有影響;因此,筆者選擇年降水量(代表降水)、經(jīng)度和緯度(代表地理位置)等3個(gè)影響因素作為輔助變量,用于回歸分析和空間插值。筆者采用的技術(shù)方法主要包括主成分分析和地理加權(quán)回歸克里格插值。

    2.3 主成分分析

    主成分分析(principal component analysis, PCA)是一種應(yīng)用廣泛的統(tǒng)計(jì)分析方法,主要用于降維和消除變量間的相關(guān)性。其通過(guò)正交變換將1組存在相關(guān)性的變量轉(zhuǎn)換為1組數(shù)量較少的、線性不相關(guān)的變量(稱之為主成分),在盡可能多地保留原來(lái)變量的信息的基礎(chǔ)上,可減少變量的數(shù)量,有效降低分析問(wèn)題的復(fù)雜性[8]。

    2.4 地理加權(quán)回歸克里格插值

    克里格插值法是一種基于地統(tǒng)計(jì)學(xué)的經(jīng)典空間插值方法,目前已在多個(gè)領(lǐng)域得到廣泛應(yīng)用。有學(xué)者提出,通過(guò)引入與目標(biāo)變量具有較好相關(guān)性的影響因素變量、構(gòu)建回歸克里格模型(regression Kriging, RK),可以有效提高克里格法的插值精度[9-11]。RK法的基本思想是建立目標(biāo)變量和影響因素之間的回歸方程,分離出趨勢(shì)項(xiàng),然后對(duì)回歸殘差(即隨機(jī)因素)采用傳統(tǒng)的普通克里格插值法(ordinary Kriging, OK)進(jìn)行插值,最后對(duì)回歸預(yù)測(cè)的趨勢(shì)項(xiàng)和回歸殘差的估計(jì)值進(jìn)行求和,從而得到目標(biāo)變量的預(yù)測(cè)值。RK法可用以下式[12]表示:

    (1)

    式中:y(pi)為位置pi處的目標(biāo)變量;βj,(pi)為回歸系數(shù);xj,(pi)為位置pi處的第j個(gè)自變量;ε(pi)為位置pi處的回歸殘差;n為自變量的個(gè)數(shù)。

    傳統(tǒng)的RK法中,回歸分析部分大都采用基于最小二乘法(ordinary least squares, OLS)的全局回歸(global regression, GR),各回歸系數(shù)為常數(shù),即為全局回歸克里格法(global regression Kriging, GRK)。由于沒(méi)有考慮變量的空間非平穩(wěn)性,因此在處理土壤、降水等具有較強(qiáng)空間非平穩(wěn)性的變量時(shí),該方法的插值精度仍受到一定的限制。

    地理加權(quán)回歸克里格法(geographically weighted regression Kriging, GWRK)法對(duì)GRK法進(jìn)行了改進(jìn),是將GRK法中基于OLS的全局回歸替換為地理加權(quán)回歸的局部加權(quán)回歸(local regression, LR),其余步驟、方法保持不變的一種空間插值方法。GWRK法的核心是地理加權(quán)回歸(geographically weighted regression, GWR)。GWR法對(duì)全局回歸模型進(jìn)行了改進(jìn),通過(guò)構(gòu)建空間權(quán)重函數(shù),對(duì)空間非平穩(wěn)性進(jìn)行了量化。在GWR中,回歸系數(shù)不再是常數(shù),而是空間位置的函數(shù),位置(ui,vi)處的回歸系數(shù)β通過(guò)下式[12]計(jì)算得到:

    (ui,vi)=[XTW(ui,vi)X]-1XTW(ui,vi)Y。

    (2)

    式中:β為回歸系數(shù)矩陣;X為自變量設(shè)計(jì)矩陣(XT為其轉(zhuǎn)置矩陣);W(ui,vi)為空間權(quán)重矩陣,由空間權(quán)重函數(shù)W(i)求得,其作用是定量衡量領(lǐng)域內(nèi)不同空間位置k(k=1,2,…,n)樣點(diǎn)的觀測(cè)值對(duì)于回歸點(diǎn)(ui,vi)回歸系數(shù)估計(jì)的影響程度[12](n為領(lǐng)域樣點(diǎn)數(shù));Y為因變量矩陣。

    已有研究表明,GWRK的插值結(jié)果能揭示被空間非平穩(wěn)性所掩蓋的一些局部變化,更真實(shí)地反映目標(biāo)變量的空間變異情況,使插值精度得到進(jìn)一步的提高[13]?;谏鲜龇治?,筆者采用GWRK法進(jìn)行四川省137Cs本底計(jì)算值的插值研究。

    2.5 模型精度的評(píng)價(jià)

    從147個(gè)站點(diǎn)數(shù)據(jù)中隨機(jī)選取30個(gè)站點(diǎn)構(gòu)建驗(yàn)證數(shù)據(jù)集。采用其余117個(gè)站點(diǎn)數(shù)據(jù)建立137Cs本底值插值模型,并求取30個(gè)驗(yàn)證站點(diǎn)的137Cs本底值預(yù)測(cè)值,通過(guò)比較驗(yàn)證站點(diǎn)的理論計(jì)算值與插值預(yù)測(cè)值之間的平均誤差(mean error, ME)、平均絕對(duì)誤差(mean absolute error,MAE)、均方根誤差(root mean square error, RMSE)[13]、平均絕對(duì)百分誤差(mean absolute relative error, MARE)和均方根相對(duì)誤差(root mean square relative error, RMSRE)等指標(biāo)評(píng)估不同插值模型的插值精度。其中,MARE和RMSRE定義如下:

    (3)

    (4)

    式中:ya,i為位置i處的137Cs本底計(jì)算值,Bq/m2;ye,i為位置i處的137Cs本底預(yù)測(cè)值,Bq/m2;n為驗(yàn)證數(shù)據(jù)集的站點(diǎn)個(gè)數(shù)。

    2.6 研究軟件平臺(tái)

    常規(guī)統(tǒng)計(jì)分析和PCA采用SPSS 23完成,基于OLS的GR和GWR回歸分析由GWR4完成,半方差函數(shù)計(jì)算通過(guò)GS+9進(jìn)行,克里格插值和各空間分布圖的繪制由ArcGIS 10.4完成。各變量空間分布圖的分辨率均為1 km×1 km。

    3 結(jié)果與分析

    3.1 137Cs本底值描述性統(tǒng)計(jì)分析

    對(duì)研究區(qū)域各站點(diǎn)的137Cs本底計(jì)算值指標(biāo)進(jìn)行了描述性統(tǒng)計(jì)分析(表1)。分析結(jié)果表明,在2個(gè)研究區(qū)域各站點(diǎn)的137Cs本底計(jì)算值在數(shù)量上存在一定的差異,H30區(qū)域和L30區(qū)域的平均值分別為1 167.42和616.39 Bq/m2,而137Cs本底計(jì)算值的變異系數(shù)分別為8.65%和7.48%,在分異程度上基本一致。2個(gè)區(qū)域的137Cs本底計(jì)算值數(shù)據(jù)均通過(guò)了J-B(Jarque-Bera)檢驗(yàn)(P>0.05),符合正態(tài)分布。

    表1 研究區(qū)域站點(diǎn)的137Cs本底計(jì)算值描述性統(tǒng)計(jì)分析和正態(tài)性檢驗(yàn)

    3.2 影響因素主成分分析

    表2為計(jì)算的研究區(qū)域137Cs本底計(jì)算值與各影響因素間的相關(guān)系數(shù),結(jié)果顯示:137Cs本底計(jì)算值與各影響因素均存在較高的相關(guān),與經(jīng)度、降水量呈極顯著正相關(guān),與緯度呈極顯著負(fù)相關(guān)。另一方面,各影響因素間也存在著不同程度的相互作用,如降水量與經(jīng)度間達(dá)到了極顯著正相關(guān)。因此,為了消除變量間的信息冗余性和多重共線性,更好地反映各變量間的真實(shí)關(guān)系,首先對(duì)降水量、經(jīng)度和緯度等3個(gè)影響因素進(jìn)行了主成分分析。

    對(duì)各影響因素變量進(jìn)行主成分分析后,根據(jù)方差累計(jì)貢獻(xiàn)率≥85%的原則提取了主成分(表3),結(jié)果表明:對(duì)于2個(gè)研究區(qū)域,前2項(xiàng)主成分(F1、F2)的方差累計(jì)貢獻(xiàn)率均超過(guò)90%,分別達(dá)到了95.39%和94.66%,表明前2項(xiàng)主成分已經(jīng)表達(dá)了原始數(shù)據(jù)中絕大部分信息,采用前2項(xiàng)主成分即可較好地替代原數(shù)據(jù)信息;因此,后續(xù)的分析中將降水量、經(jīng)度和緯度等3個(gè)原始影響因素指標(biāo)轉(zhuǎn)換為2個(gè)主成分變量F1和F2,2個(gè)主成分的空間分布如圖2所示。

    表2 137Cs本底值與各影響因素的相關(guān)系數(shù)

    注:同一變量行上方數(shù)據(jù)表示H30區(qū)域的結(jié)果,下方數(shù)據(jù)表示L30區(qū)域的結(jié)果;*表示在0.05水平上具有顯著性,** 表示在0.01水平上具有顯著性,下同. Notes: The data in the top of the same row refer to the results of H30, the data in the bottom of the same row refer to the results of L30. * refers to significance at 0.05 level, ** refers to the significance at 0.01 level. The same below.

    表3 主成分的特征值與方差貢獻(xiàn)率

    Notes:F1refers to for the first principal component,F2refers to for the second principal component, andF3refers to for the third principal component. The same below.

    (a) H30區(qū)域F1變量 F1 in H30;(b) H30區(qū)域F2變量 F2 in H30;(c) L30區(qū)域F1變量 F1 in L30;(d) L30區(qū)域F2變量 F2 in L30圖2 研究區(qū)域主成分變量空間分布圖Fig.2 Spatial distribution of principal components in the studied area

    3.3 回歸分析

    以各站點(diǎn)的2個(gè)主成分值為自變量,137Cs本底計(jì)算值為因變量,分別采用基于OLS的GR和GWR2種模型進(jìn)行了回歸分析。2種回歸分析的結(jié)果分別如表4和表5所示。

    表4 基于OLS的GR回歸計(jì)算結(jié)果

    Notes: GR stands for global regression. OLS stands for ordinary least square. VIF stands for variance inflation factor. The same below.

    表5 GWR回歸計(jì)算結(jié)果

    Notes:GWR stands for geographically weighted regression. The same below.

    方差膨脹因子(variance inflation factor,VIF)是衡量自變量間是否存在多重共線性的主要指標(biāo)。由表4的計(jì)算結(jié)果來(lái)看,對(duì)原始變量進(jìn)行主成分分析后,回歸分析各自變量的VIF值遠(yuǎn)小于10,這說(shuō)明經(jīng)過(guò)PCA處理后,各自變量間基本不存在多重共線性問(wèn)題,回歸分析的結(jié)果是可靠的。從GWR模型的計(jì)算結(jié)果來(lái)看,各解釋變量的回歸系數(shù)值不再是常數(shù),如H30區(qū)域的F1回歸系數(shù)的取值范圍介于-57.19與-34.18間,變異系數(shù)達(dá)到了33.46%,這說(shuō)明137Cs本底計(jì)算值與各解釋變量間的相互關(guān)系具有顯著的空間非平穩(wěn)性。

    對(duì)2種模型的回歸效果進(jìn)行了比較。表6的分析結(jié)果表明:GWR模型的各項(xiàng)診斷指標(biāo)均明顯優(yōu)于基于OLS的GR模型。如在H30區(qū)域,GWR模型的RSS與GR模型相比,大幅減少了57.33%;AICc值也從GR模型的674.61顯著下降到GWR模型的639.37。這說(shuō)明基于局部回歸的GWR模型能夠更好的探測(cè)和解釋變量間的空間非平穩(wěn)性,相對(duì)于GR模型,具有更好的回歸效果。

    表6 GR和GWR模型的診斷指標(biāo)

    3.4 回歸殘差的克里格插值

    正態(tài)性檢驗(yàn)的結(jié)果顯示,GR模型在2個(gè)研究區(qū)域、GWR模型在L30區(qū)域的殘差均通過(guò)了J-B檢驗(yàn),為正態(tài)分布,而GWR模型在H30區(qū)域的回歸殘差不符合正態(tài)分布;因此,首先對(duì)該殘差數(shù)據(jù)進(jìn)行5/7次方轉(zhuǎn)換,使其滿足進(jìn)行克里格插值的條件。

    在GS+9軟件中,對(duì)2種模型的回歸殘差(包括H30區(qū)域轉(zhuǎn)換后的回歸殘差)進(jìn)行半方差函數(shù)擬合。從表7的擬合結(jié)果可以看出,各半方差函數(shù)的R2均達(dá)到0.80,取得了較好的擬合效果。此外,各殘差的半方差函數(shù)的塊基比(即C0/(C0+C))大都小于25%,說(shuō)明各殘差存在較強(qiáng)的空間自相關(guān)性,適合使用克里格法進(jìn)行空間插值。

    表7 GR和GWR模型回歸殘差的半方差函數(shù)擬合結(jié)果

    根據(jù)上述半方差函數(shù)的擬合結(jié)果,采用OK法分別對(duì)2種模型的回歸殘差進(jìn)行了空間插值。根據(jù)2個(gè)主成分變量及其回歸系數(shù)、殘差插值等數(shù)據(jù)的計(jì)算結(jié)果,分別獲得了基于GRK模型和GWRK模型的137Cs本底計(jì)算值空間分布圖(圖3)。

    3.5 模型預(yù)測(cè)精度評(píng)價(jià)

    作為對(duì)照,同時(shí)采用OK法對(duì)各站點(diǎn)的137Cs本底計(jì)算值進(jìn)行了空間插值。通過(guò)構(gòu)建的驗(yàn)證數(shù)據(jù)集分別評(píng)價(jià)了OK、GRK和GWRK等模型的預(yù)測(cè)精度。各評(píng)價(jià)指標(biāo)的值越接近0,說(shuō)明相應(yīng)模型的預(yù)測(cè)精度越高。由表8的分析結(jié)果可以看出,在各研究區(qū)域,模型預(yù)測(cè)精度均為GWRK>GRK>OK。綜合考慮了降水、空間位置等多個(gè)影響因素的GRK模型和GWRK模型對(duì)137Cs本底值的預(yù)測(cè)精度較OK法均有較大提高,而進(jìn)一步探究了各影響因素空間非平穩(wěn)性的GWRK模型的預(yù)測(cè)效果是這3種模型中最好的。

    (a) H30區(qū)域GWRK結(jié)果 Result of GWRK model in H30;(b) H30區(qū)域GRK結(jié)果 Result of GRK model in H30;(c) L30區(qū)域GWRK結(jié)果 Result of GWRK model in L30;(d) L30區(qū)域GRK結(jié)果 Result of GRK model in L30圖3 四川省137Cs本底計(jì)算值空間分布Fig.3 Spatial distribution of 137Cs reference inventory of Sichuan province

    表8 不同預(yù)測(cè)方法的預(yù)測(cè)精度對(duì)比

    Notes: ME stand for mean error, MAE stand for mean absolute error, RMSE stand for root mean square error, MARE stand for mean absolute relative error, and RMSRE stand for root mean square relative error.

    同時(shí)收集了已發(fā)表文獻(xiàn)中137Cs本底值的實(shí)測(cè)數(shù)據(jù),與GWRK模型的預(yù)測(cè)結(jié)果進(jìn)行了對(duì)比分析。從表9的對(duì)比結(jié)果可以看出,GWRK模型的預(yù)測(cè)值普遍略低于實(shí)測(cè)值,差異幅度大致為10%??紤]到測(cè)定137Cs活度的γ能譜儀的測(cè)試誤差普遍在5%~10%之間,我們?nèi)∑渖舷?0%作為實(shí)測(cè)本底值和預(yù)測(cè)值是否存在差異的標(biāo)準(zhǔn)。由此可見(jiàn),雖然預(yù)測(cè)值略低于實(shí)測(cè)值,但差異不大,實(shí)測(cè)值與本文的預(yù)測(cè)結(jié)果是基本吻合的。

    3.6 四川省137Cs本底值的空間分布

    由圖3的插值結(jié)果可以看出,通過(guò)GRK和GWRK法得到的四川省137Cs本底計(jì)算值空間分布的總趨勢(shì)基本一致,均呈現(xiàn)北高南低、東高西低及高原低盆地高、盆周山地高盆中丘陵低的趨勢(shì),該預(yù)測(cè)結(jié)果與前人對(duì)全球137Cs沉降規(guī)律的認(rèn)識(shí)基本一致[7]。一些局部的、細(xì)節(jié)性的137Cs本底計(jì)算值空間分布特征在插值結(jié)果中也有所反映:在四川盆地西部邊緣出現(xiàn)了1個(gè)局部137Cs本底計(jì)算值峰值區(qū),大致位于雅安市、樂(lè)山一線,137Cs本底計(jì)算值達(dá)到了1 300~1 400 Bq/m2。出現(xiàn)這一現(xiàn)象主要是由于該區(qū)的“華西雨屏”[19]現(xiàn)象造成降水量較高,而137Cs主要來(lái)源于降水帶來(lái)的濕沉降所致。此外,由于W&H計(jì)算模型的局限——模型中按照10°的緯度帶劃分計(jì)算區(qū)域,造成北緯30°兩側(cè)區(qū)域的137Cs本底計(jì)算值存在“斷崖式”突變。在資料精度允許的條件下,今后可以通過(guò)細(xì)分基本區(qū)域,改進(jìn)模型,提高模擬結(jié)果的精度和可信度[7]。

    表9 137Cs本底值實(shí)測(cè)值與GWRK法預(yù)測(cè)值比較

    4 結(jié)論

    以地處長(zhǎng)江上游、土壤侵蝕特征具有典型性的四川省為研究對(duì)象,基于147個(gè)氣象站點(diǎn)的實(shí)測(cè)數(shù)據(jù)以及W&H模型的137Cs本底值分布模型,采用GWRK法進(jìn)行空間插值,獲得四川省137Cs本底計(jì)算值的空間分布數(shù)據(jù)。初步研究結(jié)果表明:

    1)應(yīng)用PCA法將經(jīng)緯度、降水量等3個(gè)影響因素轉(zhuǎn)換為2個(gè)主成分變量,主成分變量能表達(dá)原始數(shù)據(jù)中90%以上的信息,同時(shí)有效的消除了變量間的多重共線性,為后續(xù)的空間插值分析奠定了基礎(chǔ)。

    2)根據(jù)PCA的分析結(jié)果,通過(guò)基于OLS的GRK和GWRK模型對(duì)四川省137Cs本底計(jì)算值空間分布進(jìn)行了預(yù)測(cè),并評(píng)價(jià)了不同模型的預(yù)測(cè)精度?;貧w分析的結(jié)果表明,基于局部回歸的GWRK模型比GRK模型能更好地解釋137Cs本底計(jì)算值與各影響因素間相互關(guān)系的空間非平穩(wěn)性,具有更好的擬合效果;從最終的插值結(jié)果來(lái)看,GWRK模型綜合考慮了降水、空間位置等影響因素,并進(jìn)一步探究了各影響因素對(duì)137Cs本底值影響的空間非平穩(wěn)性,能更好地揭示影響因素和137Cs本底值間的相互關(guān)系,其預(yù)測(cè)精度最高,且與實(shí)測(cè)結(jié)果較吻合。

    3)以實(shí)測(cè)數(shù)據(jù)為基礎(chǔ)、基于GIS空間插值技術(shù)獲取的四川省137Cs本底計(jì)算值數(shù)據(jù)具有較高空間分辨率,并能比較準(zhǔn)確地表達(dá)局部區(qū)域137Cs本底值發(fā)生突變的細(xì)節(jié)信息。筆者提出的PCA+GWRK插值法可以有效的揭示137Cs本底值的空間分布規(guī)律及不同因素對(duì)其的影響作用,對(duì)運(yùn)用137Cs示蹤技術(shù)開(kāi)展侵蝕泥沙研究具有一定的借鑒作用,可以為水土流失治理、土壤退化防治及水土保持工程等工作提供有效的技術(shù)支持。

    猜你喜歡
    區(qū)域模型研究
    一半模型
    FMS與YBT相關(guān)性的實(shí)證研究
    遼代千人邑研究述論
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    視錯(cuò)覺(jué)在平面設(shè)計(jì)中的應(yīng)用與研究
    科技傳播(2019年22期)2020-01-14 03:06:54
    EMA伺服控制系統(tǒng)研究
    3D打印中的模型分割與打包
    關(guān)于四色猜想
    分區(qū)域
    亚洲精品日韩av片在线观看| .国产精品久久| 亚洲国产色片| 久久久国产成人精品二区| 首页视频小说图片口味搜索| 国产色爽女视频免费观看| 一个人看视频在线观看www免费| 国产爱豆传媒在线观看| 国产一区二区在线av高清观看| 亚洲五月天丁香| 国产成人欧美在线观看| 色综合婷婷激情| 国产精品一区二区性色av| 自拍偷自拍亚洲精品老妇| 国产精品久久久久久久电影| 亚洲,欧美,日韩| 91九色精品人成在线观看| 欧美乱色亚洲激情| 国产熟女xx| 成人高潮视频无遮挡免费网站| 看十八女毛片水多多多| 欧美最新免费一区二区三区 | 国产av一区在线观看免费| 高清在线国产一区| 男女那种视频在线观看| 9191精品国产免费久久| 国产精品国产高清国产av| 亚洲经典国产精华液单 | а√天堂www在线а√下载| 午夜福利在线观看吧| 美女cb高潮喷水在线观看| 可以在线观看的亚洲视频| 亚洲国产精品sss在线观看| 久久香蕉精品热| 欧美+日韩+精品| 欧美绝顶高潮抽搐喷水| 国产精品国产高清国产av| 欧美最新免费一区二区三区 | 黄色女人牲交| 亚洲av.av天堂| 麻豆国产97在线/欧美| 欧美日韩中文字幕国产精品一区二区三区| 我要看日韩黄色一级片| 老鸭窝网址在线观看| 国产亚洲精品久久久com| 99精品久久久久人妻精品| 久久久久国内视频| 成人无遮挡网站| 中国美女看黄片| 老鸭窝网址在线观看| 国产精华一区二区三区| xxxwww97欧美| 国产日本99.免费观看| 欧美xxxx性猛交bbbb| 99久国产av精品| 精品人妻一区二区三区麻豆 | 色综合站精品国产| 免费人成视频x8x8入口观看| 国产亚洲精品久久久com| 色视频www国产| 国产高清三级在线| 欧美激情国产日韩精品一区| 国产精品一区二区性色av| 丁香欧美五月| 国产免费av片在线观看野外av| 少妇丰满av| 免费高清视频大片| 成人性生交大片免费视频hd| 赤兔流量卡办理| 国产精品影院久久| 日本免费a在线| 首页视频小说图片口味搜索| 亚洲经典国产精华液单 | 99热精品在线国产| 大型黄色视频在线免费观看| 亚洲av电影在线进入| 很黄的视频免费| 免费av观看视频| 人人妻人人澡欧美一区二区| 两个人视频免费观看高清| 中出人妻视频一区二区| 亚洲经典国产精华液单 | 国产69精品久久久久777片| 成年女人毛片免费观看观看9| 国产单亲对白刺激| 在线看三级毛片| 蜜桃久久精品国产亚洲av| 美女被艹到高潮喷水动态| 日本熟妇午夜| 亚洲av二区三区四区| 国产精品综合久久久久久久免费| 深夜a级毛片| АⅤ资源中文在线天堂| 观看美女的网站| 18禁裸乳无遮挡免费网站照片| 免费在线观看亚洲国产| av在线老鸭窝| 日韩精品中文字幕看吧| 亚洲精品色激情综合| 在现免费观看毛片| 在线播放无遮挡| 亚洲国产精品999在线| 嫩草影院新地址| 在线看三级毛片| 熟妇人妻久久中文字幕3abv| 久久99热这里只有精品18| 桃红色精品国产亚洲av| 午夜视频国产福利| 18禁黄网站禁片免费观看直播| 俄罗斯特黄特色一大片| 校园春色视频在线观看| 国产精品亚洲一级av第二区| 国产精品爽爽va在线观看网站| 亚洲一区高清亚洲精品| 在线观看66精品国产| 美女黄网站色视频| 欧美色欧美亚洲另类二区| 久久久久久久久久成人| 欧美午夜高清在线| 国内精品久久久久久久电影| 国产精品98久久久久久宅男小说| 久久精品国产清高在天天线| 久久九九热精品免费| 成人欧美大片| 18禁黄网站禁片免费观看直播| 日本 欧美在线| 最近最新中文字幕大全电影3| 国产免费一级a男人的天堂| 麻豆国产97在线/欧美| 99久久精品国产亚洲精品| 搞女人的毛片| 在线免费观看的www视频| 精品一区二区三区视频在线| 91狼人影院| 中文字幕人成人乱码亚洲影| 国产欧美日韩一区二区精品| 欧美精品啪啪一区二区三区| 亚洲18禁久久av| 日韩 亚洲 欧美在线| 韩国av一区二区三区四区| 一级a爱片免费观看的视频| 久久午夜福利片| 久久久久久久亚洲中文字幕 | 特大巨黑吊av在线直播| .国产精品久久| 一二三四社区在线视频社区8| 午夜激情福利司机影院| 亚洲av二区三区四区| 国产精品一区二区三区四区久久| 欧美丝袜亚洲另类 | 久久精品国产自在天天线| 中亚洲国语对白在线视频| 成人午夜高清在线视频| 色尼玛亚洲综合影院| 国内精品美女久久久久久| 九九热线精品视视频播放| 91麻豆av在线| 1000部很黄的大片| 天美传媒精品一区二区| 高清在线国产一区| 好男人电影高清在线观看| 最近视频中文字幕2019在线8| 精品午夜福利视频在线观看一区| 悠悠久久av| 欧美日韩乱码在线| 久久久久精品国产欧美久久久| 色哟哟哟哟哟哟| 噜噜噜噜噜久久久久久91| 如何舔出高潮| 日韩欧美 国产精品| 俄罗斯特黄特色一大片| 久久久久久久亚洲中文字幕 | 别揉我奶头 嗯啊视频| 国产精品一区二区性色av| 一级毛片久久久久久久久女| 亚洲人成网站高清观看| 91av网一区二区| 日韩欧美精品免费久久 | 啪啪无遮挡十八禁网站| 美女被艹到高潮喷水动态| 18禁黄网站禁片免费观看直播| 久久久久九九精品影院| 特级一级黄色大片| 麻豆国产97在线/欧美| 精品午夜福利视频在线观看一区| 日本黄大片高清| 日本免费a在线| 欧美国产日韩亚洲一区| 国产三级在线视频| 欧美乱妇无乱码| 小蜜桃在线观看免费完整版高清| 亚洲综合色惰| 蜜桃亚洲精品一区二区三区| 日韩欧美 国产精品| 亚洲七黄色美女视频| 国产成人a区在线观看| 亚洲乱码一区二区免费版| 嫩草影院入口| 久久久久免费精品人妻一区二区| 亚洲av.av天堂| 老司机午夜十八禁免费视频| 亚洲精品在线美女| 欧美黑人巨大hd| 成人一区二区视频在线观看| 99久久精品热视频| 男人狂女人下面高潮的视频| 全区人妻精品视频| 岛国在线免费视频观看| 少妇的逼水好多| 国产精品av视频在线免费观看| 可以在线观看的亚洲视频| 激情在线观看视频在线高清| 噜噜噜噜噜久久久久久91| 久久久久精品国产欧美久久久| 国产精品久久久久久亚洲av鲁大| 久久久色成人| 亚洲经典国产精华液单 | 成人国产综合亚洲| 国产精品精品国产色婷婷| www.熟女人妻精品国产| 午夜福利在线在线| 久久九九热精品免费| 内射极品少妇av片p| 美女黄网站色视频| 免费黄网站久久成人精品 | 极品教师在线视频| 神马国产精品三级电影在线观看| 欧美xxxx黑人xx丫x性爽| 成人精品一区二区免费| 国产私拍福利视频在线观看| 亚洲avbb在线观看| 日本与韩国留学比较| 十八禁人妻一区二区| 99久久99久久久精品蜜桃| 国产蜜桃级精品一区二区三区| x7x7x7水蜜桃| 色哟哟哟哟哟哟| 亚洲欧美日韩卡通动漫| 色5月婷婷丁香| 国产成年人精品一区二区| 首页视频小说图片口味搜索| 精品人妻偷拍中文字幕| 一本久久中文字幕| 免费av不卡在线播放| 综合色av麻豆| 欧美一级a爱片免费观看看| 美女大奶头视频| 成年人黄色毛片网站| 国产探花极品一区二区| 亚洲国产欧美人成| 国产精品一区二区免费欧美| 午夜福利欧美成人| 少妇人妻精品综合一区二区 | 久久精品综合一区二区三区| 国产又黄又爽又无遮挡在线| 中文字幕人成人乱码亚洲影| xxxwww97欧美| 97碰自拍视频| 中国美女看黄片| 内射极品少妇av片p| 国产精品三级大全| 亚洲人成网站高清观看| 久久精品夜夜夜夜夜久久蜜豆| 午夜精品久久久久久毛片777| 黄色日韩在线| 国产高清有码在线观看视频| 美女被艹到高潮喷水动态| 亚洲综合色惰| 男人狂女人下面高潮的视频| 在线观看av片永久免费下载| 好男人在线观看高清免费视频| 国产白丝娇喘喷水9色精品| or卡值多少钱| a级一级毛片免费在线观看| 九九久久精品国产亚洲av麻豆| 在线看三级毛片| 三级国产精品欧美在线观看| 国产成人影院久久av| 亚洲最大成人中文| 免费大片18禁| 国产午夜精品论理片| 国产欧美日韩精品亚洲av| 亚洲不卡免费看| 无遮挡黄片免费观看| 麻豆国产av国片精品| 国产黄片美女视频| 国产伦人伦偷精品视频| 国产熟女xx| 国产一区二区在线av高清观看| 国产 一区 欧美 日韩| 午夜福利欧美成人| 人妻久久中文字幕网| av在线观看视频网站免费| 男人和女人高潮做爰伦理| 在线a可以看的网站| 日本a在线网址| 精品久久久久久久人妻蜜臀av| 亚洲av免费在线观看| 成人亚洲精品av一区二区| 自拍偷自拍亚洲精品老妇| 精品人妻熟女av久视频| 一进一出抽搐gif免费好疼| 18禁裸乳无遮挡免费网站照片| 又黄又爽又免费观看的视频| 制服丝袜大香蕉在线| 少妇的逼水好多| 十八禁网站免费在线| 亚洲第一电影网av| 热99在线观看视频| 亚洲国产精品久久男人天堂| 午夜a级毛片| 麻豆av噜噜一区二区三区| 99久久久亚洲精品蜜臀av| 在线观看午夜福利视频| 三级男女做爰猛烈吃奶摸视频| 日本免费a在线| 18+在线观看网站| 99久久99久久久精品蜜桃| 国产av不卡久久| 尤物成人国产欧美一区二区三区| 他把我摸到了高潮在线观看| 中文字幕免费在线视频6| 亚洲自拍偷在线| 欧美xxxx性猛交bbbb| 每晚都被弄得嗷嗷叫到高潮| 性色avwww在线观看| 欧美激情久久久久久爽电影| 久久6这里有精品| 国产精品99久久久久久久久| 国内毛片毛片毛片毛片毛片| 欧美乱色亚洲激情| 男女之事视频高清在线观看| 久久久久免费精品人妻一区二区| 国内毛片毛片毛片毛片毛片| 欧美绝顶高潮抽搐喷水| 91在线精品国自产拍蜜月| 国产成+人综合+亚洲专区| 欧美日韩中文字幕国产精品一区二区三区| 成人精品一区二区免费| 99在线人妻在线中文字幕| 午夜久久久久精精品| 免费看a级黄色片| 欧美成人一区二区免费高清观看| 午夜久久久久精精品| 99久久久亚洲精品蜜臀av| 成人特级av手机在线观看| 欧美3d第一页| 三级国产精品欧美在线观看| 深爱激情五月婷婷| 亚洲欧美精品综合久久99| 精品不卡国产一区二区三区| 床上黄色一级片| 欧美色欧美亚洲另类二区| 一个人看视频在线观看www免费| 国产精品av视频在线免费观看| 熟女人妻精品中文字幕| 在线观看舔阴道视频| 国产精品av视频在线免费观看| 精品熟女少妇八av免费久了| 麻豆一二三区av精品| 亚洲精品亚洲一区二区| 一个人观看的视频www高清免费观看| 国产综合懂色| 最好的美女福利视频网| 啪啪无遮挡十八禁网站| 日日夜夜操网爽| 精品人妻视频免费看| 美女cb高潮喷水在线观看| 亚洲最大成人中文| 国产高清视频在线观看网站| 免费黄网站久久成人精品 | 欧美色欧美亚洲另类二区| 此物有八面人人有两片| 99久久精品国产亚洲精品| 午夜免费成人在线视频| 欧美中文日本在线观看视频| 欧美bdsm另类| 国产精品久久电影中文字幕| 国产伦精品一区二区三区视频9| 国产精品久久电影中文字幕| 一进一出好大好爽视频| 99热这里只有是精品50| 精品久久久久久久久av| 亚洲中文日韩欧美视频| 免费大片18禁| 亚洲最大成人中文| 成人特级av手机在线观看| av黄色大香蕉| 婷婷精品国产亚洲av在线| 亚洲精品影视一区二区三区av| a级一级毛片免费在线观看| 宅男免费午夜| 国产单亲对白刺激| 日本黄大片高清| 久久亚洲真实| 色5月婷婷丁香| 色哟哟·www| 尤物成人国产欧美一区二区三区| 少妇人妻一区二区三区视频| 午夜精品久久久久久毛片777| 日本成人三级电影网站| 国产亚洲欧美在线一区二区| av在线老鸭窝| 日本 av在线| 亚洲成a人片在线一区二区| 两个人的视频大全免费| 天美传媒精品一区二区| 亚洲中文日韩欧美视频| 长腿黑丝高跟| 中文字幕人妻熟人妻熟丝袜美| 成人三级黄色视频| 亚洲最大成人手机在线| 两性午夜刺激爽爽歪歪视频在线观看| 精品久久国产蜜桃| 别揉我奶头 嗯啊视频| 国产一区二区三区视频了| 91久久精品国产一区二区成人| 成人欧美大片| 91av网一区二区| 久99久视频精品免费| 老司机深夜福利视频在线观看| 欧美黄色片欧美黄色片| 免费在线观看成人毛片| 婷婷亚洲欧美| 久久久久久国产a免费观看| 一本精品99久久精品77| 一区二区三区高清视频在线| 欧美日韩乱码在线| 蜜桃久久精品国产亚洲av| 日韩欧美 国产精品| 久久久久久久精品吃奶| .国产精品久久| 亚洲最大成人手机在线| 波多野结衣高清无吗| 国产成人欧美在线观看| 久久99热这里只有精品18| av在线老鸭窝| 九九热线精品视视频播放| 成人永久免费在线观看视频| 99国产精品一区二区蜜桃av| 小说图片视频综合网站| 久久久国产成人免费| 此物有八面人人有两片| 精品日产1卡2卡| 黄色女人牲交| 亚洲在线观看片| 成人特级av手机在线观看| 看黄色毛片网站| 宅男免费午夜| 精品99又大又爽又粗少妇毛片 | 少妇高潮的动态图| 久久精品国产清高在天天线| 中文字幕精品亚洲无线码一区| 亚洲第一区二区三区不卡| 俄罗斯特黄特色一大片| 精品国内亚洲2022精品成人| 久久热精品热| 亚洲aⅴ乱码一区二区在线播放| 99久久精品一区二区三区| 国产午夜精品久久久久久一区二区三区 | 精品乱码久久久久久99久播| 人人妻,人人澡人人爽秒播| 夜夜爽天天搞| 18禁黄网站禁片午夜丰满| 最近最新免费中文字幕在线| 久久久久国产精品人妻aⅴ院| 丰满人妻一区二区三区视频av| 精品久久久久久成人av| 男人舔奶头视频| 国产美女午夜福利| 男人舔奶头视频| 日韩 亚洲 欧美在线| 人妻夜夜爽99麻豆av| 最后的刺客免费高清国语| 内射极品少妇av片p| eeuss影院久久| 欧美色欧美亚洲另类二区| 精品人妻偷拍中文字幕| 国产真实乱freesex| 麻豆久久精品国产亚洲av| 久久久久久久精品吃奶| 久久精品国产亚洲av涩爱 | 国产高清有码在线观看视频| 国产成+人综合+亚洲专区| 亚洲狠狠婷婷综合久久图片| a级一级毛片免费在线观看| 床上黄色一级片| 免费在线观看日本一区| 少妇被粗大猛烈的视频| 国产精品久久久久久久久免 | а√天堂www在线а√下载| 亚洲午夜理论影院| 久久人人爽人人爽人人片va | 一区二区三区高清视频在线| 青草久久国产| 亚洲电影在线观看av| 国产成人影院久久av| 熟女电影av网| 国产麻豆成人av免费视频| 中文字幕人成人乱码亚洲影| 高潮久久久久久久久久久不卡| 最近视频中文字幕2019在线8| h日本视频在线播放| 一级作爱视频免费观看| 97人妻精品一区二区三区麻豆| 日本在线视频免费播放| 欧美日韩中文字幕国产精品一区二区三区| 国产精品一区二区三区四区免费观看 | 国产精品久久久久久亚洲av鲁大| 欧洲精品卡2卡3卡4卡5卡区| 国产三级黄色录像| 看十八女毛片水多多多| 能在线免费观看的黄片| 少妇丰满av| 脱女人内裤的视频| 性色av乱码一区二区三区2| 国产精品久久久久久久电影| 夜夜夜夜夜久久久久| 91麻豆精品激情在线观看国产| a级一级毛片免费在线观看| 最近中文字幕高清免费大全6 | 中文字幕av成人在线电影| 欧美极品一区二区三区四区| 少妇的逼好多水| 成人高潮视频无遮挡免费网站| 午夜精品在线福利| 一级a爱片免费观看的视频| 精品一区二区三区av网在线观看| 免费观看人在逋| 精品国产亚洲在线| 久久精品国产99精品国产亚洲性色| 国产欧美日韩一区二区三| 成人一区二区视频在线观看| 日韩av在线大香蕉| 国内少妇人妻偷人精品xxx网站| 蜜桃亚洲精品一区二区三区| 免费高清视频大片| 日韩精品中文字幕看吧| 亚洲人成电影免费在线| 久久亚洲精品不卡| 天堂av国产一区二区熟女人妻| 日韩av在线大香蕉| 亚洲 国产 在线| 男人狂女人下面高潮的视频| 丁香欧美五月| 在线观看av片永久免费下载| 嫩草影院新地址| 国产熟女xx| 免费看美女性在线毛片视频| 久久精品国产亚洲av涩爱 | 欧美在线一区亚洲| 一个人看的www免费观看视频| 我要搜黄色片| 国产激情偷乱视频一区二区| 国产成人福利小说| 午夜免费男女啪啪视频观看 | 欧美三级亚洲精品| 男插女下体视频免费在线播放| 免费在线观看影片大全网站| 成人性生交大片免费视频hd| 婷婷精品国产亚洲av在线| 中出人妻视频一区二区| 久久精品久久久久久噜噜老黄 | 级片在线观看| 97碰自拍视频| 91麻豆av在线| 他把我摸到了高潮在线观看| 国产探花在线观看一区二区| 国产一区二区亚洲精品在线观看| 国产精品电影一区二区三区| 精品久久久久久久久av| 18+在线观看网站| 美女xxoo啪啪120秒动态图 | 国产精品三级大全| 搡老熟女国产l中国老女人| 久久性视频一级片| 中文字幕人成人乱码亚洲影| 国产在线男女| 高清毛片免费观看视频网站| 亚洲美女黄片视频| 国产探花在线观看一区二区| 老鸭窝网址在线观看| 嫩草影视91久久| 观看免费一级毛片| 日日干狠狠操夜夜爽| 大型黄色视频在线免费观看| 一本精品99久久精品77| 亚洲av第一区精品v没综合| 1000部很黄的大片| 日本在线视频免费播放| 国产精品女同一区二区软件 | 国产精品三级大全| 欧美日韩综合久久久久久 | 国产又黄又爽又无遮挡在线| 嫩草影院新地址| 成人国产一区最新在线观看| 真实男女啪啪啪动态图| 欧美zozozo另类| 精品人妻1区二区| 人妻丰满熟妇av一区二区三区| 欧美黑人巨大hd| 午夜影院日韩av| 日日夜夜操网爽| 最后的刺客免费高清国语| 亚洲成人中文字幕在线播放| 中文字幕高清在线视频| 俺也久久电影网| 久久精品国产自在天天线| 99久久99久久久精品蜜桃| 久久午夜亚洲精品久久| 草草在线视频免费看| 欧美+日韩+精品| 可以在线观看的亚洲视频| 黄片小视频在线播放| 白带黄色成豆腐渣|