• <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ū)表層土壤重金屬分布特征研究
    很黄的视频免费| 亚洲五月婷婷丁香| 熟妇人妻久久中文字幕3abv| 可以在线观看的亚洲视频| 香蕉丝袜av| 久久久水蜜桃国产精品网| 亚洲中文字幕日韩| 亚洲 欧美一区二区三区| 精品久久久久久久久久久久久| 日韩有码中文字幕| 69av精品久久久久久| 亚洲美女视频黄频| 精品国产亚洲在线| 老司机午夜福利在线观看视频| 亚洲 欧美 日韩 在线 免费| 91老司机精品| 在线播放国产精品三级| 国产成人aa在线观看| 欧美在线一区亚洲| 国产精品日韩av在线免费观看| 天堂√8在线中文| 妹子高潮喷水视频| 国产伦在线观看视频一区| 制服人妻中文乱码| 久久久久国内视频| 成在线人永久免费视频| 男女做爰动态图高潮gif福利片| 国内少妇人妻偷人精品xxx网站 | 欧美色视频一区免费| av在线天堂中文字幕| 两个人免费观看高清视频| 一进一出抽搐gif免费好疼| 天天一区二区日本电影三级| 91大片在线观看| 99热6这里只有精品| 日韩大码丰满熟妇| 国产一区二区激情短视频| 一a级毛片在线观看| 91字幕亚洲| 变态另类成人亚洲欧美熟女| netflix在线观看网站| 丰满人妻熟妇乱又伦精品不卡| 久久香蕉激情| 一本久久中文字幕| 亚洲五月婷婷丁香| 久久香蕉国产精品| 一进一出好大好爽视频| 亚洲自拍偷在线| 又紧又爽又黄一区二区| 国产午夜精品论理片| 国产探花在线观看一区二区| 午夜福利成人在线免费观看| 亚洲成人久久爱视频| 久久性视频一级片| cao死你这个sao货| 叶爱在线成人免费视频播放| 亚洲av美国av| 99re在线观看精品视频| 午夜福利18| 日韩精品青青久久久久久| 久久久久久人人人人人| 精品一区二区三区视频在线观看免费| 日本成人三级电影网站| 又黄又爽又免费观看的视频| 在线十欧美十亚洲十日本专区| 成人国产综合亚洲| 无限看片的www在线观看| 亚洲午夜精品一区,二区,三区| 国产成人av激情在线播放| 国产人伦9x9x在线观看| 丰满人妻熟妇乱又伦精品不卡| 日本在线视频免费播放| 伦理电影免费视频| 91老司机精品| 国产激情偷乱视频一区二区| 精品久久久久久,| 日韩成人在线观看一区二区三区| 午夜久久久久精精品| 欧美国产日韩亚洲一区| 国产av一区在线观看免费| 人妻久久中文字幕网| 精品久久蜜臀av无| aaaaa片日本免费| 91国产中文字幕| 日韩大码丰满熟妇| 桃红色精品国产亚洲av| 最近最新免费中文字幕在线| 一级片免费观看大全| 久久久国产精品麻豆| 国产精品爽爽va在线观看网站| 床上黄色一级片| 岛国视频午夜一区免费看| 亚洲熟妇熟女久久| 久久久国产精品麻豆| 国产精品日韩av在线免费观看| 亚洲国产欧洲综合997久久,| 亚洲av美国av| 久久久久久久久久黄片| 久久久久久大精品| 亚洲av电影在线进入| 丝袜美腿诱惑在线| 亚洲中文av在线| 午夜激情福利司机影院| 国产激情久久老熟女| 国产精品一区二区免费欧美| 午夜久久久久精精品| 最近最新中文字幕大全电影3| 久久伊人香网站| 欧美极品一区二区三区四区| 可以免费在线观看a视频的电影网站| 亚洲欧美精品综合久久99| 成人三级做爰电影| 狂野欧美激情性xxxx| 亚洲第一欧美日韩一区二区三区| 久9热在线精品视频| 一进一出抽搐动态| 欧美一级毛片孕妇| 国产精品永久免费网站| 国产精品 欧美亚洲| 国产精品久久电影中文字幕| 免费看日本二区| 久久国产精品人妻蜜桃| 精品久久久久久久久久免费视频| 91老司机精品| netflix在线观看网站| 欧美乱色亚洲激情| 欧美不卡视频在线免费观看 | 男女下面进入的视频免费午夜| 亚洲成人国产一区在线观看| 国产精品久久久久久精品电影| 在线永久观看黄色视频| 国产av一区二区精品久久| 久久国产乱子伦精品免费另类| 午夜免费成人在线视频| 国产精品影院久久| 嫩草影视91久久| 又粗又爽又猛毛片免费看| 国内少妇人妻偷人精品xxx网站 | 三级男女做爰猛烈吃奶摸视频| 国产精品99久久99久久久不卡| 国产精品一区二区精品视频观看| 在线观看免费视频日本深夜| 国产一级毛片七仙女欲春2| 亚洲欧美精品综合一区二区三区| 国产精品久久久久久久电影 | 真人一进一出gif抽搐免费| bbb黄色大片| 精品无人区乱码1区二区| 又大又爽又粗| 亚洲一码二码三码区别大吗| 欧美性猛交黑人性爽| 亚洲18禁久久av| 亚洲精品粉嫩美女一区| 久久久精品欧美日韩精品| 国产伦人伦偷精品视频| 精品免费久久久久久久清纯| 亚洲人成77777在线视频| 这个男人来自地球电影免费观看| 99久久国产精品久久久| 精品久久久久久久久久久久久| 久久久久久亚洲精品国产蜜桃av| 少妇的丰满在线观看| 法律面前人人平等表现在哪些方面| 人妻久久中文字幕网| 欧美另类亚洲清纯唯美| 国产人伦9x9x在线观看| 九色国产91popny在线| 精品久久久久久久久久免费视频| 99热这里只有是精品50| 精品国内亚洲2022精品成人| 可以在线观看毛片的网站| 国产av麻豆久久久久久久| 成人精品一区二区免费| 欧美日本视频| 熟妇人妻久久中文字幕3abv| 两性午夜刺激爽爽歪歪视频在线观看 | 久久精品国产亚洲av高清一级| 日本一二三区视频观看| 一个人观看的视频www高清免费观看 | 欧美绝顶高潮抽搐喷水| 欧美成人午夜精品| av中文乱码字幕在线| 久久人人精品亚洲av| 国产高清激情床上av| 亚洲熟妇中文字幕五十中出| 黄色视频,在线免费观看| a级毛片在线看网站| 在线观看免费午夜福利视频| 国产成人av教育| 婷婷丁香在线五月| 久久人妻av系列| 一本大道久久a久久精品| 白带黄色成豆腐渣| 精品一区二区三区视频在线观看免费| 久久久久久九九精品二区国产 | 亚洲aⅴ乱码一区二区在线播放 | 国产单亲对白刺激| 国产欧美日韩一区二区精品| 亚洲人与动物交配视频| 午夜福利成人在线免费观看| 欧美一级a爱片免费观看看 | 天堂√8在线中文| 日韩av在线大香蕉| 99精品久久久久人妻精品| 美女扒开内裤让男人捅视频| xxx96com| 成人高潮视频无遮挡免费网站| 日本黄色视频三级网站网址| 一边摸一边做爽爽视频免费| 久久精品人妻少妇| 色在线成人网| 久久精品国产清高在天天线| 俄罗斯特黄特色一大片| 丝袜人妻中文字幕| 精品无人区乱码1区二区| 国产激情久久老熟女| 亚洲一区二区三区不卡视频| 黑人巨大精品欧美一区二区mp4| 欧美+亚洲+日韩+国产| 午夜福利在线在线| 在线观看一区二区三区| 欧美黑人欧美精品刺激| 亚洲欧美一区二区三区黑人| 国产精华一区二区三区| 欧美中文日本在线观看视频| 国产片内射在线| 日韩欧美在线乱码| 欧美精品亚洲一区二区| 久久久久久九九精品二区国产 | 国产精华一区二区三区| 亚洲国产中文字幕在线视频| 国产成人av激情在线播放| 又黄又爽又免费观看的视频| 亚洲成人久久性| av欧美777| 十八禁网站免费在线| 久久久国产精品麻豆| 免费搜索国产男女视频| 99热6这里只有精品| 午夜精品一区二区三区免费看| 香蕉久久夜色| 国产99白浆流出| 国产蜜桃级精品一区二区三区| 久久精品综合一区二区三区| 久久中文字幕一级| 后天国语完整版免费观看| 国产单亲对白刺激| 精品久久久久久久久久免费视频| 色老头精品视频在线观看| 午夜免费成人在线视频| 18禁国产床啪视频网站| 亚洲欧美一区二区三区黑人| 成年人黄色毛片网站| 中出人妻视频一区二区| 久久草成人影院| 一级毛片女人18水好多| 免费在线观看视频国产中文字幕亚洲| 日韩欧美 国产精品| 伊人久久大香线蕉亚洲五| 亚洲精品久久成人aⅴ小说| 亚洲国产欧美人成| 中亚洲国语对白在线视频| 最新在线观看一区二区三区| 久久 成人 亚洲| 久久婷婷成人综合色麻豆| 午夜福利18| 亚洲,欧美精品.| 99国产综合亚洲精品| 久久精品成人免费网站| e午夜精品久久久久久久| 中文在线观看免费www的网站 | www日本黄色视频网| 18禁国产床啪视频网站| 亚洲av第一区精品v没综合| 别揉我奶头~嗯~啊~动态视频| 日本三级黄在线观看| 9191精品国产免费久久| 国产精品亚洲一级av第二区| 久久精品人妻少妇| 日本精品一区二区三区蜜桃| 日韩欧美免费精品| 一个人观看的视频www高清免费观看 | 久久久久久久久中文| 欧美日韩国产亚洲二区| 女人爽到高潮嗷嗷叫在线视频| 亚洲专区字幕在线| 欧美精品啪啪一区二区三区| 一二三四社区在线视频社区8| 18禁美女被吸乳视频| 久久午夜亚洲精品久久| www.精华液| 在线观看舔阴道视频| 国产午夜福利久久久久久| 婷婷精品国产亚洲av在线| 欧美黑人欧美精品刺激| cao死你这个sao货| 亚洲第一电影网av| 亚洲熟女毛片儿| 国产av一区二区精品久久| 久久久久久久精品吃奶| bbb黄色大片| 亚洲精品av麻豆狂野| 身体一侧抽搐| 亚洲国产精品成人综合色| 国产av又大| 久久热在线av| 欧美不卡视频在线免费观看 | 性色av乱码一区二区三区2| 热99re8久久精品国产| 91老司机精品| 在线a可以看的网站| 午夜福利18| 亚洲熟妇熟女久久| 毛片女人毛片| 级片在线观看| 亚洲一区中文字幕在线| 日韩欧美精品v在线| 久久精品综合一区二区三区| 亚洲精华国产精华精| 成熟少妇高潮喷水视频| www日本在线高清视频| 人人妻人人澡欧美一区二区| 老汉色av国产亚洲站长工具| 亚洲成av人片免费观看| 亚洲欧美精品综合久久99| 国产成人精品久久二区二区免费| 12—13女人毛片做爰片一| 亚洲成人中文字幕在线播放| 久久中文字幕一级| 国产日本99.免费观看| 亚洲欧美日韩高清在线视频| 搞女人的毛片| 婷婷精品国产亚洲av在线| 亚洲精品粉嫩美女一区| 麻豆国产av国片精品| 色综合婷婷激情| 欧美+亚洲+日韩+国产| 99精品欧美一区二区三区四区| 国产欧美日韩一区二区精品| 精品国内亚洲2022精品成人| 国产精品亚洲美女久久久| 午夜视频精品福利| 19禁男女啪啪无遮挡网站| 久久这里只有精品19| 成人国产综合亚洲| 国产精品久久久人人做人人爽| 久久精品国产亚洲av高清一级| 亚洲av电影不卡..在线观看| 国产乱人伦免费视频| 久久久久久久久久黄片| 亚洲色图 男人天堂 中文字幕| 级片在线观看| av免费在线观看网站| 国产91精品成人一区二区三区| 两性夫妻黄色片| 人人妻人人看人人澡| 国产av在哪里看| 亚洲七黄色美女视频| 99精品久久久久人妻精品| 国产1区2区3区精品| 亚洲 国产 在线| 国产精品免费一区二区三区在线| 欧美黑人巨大hd| 午夜成年电影在线免费观看| 成年版毛片免费区| 亚洲精品久久成人aⅴ小说| 大型av网站在线播放| 美女扒开内裤让男人捅视频| 成人国产综合亚洲| 精品国产美女av久久久久小说| 91国产中文字幕| 国产久久久一区二区三区| 欧美不卡视频在线免费观看 | 欧美成人一区二区免费高清观看 | 久久久国产成人精品二区| 欧美精品亚洲一区二区| 国产高清激情床上av| 成人手机av| 国产av麻豆久久久久久久| 免费无遮挡裸体视频| 亚洲国产精品999在线| 两个人视频免费观看高清| 国产高清videossex| 日韩欧美 国产精品| 国产精品爽爽va在线观看网站| 热99re8久久精品国产| 在线国产一区二区在线| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲国产精品成人综合色| 欧美日韩亚洲国产一区二区在线观看| 欧美激情久久久久久爽电影| 午夜福利欧美成人| 色播亚洲综合网| 色尼玛亚洲综合影院| 久久婷婷成人综合色麻豆| 男女做爰动态图高潮gif福利片| 久久精品91无色码中文字幕| 精品一区二区三区av网在线观看| 国产精品久久久av美女十八| 国产亚洲av嫩草精品影院| 九色国产91popny在线| 亚洲中文字幕日韩| 哪里可以看免费的av片| 两人在一起打扑克的视频| 又黄又爽又免费观看的视频| 在线免费观看的www视频| 狂野欧美白嫩少妇大欣赏| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜影院日韩av| 国产精品 国内视频| 一本精品99久久精品77| 妹子高潮喷水视频| 亚洲18禁久久av| 国产真实乱freesex| 好看av亚洲va欧美ⅴa在| 日韩欧美一区二区三区在线观看| 午夜福利视频1000在线观看| 亚洲成av人片在线播放无| 日韩精品免费视频一区二区三区| a在线观看视频网站| 9191精品国产免费久久| 亚洲国产精品sss在线观看| 一个人观看的视频www高清免费观看 | 香蕉丝袜av| 国产亚洲欧美在线一区二区| x7x7x7水蜜桃| 777久久人妻少妇嫩草av网站| 日本精品一区二区三区蜜桃| 97超级碰碰碰精品色视频在线观看| 久久久久国产精品人妻aⅴ院| 男女床上黄色一级片免费看| 老司机午夜福利在线观看视频| 我的老师免费观看完整版| 亚洲狠狠婷婷综合久久图片| www.精华液| 欧美激情久久久久久爽电影| 欧美+亚洲+日韩+国产| 免费看十八禁软件| 欧美中文综合在线视频| 国产精品一区二区精品视频观看| 天堂√8在线中文| 婷婷丁香在线五月| 日韩欧美 国产精品| 亚洲av第一区精品v没综合| 久久久精品大字幕| 午夜免费成人在线视频| 日韩高清综合在线| 看黄色毛片网站| 久久婷婷人人爽人人干人人爱| 久久午夜亚洲精品久久| 欧美 亚洲 国产 日韩一| 精品乱码久久久久久99久播| 欧美中文日本在线观看视频| 亚洲一区二区三区色噜噜| 国内毛片毛片毛片毛片毛片| 国产片内射在线| www.自偷自拍.com| 中亚洲国语对白在线视频| 久久 成人 亚洲| 国内揄拍国产精品人妻在线| 99精品久久久久人妻精品| www日本在线高清视频| 男女床上黄色一级片免费看| 国产区一区二久久| 国产精品久久久久久久电影 | 亚洲免费av在线视频| 久久国产精品影院| 亚洲无线在线观看| 高清在线国产一区| 日韩精品免费视频一区二区三区| 国产亚洲欧美98| 久久精品人妻少妇| 国产成人aa在线观看| 精品电影一区二区在线| 日韩有码中文字幕| 亚洲一区二区三区色噜噜| 国产亚洲av嫩草精品影院| 久久香蕉精品热| a在线观看视频网站| 亚洲人成网站在线播放欧美日韩| cao死你这个sao货| 欧美 亚洲 国产 日韩一| 日韩欧美免费精品| 91在线观看av| 国产91精品成人一区二区三区| 很黄的视频免费| 久久欧美精品欧美久久欧美| 99精品久久久久人妻精品| 成人av一区二区三区在线看| 欧美zozozo另类| 成人欧美大片| 亚洲美女视频黄频| 国产日本99.免费观看| 国产99白浆流出| av视频在线观看入口| 国产激情偷乱视频一区二区| 99久久精品热视频| 国产视频内射| 999久久久精品免费观看国产| 一个人免费在线观看电影 | 国产成人av激情在线播放| 岛国视频午夜一区免费看| 最新美女视频免费是黄的| 国产不卡一卡二| 老司机午夜福利在线观看视频| 亚洲自拍偷在线| 日本一区二区免费在线视频| 国产精品98久久久久久宅男小说| 一个人免费在线观看的高清视频| 亚洲欧洲精品一区二区精品久久久| 国产三级在线视频| 欧美人与性动交α欧美精品济南到| 熟女电影av网| 欧美成人性av电影在线观看| 久久久久免费精品人妻一区二区| 女生性感内裤真人,穿戴方法视频| 级片在线观看| 午夜福利欧美成人| 成年女人毛片免费观看观看9| 色在线成人网| 女警被强在线播放| 日本黄大片高清| 最近最新中文字幕大全免费视频| 人人妻人人看人人澡| 欧美日韩中文字幕国产精品一区二区三区| 怎么达到女性高潮| 黄色女人牲交| 久久草成人影院| 美女免费视频网站| 首页视频小说图片口味搜索| 少妇粗大呻吟视频| 一本一本综合久久| 欧美日韩黄片免| 亚洲精品久久国产高清桃花| 老司机在亚洲福利影院| 狠狠狠狠99中文字幕| 一本综合久久免费| 亚洲人成77777在线视频| 人妻丰满熟妇av一区二区三区| 精品国产美女av久久久久小说| 校园春色视频在线观看| 国产在线观看jvid| 国产亚洲欧美98| 男人舔女人的私密视频| 久久香蕉激情| 人妻夜夜爽99麻豆av| 欧美乱妇无乱码| 全区人妻精品视频| 麻豆av在线久日| 老司机在亚洲福利影院| 性色av乱码一区二区三区2| netflix在线观看网站| АⅤ资源中文在线天堂| 两个人看的免费小视频| 日韩欧美一区二区三区在线观看| 99久久综合精品五月天人人| 国产亚洲精品第一综合不卡| 最新在线观看一区二区三区| 日韩欧美国产一区二区入口| 亚洲人与动物交配视频| 999久久久国产精品视频| 麻豆成人午夜福利视频| 12—13女人毛片做爰片一| 女警被强在线播放| 天天躁夜夜躁狠狠躁躁| 欧美日韩福利视频一区二区| 午夜视频精品福利| 真人一进一出gif抽搐免费| 老熟妇仑乱视频hdxx| 在线观看日韩欧美| 人妻夜夜爽99麻豆av| 免费在线观看黄色视频的| 日韩欧美 国产精品| 香蕉国产在线看| 国产午夜福利久久久久久| 黄片大片在线免费观看| 成人特级黄色片久久久久久久| 给我免费播放毛片高清在线观看| 12—13女人毛片做爰片一| 亚洲欧美日韩高清专用| 免费看美女性在线毛片视频| 女人爽到高潮嗷嗷叫在线视频| 国产精品精品国产色婷婷| 久久久久久久久中文| 露出奶头的视频| 国产亚洲精品av在线| 国产一区二区三区在线臀色熟女| 看黄色毛片网站| 舔av片在线| 无限看片的www在线观看| 亚洲九九香蕉| 女人爽到高潮嗷嗷叫在线视频| 成在线人永久免费视频| 亚洲国产日韩欧美精品在线观看 | 欧洲精品卡2卡3卡4卡5卡区| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲一区二区三区不卡视频| 国产精品精品国产色婷婷| 91九色精品人成在线观看| 制服丝袜大香蕉在线| 熟妇人妻久久中文字幕3abv| av视频在线观看入口| 熟女电影av网| 操出白浆在线播放| av天堂在线播放| 少妇的丰满在线观看| 搡老岳熟女国产| √禁漫天堂资源中文www| 国产精品久久久av美女十八| 欧美人与性动交α欧美精品济南到| 999精品在线视频| 国产麻豆成人av免费视频| 亚洲电影在线观看av| 久久久久久久精品吃奶| 国产精品国产高清国产av|