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

    山東煙臺養(yǎng)馬島南部地區(qū)土壤重金屬污染風(fēng)險及生態(tài)風(fēng)險評價

    2024-07-11 07:22:38湯世凱扈勝濤黎奕宏強萌麟于林弘張祥恒王璨
    山東國土資源 2024年6期
    關(guān)鍵詞:土壤重金屬相關(guān)性分析

    湯世凱 扈勝濤 黎奕宏 強萌麟 于林弘 張祥恒 王璨

    收稿日期:20231205;修訂日期:20240117;編輯:陶衛(wèi)衛(wèi)

    基金項目:山東省自然資源廳部省協(xié)議地質(zhì)勘查項目“山東省煙臺市養(yǎng)馬島及鄰近海域綜合地質(zhì)調(diào)查”(魯勘字〔2023〕38號);山東省第三地質(zhì)礦產(chǎn)勘查院科技創(chuàng)新項目“海灣沉積物粒度差異性分析研究——以丁字灣、芝罘灣為例”魯三地勘發(fā)〔2022〕45號

    作者簡介:湯世凱(1988—),男,山東東明人,高級工程師,主要從事海洋地質(zhì)調(diào)查研究;Email:tsktop@163.com? *通訊作者:黎奕宏(1991—),女,山東海陽人,工程師,主要從事海洋地質(zhì)調(diào)查研究;Email:553963083@qq.com

    摘要:通過對煙臺養(yǎng)馬島南部第四系分布區(qū)土壤樣品中的重金屬(Cr、Ni、Cu、Zn、As、Cd、Pb和Hg)含量測試分析,運用統(tǒng)計分析、地質(zhì)累積指數(shù)法和潛在生態(tài)風(fēng)險指數(shù)法對元素的相關(guān)性、污染程度以及生態(tài)風(fēng)險進行研究評價。結(jié)果顯示,受一定程度的人類活動影響,除Ni外,高值區(qū)主要分布在第四系中—西南部;除Cr和Cu外,元素含量平均值均超過背景值。相關(guān)性分析表明,Zn與Hg,Zn與Cu以及Cu、Cr、As相互間具較高正相關(guān)性。Hg在TR2點位具有異常高值,屬于中度污染和風(fēng)險程度重。評價各元素污染程度Hg>As>Zn>Cd>Ni>Pb>Cu>Cr,風(fēng)險程度Hg>Cd>As>Ni>Pb>Cu>Cr>Zn,綜合潛在生態(tài)風(fēng)險程度評價TR2點位為極重程度,其余為中等程度—低程度。

    關(guān)鍵詞:土壤重金屬;相關(guān)性分析;生態(tài)風(fēng)險評價;養(yǎng)馬島南部;山東煙臺

    中圖分類號:X53??? 文獻標識碼:A??? doi:10.12128/j.issn.16726979.2024.06.002

    引文格式:湯世凱,扈勝濤,黎奕宏,等.山東煙臺養(yǎng)馬島南部地區(qū)土壤重金屬污染風(fēng)險及生態(tài)風(fēng)險評價[J].山東國土資源,2024,40(6):1117. TANG Shikai, HU Shengtao, LI Yihong, et al. Evaluation of Heavy Metal Contamination and Ecological Risk of Soil in Southern Yangmadao in Yantai City in Shandong Province[J].Shandong Land and Resources,2024,40(6):1117.

    0? 引言

    海島因其狹小的面積,單一的地域結(jié)構(gòu),突顯出鮮明的生態(tài)脆弱性,極易遭受外界干擾破壞生態(tài)環(huán)境[12]。海島土壤作為生物生存的載體,在元素循環(huán)交換、維持生態(tài)環(huán)境平衡方面具有重要作用[34]。

    隨著人類對海島的開發(fā)利用,海島土壤物理、化學(xué)等性質(zhì)不可避免地受到影響。其中的土壤元素尤其是重金屬元素含量會持續(xù)積累,因其毒性較強、不易分解、沿食物鏈傳遞等特點,一定程度上損害人類的身體健康[58]。

    養(yǎng)馬島作為旅游度假區(qū),以往對該島的研究工作偏重于旅游資源的開發(fā)、評價等方面[910],對海島生態(tài)環(huán)境影響因素及變化規(guī)律等研究較少。本次工作通過對養(yǎng)馬島第四系分布區(qū)進行系統(tǒng)采樣,分析土壤中的重金屬含量,運用數(shù)據(jù)統(tǒng)計分析、地累積指數(shù)法和潛在生態(tài)風(fēng)險指數(shù)法等方法,研究養(yǎng)馬島土壤重金屬的分布特征和污染情況,以期為養(yǎng)馬島土壤污染治理、合理利用與保護土地資源,營造良好海島生態(tài)環(huán)境提供重要的地質(zhì)參考依據(jù)。

    1? 研究區(qū)概況

    養(yǎng)馬島位于山東省煙臺市牟平區(qū)北部的黃海之中,距離牟平城區(qū)9km,面積約13km2,地形北高南低,地貌類型以低丘陵為主,屬于北溫帶季風(fēng)型大陸氣候[11]。人類開發(fā)利用活動主要在第四系分布區(qū)最為劇烈,交通運輸頻繁,城鎮(zhèn)化建設(shè)程度較高。

    養(yǎng)馬島中-北部為基巖發(fā)育區(qū),面積約9km2,巖性為荊山群斜長片麻巖、斜長角閃巖、大理巖等;南部為第四系發(fā)育區(qū),面積約4km2,巖性為旭口組中細砂,本次研究對象為南部第四系分布區(qū)。

    2? 工作方法

    2.1? 樣品采集與分析

    為了查清養(yǎng)馬島南部第四系分布區(qū)土壤地質(zhì)環(huán)境質(zhì)量,按照大致均勻布設(shè)點位的原則,2023年7月在養(yǎng)馬島采集表層土壤樣品10件,采樣點位見圖1。采樣過程中首先除去表層浮土,然后將約20cm深度內(nèi)的土壤混合后裝入聚乙烯袋中密封,送往實驗室測試分析,每件土壤樣品重量約1.5kg。

    土壤化驗分析由山東省第三地質(zhì)礦產(chǎn)勘查院承擔(dān),樣品自然風(fēng)干后,進行研磨篩分,然后對土壤中的Cr、Ni、Cu、Zn、As、Cd、Pb和Hg等8種重金屬元素進行分析。Cr、Ni、Cu、Zn、As、Cd和Pb等7種重金屬元素的操作分析方法、質(zhì)量控制按照《土壤和沉積物12種金屬元素的測定王水提取-電感耦合等離子體質(zhì)譜法》(HJ803—2016)的要求進行,元素檢出限分別為2mg/kg、2mg/kg、0.5mg/kg、7mg/kg、0.6mg/kg、0.07mg/kg、2mg/kg。Hg的操作分析方法及質(zhì)量控制按照《土壤和沉積物 汞、砷、硒、鉍、銻的測定微波消解/原子熒光法》(HJ680—2013)的要求進行,檢出限為0.002mg/kg。為保證土壤樣品分析數(shù)據(jù)的準確性,按照《中國地質(zhì)調(diào)查局地質(zhì)調(diào)查技術(shù)標準·多目標區(qū)域地球化學(xué)調(diào)查規(guī)范》的要求采取重復(fù)樣品分析、對異常點重復(fù)性檢驗等方法進行嚴格的質(zhì)量監(jiān)控。

    2.2? 重金屬污染和生態(tài)風(fēng)險評價方法

    2.2.1? 地質(zhì)累積指數(shù)法

    由德國科學(xué)家Muller在1969年提出[12],后續(xù)被廣泛應(yīng)用[1314],該方法主要是參考沉積物背景值對當(dāng)前重金屬含量進行污染評價。評價公式如下:Igeo=log2[Cn/(K×Bn)],其中Cn為重金屬n的實測含量;K為因地區(qū)巖性差異引起的背景值的波動,是人為取的一個相關(guān)系數(shù),通常取K=1.5;Bn是沉積母質(zhì)中該元素的地球化學(xué)背景值,本文采用牟平區(qū)土壤地球化學(xué)背景值[15]作為參照進行計算。

    根據(jù)計算的地質(zhì)累積指數(shù)(Igeo)數(shù)值大小,參考地質(zhì)累積指數(shù)與污染程度分級對照表(表1),最終對土壤的重金屬污染程度劃分級別。

    2.2.2? 潛在生態(tài)風(fēng)險指數(shù)

    由瑞典科學(xué)家Hakanson在1980年提出[16],該方法綜合考慮了每一種重金屬元素的污染指數(shù)、毒性響應(yīng)系數(shù),最終反映的是多種元素綜合效應(yīng)。因而是目前國內(nèi)外最為常用的生態(tài)風(fēng)險評價方法[1720]。評價公式如式(1):

    Cif=Ci/Cin;Ei r=Tir·Cif;RI=∑Eir(1)

    式中:Ci為重金屬元素i的實測含量;Cin是重金屬元素i的地球化學(xué)背景值,本次評價仍采用煙臺市該元素土壤地球化學(xué)背景值參與計算;Cif為重金屬元素i的污染系數(shù);Tir為重金屬i的毒性響應(yīng)系數(shù),其中Cr、Ni、Cu、Zn、As、Cd、Pb和Hg的毒性響應(yīng)系數(shù)分別為2、5、5、1、10、30、5、40[15];Eir為重金屬i的潛在生態(tài)風(fēng)險指數(shù);RI為參與評價的所有重金屬元素的綜合潛在生態(tài)風(fēng)險指數(shù)。單元素的潛在生態(tài)風(fēng)險系數(shù)(Eir)和多元素的綜合潛在生態(tài)風(fēng)險指數(shù)(RI)分級標準見表2[21]。

    3? 結(jié)果與討論

    3.1? 重金屬含量及分布

    養(yǎng)馬島南部第四系分布區(qū)土壤重金屬含量統(tǒng)計結(jié)果(表3)表明,Cr、Ni、Cu、Zn、As、Cd、Pb和Hg的含量平均值分別為23.54mg/kg、20.18mg/kg、16.28mg/kg、59.08mg/kg、8.14mg/kg、0.12mg/kg、32.96mg/kg和0.15mg/kg。經(jīng)統(tǒng)計,Cr和Cu分別有10%和40%點位含量超過背景值,平均含量低于背景值;Ni、Zn、As、Cd、Pb和Hg分別有60%、50%、80%、50%、40%和70%點位含量超過背景值,平均含量超過背景值,其中Hg元素含量平均值超背景值倍數(shù)最高,為6倍。各重金屬元素中,As元素變異系數(shù)最小,為44%,說明在各點位之間As含量差別相對較小,分布較為均勻,Hg變異系數(shù)最大,為209%,表明Hg元素分布非常不均勻。

    對各重金屬元素含量作等值線分布圖(圖2),可以看出,養(yǎng)馬島第四系土壤重金屬Cr、Ni、Cu、Zn、As、Cd、Pb和Hg的平面分布具有如下特征:

    (1)土壤中Zn、Hg和Cd的分布規(guī)律較為相似(圖2a、圖2d和圖2g),在養(yǎng)馬島大橋以西,即第四系西南分布區(qū)出現(xiàn)含量高值區(qū),在大橋以東的第四系中、東北分布區(qū)為含量低值區(qū)。

    (2)土壤中Pb、Cu、Cr和As的分布規(guī)律較為類似(圖2b、圖2e、圖2f和圖2h),在養(yǎng)馬島大橋兩側(cè),即第四系中—西南分布區(qū)出現(xiàn)含量高值區(qū),第四系中—東北分布區(qū)出現(xiàn)含量低值區(qū)。

    (3)土壤中Ni的分布(圖2c)與其余各元素均不同,Ni僅在第四系東北部附近出現(xiàn)一個高值區(qū),其余地區(qū)為低值分布區(qū)。

    這種重金屬分布特征可能與人類生產(chǎn)活動有關(guān)。自1984年養(yǎng)馬島被列為山東省重點旅游開發(fā)區(qū)以來,人類活動主要集中在西南區(qū)域,島上的旅游業(yè)、捕撈業(yè)、建筑業(yè)蓬勃發(fā)展,與之相關(guān)的建筑管線、養(yǎng)殖浮球、農(nóng)業(yè)藥劑、生活電池等工業(yè)制品逐漸被大規(guī)模使用,產(chǎn)品廢棄后分解出來的重金屬隨著雨水、土壤孔隙擴散,從而影響土壤中重金屬的背景含量。

    3.2? 相關(guān)性分析

    定性方面,從圖2養(yǎng)馬島土壤重金屬含量分布圖中可以大致看出各元素之間的相關(guān)性,Zn和Hg之間的相關(guān)性較高,Pb、Cu和Cr之間的相關(guān)性較高,Ni分布與其余各元素不同,與其余元素相關(guān)性最低。

    定量方面,本次采用SPSS軟件對養(yǎng)馬島土壤重金屬含量進行Pearson相關(guān)分析,進一步分析各元素的同源性[2223],分析結(jié)果見表4。

    從表4中可以看出Cr分別與Cu、As之間相關(guān)系數(shù)在0.01水平上為0.842、0.930,為極高度相關(guān);Cu與As、Zn的相關(guān)系數(shù)分別在0.01和0.05水平為0.884和0.818,為極高度相關(guān);Zn與Hg相關(guān)系數(shù)在0.01水平上為0.901,為極高度相關(guān)。

    上述分析表明,Cu和As以及Zn和Hg均為親銅元素,在土壤中多具有較高的相關(guān)性,但Cr為親鐵元素和Cu也具有較高相關(guān)性,說明各元素之間除了受巖石風(fēng)化、侵蝕及水動力作用等自然過程作用影響,還受到人為的影響(例如人為排放污水,固體垃圾等),相關(guān)性越高的各元素同源性可能性也越大。Ni與其余各金屬元素之間相關(guān)性較低,同源可能性最低。

    3.3? 地質(zhì)累積指數(shù)評價

    養(yǎng)馬島土壤各重金屬的地質(zhì)累積指數(shù)(Igeo)計算結(jié)果見表5。從表中可以看出TR2點位Zn、As、Cd為偏中度污染,Hg為嚴重污染;TR3點位As、Cd、Hg為偏中度污染;TR4點位As、Pb為偏中度污染;TR5點位As為偏中度污染,Hg為中度污染;TR7點位Hg為偏中度污染;TR9點位Ni為偏中度污染;其余點位及元素為輕度或無污染。第四系分布區(qū)污染程度綜合來看,Hg均值為1.09,屬于偏中度污染,As和Zn均值分別為0.60、0.02,屬于輕度污染,其余元素為無污染,重金屬元素污染程度依次為Hg>As>Zn>Cd>Ni>Pb>Cu>Cr。

    3.4? 潛在生態(tài)風(fēng)險評價

    養(yǎng)馬島土壤中重金屬單元素及綜合潛在生態(tài)風(fēng)險指數(shù)的計算結(jié)果見表6,從計算結(jié)果來看,單元素潛在生態(tài)風(fēng)險指數(shù)(Eir),Cd在TR2、TR3、TR6和TR7點位風(fēng)險程度中等,Hg在TR2、TR3、TR4、TR5、TR7和TR9點位風(fēng)險程度依次為嚴重、較重、中等、重、較重和中等,其余點位元素風(fēng)險程度低;受Cd和Hg潛在生態(tài)風(fēng)險指數(shù)高值影響,TR2點位綜合潛在生態(tài)風(fēng)險指數(shù)(RI)最高,風(fēng)險程度極重;TR3、TR4、TR5和TR7點位綜合潛在生態(tài)風(fēng)險指數(shù)(RI)較小,風(fēng)險程度中等;其余點位綜合潛在生態(tài)風(fēng)險指數(shù)(RI)小,風(fēng)險程度低。

    對各元素風(fēng)險程度大小進行比較,依次為Hg>Cd>As>Ni>Pb>Cu>Cr>Zn,這與上述用地質(zhì)累積指數(shù)法得出的排列順序有所不同,主要原因是潛在生態(tài)風(fēng)險計算考慮了各重金屬的毒性響應(yīng)系數(shù)(Tir),結(jié)果更能反映重金屬元素的潛在生態(tài)危害。

    對各點位重金屬風(fēng)險程度大小進行比較,依次為TR2>TR5>TR3>TR7>TR4>TR9>TR6>TR1>TR10>TR8,綜合潛在生態(tài)風(fēng)險指數(shù)除TR2點位偏高外,其余點位均小于300,綜合來看,養(yǎng)馬島土壤重金屬潛在生態(tài)風(fēng)險程度整體屬于中等程度—低程度。

    4? 結(jié)論

    本次通過對養(yǎng)馬島第四系土壤的重金屬含量分析研究,得出如下結(jié)論:

    (1)Cr、Ni、Cu、Zn、As、Cd、Pb和Hg的含量平均值分別為23.54mg/kg、20.18mg/kg、16.28mg/kg、59.08mg/kg、8.14mg/kg、0.12mg/kg、32.96mg/kg和0.15mg/kg。除Cr和Cu外,其余元素含量受一定程度的人類活動影響,均超過了該元素牟平地區(qū)背景值。

    (2)空間分布上,Ni元素與其余元素不同,高值區(qū)分布在第四系的東北部,其余元素高值區(qū)均分布在第四系的中—西南部;相關(guān)性方面,Ni與其余元素相關(guān)性均較低,Zn和Hg以及Zn和Cu之間相關(guān)性較高,Cu、Cr、As互相之間相關(guān)性較高。

    (3)TR2點位Hg異常高值,Hg元素為偏中度污染,元素風(fēng)險程度重。地質(zhì)累積指數(shù)法評價各元素污染程度為Hg>As>Zn>Cd>Ni>Pb>Cu>Cr,潛在生態(tài)風(fēng)險法評價各元素風(fēng)險程度為Hg>Cd>As>Ni>Pb>Cu>Cr>Zn,綜合潛在生態(tài)風(fēng)險程度除TR2點位為極重外,其余點位為中等程度—低程度。

    參考文獻:

    [1]? 趙鳴,石洪華,池源,等.廟島群島南五島土壤性質(zhì)空間分布及其影響因子[J].海洋環(huán)境科學(xué),2017,36(6):884891.

    [2]? 王利英,程雄,黃文鑫,等.廣東內(nèi)伶仃島土壤調(diào)查與污染評價[J].廣東農(nóng)業(yè)科學(xué),2017,44(3):101107.

    [3]? 魏艷艷,張凱迪,徐良,等.浙江省無居民海島土壤重金屬污染與植物多樣性的關(guān)系[J].土壤通報,2019,50(5):12261232.

    [4]? 許昊,郭振,張志衛(wèi),等.有居民海島土壤重金屬空間異質(zhì)性及其入海風(fēng)險分析:以田橫島為例[J].海洋環(huán)境科學(xué),2023,42(5):704712.

    [5]? 劉華峰.章丘區(qū)北部土壤重金屬地球化學(xué)特征及潛在生態(tài)風(fēng)險[J].山東國土資源,2020,36(9):5057.

    [6]? 湯世凱,張杰,于曉靜,等.山東丁字灣海域沉積物重金屬含量、分布及與粒徑之間的關(guān)系研究[J].現(xiàn)代地質(zhì),2020,34(5):928935.

    [7]? 王美娥,彭馳,陳衛(wèi)平.寧夏干旱地區(qū)工業(yè)區(qū)對農(nóng)田土壤重金屬累積的影響[J].環(huán)境科學(xué),2016,37(9):35323539.

    [8]? 李金哲,劉寧強,龔慶杰,等.廣東汕頭市內(nèi)海灣沉積物重金屬環(huán)境質(zhì)量調(diào)查與評價[J].現(xiàn)代地質(zhì),2021,35(5):14411449.

    [9]? 林龍飛,彭濛萌.論海島型旅游非優(yōu)區(qū)資源評價:以煙臺養(yǎng)馬島和芝罘島為例[J].山東工商學(xué)院學(xué)報,2014,28(6):2329.

    [10]? 劉真真,張廣海.基于博弈論的中國有居民海島旅游開發(fā)決策研究:以山東煙臺養(yǎng)馬島為例[J].海洋開發(fā)與管理,2016,33(7):38.

    [11]? 山東省海洋與漁業(yè)廳.山東海島[M].北京:海洋出版社,2010:224228.

    [12]? MULLER G.Index of geoaccumulation in sediments of the Rhine River[J].Geojournal,1969,2(2):108118.

    [13]? 顧效源,孔祥淮,王偉,等.山東丁字灣表層沉積物重金屬分布及污染評價[J].海洋地質(zhì)前沿,2019,35(3):1321.

    [14]? 康桂玲,徐佳,平艷麗,等.魯西萊陽茌梨產(chǎn)地小流域沉積物重金屬生態(tài)風(fēng)險評價及來源分析[J].山東國土資源,2023,39(1):3039.

    [15]? 代杰瑞,龐緒貴.山東省縣(區(qū))級土壤地球化學(xué)基準值與背景值[M].北京:海洋出版社,2019:396397.

    [16]? HAKANSON L.An ecological risk index for aquatic pollution control. A sedimentological approach[J].Water research,1980,14(8):9751001.

    [17]? 胡慶海,王學(xué)求,韓志軒,等.京津冀地區(qū)永清縣土壤重金屬地球化學(xué)特征及綠色食品產(chǎn)地的土壤質(zhì)量評價[J].現(xiàn)代地質(zhì),2023,37(3):778789.

    [18]? 張朋朋,胡蕾,王小丹,等.長島北四島海洋牧場海域環(huán)境適宜性分析[J].山東國土資源,2022,38(5):5462.

    [19]? 李蘋,黃勇,林赟,等.北京市懷柔區(qū)土壤重金屬的分布特征、來源分析及風(fēng)險評價[J].現(xiàn)代地質(zhì),2018,32(1):8694.

    [20]? 丁婷婷,李強,杜士林,等.沙潁河流域水環(huán)境重金屬污染特征及生態(tài)風(fēng)險評價[J].環(huán)境化學(xué),2019,38(10):23862401.

    [21]? 姚春卉,張春榮,李少勇,等.膠州灣沿岸土壤重金屬元素分布特征及其生態(tài)風(fēng)險評價[J].中國科技論文,2021,16(1):112120.

    [22]? 楊振宇,廖超林,鄒炎,等.湘東北典型河源區(qū)土壤重金屬分布特征、來源解析及潛在生態(tài)風(fēng)險評價[J].環(huán)境化學(xué),2023,44(9):52885298.

    [23]? 張又文,韓建華,涂棋,等.天津市郊農(nóng)田土壤重金屬積累特征及評價[J].生態(tài)與農(nóng)村環(huán)境學(xué)報,2019,35(11):14451452.

    Evaluation of Heavy Metal Contamination and Ecological Risk of Soil in Southern Yangmadao in Yantai City in Shandong Province

    TANG Shikai1,2,HU Shengtao1,2,LI Yihong1,2,QIANG Menglin1,2,YU Linhong1,2,ZHANG Xiangheng3,WANG Can4

    (1.Shandong Marine Geological Exploration Institute, Shandong Yantai 264004,China;2.No.3 Exploration Institute of Geology and Mineral Resources, Shandong Yantai 264004,China;3.Yantai Ludong Surveying and Mappinging Limited Corporation,Shandong Yantai 264004,China;4.Shandong Geophysical and Geochemical Exploration Institute, Shandong Ji'nan 250000, China)

    Abstracr:Based on content testing analysis of heavy metal(Cr, Ni, Cu, Zn, As, Cd, Pb and Hg) of soil in Quaternary distribution area in southern Yangmadao in Yantai city, by using statistical analysis, geological accumulation index method and potential ecological risk index method, the correlation, pollution level and ecological risk of elements have been studied and evaluated. It is showed that influenced by human activities to a certain extent, except for Ni, high value areas are mainly distributed in the middle to southwest of the Quaternary system. Except for Cr and Cu, the average element content exceeds the background value. Correlation analysis shows that there is a high positive correlation between Zn and Hg, Zn and Cu, as well as between Cu, Cr and As. Hg has an abnormally high value at TR2 point, indicating moderate pollution and a high level of risk. The pollution level of each element is evaluated as Hg>As>Zn>Cd>Ni>Pb>Cu>Cr, with a risk level of Hg>Cd>As>Ni>Pb>Cu>Cr>Zn. Based on comprehensive evaluation of the potential ecological risk level, TR2 points are classified as extremely severe, while the rest are classified as moderate to low.

    Key words:Heavy metals in soil; correlation analysis; ecological risk assessment; southern Yangmadao; Yantai city in Shandong province

    猜你喜歡
    土壤重金屬相關(guān)性分析
    有機肥對土壤重金屬生物有效性影響研究進展
    濱州市城區(qū)苔蘚植物主要重金屬含量的調(diào)查與分析
    人民幣匯率變動與中國入境旅游相關(guān)性分析(2002—2016)
    上市公司財務(wù)指標與股票價格的相關(guān)性實證分析
    時代金融(2016年29期)2016-12-05 16:15:54
    淘寶星店成長中的粉絲力量
    中國市場(2016年40期)2016-11-28 03:35:20
    中國城市化與經(jīng)濟發(fā)展水平關(guān)系研究
    商(2016年33期)2016-11-24 22:04:19
    我國物流企業(yè)規(guī)模與效益的相關(guān)性分析
    商(2016年22期)2016-07-08 21:59:09
    基于Kriging插值的礦區(qū)周邊土壤重金屬空間分布規(guī)律研究
    科技視界(2016年13期)2016-06-13 08:06:36
    貴州擺紀磷石膏尾礦庫周邊土壤重金屬污染評價
    某地區(qū)表層土壤重金屬分布特征研究
    丰满人妻熟妇乱又伦精品不卡| 精品一区二区三区视频在线观看免费 | 久久久久久亚洲精品国产蜜桃av| 老司机午夜福利在线观看视频 | 黄色毛片三级朝国网站| 一本综合久久免费| 亚洲av成人一区二区三| 蜜桃在线观看..| 妹子高潮喷水视频| 亚洲欧美日韩高清在线视频 | 无限看片的www在线观看| 在线看a的网站| 国产深夜福利视频在线观看| 黄网站色视频无遮挡免费观看| 亚洲欧洲日产国产| 操出白浆在线播放| 美女高潮到喷水免费观看| av不卡在线播放| 黄片大片在线免费观看| 91精品三级在线观看| 亚洲精品在线美女| 色综合欧美亚洲国产小说| 成人黄色视频免费在线看| 久久亚洲真实| 肉色欧美久久久久久久蜜桃| 黄色视频不卡| 亚洲久久久国产精品| 亚洲 国产 在线| 激情视频va一区二区三区| 一级黄色大片毛片| 啦啦啦在线免费观看视频4| 日韩有码中文字幕| 午夜视频精品福利| 亚洲avbb在线观看| 99香蕉大伊视频| 热99久久久久精品小说推荐| 99精品在免费线老司机午夜| 大香蕉久久成人网| 国产精品久久久人人做人人爽| 精品福利观看| 亚洲专区字幕在线| 欧美亚洲 丝袜 人妻 在线| 蜜桃在线观看..| 免费黄频网站在线观看国产| 黄色 视频免费看| 少妇猛男粗大的猛烈进出视频| 久久这里只有精品19| 黄色视频不卡| 欧美+亚洲+日韩+国产| 女性被躁到高潮视频| 久久精品国产99精品国产亚洲性色 | 日韩免费高清中文字幕av| 美女主播在线视频| av在线播放免费不卡| 国产成人精品久久二区二区免费| 亚洲成a人片在线一区二区| 亚洲国产精品一区二区三区在线| 99re在线观看精品视频| 在线观看舔阴道视频| 高清在线国产一区| 欧美在线一区亚洲| 亚洲性夜色夜夜综合| 婷婷成人精品国产| 久久久国产一区二区| 欧美精品人与动牲交sv欧美| 国产野战对白在线观看| 成年人午夜在线观看视频| 男女边摸边吃奶| 乱人伦中国视频| 脱女人内裤的视频| 桃花免费在线播放| 狠狠狠狠99中文字幕| videosex国产| 国产午夜精品久久久久久| 建设人人有责人人尽责人人享有的| 窝窝影院91人妻| 久久久精品区二区三区| 日本vs欧美在线观看视频| 视频区欧美日本亚洲| 久热爱精品视频在线9| 如日韩欧美国产精品一区二区三区| 一区福利在线观看| 亚洲中文日韩欧美视频| 后天国语完整版免费观看| 精品人妻熟女毛片av久久网站| 好男人电影高清在线观看| 国产亚洲午夜精品一区二区久久| 一个人免费在线观看的高清视频| 99在线人妻在线中文字幕 | 国产在视频线精品| 欧美日韩黄片免| 日韩 欧美 亚洲 中文字幕| 国产午夜精品久久久久久| 国产成人影院久久av| e午夜精品久久久久久久| 国产精品一区二区免费欧美| 在线观看免费高清a一片| 9热在线视频观看99| 每晚都被弄得嗷嗷叫到高潮| 国产精品亚洲av一区麻豆| 激情视频va一区二区三区| 男男h啪啪无遮挡| 日本vs欧美在线观看视频| 成人国产av品久久久| 在线观看免费高清a一片| 在线十欧美十亚洲十日本专区| 亚洲三区欧美一区| 黑人巨大精品欧美一区二区蜜桃| 国产日韩欧美在线精品| 国产精品久久久人人做人人爽| av有码第一页| 欧美精品一区二区大全| 后天国语完整版免费观看| 免费久久久久久久精品成人欧美视频| 国产精品一区二区在线不卡| 欧美日韩亚洲高清精品| 肉色欧美久久久久久久蜜桃| 精品第一国产精品| 老熟女久久久| 人成视频在线观看免费观看| 精品一区二区三区av网在线观看 | 男女下面插进去视频免费观看| 国产男女内射视频| 婷婷成人精品国产| 亚洲av第一区精品v没综合| av网站在线播放免费| 777久久人妻少妇嫩草av网站| 亚洲精品一二三| 成人手机av| 脱女人内裤的视频| 天天影视国产精品| 91国产中文字幕| 国产精品电影一区二区三区 | 国产精品自产拍在线观看55亚洲 | 欧美精品啪啪一区二区三区| 国产伦理片在线播放av一区| 亚洲人成伊人成综合网2020| 亚洲色图av天堂| 免费在线观看完整版高清| 亚洲天堂av无毛| 午夜福利视频精品| 99re在线观看精品视频| 69精品国产乱码久久久| 曰老女人黄片| 国产欧美日韩精品亚洲av| 99久久国产精品久久久| 下体分泌物呈黄色| 欧美激情极品国产一区二区三区| 热99re8久久精品国产| 9色porny在线观看| 国产精品麻豆人妻色哟哟久久| 最近最新免费中文字幕在线| 久久 成人 亚洲| 国产亚洲精品第一综合不卡| 久久久久国内视频| 在线观看免费高清a一片| 国产精品麻豆人妻色哟哟久久| 久久免费观看电影| 一区二区三区国产精品乱码| 午夜福利乱码中文字幕| 国产色视频综合| 国产av国产精品国产| av天堂在线播放| 亚洲国产精品一区二区三区在线| 窝窝影院91人妻| 国产在线观看jvid| 狠狠狠狠99中文字幕| 欧美大码av| 丝瓜视频免费看黄片| 男女之事视频高清在线观看| 天堂8中文在线网| 亚洲少妇的诱惑av| 欧美亚洲 丝袜 人妻 在线| 99国产精品99久久久久| 久热这里只有精品99| 日本av免费视频播放| 男男h啪啪无遮挡| 成年动漫av网址| 久久九九热精品免费| 亚洲成人手机| 亚洲,欧美精品.| 99久久国产精品久久久| 国产精品亚洲一级av第二区| 国产精品av久久久久免费| 高清在线国产一区| 波多野结衣av一区二区av| 亚洲精华国产精华精| 高潮久久久久久久久久久不卡| 精品亚洲乱码少妇综合久久| 国产色视频综合| 久久精品亚洲精品国产色婷小说| 国产一区二区三区综合在线观看| 欧美国产精品一级二级三级| 亚洲精品在线观看二区| 十八禁人妻一区二区| 国产免费视频播放在线视频| 少妇被粗大的猛进出69影院| 国产成人免费观看mmmm| 欧美日韩黄片免| 免费女性裸体啪啪无遮挡网站| 正在播放国产对白刺激| 国产日韩欧美亚洲二区| 欧美成狂野欧美在线观看| 免费在线观看日本一区| 老汉色∧v一级毛片| a级片在线免费高清观看视频| 青草久久国产| 高潮久久久久久久久久久不卡| 欧美+亚洲+日韩+国产| 天天添夜夜摸| 老司机午夜福利在线观看视频 | 91麻豆精品激情在线观看国产 | av欧美777| 日本a在线网址| √禁漫天堂资源中文www| 午夜两性在线视频| 久久久水蜜桃国产精品网| 国产色视频综合| 午夜日韩欧美国产| 亚洲国产毛片av蜜桃av| 女性生殖器流出的白浆| 精品一区二区三区四区五区乱码| 色播在线永久视频| 精品亚洲乱码少妇综合久久| 国产精品自产拍在线观看55亚洲 | 亚洲色图 男人天堂 中文字幕| 欧美大码av| 在线观看66精品国产| 80岁老熟妇乱子伦牲交| 99riav亚洲国产免费| 国产人伦9x9x在线观看| 下体分泌物呈黄色| 精品一品国产午夜福利视频| 高清黄色对白视频在线免费看| 欧美一级毛片孕妇| 一区二区三区国产精品乱码| 成人国产一区最新在线观看| 在线天堂中文资源库| 国产av又大| svipshipincom国产片| a级毛片黄视频| 国产又色又爽无遮挡免费看| 人成视频在线观看免费观看| 亚洲第一欧美日韩一区二区三区 | 国产亚洲精品第一综合不卡| 久久中文字幕一级| 最新的欧美精品一区二区| 91大片在线观看| 大片免费播放器 马上看| 一本久久中文字幕| 亚洲在线观看片| 国产爱豆传媒在线观看| 九色成人免费人妻av| 99riav亚洲国产免费| 亚洲国产精品久久男人天堂| 又紧又爽又黄一区二区| 色噜噜av男人的天堂激情| 日本一二三区视频观看| 一卡2卡三卡四卡精品乱码亚洲| 91在线观看av| 久久国产乱子伦精品免费另类| 欧美成人性av电影在线观看| 亚洲男人的天堂狠狠| 国产精品久久久久久久电影 | 色综合亚洲欧美另类图片| 午夜福利视频1000在线观看| 午夜成年电影在线免费观看| 久久精品亚洲精品国产色婷小说| 久久人人精品亚洲av| 搡老妇女老女人老熟妇| 亚洲人成伊人成综合网2020| 一区二区三区国产精品乱码| 日本 av在线| 亚洲五月天丁香| 国产综合懂色| 一进一出抽搐gif免费好疼| 中文字幕熟女人妻在线| 国产真实乱freesex| 国产毛片a区久久久久| 99久久精品热视频| 亚洲五月天丁香| 人妻丰满熟妇av一区二区三区| 97碰自拍视频| 久久久久久大精品| 51午夜福利影视在线观看| 最近在线观看免费完整版| 欧美日韩瑟瑟在线播放| 国产主播在线观看一区二区| 日韩高清综合在线| 久久人妻av系列| 欧美+亚洲+日韩+国产| 免费搜索国产男女视频| 国产精品一区二区免费欧美| 日韩欧美 国产精品| 午夜福利免费观看在线| 久久久久国内视频| 亚洲av中文字字幕乱码综合| 日本与韩国留学比较| 男插女下体视频免费在线播放| 亚洲第一电影网av| 欧美日韩黄片免| 国产三级黄色录像| 久9热在线精品视频| 麻豆一二三区av精品| 最好的美女福利视频网| а√天堂www在线а√下载| 中出人妻视频一区二区| 岛国在线观看网站| 亚洲激情在线av| 国产精品美女特级片免费视频播放器 | 三级男女做爰猛烈吃奶摸视频| 中文字幕人妻丝袜一区二区| 美女 人体艺术 gogo| 国产午夜精品论理片| 亚洲激情在线av| 国产精品乱码一区二三区的特点| 日日摸夜夜添夜夜添小说| 国产精品av视频在线免费观看| 成人无遮挡网站| 国产aⅴ精品一区二区三区波| 午夜福利在线观看免费完整高清在 | 在线观看一区二区三区| 久久久精品欧美日韩精品| 99热6这里只有精品| 国产av麻豆久久久久久久| 12—13女人毛片做爰片一| 舔av片在线| www.精华液| 女同久久另类99精品国产91| 黄片小视频在线播放| 热99re8久久精品国产| 成人高潮视频无遮挡免费网站| 国产精品av视频在线免费观看| tocl精华| 熟妇人妻久久中文字幕3abv| 他把我摸到了高潮在线观看| 久9热在线精品视频| 亚洲自偷自拍图片 自拍| 国产欧美日韩一区二区三| 成人性生交大片免费视频hd| 黑人操中国人逼视频| 国产蜜桃级精品一区二区三区| 可以在线观看的亚洲视频| 日韩国内少妇激情av| 男人舔女人下体高潮全视频| 免费看日本二区| 99国产极品粉嫩在线观看| 日韩高清综合在线| 真人做人爱边吃奶动态| 欧美成狂野欧美在线观看| 精品久久久久久久久久久久久| 在线a可以看的网站| 久久精品影院6| 首页视频小说图片口味搜索| 国产精品1区2区在线观看.| 色吧在线观看| 桃红色精品国产亚洲av| 久久精品国产综合久久久| 国产蜜桃级精品一区二区三区| 欧美日韩福利视频一区二区| 日韩国内少妇激情av| 中文字幕人妻丝袜一区二区| 国产单亲对白刺激| 国产午夜福利久久久久久| 黄色视频,在线免费观看| 天堂影院成人在线观看| 日韩免费av在线播放| 夜夜爽天天搞| 国产探花在线观看一区二区| 少妇丰满av| 在线免费观看不下载黄p国产 | 一区福利在线观看| 免费在线观看成人毛片| 在线视频色国产色| 亚洲国产中文字幕在线视频| 午夜免费观看网址| 人妻丰满熟妇av一区二区三区| 亚洲熟妇中文字幕五十中出| 国产1区2区3区精品| 视频区欧美日本亚洲| 午夜日韩欧美国产| 国产黄色小视频在线观看| 中文资源天堂在线| 国产黄色小视频在线观看| 亚洲va日本ⅴa欧美va伊人久久| 精品免费久久久久久久清纯| 国产精品99久久99久久久不卡| 欧美xxxx黑人xx丫x性爽| 黄色成人免费大全| 午夜精品在线福利| 两个人的视频大全免费| 99精品久久久久人妻精品| 日本黄大片高清| 精品欧美国产一区二区三| 成人av一区二区三区在线看| 狂野欧美白嫩少妇大欣赏| 国产乱人视频| aaaaa片日本免费| 亚洲成人精品中文字幕电影| 特大巨黑吊av在线直播| 免费av毛片视频| 蜜桃久久精品国产亚洲av| 色综合亚洲欧美另类图片| 久久热在线av| 琪琪午夜伦伦电影理论片6080| 午夜日韩欧美国产| 91在线精品国自产拍蜜月 | 国产蜜桃级精品一区二区三区| 午夜久久久久精精品| 欧美黄色片欧美黄色片| 一级毛片高清免费大全| 国产伦人伦偷精品视频| 欧美日韩精品网址| 高清毛片免费观看视频网站| 1000部很黄的大片| 欧美色欧美亚洲另类二区| 国产成人系列免费观看| 日韩欧美在线二视频| 宅男免费午夜| 免费一级毛片在线播放高清视频| 日韩欧美免费精品| 搡老妇女老女人老熟妇| 老司机在亚洲福利影院| 高清毛片免费观看视频网站| 午夜两性在线视频| 亚洲国产色片| 婷婷精品国产亚洲av在线| 亚洲成a人片在线一区二区| 美女大奶头视频| 国产69精品久久久久777片 | www.熟女人妻精品国产| 日本精品一区二区三区蜜桃| 亚洲av片天天在线观看| 搡老妇女老女人老熟妇| 亚洲国产精品合色在线| 香蕉丝袜av| 亚洲国产精品合色在线| 午夜福利视频1000在线观看| 欧美+亚洲+日韩+国产| 女人被狂操c到高潮| 欧美在线黄色| 亚洲天堂国产精品一区在线| 成年女人永久免费观看视频| 国产精品乱码一区二三区的特点| 99riav亚洲国产免费| 两个人的视频大全免费| 两个人视频免费观看高清| 色老头精品视频在线观看| a级毛片在线看网站| 亚洲熟女毛片儿| 最近在线观看免费完整版| 神马国产精品三级电影在线观看| 精品日产1卡2卡| 九九久久精品国产亚洲av麻豆 | 国产精品久久久人人做人人爽| 淫妇啪啪啪对白视频| 午夜福利免费观看在线| 成人国产综合亚洲| 一本一本综合久久| 欧美日本亚洲视频在线播放| 国产精品 欧美亚洲| 岛国视频午夜一区免费看| 小蜜桃在线观看免费完整版高清| 亚洲欧美日韩卡通动漫| 国产午夜福利久久久久久| 中文字幕精品亚洲无线码一区| 亚洲美女黄片视频| 日韩欧美国产在线观看| 欧美性猛交黑人性爽| 国产精品99久久久久久久久| 精品久久久久久久末码| 日韩免费av在线播放| www国产在线视频色| 丰满的人妻完整版| 一级毛片精品| 亚洲av熟女| 中文亚洲av片在线观看爽| 一级作爱视频免费观看| 亚洲自拍偷在线| 亚洲男人的天堂狠狠| 天堂动漫精品| 欧美乱色亚洲激情| 男插女下体视频免费在线播放| 成人国产一区最新在线观看| 校园春色视频在线观看| 99在线人妻在线中文字幕| 熟妇人妻久久中文字幕3abv| 男插女下体视频免费在线播放| 男人舔奶头视频| 美女扒开内裤让男人捅视频| 亚洲18禁久久av| 国产伦在线观看视频一区| 免费看十八禁软件| 国产精品影院久久| 美女扒开内裤让男人捅视频| 床上黄色一级片| 不卡av一区二区三区| 亚洲精品久久国产高清桃花| 亚洲专区中文字幕在线| 亚洲国产中文字幕在线视频| 久久香蕉国产精品| 手机成人av网站| 国产淫片久久久久久久久 | 不卡av一区二区三区| 桃红色精品国产亚洲av| 午夜影院日韩av| 日韩欧美国产在线观看| 精品久久久久久久毛片微露脸| 久久久久国产精品人妻aⅴ院| 欧美黑人巨大hd| 日韩免费av在线播放| 欧美极品一区二区三区四区| 国产黄a三级三级三级人| netflix在线观看网站| 亚洲黑人精品在线| 18禁观看日本| 亚洲精品乱码久久久v下载方式 | 免费在线观看亚洲国产| 国产高清激情床上av| 免费一级毛片在线播放高清视频| 免费人成视频x8x8入口观看| 欧美一区二区精品小视频在线| 特大巨黑吊av在线直播| 变态另类丝袜制服| 老熟妇乱子伦视频在线观看| 国产成人aa在线观看| www日本黄色视频网| 免费观看精品视频网站| 18禁美女被吸乳视频| 亚洲电影在线观看av| 欧美日韩一级在线毛片| 亚洲av成人不卡在线观看播放网| 免费在线观看成人毛片| 日本与韩国留学比较| 亚洲人成伊人成综合网2020| 伊人久久大香线蕉亚洲五| 免费观看的影片在线观看| 全区人妻精品视频| 毛片女人毛片| 一进一出抽搐gif免费好疼| 三级毛片av免费| 国产野战对白在线观看| 99久国产av精品| 久久这里只有精品中国| 制服人妻中文乱码| 国产欧美日韩精品亚洲av| 国产激情偷乱视频一区二区| 舔av片在线| 天天一区二区日本电影三级| 亚洲黑人精品在线| 久久香蕉精品热| 毛片女人毛片| 精品电影一区二区在线| 久久这里只有精品19| 久久国产乱子伦精品免费另类| 亚洲欧美精品综合久久99| 黄色女人牲交| 久久久精品大字幕| 亚洲精品国产精品久久久不卡| netflix在线观看网站| 欧美日韩乱码在线| 曰老女人黄片| 国产亚洲精品久久久com| 成人国产一区最新在线观看| 国产成人欧美在线观看| 亚洲精品色激情综合| 日本熟妇午夜| 最近在线观看免费完整版| 午夜a级毛片| 亚洲av片天天在线观看| 禁无遮挡网站| 99re在线观看精品视频| 精品久久久久久久人妻蜜臀av| 18美女黄网站色大片免费观看| 国产男靠女视频免费网站| 精品久久久久久久久久免费视频| 亚洲欧美日韩卡通动漫| 国产三级中文精品| 久久国产精品人妻蜜桃| 成年女人永久免费观看视频| 婷婷精品国产亚洲av| 精品无人区乱码1区二区| 国产欧美日韩一区二区三| 国产伦精品一区二区三区四那| 欧美xxxx黑人xx丫x性爽| 日本与韩国留学比较| 怎么达到女性高潮| 亚洲五月天丁香| 色在线成人网| 日日摸夜夜添夜夜添小说| 一边摸一边抽搐一进一小说| 国产三级在线视频| 99热精品在线国产| 亚洲av熟女| 国产乱人视频| 日韩欧美国产在线观看| 亚洲欧美一区二区三区黑人| 91字幕亚洲| 午夜久久久久精精品| 亚洲七黄色美女视频| ponron亚洲| 亚洲第一电影网av| 国产黄色小视频在线观看| 大型黄色视频在线免费观看| 观看免费一级毛片| 国产成人精品久久二区二区91| 91字幕亚洲| 制服人妻中文乱码| 亚洲欧美日韩卡通动漫| av天堂在线播放| 免费观看精品视频网站| 国模一区二区三区四区视频 | www.熟女人妻精品国产| 精品免费久久久久久久清纯| 男人舔女人的私密视频|