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

    基于因子分析的不同采煤地表沉陷年限土壤質(zhì)量評(píng)價(jià)

    2016-12-12 09:25:58楊春艷王丹丹丁麗麗
    地下水 2016年6期
    關(guān)鍵詞:土壤水分年限速效

    朱 莉,孫 超,楊春艷,王丹丹,丁麗麗

    (西安科技大學(xué)地質(zhì)與環(huán)境學(xué)院, 陜西 西安 710054)

    ?

    基于因子分析的不同采煤地表沉陷年限土壤質(zhì)量評(píng)價(jià)

    朱 莉,孫 超,楊春艷,王丹丹,丁麗麗

    (西安科技大學(xué)地質(zhì)與環(huán)境學(xué)院, 陜西 西安 710054)

    基于對(duì)風(fēng)蝕水蝕交錯(cuò)地帶的榆神府采煤礦區(qū)地表沉陷地表層土壤理化生指標(biāo)的測(cè)量,選擇18項(xiàng)反映土壤理化性質(zhì)的定量指標(biāo),采用因子分析法對(duì)不同沉陷年限下的土壤質(zhì)量進(jìn)行綜合評(píng)價(jià)。結(jié)果表明:18項(xiàng)土壤指標(biāo)可轉(zhuǎn)化為5個(gè)主因子:土壤水分因子、磷素因子、土壤速效養(yǎng)分因子、全鉀因子和氮素因子;不同沉陷年限5個(gè)公因子差異明顯:未沉陷地土壤養(yǎng)分因子和磷素因子指標(biāo)較好,地表沉陷5 a土壤水分因子和氮素因子指標(biāo)較好,而全鉀因子在沉陷10 a呈現(xiàn)出較高水平;綜合因子分析評(píng)價(jià)得出采煤地表沉陷1~2 a土壤質(zhì)量較差,隨著時(shí)間推移,經(jīng)過土壤的自然演替土壤質(zhì)量有轉(zhuǎn)好的趨勢(shì),但土壤質(zhì)量最好的仍為未沉陷地,說明榆神府采煤沉陷區(qū)地表土壤經(jīng)過10 a自然演替土壤質(zhì)量仍未恢復(fù)到未沉陷區(qū)狀態(tài)。

    榆神府礦區(qū);采煤沉陷;土壤理化性質(zhì)評(píng)價(jià);因子分析

    煤炭作為一次能源在我國(guó)社會(huì)經(jīng)濟(jì)發(fā)展中處于重要地位,但煤炭的大量開采給當(dāng)?shù)卦斐闪藝?yán)重的生態(tài)和社會(huì)問題[1]。因采煤引起土地資源的占用與破壞、水土流失、采空區(qū)地裂縫與地面沉陷等生態(tài)環(huán)境問題日趨嚴(yán)重[2],其中,采煤沉陷對(duì)生態(tài)環(huán)境損害嚴(yán)重且影響范圍廣,已成為人們共同關(guān)注的焦點(diǎn)。

    煤炭在井工開采完成后, 上覆巖層失去支撐而使其應(yīng)力平衡狀態(tài)被破壞, 導(dǎo)致地表沉陷的發(fā)生[3]。采煤沉陷首先會(huì)破壞建筑設(shè)施,造成人類生命及財(cái)產(chǎn)損失[4];其次,采煤沉陷造成地表變形,阻斷天然水系,改變地表徑流和地下水循環(huán)[5];采煤沉陷后,土壤理化性質(zhì)發(fā)生變化,土壤肥力下降,影響植物生長(zhǎng),導(dǎo)致水土流失和土地沙漠化等現(xiàn)象的發(fā)生,嚴(yán)重影響生態(tài)環(huán)境的平衡與穩(wěn)定[6]。

    目前,國(guó)內(nèi)眾多學(xué)者對(duì)榆神府礦區(qū)采煤沉陷后植被群落特征和土壤理化生性質(zhì)變化進(jìn)行了大量的研究,但對(duì)采煤地表沉陷后土壤理化生性質(zhì)隨沉陷年限增加的定量評(píng)價(jià)研究較為薄弱。本文采取因子分析法,對(duì)不同沉陷年限土壤理化生性質(zhì)進(jìn)行系統(tǒng)分析研究,旨在準(zhǔn)確評(píng)價(jià)不同沉陷年限對(duì)土壤理化生性質(zhì)的影響,以期為礦區(qū)生態(tài)恢復(fù)與重建提供理論依據(jù)。

    1 研究區(qū)概況

    研究區(qū)位于陜西省榆林市神木縣北部,屬于黃土高原北部與毛烏素沙地東南邊緣的過渡帶,平均海拔為1 120~1 280 m;年平均降水量為415.0 mm,且降水多集中在 7~8 月的暴雨期;多年平均蒸發(fā)量為 1 788.4 mm,是降水量的4~5倍。研究區(qū)土壤以黃土、紅土、風(fēng)沙土為主,土壤貧瘠,水分流失快。植被帶屬于草原與森林草原的過渡帶,主要植被類型有草原群落(本氏針茅Stipabungeana、百里香Thymus mongolicus等)、低矮灌木(中間錦雞兒Caraganaintermedia、沙棘Hippophaerhamnoides等)、沙生植被(黑沙蒿Artemisia ordosica、沙柳Salix psammophila等)。

    2 研究方法

    2.1 樣地選擇

    本文主要研究大柳塔礦區(qū)覆沙地。首先收集研究區(qū)內(nèi)主要礦井資料確定采空區(qū)范圍,野外通過GPS定點(diǎn)確定沉陷區(qū)邊界或地裂縫分布地點(diǎn),最后確定地表沉陷1 a、2 a、5 a和10 a工作面。每種沉陷年限選擇7塊下墊面狀況(海拔、坡形等)比較相近的沉陷地作為研究樣地,同時(shí)每塊樣地選擇一塊鄰近下墊面狀況一致且沒有受到沉陷影響的坡地作為對(duì)照樣地。由于土壤理化性質(zhì)變化最強(qiáng)烈和植物根系吸收水分最活躍的部位主要集中在0~60 cm土層,因此每塊樣地隨機(jī)選擇5個(gè)點(diǎn),根據(jù)不同土層深度(0~10 cm、10~20 cm、20~30 cm、30~40 cm、40~50 cm、50~60 cm)分別取樣。

    2.2 試驗(yàn)方法

    在樣地內(nèi)每個(gè)取樣點(diǎn)分別挖取一個(gè)面積約為1 m×1 m、深為1 m 的土壤剖面,依據(jù)不同的土壤深度用標(biāo)準(zhǔn)環(huán)刀(100 cm3)分層取樣用于土壤容重的測(cè)定;用20 cm×10 cm鋁制飯盒分層取原狀土用于測(cè)定土壤機(jī)械組成;用水分鉆分別取不同土層土樣放入鋁盒用于分析土壤水分含量;同時(shí)在土壤剖面分3個(gè)點(diǎn)分層取適量土壤樣品,帶回室內(nèi)自然風(fēng)干,進(jìn)行土壤理化指標(biāo)的測(cè)定與分析;用土壤硬度計(jì)在土壤剖面內(nèi)分層直接測(cè)定各土層土壤硬度。

    土壤理化生指標(biāo)測(cè)定:土壤含水量采用烘干法;土壤孔隙度采用土壤密度換算法;土壤容重采用標(biāo)準(zhǔn)環(huán)刀法;有機(jī)質(zhì)測(cè)定采用容量法;全氮測(cè)定采用開氏法;速效氮測(cè)定采用堿解擴(kuò)散法;全磷測(cè)定采用酸溶-鉬銻抗比色法;速效磷測(cè)定采用NaHCO3浸提-鉬銻抗比色法;全鉀測(cè)定采用HF消解-火焰光度計(jì)法;速效鉀測(cè)定采用NH4COOH浸提-火焰光度計(jì)法;pH用電位測(cè)定法;土壤脲酶采用比色法;蔗糖酶采用 3,5-二硝基水楊酸比色法;細(xì)菌、放線菌、真菌采用平板菌落測(cè)定法[7]。

    2.3 數(shù)據(jù)分析

    論文所有數(shù)據(jù)采用SPSS20.0軟件對(duì)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)檢驗(yàn),結(jié)果用平均值±標(biāo)準(zhǔn)偏差(SD)表示,用單因素方差分析(one-way ANOVA)和最小顯著差異法(LSD)比較不同數(shù)據(jù)組間的差異,差異顯著性水平設(shè)定為α=0.05;同時(shí)采用SPSS20.0和Excel軟件對(duì)原始數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化,并對(duì)土壤研究指標(biāo)進(jìn)行因子分析和土壤質(zhì)量評(píng)價(jià)。

    3 結(jié)果與討論

    3.1 評(píng)價(jià)指標(biāo)體系的建立

    土壤質(zhì)量評(píng)價(jià)指標(biāo)以能較顯著影響土壤生產(chǎn)力的土壤養(yǎng)分為主[8]。一般來說, 評(píng)價(jià)土壤質(zhì)量需要土壤物理、化學(xué)和生物指標(biāo)[9],同時(shí)大量研究表明通過因子分析的數(shù)理統(tǒng)計(jì)方法能夠客觀準(zhǔn)確地解釋土壤屬性的變異性,因此本研究基于土壤水分(X1)、孔隙度(X2)、容重(X3)、硬度(X4)、有機(jī)質(zhì)(X5)、全氮(X6)、速效氮(X7)、全磷(X8)、速效磷(X9)、全鉀(X10)、速效鉀(X11)、pH(X12)、脲酶(X13)、磷酸酶(X14)、蔗糖酶(X15)、細(xì)菌(X16)、放線菌(X17)、真菌(X18)共18項(xiàng)土壤理化生指標(biāo)的測(cè)量(表1),運(yùn)用因子分析法對(duì)榆神府礦區(qū)不同采煤沉陷年限土壤的理化生性質(zhì)進(jìn)行評(píng)價(jià),探討影響該區(qū)土壤質(zhì)量的驅(qū)動(dòng)因子。

    表1 不同采煤沉陷年限土壤理化生指標(biāo)分析結(jié)果

    注:不同小寫字母表示不同地表沉陷年限差異顯著(0.05水平)。

    3.2 不同采煤沉陷年限下土壤質(zhì)量的分析比較

    根據(jù)選取的18項(xiàng)土壤指標(biāo),基于因子分析方法進(jìn)行不同采煤沉陷年限下土壤質(zhì)量的驅(qū)動(dòng)因子分析及土壤質(zhì)量評(píng)價(jià),主要步驟為:(1)將選取的指標(biāo)變量建立原始數(shù)據(jù)矩陣;(2)原始數(shù)據(jù)矩陣的標(biāo)準(zhǔn)化;(3)依據(jù)特征值≥1的原則提取公因子,確定公因子數(shù),并運(yùn)用Varimax旋轉(zhuǎn)法對(duì)初始因子進(jìn)行旋轉(zhuǎn),以便給予明確的因子解釋;(4)將公因子作為自變量進(jìn)行回歸分析,計(jì)算各處理的因子得分值,并以因子貢獻(xiàn)率為權(quán)重計(jì)算每個(gè)處理土壤質(zhì)量的綜合得分及排名[10]。

    3.2.1 不同采煤沉陷年限土壤質(zhì)量的因子分析

    因子分析首先是公因子的提取,公因子提取的原則為特征值累計(jì)貢獻(xiàn)率≥80%以及結(jié)合旋轉(zhuǎn)后因子載荷矩陣中變量不出現(xiàn)丟失的前m個(gè)主因子。結(jié)果表明,5個(gè)公因子對(duì)樣本方差的貢獻(xiàn)率為81.3%,即5個(gè)公因子包含了原始數(shù)據(jù)信息總量的81.3%,可以用這5個(gè)公因子代表原來的18個(gè)評(píng)價(jià)指標(biāo),這5個(gè)公因子特征值、貢獻(xiàn)率、累計(jì)貢獻(xiàn)率及通過方差最大化旋轉(zhuǎn)后的因子載荷矩陣見表2。

    表2 旋轉(zhuǎn)成份矩陣

    表2中旋轉(zhuǎn)因子載荷矩陣顯示, 第1個(gè)主因子上,具有高載荷因子的指標(biāo)有土壤水分、孔隙度、容重、硬度、pH、磷酸酶、細(xì)菌、放線菌、真菌,反映了土壤保持水分能力的大小和生物活性情況。其中土壤水分含量是反映土壤質(zhì)量的重要指標(biāo),對(duì)土壤物理、化學(xué)和生物性質(zhì)有重要影響,且土壤水分對(duì)土壤生物活性具有決定性影響,因此第1主因子可稱為土壤水分因子(F1);第2個(gè)主因子上,高載荷因子有全磷、速效磷、蔗糖酶,其中全磷、速效磷反映了土壤中磷的儲(chǔ)量,可見第2主因子可稱為磷素因子(F2);在因子3上,高載荷因子有有機(jī)質(zhì)、速效氮、速效鉀,速效氮、速效鉀,都反映了土壤提供速效養(yǎng)分能力狀況,因此第3主因子可稱為土壤速效養(yǎng)分因子(F3);全鉀在第4因子上有較高的載荷,因此,將第4主因子命名為全鉀因子(F4);因子5上高載荷因子是全氮和脲酶,其都反映了土壤全氮養(yǎng)分貯存與轉(zhuǎn)化能力大小,故第5主因子命名為氮素因子(F5)。

    表3 因子得分系數(shù)矩陣

    3.2.2 不同采煤沉陷年限土壤質(zhì)量的得分及排名

    用線性回歸法計(jì)算各處理的因子得分,因子得分回歸函數(shù)表達(dá)式為:

    Fi=biX

    (1)

    式中:Fi是各公因子的得分、X是不同因子、bi是因子得分系數(shù)矩陣的第i列向量(見表 3)。

    計(jì)算可得出不同采煤沉陷年限土壤質(zhì)量綜合得分(表4),表明不同沉陷年限下土壤理化性質(zhì)在各主分量上存在很大的差異。其中土壤水分因子(F1)中,沉陷1~2 a土壤水分因子損害較為嚴(yán)重,主要原因是采煤地表沉陷初期1~2 a,土壤原始結(jié)構(gòu)遭到破壞而使土壤孔隙度增加,表層土壤水分蒸發(fā)迅速,同時(shí)引起的土壤含水量的降低影響了與水分關(guān)系較大的土壤酶和微生物活性;但由于研究區(qū)處于覆沙區(qū)地表裂隙彌合較快,土壤水分和通氣性在采煤地表沉陷后5 a左右達(dá)到較優(yōu)組合,水分決定的因子中某些土壤指標(biāo)性質(zhì)比未沉陷區(qū)有所改善,但隨著沉陷年限的增加土壤孔隙度不斷增加,土壤水分和通氣性較優(yōu)組合打破,土壤水分相關(guān)因子接近未沉陷區(qū)。磷素因子(F2)中,由于采煤沉陷易造成土壤中磷元素轉(zhuǎn)化受阻,同時(shí)有效磷的流失率也增加,因此未沉陷地磷因子水平較高;土壤速效養(yǎng)分因子(F3)中,在沉陷初期由于地表裂縫土壤速效養(yǎng)分易垂直溶淋至較深的土層處,且沉陷初期土壤水分減小使得地區(qū)植被退化,地區(qū)風(fēng)蝕作用加強(qiáng)而使土壤養(yǎng)分含量較為豐富的黏粒減小,土壤養(yǎng)分迅速降低,因此沉陷初期1~2 a土壤養(yǎng)分因子最差;而隨著沉陷年限的增加,沉陷區(qū)土壤物理結(jié)構(gòu)和地上植被不斷正向演替,使得研究區(qū)土壤速效養(yǎng)分在沉陷后期得到一定程度的恢復(fù),但經(jīng)過10 a恢復(fù)土壤養(yǎng)分因子仍未恢復(fù)至未沉陷狀態(tài)。全鉀因子(F4)中,采煤塌陷使土壤孔隙度提高進(jìn)而加強(qiáng)了土壤好氧微生物的活性,使得鉀元素釋放速率提高,因此地表沉陷5~10 a研究區(qū)鉀元素含量較高。氮素因子(F5)中,由于采煤沉陷后土壤氮元素流失率增加,土壤氮元素迅速下降,但因地表沉陷后5 a左右土壤水分和通氣性的較優(yōu)組合,塌陷裂縫消失,土壤氮元素流失減弱,因此采煤地表沉陷后5年土壤氮可利用率最高。

    把各公因子的特征值貢獻(xiàn)率作為權(quán)數(shù)進(jìn)行加權(quán)求和可得到綜合評(píng)價(jià)指標(biāo)值(表4)。

    (2)

    式中:Vi/p為 SPSS 軟件中表“Total Variance Explained”下“Rotation Sums of SquaredLoadings(旋轉(zhuǎn)后因子對(duì) X 的方差) ”欄的“% of Variance”,F(xiàn)i為各公因子的得分。

    由表4可得出,土壤綜合質(zhì)量最高的是未沉陷區(qū),未沉陷區(qū)土壤受破壞程度較小,土壤熟化程度較高且土壤中酶活性強(qiáng)、微生物的種類及數(shù)量多,使得未沉陷區(qū)土壤磷素因子和土壤速效養(yǎng)分因子含量較高,可為植物生長(zhǎng)提供足夠的養(yǎng)分,綜合得出未發(fā)生采煤沉陷區(qū)土壤質(zhì)量較好;采煤沉陷1~2 a時(shí),原來穩(wěn)定的性質(zhì)受到破壞,土壤中水分因子受影響最為明顯,并迅速降至最低,土壤中各類速效養(yǎng)分的供應(yīng)也受到阻斷,土壤質(zhì)量較差;隨著沉陷時(shí)間的增長(zhǎng)覆沙區(qū)地表裂隙逐漸彌合,土壤的各項(xiàng)理化性質(zhì)趨于穩(wěn)定,最突出的表現(xiàn)是土壤中酶的活性增強(qiáng),大量微生物開始活躍,有益于植物的生長(zhǎng)及土質(zhì)的改變,土壤質(zhì)量有變好的趨勢(shì),同時(shí)由于土壤水分和通氣性在采煤地表沉陷后5 a左右達(dá)到較優(yōu)組合,一定時(shí)間內(nèi)隨時(shí)間的增長(zhǎng)對(duì)沉陷地土壤質(zhì)量改善有較明顯的作用,但隨著土壤水分和孔隙度較優(yōu)組合的破壞,土壤各項(xiàng)指標(biāo)變化幅度將減弱甚至后退,自然恢復(fù)10年亦未恢復(fù)到到未沉陷區(qū)的土壤各項(xiàng)指標(biāo)的狀態(tài)。

    表4 各因子得分及不同采煤沉陷年限土壤質(zhì)量綜合得分

    4 結(jié)語

    通過對(duì)不同沉陷年限條件下的土壤理化生指標(biāo)進(jìn)行因子分析,可以客觀、清晰、 準(zhǔn)確地分析不同采煤沉陷年限下地表土壤質(zhì)量狀況。18項(xiàng)土壤研究指標(biāo)可轉(zhuǎn)化為5個(gè)主因子:土壤水分因子、磷素因子、土壤速效養(yǎng)分因子、全鉀因子和氮素因子;未沉陷地土壤養(yǎng)分因子和磷素因子指標(biāo)較好,地表沉陷5年土壤水分因子和氮素因子指標(biāo)較好,而全鉀因子在沉陷10 a呈現(xiàn)出較高水平;綜合因子分析評(píng)價(jià)得出采煤沉陷1~2 a時(shí),土壤原始結(jié)構(gòu)遭到破壞,土壤中水分相關(guān)因子受影響最為明顯,土壤酶活性和微生物數(shù)量降低,土壤質(zhì)量較差。地表沉陷后5年由于土壤水分和通氣性達(dá)到較優(yōu)組合,一定時(shí)間內(nèi)沉陷地土壤質(zhì)量的改善作用較明顯。但土壤綜合質(zhì)量最高的是未沉陷區(qū),土壤的各項(xiàng)理化生性質(zhì)較為穩(wěn)定,土壤養(yǎng)分效應(yīng)最為明顯,為地上植被生長(zhǎng)提供養(yǎng)分的能力充足。研究表明榆神府采煤沉陷區(qū)地表土壤隨著時(shí)間推移進(jìn)行土壤演替,土壤質(zhì)量有轉(zhuǎn)好的趨勢(shì),但經(jīng)過10 a自然演替土壤各項(xiàng)指標(biāo)仍未恢復(fù)到未沉陷區(qū)狀態(tài)。

    [1]張延旭,畢銀麗,陳書琳,等.半干旱風(fēng)沙區(qū)采煤后裂縫發(fā)育對(duì)土壤水分的影響[J].環(huán)境科學(xué)與技術(shù).2015,38(3):11-14.

    [2]Zhou D W, Wu K, Cheng G L, et al. Mechanism of mining subsidence in coal mining area with thick alluvium soil in China[J]. Arabian Journal of Geosciences, 2015, 8(4): 1855-1867.

    [3]黎煒,陳龍乾,周天建,等.我國(guó)采煤沉陷地土壤質(zhì)量研究進(jìn)展[J].煤炭科學(xué)技術(shù).2011,39(5):125-128.

    [4]王燕.煤炭開采對(duì)生態(tài)環(huán)境的影響及治理對(duì)策[J].煤炭科學(xué)技術(shù).2009,37(12):125-128.

    [5]王楨,高雪,高伯賢.榆神礦區(qū)小壕兔區(qū)采煤對(duì)地下水環(huán)境影響分析[J].地下水.2014,36(6):114-116.

    [6]Wright I A, McCarthy B, Belmer N, et al. Subsidence from an underground coal mine and mine wastewater discharge causing water pollution and degradation of aquatic ecosystems[J]. Water, Air, & Soil Pollution, 2015, 226(10): 1-14.

    [7]王琦,全占軍,韓煜,等.采煤塌陷對(duì)風(fēng)沙區(qū)土壤性質(zhì)的影響[J].中國(guó)水土保持科學(xué).2013,11(6):110-118.

    [8]劉曉冰,刑寶山.土壤質(zhì)量及其評(píng)價(jià)指標(biāo)[J].農(nóng)業(yè)系統(tǒng)科學(xué)與綜合研究.2002,18(2): 109-111.

    [9]劉世梁,傅伯杰,劉國(guó)華,等.我國(guó)土壤質(zhì)量及其評(píng)價(jià)研究的進(jìn)展[J].土壤通報(bào).2006,37(1): 137-143.

    [10]陳琳,趙陟峰,趙廷寧,等.基于因子與聚類分析的不同林分土壤理化性質(zhì)評(píng)價(jià)[J].水土保持研究.2011,18(5): 192-196.

    Evaluation of Soil Quality in Different Years of Coal Mining Subsidence Based on Factor Analysis

    ZHU Li,SUN Chao,YANG Chun-yan,WANG Dan-dan,DING Li-li

    (College of Geology & Environment, Xi’an University of Science and Technology, Xi’an 710054, Shaanxi)

    Based on the investigation of soil physical-chemical-biological properties on Yushenfu coal mining subsidence mining area which be in wind-water erosion crisscross region, 18 soil properties were used to evaluatesoil quality in different years of coal mining subsidence using factor analysis method. The results showed that 18 soil propertiescan be classified into 5 main factors, namely soil moisture factor, phosphorus factor, soil nutrient factor, total potassium, and nitrogen factor. These 5 factorstook significantly different in different years of coal mining subsidence: un-subsidencearea, 5 years subsidence area, and 10 years subsidence area displayed higher levels of soil phosphorus and nutrient factors, moisture and nitrogen factors, and total potassium than others plotsrespectively. Evaluationof comprehensive factor analysis showed that soil quality of 1~2 years subsidence area is the worst, then the soil quality (physical-chemical-biological properties) gradually improved through natural soil succession with increasing years of subsidence. The best quality, however, still appeared at un-subsidence area. These results reveal that soil quality of Yushenfu coal mining subsidence area has not been restored tooriginal levels after 10 years natural succession.

    Yushenfu mining area;coal mining subsidence;evaluation of soil physical and chemical properties;factor analysis

    2016-05-18

    2015年國(guó)家級(jí)大學(xué)生創(chuàng)新創(chuàng)業(yè)訓(xùn)練計(jì)劃(201510704194)

    朱莉(1994-),女,寧夏賀蘭人,主攻方向:礦山環(huán)境修復(fù)。

    TE122.2+3

    A

    1004-1184(2016)06-0157-04

    猜你喜歡
    土壤水分年限速效
    影響種公牛使用年限的幾個(gè)因素與解決辦法
    不同殺螨劑對(duì)柑桔紅蜘蛛速效與持效性防治田間試驗(yàn)
    不同產(chǎn)地、生長(zhǎng)年限銀杏葉總多酚含量比較
    中成藥(2017年6期)2017-06-13 07:30:35
    西藏高原土壤水分遙感監(jiān)測(cè)方法研究
    土壤速效鉀測(cè)定方法的影響因素
    綠色科技(2016年16期)2016-10-11 06:53:51
    體外發(fā)酵法評(píng)定不同茬次和生長(zhǎng)年限苜蓿的營(yíng)養(yǎng)價(jià)值
    施用不同有機(jī)物料對(duì)山地果園紅壤速效氮磷鉀含量的影響
    豫中植煙區(qū)土壤速效養(yǎng)分的空間變異分析
    不同覆蓋措施對(duì)棗園土壤水分和溫度的影響
    植被覆蓋區(qū)土壤水分反演研究——以北京市為例
    欧美不卡视频在线免费观看| 成人无遮挡网站| 干丝袜人妻中文字幕| 可以在线观看毛片的网站| 黄色日韩在线| 精品国产三级普通话版| 亚洲欧美日韩卡通动漫| 人妻久久中文字幕网| 亚洲中文日韩欧美视频| 97碰自拍视频| 露出奶头的视频| 三级毛片av免费| 亚洲在线自拍视频| 成人特级黄色片久久久久久久| 美女大奶头视频| 成人漫画全彩无遮挡| 国产精品一区www在线观看| 色视频www国产| 国产精品久久电影中文字幕| 99热这里只有是精品在线观看| 亚洲av免费高清在线观看| videossex国产| 久久久久国内视频| 99在线视频只有这里精品首页| 桃色一区二区三区在线观看| 欧美极品一区二区三区四区| 精品久久久噜噜| 麻豆一二三区av精品| 91狼人影院| 欧美人与善性xxx| 国产爱豆传媒在线观看| 亚洲精品一区av在线观看| 精品国产三级普通话版| 丰满的人妻完整版| 长腿黑丝高跟| 亚洲最大成人av| 99久久无色码亚洲精品果冻| 欧美成人a在线观看| 最后的刺客免费高清国语| 天美传媒精品一区二区| 国产人妻一区二区三区在| 亚洲精品国产av成人精品 | 久久久久久久久久久丰满| 精品久久久久久久久久免费视频| 国产单亲对白刺激| 亚洲久久久久久中文字幕| 国产在视频线在精品| 成人亚洲精品av一区二区| 国产乱人偷精品视频| 成人一区二区视频在线观看| .国产精品久久| 不卡一级毛片| 99国产极品粉嫩在线观看| 久久中文看片网| 一级毛片电影观看 | 国产69精品久久久久777片| 久久久久国内视频| 简卡轻食公司| 国产日本99.免费观看| 一级av片app| 午夜福利成人在线免费观看| 久久久色成人| 精品一区二区三区视频在线| 欧美性猛交╳xxx乱大交人| 免费大片18禁| 日韩成人av中文字幕在线观看 | 少妇高潮的动态图| 精品不卡国产一区二区三区| 又粗又爽又猛毛片免费看| 黄片wwwwww| 哪里可以看免费的av片| 男女之事视频高清在线观看| 综合色丁香网| 午夜精品在线福利| 女的被弄到高潮叫床怎么办| 国产色爽女视频免费观看| 午夜日韩欧美国产| 看十八女毛片水多多多| 国产免费一级a男人的天堂| 两个人的视频大全免费| 嫩草影视91久久| 精品一区二区免费观看| 99热网站在线观看| 亚洲内射少妇av| 亚洲国产欧美人成| a级毛片a级免费在线| 激情 狠狠 欧美| 精品免费久久久久久久清纯| 97在线视频观看| 中文字幕熟女人妻在线| 嫩草影院精品99| 亚洲av成人av| АⅤ资源中文在线天堂| 国产真实伦视频高清在线观看| 亚洲欧美成人综合另类久久久 | 91久久精品电影网| 精品免费久久久久久久清纯| 国产 一区 欧美 日韩| 精品一区二区三区人妻视频| 亚洲熟妇中文字幕五十中出| 桃色一区二区三区在线观看| 日本成人三级电影网站| 97超视频在线观看视频| 麻豆成人午夜福利视频| 国产成人一区二区在线| 菩萨蛮人人尽说江南好唐韦庄 | 久久中文看片网| 精品欧美国产一区二区三| 亚洲丝袜综合中文字幕| 男人舔奶头视频| 国产精品精品国产色婷婷| 久久久国产成人精品二区| 亚洲国产精品合色在线| 国产伦精品一区二区三区视频9| 亚洲精品国产av成人精品 | 综合色丁香网| 真实男女啪啪啪动态图| 日本五十路高清| 婷婷亚洲欧美| 国产成人91sexporn| 亚洲四区av| 亚洲丝袜综合中文字幕| 午夜a级毛片| 99riav亚洲国产免费| 淫妇啪啪啪对白视频| 日韩欧美 国产精品| 亚洲无线在线观看| 亚洲精品乱码久久久v下载方式| 国产精品亚洲一级av第二区| 久久人人爽人人爽人人片va| 久久精品国产99精品国产亚洲性色| 日韩,欧美,国产一区二区三区 | 成人三级黄色视频| 亚洲av免费高清在线观看| 亚洲av中文字字幕乱码综合| 日韩精品有码人妻一区| 别揉我奶头~嗯~啊~动态视频| 中国美白少妇内射xxxbb| eeuss影院久久| 97热精品久久久久久| 国产成人a区在线观看| 欧美+亚洲+日韩+国产| 我要搜黄色片| 亚洲在线观看片| 久久久午夜欧美精品| 岛国在线免费视频观看| 免费人成在线观看视频色| 久久欧美精品欧美久久欧美| 一个人观看的视频www高清免费观看| 国产在视频线在精品| 成人av一区二区三区在线看| 国产高清有码在线观看视频| 最近最新中文字幕大全电影3| 成年女人看的毛片在线观看| 日韩欧美免费精品| 国产久久久一区二区三区| 欧美日韩一区二区视频在线观看视频在线 | 免费av不卡在线播放| 国产高清视频在线观看网站| 亚洲av不卡在线观看| 你懂的网址亚洲精品在线观看 | 日韩三级伦理在线观看| 美女 人体艺术 gogo| 国产午夜精品论理片| 久久精品国产亚洲av香蕉五月| 国产精品嫩草影院av在线观看| 国产亚洲精品综合一区在线观看| 国产精品免费一区二区三区在线| 性欧美人与动物交配| 国产老妇女一区| 精品久久久久久久久av| 又爽又黄无遮挡网站| 国产精品亚洲一级av第二区| 一个人免费在线观看电影| 久久精品国产鲁丝片午夜精品| 91久久精品电影网| 欧美xxxx黑人xx丫x性爽| 日韩大尺度精品在线看网址| 2021天堂中文幕一二区在线观| 亚洲精品影视一区二区三区av| 大又大粗又爽又黄少妇毛片口| 村上凉子中文字幕在线| 最近在线观看免费完整版| 亚洲性久久影院| 久久久成人免费电影| 午夜爱爱视频在线播放| 国产高清不卡午夜福利| 久久综合国产亚洲精品| 免费观看精品视频网站| 日本 av在线| 最近在线观看免费完整版| 特级一级黄色大片| 女生性感内裤真人,穿戴方法视频| 一级av片app| 欧美三级亚洲精品| 欧美激情在线99| 99热精品在线国产| 一进一出好大好爽视频| 久久久国产成人精品二区| 中文亚洲av片在线观看爽| 露出奶头的视频| 亚洲国产精品国产精品| 日韩一本色道免费dvd| 日日摸夜夜添夜夜爱| 欧美最黄视频在线播放免费| eeuss影院久久| 国产真实乱freesex| 内射极品少妇av片p| 99热只有精品国产| 亚洲久久久久久中文字幕| 午夜福利高清视频| 午夜免费激情av| 亚洲精华国产精华液的使用体验 | 国产精品久久视频播放| 搞女人的毛片| 国产精品日韩av在线免费观看| 色av中文字幕| 亚洲精品在线观看二区| 国产欧美日韩一区二区精品| 欧美激情久久久久久爽电影| 国产一区亚洲一区在线观看| 国产精品久久电影中文字幕| 国产一区二区三区av在线 | 午夜精品一区二区三区免费看| 午夜影院日韩av| 国产亚洲av嫩草精品影院| 伦精品一区二区三区| 亚洲自偷自拍三级| 一卡2卡三卡四卡精品乱码亚洲| 日韩欧美 国产精品| 卡戴珊不雅视频在线播放| 国产伦精品一区二区三区四那| 国产一区二区在线观看日韩| 色视频www国产| 亚洲中文字幕一区二区三区有码在线看| 超碰av人人做人人爽久久| 国产在线精品亚洲第一网站| 久久这里只有精品中国| 国产精品久久久久久久久免| 此物有八面人人有两片| 亚洲精品国产av成人精品 | 一本久久中文字幕| 又黄又爽又免费观看的视频| 成人漫画全彩无遮挡| 国产高清三级在线| www.色视频.com| 97超级碰碰碰精品色视频在线观看| 日本黄色视频三级网站网址| 久久韩国三级中文字幕| 免费看光身美女| 成年女人永久免费观看视频| 又黄又爽又刺激的免费视频.| 91av网一区二区| 久99久视频精品免费| 国产精品1区2区在线观看.| 我要看日韩黄色一级片| 国产伦精品一区二区三区四那| 丰满的人妻完整版| 中文字幕免费在线视频6| 欧美+亚洲+日韩+国产| 一本久久中文字幕| 一个人看视频在线观看www免费| 少妇的逼好多水| 国产乱人偷精品视频| 国产精品无大码| 亚洲欧美日韩高清专用| 在线天堂最新版资源| 最近中文字幕高清免费大全6| 国产男靠女视频免费网站| 国产亚洲精品久久久com| 黄色配什么色好看| 在线观看美女被高潮喷水网站| 黄色欧美视频在线观看| 99热只有精品国产| 乱码一卡2卡4卡精品| 欧美日韩精品成人综合77777| 久久精品夜色国产| 欧美一区二区精品小视频在线| 亚洲精品久久国产高清桃花| 三级经典国产精品| 99国产精品一区二区蜜桃av| 伊人久久精品亚洲午夜| 日韩欧美免费精品| 麻豆一二三区av精品| 97超视频在线观看视频| 又爽又黄a免费视频| 亚洲精品成人久久久久久| 99在线视频只有这里精品首页| 波野结衣二区三区在线| 亚洲欧美日韩卡通动漫| 日韩,欧美,国产一区二区三区 | www.色视频.com| 国产精品久久久久久久久免| 少妇人妻一区二区三区视频| 国产精品一二三区在线看| 人妻夜夜爽99麻豆av| 欧美成人a在线观看| 国产免费男女视频| 俄罗斯特黄特色一大片| 欧美日本视频| 亚洲精品国产成人久久av| 国产色爽女视频免费观看| 日本 av在线| 久久草成人影院| 亚洲精品亚洲一区二区| 成人性生交大片免费视频hd| 色哟哟哟哟哟哟| 插阴视频在线观看视频| 变态另类成人亚洲欧美熟女| 18禁在线无遮挡免费观看视频 | 嫩草影院入口| 国产黄片美女视频| 国产成人a∨麻豆精品| 伊人久久精品亚洲午夜| 日本爱情动作片www.在线观看 | 性插视频无遮挡在线免费观看| 国产伦在线观看视频一区| 女生性感内裤真人,穿戴方法视频| 亚洲人成网站高清观看| 亚洲18禁久久av| 直男gayav资源| av国产免费在线观看| 色视频www国产| 神马国产精品三级电影在线观看| 少妇人妻精品综合一区二区 | 美女大奶头视频| 韩国av在线不卡| 国产乱人视频| 最近的中文字幕免费完整| 色哟哟哟哟哟哟| 午夜福利视频1000在线观看| 精品一区二区三区视频在线观看免费| 日韩欧美精品免费久久| 中文亚洲av片在线观看爽| videossex国产| 国产黄片美女视频| 免费搜索国产男女视频| 伦理电影大哥的女人| 综合色av麻豆| 欧美最黄视频在线播放免费| 午夜免费激情av| 久久久久国产精品人妻aⅴ院| 欧美性猛交╳xxx乱大交人| 中文字幕av在线有码专区| 天堂√8在线中文| 国产高清三级在线| av免费在线看不卡| 男女那种视频在线观看| 美女黄网站色视频| 国产乱人偷精品视频| 久久精品国产亚洲av香蕉五月| 精品午夜福利视频在线观看一区| 少妇被粗大猛烈的视频| 亚洲国产精品久久男人天堂| 日韩制服骚丝袜av| 直男gayav资源| avwww免费| 晚上一个人看的免费电影| 国产精品无大码| 麻豆av噜噜一区二区三区| 亚洲av中文av极速乱| 99在线视频只有这里精品首页| 国内精品美女久久久久久| 亚洲精品久久国产高清桃花| 中文亚洲av片在线观看爽| 插阴视频在线观看视频| 国产精品久久久久久久电影| 老司机影院成人| 午夜福利在线观看免费完整高清在 | 美女内射精品一级片tv| 变态另类成人亚洲欧美熟女| 欧美xxxx黑人xx丫x性爽| 精品福利观看| 在线观看午夜福利视频| 国产视频一区二区在线看| 99久久无色码亚洲精品果冻| av天堂中文字幕网| 美女cb高潮喷水在线观看| 高清午夜精品一区二区三区 | 久久久精品大字幕| 深夜精品福利| 99热精品在线国产| 国产成人福利小说| 亚洲欧美成人精品一区二区| 搡老妇女老女人老熟妇| 国内精品一区二区在线观看| 天天躁日日操中文字幕| 一级黄色大片毛片| 国产精品嫩草影院av在线观看| 国产午夜福利久久久久久| 亚洲四区av| 免费在线观看影片大全网站| 亚洲成人久久爱视频| 欧美+亚洲+日韩+国产| 最近的中文字幕免费完整| 日韩一本色道免费dvd| 精品无人区乱码1区二区| 日韩精品青青久久久久久| 狂野欧美白嫩少妇大欣赏| 欧美+亚洲+日韩+国产| 久久精品国产亚洲网站| 午夜福利18| 成年女人毛片免费观看观看9| 久久精品人妻少妇| 色视频www国产| 色在线成人网| 一卡2卡三卡四卡精品乱码亚洲| 草草在线视频免费看| 亚洲图色成人| 国产私拍福利视频在线观看| 欧美一级a爱片免费观看看| 国产成人福利小说| 亚洲精华国产精华液的使用体验 | 成人二区视频| 在线观看一区二区三区| 久久精品国产亚洲av天美| 国产私拍福利视频在线观看| 亚洲欧美日韩高清在线视频| 国产精品日韩av在线免费观看| 国产精品不卡视频一区二区| 我要搜黄色片| 波野结衣二区三区在线| 毛片女人毛片| 久久欧美精品欧美久久欧美| 又爽又黄无遮挡网站| 久久精品国产自在天天线| 高清毛片免费看| 99热只有精品国产| 日日摸夜夜添夜夜添小说| 亚洲性夜色夜夜综合| av在线亚洲专区| 成人午夜高清在线视频| 国产老妇女一区| 欧美色视频一区免费| 日韩人妻高清精品专区| 欧美3d第一页| 精品一区二区免费观看| 午夜精品一区二区三区免费看| 久久久精品94久久精品| 日韩成人伦理影院| 丝袜美腿在线中文| 国产极品精品免费视频能看的| 亚洲四区av| 又粗又爽又猛毛片免费看| 在线看三级毛片| 日韩欧美精品v在线| 亚洲美女搞黄在线观看 | 成人精品一区二区免费| 乱人视频在线观看| 国内精品一区二区在线观看| 日本一二三区视频观看| 国产高清视频在线播放一区| 亚洲一区高清亚洲精品| avwww免费| 中文在线观看免费www的网站| 久久精品91蜜桃| 日本一二三区视频观看| 日韩一本色道免费dvd| 久久久国产成人精品二区| 久久精品国产自在天天线| 老师上课跳d突然被开到最大视频| 在线免费观看的www视频| 成人综合一区亚洲| 深爱激情五月婷婷| 国产精品国产高清国产av| 毛片一级片免费看久久久久| 全区人妻精品视频| 国产麻豆成人av免费视频| 2021天堂中文幕一二区在线观| 精品午夜福利视频在线观看一区| 久久久久久久久大av| 国产欧美日韩精品一区二区| 少妇裸体淫交视频免费看高清| 老熟妇乱子伦视频在线观看| 毛片女人毛片| 日韩精品青青久久久久久| 欧美成人免费av一区二区三区| 久久99热这里只有精品18| 淫妇啪啪啪对白视频| 精品午夜福利在线看| 亚洲人成网站在线观看播放| 欧美zozozo另类| av福利片在线观看| 亚洲欧美日韩高清专用| 国产一区二区三区av在线 | 人人妻人人澡人人爽人人夜夜 | 久久久久国产精品人妻aⅴ院| 久久久久国内视频| 在线观看av片永久免费下载| 成年女人毛片免费观看观看9| 我的女老师完整版在线观看| 欧美最黄视频在线播放免费| 日日干狠狠操夜夜爽| 女同久久另类99精品国产91| 国产精品美女特级片免费视频播放器| 成人三级黄色视频| 中文在线观看免费www的网站| 久久久久久久午夜电影| 国产探花极品一区二区| 午夜免费激情av| 日韩人妻高清精品专区| 成年女人毛片免费观看观看9| 国产私拍福利视频在线观看| 麻豆一二三区av精品| 欧美激情在线99| av.在线天堂| 人妻少妇偷人精品九色| 色噜噜av男人的天堂激情| 午夜久久久久精精品| 露出奶头的视频| 国产毛片a区久久久久| 男人的好看免费观看在线视频| 亚洲激情五月婷婷啪啪| 成人亚洲精品av一区二区| 国产亚洲91精品色在线| 成人永久免费在线观看视频| 午夜老司机福利剧场| 老熟妇仑乱视频hdxx| 亚洲一区高清亚洲精品| 日韩三级伦理在线观看| 亚洲国产日韩欧美精品在线观看| 最后的刺客免费高清国语| 精品人妻偷拍中文字幕| 老司机福利观看| 久久人人爽人人爽人人片va| 欧美人与善性xxx| 日本一二三区视频观看| av天堂在线播放| 日本三级黄在线观看| 国产成人福利小说| 久久鲁丝午夜福利片| 久久久色成人| 欧美日韩在线观看h| 99久国产av精品| 91av网一区二区| 中文字幕免费在线视频6| 天堂网av新在线| 色哟哟哟哟哟哟| 日日摸夜夜添夜夜爱| 伦理电影大哥的女人| 六月丁香七月| 亚洲成人av在线免费| 午夜亚洲福利在线播放| 看免费成人av毛片| 18禁在线无遮挡免费观看视频 | 综合色av麻豆| 国内久久婷婷六月综合欲色啪| 少妇被粗大猛烈的视频| 直男gayav资源| 国内揄拍国产精品人妻在线| 可以在线观看的亚洲视频| 夜夜看夜夜爽夜夜摸| 欧美激情在线99| 午夜视频国产福利| АⅤ资源中文在线天堂| 日本熟妇午夜| 精品人妻熟女av久视频| avwww免费| 国产伦精品一区二区三区四那| 一个人观看的视频www高清免费观看| 亚洲欧美精品自产自拍| 特大巨黑吊av在线直播| 欧美精品国产亚洲| 99久久九九国产精品国产免费| 午夜精品在线福利| 日本色播在线视频| 日日摸夜夜添夜夜添小说| 69av精品久久久久久| 成人毛片a级毛片在线播放| 中文字幕免费在线视频6| 欧美性猛交黑人性爽| 人人妻人人澡人人爽人人夜夜 | 我的老师免费观看完整版| 国产单亲对白刺激| 精品午夜福利在线看| 少妇高潮的动态图| 狂野欧美白嫩少妇大欣赏| 婷婷六月久久综合丁香| 日韩成人伦理影院| 亚洲精品国产av成人精品 | 成人美女网站在线观看视频| 亚洲av第一区精品v没综合| 国产精品一区二区性色av| 欧美高清性xxxxhd video| 有码 亚洲区| 亚洲色图av天堂| 久久精品国产亚洲网站| 久久精品国产清高在天天线| 两性午夜刺激爽爽歪歪视频在线观看| 久久久色成人| 男女做爰动态图高潮gif福利片| 精品久久久久久久人妻蜜臀av| 18禁裸乳无遮挡免费网站照片| av天堂在线播放| 国内精品久久久久精免费| 天美传媒精品一区二区| 我的老师免费观看完整版| 国内久久婷婷六月综合欲色啪| 国产一区二区三区在线臀色熟女| 国产午夜福利久久久久久| 99久国产av精品| 亚洲精品国产成人久久av| 精品不卡国产一区二区三区| 中出人妻视频一区二区| 亚洲内射少妇av| 插逼视频在线观看| 久久久国产成人精品二区| 亚洲国产欧美人成| 免费av不卡在线播放| 男人狂女人下面高潮的视频| 五月玫瑰六月丁香| 免费av不卡在线播放| 国产视频内射| 国产黄a三级三级三级人| 精品福利观看|