• <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ū)土壤水分反演研究——以北京市為例
    色哟哟·www| av在线播放精品| 男男h啪啪无遮挡| 国产视频内射| 成人午夜精彩视频在线观看| 制服丝袜香蕉在线| 亚洲一级一片aⅴ在线观看| 校园人妻丝袜中文字幕| 男人爽女人下面视频在线观看| 伊人久久国产一区二区| 欧美日韩综合久久久久久| 日本午夜av视频| 国产一区二区三区综合在线观看 | 久热这里只有精品99| 高清视频免费观看一区二区| .国产精品久久| 亚洲精品aⅴ在线观看| 两个人的视频大全免费| 五月伊人婷婷丁香| 国产精品偷伦视频观看了| .国产精品久久| 成人亚洲精品一区在线观看| 国产探花极品一区二区| 国产熟女欧美一区二区| 男人舔奶头视频| 少妇人妻 视频| 亚洲av福利一区| 日韩中文字幕视频在线看片| 国产一区亚洲一区在线观看| 国产男女内射视频| av又黄又爽大尺度在线免费看| 精华霜和精华液先用哪个| 久久久久精品久久久久真实原创| 国产欧美日韩综合在线一区二区 | 国产亚洲午夜精品一区二区久久| 久久免费观看电影| 日本黄色日本黄色录像| 又爽又黄a免费视频| 久久久久网色| 简卡轻食公司| 亚洲国产欧美日韩在线播放 | 精品人妻偷拍中文字幕| 成人亚洲欧美一区二区av| 国产免费又黄又爽又色| 国产淫片久久久久久久久| 国产精品国产三级国产av玫瑰| 人人妻人人澡人人爽人人夜夜| 青春草国产在线视频| 精品一区在线观看国产| 亚洲欧美一区二区三区黑人 | 国产精品一二三区在线看| 国产综合精华液| 视频区图区小说| 伊人亚洲综合成人网| 精品国产一区二区三区久久久樱花| 成人亚洲欧美一区二区av| 中文资源天堂在线| 亚洲在久久综合| 久久毛片免费看一区二区三区| 99久久精品一区二区三区| 国产av一区二区精品久久| 亚洲av成人精品一区久久| 99久久精品一区二区三区| 亚洲丝袜综合中文字幕| av视频免费观看在线观看| 日韩成人伦理影院| 亚洲,欧美,日韩| 日韩成人伦理影院| 天堂中文最新版在线下载| 日韩精品有码人妻一区| 中文字幕精品免费在线观看视频 | 亚洲成人一二三区av| 日韩成人伦理影院| 26uuu在线亚洲综合色| 少妇裸体淫交视频免费看高清| 最近手机中文字幕大全| 免费在线观看成人毛片| 嘟嘟电影网在线观看| 亚洲精品456在线播放app| 国产乱人偷精品视频| 美女中出高潮动态图| 精品视频人人做人人爽| 欧美变态另类bdsm刘玥| 国产精品国产av在线观看| 亚洲在久久综合| 久久女婷五月综合色啪小说| 久久女婷五月综合色啪小说| 国产精品一区二区性色av| 亚洲精品日韩在线中文字幕| 午夜av观看不卡| av又黄又爽大尺度在线免费看| 91久久精品国产一区二区成人| 成年人午夜在线观看视频| 精品久久久噜噜| 天天操日日干夜夜撸| 中国三级夫妇交换| 乱系列少妇在线播放| 久久这里有精品视频免费| 亚洲精品自拍成人| 亚洲,一卡二卡三卡| 成年av动漫网址| 日韩熟女老妇一区二区性免费视频| 26uuu在线亚洲综合色| 午夜福利视频精品| 亚洲av成人精品一区久久| 成人毛片60女人毛片免费| 欧美精品国产亚洲| 国产日韩欧美视频二区| 少妇的逼水好多| 国产精品国产三级专区第一集| 一级毛片久久久久久久久女| 亚洲不卡免费看| 久久久精品免费免费高清| 久久久久久久久久久久大奶| 国产综合精华液| 国产精品久久久久久精品古装| 国产精品99久久久久久久久| 少妇被粗大猛烈的视频| 插阴视频在线观看视频| 亚洲人成网站在线观看播放| 久久精品国产自在天天线| 毛片一级片免费看久久久久| 亚洲av成人精品一区久久| 中文天堂在线官网| 伦理电影免费视频| 国内精品宾馆在线| 亚洲精品456在线播放app| 亚洲国产毛片av蜜桃av| 成人免费观看视频高清| 交换朋友夫妻互换小说| 哪个播放器可以免费观看大片| 国产精品国产三级国产av玫瑰| 成人18禁高潮啪啪吃奶动态图 | 亚洲欧洲日产国产| 国产在视频线精品| 亚洲国产精品国产精品| 国产在线免费精品| 欧美日韩在线观看h| 国产一区亚洲一区在线观看| 欧美日韩国产mv在线观看视频| 三级国产精品片| 中文字幕人妻熟人妻熟丝袜美| kizo精华| 精品人妻偷拍中文字幕| 大码成人一级视频| 国产精品久久久久成人av| 亚洲美女搞黄在线观看| 精品一品国产午夜福利视频| 成年女人在线观看亚洲视频| 尾随美女入室| 伦理电影大哥的女人| 高清黄色对白视频在线免费看 | 国产爽快片一区二区三区| 熟妇人妻不卡中文字幕| 久热这里只有精品99| 精品国产一区二区三区久久久樱花| 蜜臀久久99精品久久宅男| 自拍欧美九色日韩亚洲蝌蚪91 | 看十八女毛片水多多多| 亚洲激情五月婷婷啪啪| 日韩av免费高清视频| 热re99久久精品国产66热6| 国产成人freesex在线| www.色视频.com| 丝袜脚勾引网站| 在线播放无遮挡| 成人国产麻豆网| 女性被躁到高潮视频| 国产av一区二区精品久久| 日本色播在线视频| 五月天丁香电影| 日韩在线高清观看一区二区三区| 久久精品久久久久久噜噜老黄| 大香蕉97超碰在线| www.av在线官网国产| 成人国产av品久久久| 晚上一个人看的免费电影| av天堂中文字幕网| 国产精品欧美亚洲77777| 嘟嘟电影网在线观看| 高清午夜精品一区二区三区| av天堂久久9| 三上悠亚av全集在线观看 | 三级国产精品片| 成人亚洲精品一区在线观看| 亚洲av电影在线观看一区二区三区| 亚洲人成网站在线播| 成年美女黄网站色视频大全免费 | 亚洲va在线va天堂va国产| 高清在线视频一区二区三区| 亚洲欧美日韩东京热| 免费av不卡在线播放| 日韩欧美 国产精品| 国产一区二区三区av在线| 成人影院久久| 嫩草影院入口| 国产精品99久久99久久久不卡 | 亚洲色图综合在线观看| 午夜激情福利司机影院| 免费少妇av软件| 美女大奶头黄色视频| 欧美另类一区| 高清欧美精品videossex| 亚洲国产色片| 国产无遮挡羞羞视频在线观看| 中文字幕人妻熟人妻熟丝袜美| 午夜免费观看性视频| 九草在线视频观看| 自拍偷自拍亚洲精品老妇| 女人精品久久久久毛片| 精品少妇久久久久久888优播| 国产美女午夜福利| 97在线人人人人妻| 国产成人一区二区在线| 亚洲国产精品一区二区三区在线| 日本wwww免费看| av在线播放精品| 国产精品久久久久久久久免| 两个人免费观看高清视频 | 十八禁高潮呻吟视频 | 国产欧美日韩综合在线一区二区 | 久久99热6这里只有精品| 免费黄频网站在线观看国产| 亚洲成人av在线免费| 激情五月婷婷亚洲| 日本欧美视频一区| 亚洲精品日本国产第一区| 日韩精品免费视频一区二区三区 | 欧美激情国产日韩精品一区| 日本猛色少妇xxxxx猛交久久| 成年av动漫网址| 亚洲精品乱码久久久久久按摩| 精品少妇内射三级| 国产一区二区三区综合在线观看 | 爱豆传媒免费全集在线观看| 亚洲精品456在线播放app| 国产成人精品无人区| 一二三四中文在线观看免费高清| 久久国产精品男人的天堂亚洲 | 交换朋友夫妻互换小说| 最近手机中文字幕大全| 国产一区二区三区综合在线观看 | 国产精品国产三级国产av玫瑰| 超碰97精品在线观看| 永久免费av网站大全| 免费观看无遮挡的男女| 日日摸夜夜添夜夜添av毛片| 亚洲内射少妇av| 亚洲不卡免费看| 欧美亚洲 丝袜 人妻 在线| 一区二区三区乱码不卡18| 久久久欧美国产精品| 菩萨蛮人人尽说江南好唐韦庄| 建设人人有责人人尽责人人享有的| 国产av码专区亚洲av| 嫩草影院新地址| 国产淫语在线视频| 日韩大片免费观看网站| 高清毛片免费看| 欧美亚洲 丝袜 人妻 在线| 99久久精品国产国产毛片| 欧美激情极品国产一区二区三区 | 一本—道久久a久久精品蜜桃钙片| 成人二区视频| 免费观看在线日韩| 色哟哟·www| 欧美成人精品欧美一级黄| 高清视频免费观看一区二区| 一级爰片在线观看| 中文字幕制服av| 寂寞人妻少妇视频99o| 一级毛片电影观看| 少妇熟女欧美另类| 亚洲性久久影院| 一区二区av电影网| 亚洲精品第二区| 五月开心婷婷网| 国产男女内射视频| 内地一区二区视频在线| 99精国产麻豆久久婷婷| xxx大片免费视频| 丝袜脚勾引网站| 两个人免费观看高清视频 | 国产在视频线精品| 夜夜爽夜夜爽视频| 日日摸夜夜添夜夜添av毛片| 哪个播放器可以免费观看大片| 好男人视频免费观看在线| 国产日韩欧美在线精品| 最新的欧美精品一区二区| 毛片一级片免费看久久久久| 两个人免费观看高清视频 | 狠狠精品人妻久久久久久综合| 日本-黄色视频高清免费观看| 久久综合国产亚洲精品| 韩国高清视频一区二区三区| 大陆偷拍与自拍| a级毛片免费高清观看在线播放| 免费黄色在线免费观看| 久久久国产精品麻豆| 有码 亚洲区| 天堂俺去俺来也www色官网| 国产午夜精品一二区理论片| 亚洲熟女精品中文字幕| 国产黄片视频在线免费观看| 一个人免费看片子| 交换朋友夫妻互换小说| 日日撸夜夜添| 狂野欧美激情性bbbbbb| 亚洲欧美日韩东京热| 国产精品人妻久久久影院| 少妇 在线观看| 亚洲性久久影院| 欧美日韩在线观看h| 国产美女午夜福利| 99视频精品全部免费 在线| 一级二级三级毛片免费看| 新久久久久国产一级毛片| av播播在线观看一区| 精品卡一卡二卡四卡免费| 久久这里有精品视频免费| 欧美日韩在线观看h| 老司机亚洲免费影院| 人妻制服诱惑在线中文字幕| 日本av免费视频播放| 春色校园在线视频观看| 在线精品无人区一区二区三| 97超视频在线观看视频| 国产 精品1| 亚洲在久久综合| 有码 亚洲区| av女优亚洲男人天堂| 在线精品无人区一区二区三| 这个男人来自地球电影免费观看 | 国产精品人妻久久久影院| 国模一区二区三区四区视频| 免费黄频网站在线观看国产| 亚洲av中文av极速乱| 色视频在线一区二区三区| 免费黄频网站在线观看国产| 观看av在线不卡| 亚洲av免费高清在线观看| 岛国毛片在线播放| 99热国产这里只有精品6| 国产在视频线精品| 免费黄频网站在线观看国产| 欧美精品人与动牲交sv欧美| 日本爱情动作片www.在线观看| 国产淫语在线视频| 五月伊人婷婷丁香| 免费观看无遮挡的男女| 人人妻人人看人人澡| 少妇被粗大猛烈的视频| 亚洲欧美一区二区三区国产| 国产一区二区三区综合在线观看 | 国产一区二区三区综合在线观看 | 夫妻性生交免费视频一级片| 亚洲欧洲精品一区二区精品久久久 | 乱系列少妇在线播放| 少妇裸体淫交视频免费看高清| 青青草视频在线视频观看| 国产伦精品一区二区三区四那| 一级毛片 在线播放| 国产69精品久久久久777片| 色网站视频免费| 国产日韩一区二区三区精品不卡 | 久久久久人妻精品一区果冻| 新久久久久国产一级毛片| 在线观看免费日韩欧美大片 | 99热这里只有是精品在线观看| 亚洲无线观看免费| 亚洲精品乱码久久久久久按摩| 在线观看免费日韩欧美大片 | 噜噜噜噜噜久久久久久91| 国产在线一区二区三区精| 日本黄大片高清| 91精品伊人久久大香线蕉| 久久久a久久爽久久v久久| 久久免费观看电影| 亚洲欧洲精品一区二区精品久久久 | 成年女人在线观看亚洲视频| av福利片在线| 久久久久精品性色| 亚洲av在线观看美女高潮| 肉色欧美久久久久久久蜜桃| 色视频在线一区二区三区| 日韩电影二区| 极品人妻少妇av视频| 亚洲国产欧美日韩在线播放 | 欧美精品亚洲一区二区| 亚洲无线观看免费| 九色成人免费人妻av| 美女视频免费永久观看网站| 男女啪啪激烈高潮av片| 亚洲成人手机| 在线播放无遮挡| 亚洲内射少妇av| 国产精品麻豆人妻色哟哟久久| 搡女人真爽免费视频火全软件| 91精品国产九色| 亚洲激情五月婷婷啪啪| 国语对白做爰xxxⅹ性视频网站| 少妇的逼水好多| 少妇被粗大猛烈的视频| 妹子高潮喷水视频| 大话2 男鬼变身卡| 久久6这里有精品| 观看免费一级毛片| av免费观看日本| 熟女人妻精品中文字幕| 成人毛片60女人毛片免费| 婷婷色综合www| 欧美日本中文国产一区发布| 偷拍熟女少妇极品色| 国产免费又黄又爽又色| 日本av免费视频播放| 最近的中文字幕免费完整| 久久99一区二区三区| 精品少妇内射三级| 看免费成人av毛片| 搡老乐熟女国产| 久久狼人影院| 国产熟女午夜一区二区三区 | 免费看日本二区| 国产午夜精品久久久久久一区二区三区| 国产精品99久久99久久久不卡 | 免费看不卡的av| 精品久久国产蜜桃| 中文字幕人妻熟人妻熟丝袜美| 欧美日韩精品成人综合77777| 青春草国产在线视频| 我的老师免费观看完整版| 中文字幕av电影在线播放| 日韩成人伦理影院| 免费久久久久久久精品成人欧美视频 | a级片在线免费高清观看视频| 大码成人一级视频| 中文在线观看免费www的网站| 如日韩欧美国产精品一区二区三区 | 亚洲av国产av综合av卡| 美女内射精品一级片tv| 国产在线免费精品| 国产av一区二区精品久久| 看免费成人av毛片| 一区二区av电影网| 久久久久久人妻| 观看av在线不卡| 成人二区视频| 免费高清在线观看视频在线观看| 亚洲精品国产av成人精品| 国产一区二区三区综合在线观看 | 日本免费在线观看一区| 国产成人免费无遮挡视频| 成人免费观看视频高清| 夫妻午夜视频| av在线老鸭窝| 最新中文字幕久久久久| 中文字幕免费在线视频6| av免费观看日本| 男人狂女人下面高潮的视频| 日本av免费视频播放| 韩国高清视频一区二区三区| 婷婷色麻豆天堂久久| 精品人妻熟女av久视频| 少妇的逼水好多| 精品国产乱码久久久久久小说| 亚洲国产精品成人久久小说| 亚洲成人av在线免费| 男女边吃奶边做爰视频| 91成人精品电影| 最新中文字幕久久久久| 久热久热在线精品观看| 少妇人妻一区二区三区视频| 亚洲成人一二三区av| 亚洲av日韩在线播放| 国产精品久久久久成人av| 国产成人freesex在线| 久久久久久久精品精品| 高清视频免费观看一区二区| 久久影院123| 亚洲精华国产精华液的使用体验| 午夜激情久久久久久久| 亚洲三级黄色毛片| 亚洲内射少妇av| 青春草视频在线免费观看| av又黄又爽大尺度在线免费看| tube8黄色片| 啦啦啦在线观看免费高清www| 国产淫语在线视频| 午夜久久久在线观看| 成人综合一区亚洲| 97超视频在线观看视频| 免费黄频网站在线观看国产| 日本与韩国留学比较| 一级毛片我不卡| 黑丝袜美女国产一区| 91精品国产国语对白视频| 精品酒店卫生间| 久久 成人 亚洲| 97超碰精品成人国产| 国内揄拍国产精品人妻在线| 国产69精品久久久久777片| 亚洲真实伦在线观看| 午夜老司机福利剧场| 两个人免费观看高清视频 | 中文字幕久久专区| 色婷婷av一区二区三区视频| 亚洲美女搞黄在线观看| 亚洲熟女精品中文字幕| 男人狂女人下面高潮的视频| 深夜a级毛片| 国产无遮挡羞羞视频在线观看| av播播在线观看一区| 另类精品久久| 欧美亚洲 丝袜 人妻 在线| 狂野欧美激情性bbbbbb| 一级,二级,三级黄色视频| 国产美女午夜福利| 日韩av在线免费看完整版不卡| 97精品久久久久久久久久精品| 22中文网久久字幕| 国产美女午夜福利| 亚洲精品自拍成人| 国产精品国产三级国产专区5o| 国产亚洲5aaaaa淫片| 丰满饥渴人妻一区二区三| 亚洲欧美精品专区久久| 国产一级毛片在线| 五月伊人婷婷丁香| 免费黄网站久久成人精品| 免费久久久久久久精品成人欧美视频 | 五月玫瑰六月丁香| 亚洲,欧美,日韩| 一本大道久久a久久精品| 国产精品偷伦视频观看了| 丁香六月天网| 精品一区在线观看国产| 国产男女超爽视频在线观看| 麻豆乱淫一区二区| 亚洲av男天堂| kizo精华| 日韩免费高清中文字幕av| 国产精品福利在线免费观看| tube8黄色片| 视频中文字幕在线观看| 制服丝袜香蕉在线| 欧美 日韩 精品 国产| 午夜激情久久久久久久| 大香蕉久久网| 狂野欧美白嫩少妇大欣赏| 精品国产一区二区久久| 欧美变态另类bdsm刘玥| 午夜91福利影院| 亚洲真实伦在线观看| 久久久国产欧美日韩av| 日本av手机在线免费观看| av天堂中文字幕网| a级毛色黄片| 免费不卡的大黄色大毛片视频在线观看| 中文资源天堂在线| 欧美 亚洲 国产 日韩一| 国产一区二区三区av在线| 午夜影院在线不卡| 国产高清有码在线观看视频| 国内少妇人妻偷人精品xxx网站| 自线自在国产av| 国产视频首页在线观看| 高清午夜精品一区二区三区| 免费看日本二区| 一边亲一边摸免费视频| 久久久久久久久久久免费av| 久久人人爽人人片av| 桃花免费在线播放| 久久av网站| 亚洲四区av| 少妇被粗大猛烈的视频| 赤兔流量卡办理| 亚洲综合精品二区| 久久国产乱子免费精品| 中文字幕久久专区| 美女福利国产在线| 最新的欧美精品一区二区| 丝袜喷水一区| 亚洲欧美日韩卡通动漫| 亚洲精品视频女| 涩涩av久久男人的天堂| 日韩成人av中文字幕在线观看| 天天操日日干夜夜撸| 老熟女久久久| 9色porny在线观看| 大话2 男鬼变身卡| 自线自在国产av| 丰满迷人的少妇在线观看| 寂寞人妻少妇视频99o| 99热这里只有精品一区| 国产精品一区二区在线不卡| 欧美xxⅹ黑人| 免费黄网站久久成人精品| 中国国产av一级| 亚洲av日韩在线播放| 蜜桃久久精品国产亚洲av| 久久影院123| 国产精品一区二区性色av| 丝瓜视频免费看黄片| 丰满迷人的少妇在线观看| 国产成人精品无人区| 国产在线一区二区三区精| 综合色丁香网| 日韩人妻高清精品专区| 免费黄频网站在线观看国产| 亚洲色图综合在线观看| 婷婷色麻豆天堂久久| 免费少妇av软件| 好男人视频免费观看在线| 欧美日韩视频精品一区|