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

    縣南溝流域農(nóng)業(yè)產(chǎn)業(yè)-資源系統(tǒng)耦合過程模型的建立與分析

    2017-09-12 07:42:09劉罕奇王繼軍
    水土保持研究 2017年5期
    關(guān)鍵詞:農(nóng)業(yè)產(chǎn)業(yè)流域耦合

    劉罕奇, 王繼軍

    (1.中國科學(xué)院 水利部 水土保持研究所, 陜西 楊凌 712100;2.西北農(nóng)林科技大學(xué), 陜西 楊凌 712100; 3.中國科學(xué)院大學(xué), 北京 100049)

    縣南溝流域農(nóng)業(yè)產(chǎn)業(yè)-資源系統(tǒng)耦合過程模型的建立與分析

    劉罕奇1,3, 王繼軍1,2

    (1.中國科學(xué)院 水利部 水土保持研究所, 陜西 楊凌 712100;2.西北農(nóng)林科技大學(xué), 陜西 楊凌 712100; 3.中國科學(xué)院大學(xué), 北京 100049)

    針對黃土丘陵區(qū)農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)演變軌跡發(fā)生變化的現(xiàn)實,參照農(nóng)業(yè)生態(tài)經(jīng)濟(jì)系統(tǒng)耦合過程模型,建立農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)耦合過程模型。并對模型參數(shù)(k1,k2,k3,k4)的變化過程及其代表的生態(tài)經(jīng)濟(jì)意義進(jìn)行重點討論,在此基礎(chǔ)上對縣南溝流域進(jìn)行了實證。結(jié)果表明:不同參數(shù)k值組合,形成了不同的農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)互動關(guān)系。退耕后,縣南溝流域資源所依賴的環(huán)境系統(tǒng)的提升空間現(xiàn)已達(dá)到極限,產(chǎn)業(yè)發(fā)展未建立在對資源的合理、有效利用上?,F(xiàn)階段資源系統(tǒng)發(fā)展態(tài)勢與資源系統(tǒng)演化速度和產(chǎn)業(yè)系統(tǒng)演化速度都成負(fù)相關(guān),產(chǎn)業(yè)系統(tǒng)發(fā)展態(tài)勢與資源系統(tǒng)演化速度和產(chǎn)業(yè)系統(tǒng)演化速度都成正相關(guān)。將碳匯產(chǎn)業(yè)融入系統(tǒng)后,可使資源系統(tǒng)演化速度提升2.56倍,若通過“輪牧+舍飼”養(yǎng)殖模式進(jìn)一步改良,可使資源系統(tǒng)演化速度提升4.66倍,使得資源得到合理、高效利用,加速產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化升級。

    農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng); 耦合過程模型; 碳匯; 半輪牧半舍飼; 縣南溝流域

    退耕還林還草工程的實施,改變了退耕區(qū)農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)演變的原軌跡,農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)發(fā)展未能與農(nóng)業(yè)資源系統(tǒng)相適應(yīng),區(qū)域農(nóng)業(yè)產(chǎn)業(yè)與資源(量)的一致性程度較弱,農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)與農(nóng)業(yè)資源系統(tǒng)發(fā)展處于相悖態(tài)勢[1-4],而系統(tǒng)的演替方向和速度直接影響區(qū)域的可持續(xù)發(fā)展。為解決系統(tǒng)發(fā)展的相悖態(tài)勢,協(xié)調(diào)農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)和農(nóng)業(yè)資源系統(tǒng)的關(guān)系,就要對農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)進(jìn)行優(yōu)化耦合,其關(guān)鍵在于弄清農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)的耦合過程和演變規(guī)律。耦合過程模型從定量的角度可以較好地闡明農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)與農(nóng)業(yè)資源系統(tǒng)的相互作用關(guān)系,判定農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)各階段的耦合態(tài)勢,從而為系統(tǒng)結(jié)構(gòu)優(yōu)化提供依據(jù)[5-6]。從現(xiàn)有研究來看,大部分學(xué)者主要對系統(tǒng)耦合模式、系統(tǒng)耦合態(tài)勢、系統(tǒng)耦合效應(yīng)方面進(jìn)行了分析,系統(tǒng)耦合過程模型研究相對滯后,且在已有的農(nóng)業(yè)生態(tài)經(jīng)濟(jì)耦合過程模型的研究中,未能夠?qū)⑥r(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)單獨考慮,而農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)作為農(nóng)業(yè)生態(tài)經(jīng)濟(jì)系統(tǒng)的核心,其演變過程更能夠直接地說明系統(tǒng)發(fā)展的耦合態(tài)勢[5,7-18]。為充分開發(fā)利用農(nóng)業(yè)資源,協(xié)調(diào)產(chǎn)業(yè)發(fā)展與資源利用的矛盾,必須分析農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)的耦合關(guān)系,建立農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)耦合過程模型,從而判定系統(tǒng)發(fā)展的耦合態(tài)勢,這對于黃土丘陵區(qū)可持續(xù)發(fā)展方案的制定具有重要意義。

    縣南溝流域是黃土高原丘陵區(qū)退耕還林還草工程實施的典型代表流域,退耕還林還草工程實施后,農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)的耦合關(guān)系發(fā)生顯著變化,林草資源量過剩,而依賴資源發(fā)展的產(chǎn)業(yè)規(guī)模過小且發(fā)展趨于穩(wěn)定。在此,通過應(yīng)用農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)與農(nóng)業(yè)資源系統(tǒng)耦合過程模型對縣南溝流域發(fā)展耦合關(guān)系進(jìn)行分析,旨在為協(xié)調(diào)好產(chǎn)業(yè)與資源的關(guān)系,為實現(xiàn)農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)良性循環(huán)提供依據(jù)。

    1 研究區(qū)概況與資料來源

    1.1 研究區(qū)概況

    縣南溝流域位于安塞縣沿河灣鎮(zhèn)(109°12′12″—109°22′12″E,36°41′24″—36°46′12″N),流域面積50.64 km2,屬典型暖溫帶干旱半干旱氣候,年平均氣溫8.8℃,年降水量500~550 mm,降雨年際差異大,且年內(nèi)分配不均,60%~80%降雨集中在7—9月。該流域包括5個行政村(磚窯溝、方家河、畔坡山、崖窯、寨子灣)和1個自然村(何塌),2015年有農(nóng)戶544戶,人口2 534人。自退耕以來,流域的土地利用結(jié)構(gòu)發(fā)生了顯著的變化,耕地∶林地∶草地面積比由1998年的1∶0.38∶5.48變?yōu)?015年的1∶7.3∶11.6。農(nóng)戶收入主要來自于種植業(yè)、果業(yè)、畜牧業(yè)、工副業(yè)等,種植業(yè)主要類型包括糧食作物、經(jīng)濟(jì)作物,人均純收入總體呈上升趨勢,由1999年的1 300多元/人增長至2015年的6 645元/人,農(nóng)業(yè)總收入為858萬元。

    1.2 資料來源

    1999年縣南溝流域成為退耕還林還草工程試點單元,開始實施退耕還林還草工程。因此,以1998年為基礎(chǔ)年,選擇1998—2015年為研究時段,探討退耕還林還草工程實施后農(nóng)業(yè)產(chǎn)業(yè)與農(nóng)業(yè)資源系統(tǒng)的變化情況。本研究所用到的數(shù)據(jù)來源于1998—2015年課題組的調(diào)研資料以及1998—2015年《安塞縣統(tǒng)計年鑒》。文中的圖表分別用Word2010,Excel 2010制作,數(shù)據(jù)用Excel 2010,MATLAB 7.0和SPSS 21.0處理。

    2 農(nóng)業(yè)產(chǎn)業(yè)-資源系統(tǒng)耦合過程模型設(shè)計與分析

    2.1 模型設(shè)計

    2.1.1 基礎(chǔ)模型設(shè)計 為使農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)及其互動過程定量化,借鑒農(nóng)業(yè)生態(tài)經(jīng)濟(jì)耦合過程模型[6],建立農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)耦合過程模型。用ai,ar分別表示農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)與農(nóng)業(yè)資源系統(tǒng),ai(t),ar(t)分別表示時間t時農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)與農(nóng)業(yè)資源系統(tǒng)總量,則可用dai(t)/dt、dar(t)/dt分別表示t時農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)與農(nóng)業(yè)資源系統(tǒng)的演化速度。按照貝塔蘭菲的一般系統(tǒng)理論[19],農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)演化方程的形式為:

    (1)

    (2)

    式中:k1,k2分別是農(nóng)業(yè)資源系統(tǒng)與農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)對資源系統(tǒng)演化過程的影響系數(shù),是資源系統(tǒng)與產(chǎn)業(yè)系統(tǒng)狀態(tài)表征量對資源系統(tǒng)演變速度的貢獻(xiàn)度;k3,k4分別為農(nóng)業(yè)資源系統(tǒng)與農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)對產(chǎn)業(yè)系統(tǒng)演化過程的影響系數(shù),是資源系統(tǒng)與產(chǎn)業(yè)系統(tǒng)狀態(tài)表征量對產(chǎn)業(yè)系統(tǒng)演變速度的貢獻(xiàn)度,它揭示了農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)與農(nóng)業(yè)資源系統(tǒng)的互動方向與強(qiáng)度,在特定區(qū)域k1,k2,k3,k4都是常數(shù)。通過已知的k1,k2,k3、k4計算出當(dāng)年的農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)與農(nóng)業(yè)資源系統(tǒng)綜合值,則可判定其耦合態(tài)勢。

    2.1.2 擴(kuò)展模型 由于農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)在它的演變過程會受到農(nóng)業(yè)系統(tǒng)外力的作用,從而改變其演變軌跡。因此要重新探討新背景下系統(tǒng)的演變軌跡,對耦合過程基礎(chǔ)模型進(jìn)行修正。以公式(1),(2)為基礎(chǔ),假設(shè)e1,e2分別為外力對農(nóng)業(yè)資源系統(tǒng)與農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)影響后導(dǎo)致其演變軌跡發(fā)生變化后的調(diào)整系數(shù),通過對農(nóng)業(yè)資源系統(tǒng)和農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)演變軌跡的修正,可形成公式(3),(4)。

    (3)

    (4)

    外力對農(nóng)業(yè)資源系統(tǒng)與農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)正向影響時,e1,e2取正,逆向影響時,e1,e2取負(fù)。e1,e2主要根據(jù)外力對系統(tǒng)的演變軌跡的偏離程度來確定,可用以下公式求得:

    (5)

    式中:xi(i=1,2,…,m)表示農(nóng)業(yè)資源系統(tǒng)表征量;yj

    (j=1,2,…,n)表示農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)表征量;pi,qj分別為對應(yīng)的權(quán)重;pixi,qjyj分別為外力引起資源系統(tǒng)與產(chǎn)業(yè)系統(tǒng)偏離后的增量。

    2.2 系統(tǒng)耦合過程分析

    農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)耦合過程模型由基礎(chǔ)模型和擴(kuò)展模型共同構(gòu)成。系統(tǒng)發(fā)展過程中,農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)與農(nóng)業(yè)資源系統(tǒng)的互動關(guān)系由系數(shù)k值直接反映。通過對k值正負(fù)的分析,表明相應(yīng)系統(tǒng)發(fā)展對演化速度的促進(jìn)或抑制作用,反映系統(tǒng)間的耦合態(tài)勢(表1)。

    為了更好地說明系統(tǒng)發(fā)展的耦合態(tài)勢,基于表1的分析,分別從資源系統(tǒng)發(fā)展、產(chǎn)業(yè)系統(tǒng)發(fā)展、農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)發(fā)展角度對系統(tǒng)耦合態(tài)勢進(jìn)行劃分。

    (1) 基于資源系統(tǒng)發(fā)展角度。在k1,k2都為正的情境下,無論是產(chǎn)業(yè)的發(fā)展還是資源的建設(shè),對資源系統(tǒng)的演化速度都起促進(jìn)作用。在此階段為滿足資源系統(tǒng)建設(shè)、修復(fù)的需求,首先考慮資源的發(fā)展,可適當(dāng)降低經(jīng)濟(jì)發(fā)展速度,即產(chǎn)業(yè)的發(fā)展和資源的發(fā)展始終圍繞著資源所依賴的生態(tài)系統(tǒng)的修復(fù)、建設(shè)、改良。

    表1 k1,k2,k3,k4值分析

    注:k值的正負(fù)用“+/-”表示,“-”表示該種耦合態(tài)勢基本不存在。

    在k1>0,k2<0的情境下,資源的開發(fā)、利用對資源所依賴的環(huán)境以及資源的演變依然起促進(jìn)作用,資源還有進(jìn)一步增長的空間。若|k1|>|k2|,產(chǎn)業(yè)的發(fā)展對資源已經(jīng)起到了一個局部破壞作用,建立在資源基礎(chǔ)上的產(chǎn)業(yè)對資源的開發(fā)利用沒有超越度;若|k1|<|k2|,產(chǎn)業(yè)的發(fā)展已經(jīng)超越了資源的發(fā)展速度,這種情況下產(chǎn)業(yè)的發(fā)展就對資源是一種破壞。整體上來說,產(chǎn)業(yè)在這個階段對資源是一種局部破壞,處于可恢復(fù)階段。

    在k1<0的情境下,資源發(fā)展的空間已經(jīng)達(dá)到了極值,資源量趨于穩(wěn)定,資源發(fā)展對資源系統(tǒng)的演化速度開始起抑制作用,產(chǎn)業(yè)發(fā)展應(yīng)圍繞對資源的開發(fā)、利用調(diào)整發(fā)展方向和規(guī)模。

    (2) 基于產(chǎn)業(yè)系統(tǒng)發(fā)展角度。在k4>0的情境下,產(chǎn)業(yè)的發(fā)展對產(chǎn)業(yè)系統(tǒng)的演化速度起促進(jìn)作用。若k3>0,產(chǎn)業(yè)系統(tǒng)與資源系統(tǒng)相互支撐,產(chǎn)業(yè)發(fā)展建立在對資源的開發(fā)、利用上;若k3<0,資源支撐產(chǎn)業(yè)的發(fā)展較弱,產(chǎn)業(yè)系統(tǒng)的發(fā)展方向和規(guī)模需要進(jìn)一步調(diào)整。此階段以產(chǎn)業(yè)發(fā)展為主,保證在資源可利用閾值范圍內(nèi),使得產(chǎn)業(yè)系統(tǒng)發(fā)展盡快進(jìn)入成熟期。在k4<0的情境下,產(chǎn)業(yè)的發(fā)展對產(chǎn)業(yè)系統(tǒng)演化速度起抑制作用。產(chǎn)業(yè)間競爭較大,產(chǎn)業(yè)系統(tǒng)需根據(jù)產(chǎn)業(yè)與產(chǎn)業(yè)、產(chǎn)業(yè)與資源的互動關(guān)系調(diào)整發(fā)展方向和規(guī)模。

    (3) 基于產(chǎn)業(yè)系統(tǒng)與資源系統(tǒng)的耦合角度。產(chǎn)業(yè)系統(tǒng)與資源系統(tǒng)的互動關(guān)系由k2,k3反映,當(dāng)都是正的時,互相促進(jìn),當(dāng)都是負(fù)的時,互相抑制。除此之外,根據(jù)農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)與農(nóng)業(yè)資源系統(tǒng)各自的發(fā)展情況并結(jié)合產(chǎn)業(yè)與資源的結(jié)合程度,出現(xiàn)兩種情況:① 資源系統(tǒng)對產(chǎn)業(yè)系統(tǒng)演化速度起抑制作用,產(chǎn)業(yè)系統(tǒng)的發(fā)展對資源系統(tǒng)演化速度起到促進(jìn)作用。在此發(fā)展態(tài)勢下,資源系統(tǒng)與產(chǎn)業(yè)系統(tǒng)發(fā)展類型和方向已經(jīng)確定,但與資源利用相關(guān)的產(chǎn)業(yè)發(fā)展剛起步或規(guī)模較小。② 資源系統(tǒng)對產(chǎn)業(yè)系統(tǒng)演化速度起促進(jìn)作用,產(chǎn)業(yè)系統(tǒng)的發(fā)展對資源系統(tǒng)演化速度起到抑制作用。在此發(fā)展態(tài)勢下,資源利用型產(chǎn)業(yè)得到發(fā)展,而產(chǎn)業(yè)發(fā)展造成資源鏈網(wǎng)的局部破壞,使得資源系統(tǒng)發(fā)展受到抑制。無論是在資源系統(tǒng)建設(shè)初期產(chǎn)業(yè)系統(tǒng)對資源系統(tǒng)的過度利用,還是資源系統(tǒng)發(fā)展成熟期產(chǎn)業(yè)系統(tǒng)對資源系統(tǒng)開發(fā)利用不足的情況,都需調(diào)節(jié)資源系統(tǒng)與產(chǎn)業(yè)系統(tǒng)的互動關(guān)系,使資源利用發(fā)展的產(chǎn)業(yè)不斷成長。

    3 縣南溝流域耦合過程分析

    縣南溝流域自退耕以來,流域生態(tài)林面積由退耕前的149.6 hm2增加到2015年的1 370.8 hm2,碳匯潛在經(jīng)濟(jì)價值達(dá)到476.3萬元,是退耕前的8.5倍,碳匯發(fā)展?jié)摿薮?。考慮到碳匯產(chǎn)業(yè)可能會是流域產(chǎn)業(yè)鏈網(wǎng)結(jié)構(gòu)優(yōu)化的路徑,本研究將引入碳匯產(chǎn)業(yè),對流域農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)耦合過程進(jìn)行綜合評估。

    3.1.1 權(quán)重確定 各項指標(biāo)在綜合評價指標(biāo)體系中的重要程度用權(quán)重來表示。為了使其具有可比性、更加符合農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)演變與實際情況,通過德爾菲法(Dalphi)確定各個指標(biāo)權(quán)重賦值,其值介于0—1之間(表2)[5,20]。

    表2 縣南溝流域農(nóng)業(yè)產(chǎn)業(yè)-資源系統(tǒng)綜合評價指標(biāo)體系與權(quán)重

    3.1.2 數(shù)據(jù)標(biāo)準(zhǔn)化 為了消除數(shù)據(jù)的數(shù)量級以及量綱的不同而造成的影響,對數(shù)據(jù)用極差標(biāo)準(zhǔn)化的方法進(jìn)行標(biāo)準(zhǔn)化處理[14]。通過公式(4),(5) 得出ar,ai的值(表3)。為更好地體現(xiàn)系統(tǒng)發(fā)展過程,在對農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)和農(nóng)業(yè)資源系統(tǒng)綜合指數(shù)計算的基礎(chǔ)上,對其進(jìn)行非線性擬合(圖1)。

    圖11998-2015年縣南溝流域農(nóng)業(yè)產(chǎn)業(yè)-資源

    系統(tǒng)綜合指數(shù)變化情況

    由圖1可以看出,在1998—2015年,縣南溝流域農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)和農(nóng)業(yè)資源系統(tǒng)綜合指數(shù)發(fā)展曲線呈S型,農(nóng)業(yè)資源系統(tǒng)發(fā)展逐漸趨于穩(wěn)定,農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)發(fā)展處于上升態(tài)勢。總體上,縣南溝流域農(nóng)業(yè)產(chǎn)業(yè)與資源狀況明顯改善,產(chǎn)業(yè)系統(tǒng)與資源系統(tǒng)綜合指數(shù)明顯提高。2015年資源系統(tǒng)綜合指數(shù)較1998年升高了23.21%,產(chǎn)業(yè)系統(tǒng)綜合指數(shù)升高了88.20%。1999—2003年由于退耕還林工程的實施,農(nóng)業(yè)資源系統(tǒng)變化劇烈,資源系統(tǒng)綜合指數(shù)明顯提高,年平均遞增7.57%,隨后生態(tài)環(huán)境重建基本完成,資源系統(tǒng)發(fā)展態(tài)勢平穩(wěn);2000—2007年農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)受農(nóng)業(yè)資源系統(tǒng)影響,產(chǎn)業(yè)系統(tǒng)綜合指數(shù)呈上升發(fā)展態(tài)勢,年平均遞增9.98%,隨后農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)與農(nóng)業(yè)資源系統(tǒng)相互調(diào)節(jié),產(chǎn)業(yè)系統(tǒng)發(fā)展緩慢上升。

    3.2 k1,k2,k3,k4的確定及分析

    根據(jù)縣南溝流域1998—2015年各年份的ar(t),ai(t)值及各自擬合曲線,得出Var,Vai(表3),按照公式(1),(2)回歸得出k1,k2,k3,k4。

    表3 縣南溝流域1998-2015年農(nóng)業(yè)產(chǎn)業(yè)-資源系統(tǒng)互動過程相關(guān)指標(biāo)值

    根據(jù)縣南溝流域1998-2015年農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)和農(nóng)業(yè)資源系統(tǒng)的演變過程,并結(jié)合產(chǎn)業(yè)系統(tǒng)與資源系統(tǒng)的耦合關(guān)系分析,可將農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)的發(fā)展劃分為三個階段[21]:資源修復(fù)階段(1999—2003)、產(chǎn)業(yè)適應(yīng)階段(2004—2010)、系統(tǒng)相悖階段(2011—2015年),分別對k值進(jìn)行計算(表4),對各階段耦合態(tài)勢進(jìn)行判定。

    表4 k1,k2,k3,k4及其標(biāo)準(zhǔn)誤

    由表4可以看出,k1,k2,k3,k4的值在系統(tǒng)不同發(fā)展階段取值有所差異。1998—2003年流域?qū)嵤┩烁M(jìn)行資源修復(fù)過程,此階段以資源系統(tǒng)所依賴的環(huán)境建設(shè)為主,基于資源利用的產(chǎn)業(yè)發(fā)展受到抑制(k1>0,k3<0)。2004—2010年資源系統(tǒng)從低層次的水平經(jīng)過生態(tài)建設(shè)不斷的提升,k1值有所上升,產(chǎn)業(yè)逐漸與資源相適應(yīng),產(chǎn)業(yè)鏈與資源相關(guān)度有所提升,產(chǎn)業(yè)的發(fā)展加大了對有限資源的利用,造成資源鏈網(wǎng)結(jié)局部破壞,從而對資源系統(tǒng)的演化起負(fù)面作用(k2<0);與此同時,產(chǎn)業(yè)發(fā)展脫離對現(xiàn)有資源利用的局限,各產(chǎn)業(yè)在新環(huán)境下逐漸趨于穩(wěn)定,部分產(chǎn)業(yè)為爭奪資源的利用而出現(xiàn)競爭,各產(chǎn)業(yè)間發(fā)展速度受到抑制(k4<0)。2011—2015年資源系統(tǒng)在現(xiàn)有土地利用結(jié)構(gòu)上已經(jīng)處于穩(wěn)定發(fā)展階段,而現(xiàn)有資源結(jié)構(gòu)未達(dá)到優(yōu)化閾值且產(chǎn)業(yè)發(fā)展未建立在對資源的合理、有效利用上,由于沒有得到及時有效的調(diào)整,產(chǎn)業(yè)與資源的發(fā)展都受到了抑制(k1<0,k3<0);而在系統(tǒng)融入碳匯產(chǎn)業(yè)之后,產(chǎn)業(yè)鏈與資源相關(guān)度提升,對產(chǎn)業(yè)演化速度起促進(jìn)作用(k4>0)。從農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)演化過程的整體來看,1998—2003年的資源修復(fù)過程使得產(chǎn)業(yè)系統(tǒng)發(fā)展受阻,但是其修復(fù)的資源為第二階段的產(chǎn)業(yè)發(fā)展奠定了良好的基礎(chǔ);2004—2010年經(jīng)濟(jì)適應(yīng)階段,產(chǎn)業(yè)結(jié)構(gòu)根據(jù)資源發(fā)展?fàn)顩r進(jìn)行調(diào)整,使產(chǎn)業(yè)發(fā)展與資源建設(shè)相適應(yīng);2011—2015年是對第二階段發(fā)展的穩(wěn)固,但仍需對產(chǎn)業(yè)結(jié)構(gòu)進(jìn)行調(diào)整,使得產(chǎn)業(yè)系統(tǒng)與資源系統(tǒng)協(xié)同發(fā)展,現(xiàn)階段耦合公式(3),(4)變?yōu)椋?/p>

    (6)

    (7)

    3.3 優(yōu)化情境下e1,e2的確定及分析

    對于縣南溝流域,盡管已經(jīng)計算了碳匯的潛在價值,但必定碳匯產(chǎn)業(yè)尚未成為現(xiàn)實,在上述分析過程中,并未考慮碳匯產(chǎn)業(yè)發(fā)展與其他相關(guān)產(chǎn)業(yè)或要素的關(guān)系。在此,將碳匯產(chǎn)業(yè)融入系統(tǒng)發(fā)展,從量變角度探討各產(chǎn)業(yè)之間、產(chǎn)業(yè)與資源的結(jié)構(gòu)變化,分析農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)耦合發(fā)展態(tài)勢。

    (1) 原系統(tǒng)融入碳匯產(chǎn)業(yè)情景。從系統(tǒng)各要素考慮,發(fā)展碳匯產(chǎn)業(yè)將會提高人均收入(y8)并增加產(chǎn)業(yè)鏈與資源量相關(guān)度(y1)。資源系統(tǒng)基本不受影響,產(chǎn)業(yè)系統(tǒng)改變了原演變軌跡,e2大于0。公式(5)變?yōu)椋?/p>

    (8)

    在融入碳匯產(chǎn)業(yè)后,縣南溝流域經(jīng)濟(jì)純收入可增加1 243.6萬元,人均純收入增加4 924.4元,根據(jù)表1指標(biāo),農(nóng)業(yè)產(chǎn)業(yè)鏈與資源量相關(guān)度將增加1個標(biāo)度,前者在上述分析時已考慮,后者需要補充修訂。由計算得到e2=0.1294,現(xiàn)階段公式(6),(7)為:

    (9)

    (10)

    此時,產(chǎn)業(yè)發(fā)展對產(chǎn)業(yè)系統(tǒng)演化起到促進(jìn)作用,同時也使資源系統(tǒng)演化速度提升為:Var=0.0238,是融入碳匯產(chǎn)業(yè)前的3.56倍。可以看出,碳匯產(chǎn)業(yè)的融入可以改善農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)和農(nóng)業(yè)資源系統(tǒng)的相悖態(tài)勢,對于現(xiàn)階段農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)的發(fā)展起到了促進(jìn)作用,資源利用更加合理有效,加速了產(chǎn)業(yè)優(yōu)化升級,使農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)和農(nóng)業(yè)資源系統(tǒng)的關(guān)系更加緊密。

    (2) 碳匯產(chǎn)業(yè)+畜牧業(yè)“輪牧+舍飼”情景。退耕后,林草資源雖還沒有達(dá)到最大化,但受禁牧政策的影響使得現(xiàn)有林草資源沒有得到合理利用。在融入碳匯產(chǎn)業(yè)的基礎(chǔ)上,針對縣南溝流域退耕政策實施使得流域林草資源富足的現(xiàn)狀,采納“輪牧+舍飼”的養(yǎng)殖政策,在保持現(xiàn)有生態(tài)環(huán)境的基礎(chǔ)上,使得林草資源得到有序、合理利用,產(chǎn)業(yè)發(fā)展得到優(yōu)化升級。此時系統(tǒng)演變軌跡發(fā)生了變化,其中產(chǎn)業(yè)系統(tǒng)中變化最大的指標(biāo)是農(nóng)業(yè)產(chǎn)業(yè)鏈與資源量相關(guān)度(y1)、農(nóng)產(chǎn)品商品率(y6)、人均純收入(y8),這些指標(biāo)都會提高,即產(chǎn)業(yè)系統(tǒng)改變了原演變軌跡,e2大于0。公式(5)變?yōu)椋?/p>

    (11)

    縣南溝流域2015年有草地2 420.43 hm2,林地1 589.27 hm2,其中生態(tài)林1 370.8 hm2。按照“輪牧+舍飼”管理方式,冬春舍飼,夏秋輪牧,極限放牧量按1 000只算(歷史最大1 070只),為使林草資源穩(wěn)固、持續(xù)發(fā)展,假設(shè)600只羊可以實現(xiàn)系統(tǒng)良性互動。每只羊按純收入為431元,可對流域增加收入25.86萬元,人均純收入增加102.05元。根據(jù)表1指標(biāo),農(nóng)業(yè)產(chǎn)業(yè)鏈與資源量相關(guān)度、農(nóng)產(chǎn)品商品率都將增加1個標(biāo)度,人均純收入由于變化不大所以保持不變。通過計算得到e2=0.0562,公式(6),(7)變?yōu)椋?/p>

    (12)

    (13)

    此時,產(chǎn)業(yè)發(fā)展對產(chǎn)業(yè)系統(tǒng)演化進(jìn)一步促進(jìn),資源系統(tǒng)演化速度提升為:Var=0.0312,是產(chǎn)業(yè)系統(tǒng)改良前的4.66倍,是單純?nèi)谌胩紖R產(chǎn)業(yè)的1.31倍??梢钥闯?,在融入碳匯產(chǎn)業(yè)并結(jié)合畜牧業(yè)“輪牧+舍飼”制度的實施,不僅產(chǎn)業(yè)系統(tǒng)結(jié)構(gòu)得到優(yōu)化改良,也使資源系統(tǒng)得到較好的開發(fā)利用,有效地促進(jìn)了產(chǎn)業(yè)系統(tǒng)與資源系統(tǒng)的協(xié)同發(fā)展。

    綜合來看,現(xiàn)有資源利用模式與產(chǎn)業(yè)態(tài)勢已形成了農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)相悖態(tài)勢,無論是對農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)單純地融入碳匯產(chǎn)業(yè)還是碳匯產(chǎn)業(yè)結(jié)合“輪牧+舍飼”養(yǎng)殖模式的優(yōu)化方案,都有利于產(chǎn)業(yè)系統(tǒng)發(fā)展并促進(jìn)資源系統(tǒng)的演化。

    4 結(jié)論與討論

    通過構(gòu)建農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)耦合過程模型,使得系統(tǒng)間耦合過程得以量化,通過對k值可能組合方式下所表征的生態(tài)經(jīng)濟(jì)意義的分析,明確了農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)和農(nóng)業(yè)資源系統(tǒng)的互動過程,為產(chǎn)業(yè)系統(tǒng)與資源系統(tǒng)優(yōu)化發(fā)展方向和規(guī)模提供依據(jù)。

    退耕后,縣南溝流域農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)經(jīng)過三個發(fā)展階段,資源所依賴的環(huán)境系統(tǒng)的提升空間已達(dá)到極限,產(chǎn)業(yè)發(fā)展未建立在對資源的合理、有效利用上。現(xiàn)階段農(nóng)業(yè)資源系統(tǒng)發(fā)展態(tài)勢與資源系統(tǒng)演化速度和產(chǎn)業(yè)系統(tǒng)演化速度都成負(fù)向相關(guān),農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)發(fā)展態(tài)勢與資源系統(tǒng)演化速度和產(chǎn)業(yè)系統(tǒng)演化速度都成正向相關(guān)。通過對縣南溝流域耦合過程優(yōu)化情景的模擬,得出通過融入碳匯產(chǎn)業(yè)可使農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)發(fā)展和農(nóng)業(yè)資源系統(tǒng)發(fā)展相互促進(jìn),資源系統(tǒng)演化速度提升2.56倍,且若通過“輪牧+舍飼”養(yǎng)殖模式進(jìn)一步優(yōu)化,可使農(nóng)業(yè)產(chǎn)業(yè)系統(tǒng)發(fā)展速度在融入碳匯產(chǎn)業(yè)的基礎(chǔ)上提升0.31倍,提升了產(chǎn)業(yè)發(fā)展與資源利用的相關(guān)度,擴(kuò)大了產(chǎn)業(yè)發(fā)展的空間。

    本研究所建農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)耦合過程模型比較成功地解釋了縣南溝流域的發(fā)展?fàn)顩r,但仍需在今后的實踐和應(yīng)用中對其進(jìn)一步完善,并對模型參數(shù)不斷地驗證和修訂,使模型更符合黃土丘陵區(qū)退耕區(qū)域的發(fā)展?fàn)顩r。通過模型對縣南溝流域的分析,證實退耕還林還草工程實施過程及方案有待于根據(jù)黃土丘陵區(qū)實際發(fā)展?fàn)顩r進(jìn)行優(yōu)化,促進(jìn)碳匯產(chǎn)業(yè)的發(fā)展,使農(nóng)業(yè)資源得到合理、高效利用,加速產(chǎn)業(yè)優(yōu)化升級,實現(xiàn)農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)的可持續(xù)發(fā)展。

    [1] 王繼軍,謝永生,盧宗凡,等.退耕還林還草下生態(tài)農(nóng)業(yè)發(fā)展模式初探[J].水土保持學(xué)報,2004,18(1):134-137.

    [2] 盧宗凡,梁一民,劉國彬.中國黃土高原生態(tài)農(nóng)業(yè)[M].西安:陜西科學(xué)技術(shù)出版社,1997.

    [3] 王繼軍.黃土丘陵區(qū)紙坊溝流域農(nóng)業(yè)生態(tài)經(jīng)濟(jì)系統(tǒng)耦合過程分析[J].應(yīng)用生態(tài)學(xué)報,2009,20(11):2723-2729.

    [4] 劉耀彬,陳斐,李仁東.區(qū)域城市化與生態(tài)環(huán)境耦合發(fā)展模擬及調(diào)控策略:以江蘇省為例[J].地理研究,2007,26(1):187-196.

    [5] 王繼軍,姜志德,連坡,等.70年來陜西省紙坊溝流域農(nóng)業(yè)生態(tài)經(jīng)濟(jì)系統(tǒng)耦合態(tài)勢[J].生態(tài)學(xué)報,2009,29(9):5130-5137.

    [6] 王繼軍,郭滿才,姜志德,等.農(nóng)業(yè)生態(tài)經(jīng)濟(jì)系統(tǒng)耦合過程模型的建立及應(yīng)用[J].生態(tài)學(xué)報,2010,30(9):2371-2378.

    [7] 劉興元,王鎖民,郭正剛.半干旱地區(qū)農(nóng)業(yè)資源的復(fù)合經(jīng)營模式及生態(tài)經(jīng)濟(jì)耦合效應(yīng)研究[J].自然資源學(xué)報,2004,19(5):624-631.

    [8] 周立華,樊勝岳,王濤.黑河流域生態(tài)經(jīng)濟(jì)系統(tǒng)分析與耦合發(fā)展模式[J].干旱區(qū)資源與環(huán)境,2005,19(5):67-72.

    [9] 周立華,王濤,樊勝岳,等.內(nèi)陸河流域的生態(tài)經(jīng)濟(jì)問題與協(xié)調(diào)發(fā)展模式:以黑河流域為例[J].中國軟科學(xué),2005(1):114-119.

    [10] 曹洪華.生態(tài)文明視角下流域生態(tài)—經(jīng)濟(jì)系統(tǒng)耦合模式研究:以洱海流域為例[D].長春:東北師范大學(xué),2014.

    [11] 邵權(quán)熙.當(dāng)代中國林業(yè)生態(tài)經(jīng)濟(jì)社會耦合系統(tǒng)及耦合模式研究[D].北京:北京林業(yè)大學(xué),2008.

    [12] 李崇明,丁烈云.小城鎮(zhèn)資源環(huán)境與社會經(jīng)濟(jì)協(xié)調(diào)發(fā)展評價模型及應(yīng)用研究[J].系統(tǒng)工程理論與實踐,2004,24(11):134-139.

    [13] 梁紅梅,劉衛(wèi)東,林育欣,等.土地利用效益的耦合模型及其應(yīng)用[J].浙江大學(xué)學(xué)報:農(nóng)業(yè)與生命科學(xué)版,2008,34(2):230-236.

    [14] 喬標(biāo),方創(chuàng)琳.城市化與生態(tài)環(huán)境協(xié)調(diào)發(fā)展的動態(tài)耦合模型及其在干旱區(qū)的應(yīng)用[J].生態(tài)學(xué)報,2005,25(11):3003-3009.

    [15] 劉佳,王繼軍.黃土丘陵區(qū)紙坊溝流域農(nóng)業(yè)生態(tài)經(jīng)濟(jì)系統(tǒng)耦合態(tài)勢[J].應(yīng)用生態(tài)學(xué)報,2010,21(6):1511-1517.

    [16] 董孝斌,高旺盛,嚴(yán)茂超.基于能值理論的農(nóng)牧交錯帶兩個典型縣域生態(tài)經(jīng)濟(jì)系統(tǒng)的耦合效應(yīng)分析[J].農(nóng)業(yè)工程學(xué)報,2005,21(11):1-6.

    [17] 盧良恕,梅方權(quán),許世衛(wèi).中國解決溫飽問題的成就及其發(fā)展戰(zhàn)略[J].農(nóng)業(yè)現(xiàn)代化研究,1990,11(4):2-25.

    [18] 夏自蘭,王繼軍.基于水土保持下的紙坊溝流域農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)耦合效應(yīng)評價[J].自然資源學(xué)報,2011,26(10):1647-1657.

    [19] Bertalanffy L V. Clean coal technologies: A status report[J]. Electrical World, 1992(2):37-42.

    [20] 李芬,王繼軍.黃土丘陵區(qū)紙坊溝流域近70年農(nóng)業(yè)生態(tài)安全評價[J].生態(tài)學(xué)報,2008,28(5):2380-2388.

    [21] 劉罕奇,王繼軍.黃土丘陵區(qū)縣南溝流域農(nóng)業(yè)產(chǎn)業(yè)—資源系統(tǒng)鏈網(wǎng)結(jié)構(gòu)分析[J].水土保持研究,2017,24(2):207-212.

    EstablishmentandAnalysisofAgriculturalIndustry-ResourceSystemCouplingProcessModelinXiannangouWatershed

    LIU Hanqi1,3, WANG Jijun1,2

    (1.InstituteofSoilandWaterConversation,ChineseAcademyofSciencesandMinistryofWaterResources,Yangling,Shaanxi712100,China; 2.NorthwestA&FUniversity,Yangling,Shaanxi712100,China; 3.UniversityofChineseAcademyofSciences,Beijing100049,China)

    In view of the development of agricultural industry-resource system in Loess Hilly region, we established a coupling process model of agricultural industry-resource system following the agricultural ecological-economic coupling process model, and we analyzed the coupling situation of system. We focused on the ecological and economic significance of model parameters (k1,k2,k3,k4) and conducted an empirical study on the Xiannangou watershed. The results show that different parameter combination forms the different agricultural industry-resource system interactive relationship. After returning farmlands to the natural lands, the space of resources that relied on environmental systems has reached the limit now, and the development of industry is not based on the rational and efficient use of resources in Xiannangou watershed. The development trend of resource system is negatively related to the evolution speed of the resource system and the industrial system. The development trend of industry system is positively related to the evolution speed of the resource system and the industrial system. The account in of carbon sink can make the optimized system evolution speed increase by 2.56 times. On the basis of that, the agricultural system can be promoted by 4.66 times through rotational grazing and house feeding mode, which makes the resources be rationally and efficiently used and accelerates the optimization and upgrade of industrial structure.

    agricultural industry-resource system; coupling process model; carbon sequestration; rotational grazing and house feeding; Xiannangou watershed

    2017-02-05

    :2017-03-16

    國家自然科學(xué)基金(41571515);國家重點研發(fā)計劃課題(2016YFC0503702);“十三五”國家重點研發(fā)計劃項目(2016YFC0501707)

    劉罕奇(1990—),男,山東淄博人,碩士研究生,研究方向:水土保持監(jiān)測與評價。E-mail:liuhq@nwsuaf.edu.cn

    王繼軍(1964—),男,陜西渭南人,研究員,研究方向:生態(tài)經(jīng)濟(jì)。E-mail:jjwang@ms.iswc.ac.cn

    F062.2

    :A

    :1005-3409(2017)05-0259-07

    猜你喜歡
    農(nóng)業(yè)產(chǎn)業(yè)流域耦合
    壓油溝小流域
    非Lipschitz條件下超前帶跳倒向耦合隨機(jī)微分方程的Wong-Zakai逼近
    堡子溝流域綜合治理
    羅堰小流域
    以“黨建”推進(jìn)特色農(nóng)業(yè)產(chǎn)業(yè)發(fā)展
    甘肅省農(nóng)業(yè)產(chǎn)業(yè)聯(lián)合會會員單位
    打造智慧流域的思路及構(gòu)想——以討賴河流域為例
    甘肅省農(nóng)業(yè)產(chǎn)業(yè)聯(lián)合會會員單位
    農(nóng)業(yè)產(chǎn)業(yè)的改革與創(chuàng)新
    基于“殼-固”耦合方法模擬焊接裝配
    大型鑄鍛件(2015年5期)2015-12-16 11:43:20
    成人特级黄色片久久久久久久| 麻豆国产97在线/欧美 | 超碰成人久久| 国产亚洲欧美98| а√天堂www在线а√下载| 亚洲精品一区av在线观看| 久久久久免费精品人妻一区二区| 婷婷丁香在线五月| 色播亚洲综合网| 18禁美女被吸乳视频| 一进一出抽搐gif免费好疼| 日日摸夜夜添夜夜添小说| 国产69精品久久久久777片 | 亚洲精品色激情综合| 12—13女人毛片做爰片一| 中文资源天堂在线| 欧洲精品卡2卡3卡4卡5卡区| 日本一二三区视频观看| 女生性感内裤真人,穿戴方法视频| 国产av在哪里看| 亚洲av五月六月丁香网| 国产精品一区二区免费欧美| 18禁裸乳无遮挡免费网站照片| 亚洲18禁久久av| 国产三级中文精品| 欧美日韩亚洲国产一区二区在线观看| 欧美午夜高清在线| 777久久人妻少妇嫩草av网站| 久久久精品大字幕| 欧美3d第一页| 十八禁人妻一区二区| www.www免费av| 天天躁狠狠躁夜夜躁狠狠躁| 黑人欧美特级aaaaaa片| 亚洲国产欧美一区二区综合| 国产蜜桃级精品一区二区三区| 少妇被粗大的猛进出69影院| АⅤ资源中文在线天堂| 日韩欧美免费精品| a在线观看视频网站| 十八禁人妻一区二区| 午夜老司机福利片| 一本精品99久久精品77| 亚洲国产日韩欧美精品在线观看 | 在线观看免费日韩欧美大片| 亚洲色图 男人天堂 中文字幕| 日韩av在线大香蕉| 麻豆国产97在线/欧美 | 日本精品一区二区三区蜜桃| 亚洲欧美日韩无卡精品| 亚洲五月天丁香| 久久亚洲精品不卡| 悠悠久久av| 亚洲一区二区三区不卡视频| 成人亚洲精品av一区二区| 毛片女人毛片| 久久久久久人人人人人| 精品一区二区三区av网在线观看| 三级毛片av免费| 99国产精品一区二区蜜桃av| 五月玫瑰六月丁香| 99国产极品粉嫩在线观看| 中文字幕人成人乱码亚洲影| 十八禁网站免费在线| 欧美日韩亚洲综合一区二区三区_| 黄片大片在线免费观看| 三级毛片av免费| 亚洲一区高清亚洲精品| 色综合欧美亚洲国产小说| 女人被狂操c到高潮| 国内久久婷婷六月综合欲色啪| 最新在线观看一区二区三区| 好男人在线观看高清免费视频| 很黄的视频免费| 中出人妻视频一区二区| 精品久久久久久成人av| 亚洲第一电影网av| 九色国产91popny在线| 熟女少妇亚洲综合色aaa.| 国产成人欧美在线观看| 中文字幕精品亚洲无线码一区| 欧美色视频一区免费| 一进一出抽搐动态| 国内精品久久久久久久电影| 国产熟女午夜一区二区三区| 18美女黄网站色大片免费观看| 岛国在线免费视频观看| 叶爱在线成人免费视频播放| 国产高清视频在线观看网站| 欧美黑人精品巨大| 巨乳人妻的诱惑在线观看| 国产精品久久久久久亚洲av鲁大| 精品电影一区二区在线| 色哟哟哟哟哟哟| 国产亚洲av高清不卡| 国产精品一区二区精品视频观看| 国产免费男女视频| 国产激情偷乱视频一区二区| 夜夜躁狠狠躁天天躁| 午夜久久久久精精品| 老司机午夜福利在线观看视频| 99热6这里只有精品| 成人手机av| 俺也久久电影网| 1024视频免费在线观看| 18禁美女被吸乳视频| 听说在线观看完整版免费高清| 九色成人免费人妻av| 1024香蕉在线观看| 免费av毛片视频| 亚洲精品久久成人aⅴ小说| 天天躁狠狠躁夜夜躁狠狠躁| 午夜免费成人在线视频| 久久精品人妻少妇| 日韩欧美国产一区二区入口| 精品电影一区二区在线| 少妇粗大呻吟视频| 成人av一区二区三区在线看| 久久久久久九九精品二区国产 | 黄色视频不卡| 五月玫瑰六月丁香| 五月玫瑰六月丁香| 99精品欧美一区二区三区四区| 亚洲,欧美精品.| 国产精品一及| 久久99热这里只有精品18| 婷婷六月久久综合丁香| 亚洲片人在线观看| 在线永久观看黄色视频| 国产av一区二区精品久久| 国产91精品成人一区二区三区| 日韩精品免费视频一区二区三区| 亚洲av美国av| 日本精品一区二区三区蜜桃| 99精品欧美一区二区三区四区| 国产私拍福利视频在线观看| 香蕉av资源在线| 91麻豆精品激情在线观看国产| 国产欧美日韩一区二区三| 他把我摸到了高潮在线观看| 午夜免费观看网址| 亚洲成人国产一区在线观看| 1024香蕉在线观看| 精品乱码久久久久久99久播| 国产精品一及| 宅男免费午夜| 99精品欧美一区二区三区四区| 亚洲男人天堂网一区| 天堂动漫精品| 久久精品影院6| www.999成人在线观看| 亚洲18禁久久av| 亚洲av成人av| 神马国产精品三级电影在线观看 | 国产一区二区激情短视频| 一个人免费在线观看的高清视频| 国产精品自产拍在线观看55亚洲| 91av网站免费观看| 欧美日韩乱码在线| 久久精品国产亚洲av高清一级| 亚洲精品中文字幕一二三四区| 51午夜福利影视在线观看| 69av精品久久久久久| 中国美女看黄片| 男插女下体视频免费在线播放| 给我免费播放毛片高清在线观看| 久久午夜亚洲精品久久| 精品欧美国产一区二区三| 久久九九热精品免费| 国产成人一区二区三区免费视频网站| 国产aⅴ精品一区二区三区波| 一本综合久久免费| 久久欧美精品欧美久久欧美| 国产亚洲欧美在线一区二区| 99在线视频只有这里精品首页| 色综合欧美亚洲国产小说| 91在线观看av| 日本免费一区二区三区高清不卡| av国产免费在线观看| 91成年电影在线观看| 久久天躁狠狠躁夜夜2o2o| 精品欧美国产一区二区三| 亚洲天堂国产精品一区在线| 成人高潮视频无遮挡免费网站| 三级男女做爰猛烈吃奶摸视频| 欧美av亚洲av综合av国产av| 女生性感内裤真人,穿戴方法视频| 99久久无色码亚洲精品果冻| 免费在线观看亚洲国产| 亚洲欧洲精品一区二区精品久久久| 夜夜躁狠狠躁天天躁| 超碰成人久久| 夜夜看夜夜爽夜夜摸| 午夜a级毛片| 19禁男女啪啪无遮挡网站| 18禁美女被吸乳视频| 制服诱惑二区| 亚洲欧美激情综合另类| 一本精品99久久精品77| 免费一级毛片在线播放高清视频| 精品福利观看| 又粗又爽又猛毛片免费看| 婷婷六月久久综合丁香| 国产午夜精品久久久久久| 18美女黄网站色大片免费观看| 一级作爱视频免费观看| 免费在线观看亚洲国产| 国产精品久久电影中文字幕| 久久热在线av| 一本大道久久a久久精品| 久久久久久人人人人人| 一区福利在线观看| av有码第一页| 桃红色精品国产亚洲av| 热99re8久久精品国产| 熟女少妇亚洲综合色aaa.| 中文字幕久久专区| 精品不卡国产一区二区三区| 国产精品久久久av美女十八| 国产午夜精品久久久久久| 成年人黄色毛片网站| 国产精品 欧美亚洲| 99国产精品一区二区三区| 成人特级黄色片久久久久久久| av国产免费在线观看| 国产一区二区三区在线臀色熟女| 一本久久中文字幕| 久久伊人香网站| 国产精品亚洲美女久久久| 国产精品久久久久久亚洲av鲁大| 中文字幕人妻丝袜一区二区| 日本 欧美在线| www国产在线视频色| 国产免费男女视频| 日韩欧美精品v在线| 亚洲无线在线观看| 国产欧美日韩一区二区精品| 国产精品一区二区免费欧美| 欧美日韩精品网址| 国产三级在线视频| 91成年电影在线观看| 99riav亚洲国产免费| 久久久久久亚洲精品国产蜜桃av| 黑人欧美特级aaaaaa片| 亚洲中文字幕日韩| 亚洲美女视频黄频| 18美女黄网站色大片免费观看| 午夜老司机福利片| 最新美女视频免费是黄的| 国产午夜精品论理片| 亚洲av成人不卡在线观看播放网| 一进一出好大好爽视频| 久久久久久大精品| 一本久久中文字幕| 国产精品亚洲av一区麻豆| 亚洲av中文字字幕乱码综合| 夜夜看夜夜爽夜夜摸| 中文字幕人妻丝袜一区二区| 久久久久国内视频| 国产成人啪精品午夜网站| 国产免费男女视频| 此物有八面人人有两片| 久久精品影院6| 成年人黄色毛片网站| 国产激情久久老熟女| 欧美高清成人免费视频www| 天天躁夜夜躁狠狠躁躁| 国产精品九九99| 成人av在线播放网站| 看黄色毛片网站| 亚洲第一欧美日韩一区二区三区| 窝窝影院91人妻| av免费在线观看网站| 国产精品久久久久久亚洲av鲁大| 国产真人三级小视频在线观看| 女生性感内裤真人,穿戴方法视频| 毛片女人毛片| 一本精品99久久精品77| 三级国产精品欧美在线观看 | 亚洲中文字幕日韩| 99re在线观看精品视频| 亚洲欧美一区二区三区黑人| 国产v大片淫在线免费观看| 一区二区三区国产精品乱码| 国产黄片美女视频| 亚洲欧美日韩无卡精品| 成在线人永久免费视频| 久99久视频精品免费| 久久午夜综合久久蜜桃| 婷婷精品国产亚洲av| av国产免费在线观看| 嫩草影院精品99| 欧美人与性动交α欧美精品济南到| 观看免费一级毛片| 国产精品国产高清国产av| 2021天堂中文幕一二区在线观| 久久久国产精品麻豆| 黄色片一级片一级黄色片| 精品一区二区三区av网在线观看| 亚洲av五月六月丁香网| 国产成人系列免费观看| 国产精品精品国产色婷婷| 日本一区二区免费在线视频| 亚洲欧美日韩高清在线视频| 欧美成人午夜精品| 少妇被粗大的猛进出69影院| 亚洲精品中文字幕一二三四区| 在线看三级毛片| 色综合婷婷激情| 久久久久久人人人人人| 精品免费久久久久久久清纯| 亚洲欧美激情综合另类| 一本大道久久a久久精品| 亚洲成人精品中文字幕电影| 久久久久久国产a免费观看| 香蕉国产在线看| 观看免费一级毛片| 亚洲全国av大片| 国产黄色小视频在线观看| 久久国产精品人妻蜜桃| av中文乱码字幕在线| 最近最新免费中文字幕在线| 免费一级毛片在线播放高清视频| 欧美在线黄色| 丝袜人妻中文字幕| netflix在线观看网站| 村上凉子中文字幕在线| 中文字幕熟女人妻在线| 日本一区二区免费在线视频| 成人一区二区视频在线观看| av天堂在线播放| 久久久国产精品麻豆| 18禁黄网站禁片午夜丰满| 亚洲av中文字字幕乱码综合| 三级国产精品欧美在线观看 | 久久欧美精品欧美久久欧美| 亚洲精品国产精品久久久不卡| 黄片大片在线免费观看| 精品久久久久久久久久久久久| www日本在线高清视频| 日韩精品青青久久久久久| 变态另类丝袜制服| 亚洲人成电影免费在线| 精品国内亚洲2022精品成人| 日本黄色视频三级网站网址| 国产欧美日韩一区二区三| 中亚洲国语对白在线视频| 蜜桃久久精品国产亚洲av| 又黄又爽又免费观看的视频| 精华霜和精华液先用哪个| 精品免费久久久久久久清纯| 久久久国产成人精品二区| 91在线观看av| av天堂在线播放| 我要搜黄色片| 国产97色在线日韩免费| 欧美乱码精品一区二区三区| 男人舔奶头视频| 亚洲无线在线观看| 午夜激情福利司机影院| 精品久久久久久久人妻蜜臀av| 一进一出抽搐gif免费好疼| 国产一区二区在线av高清观看| 好男人在线观看高清免费视频| 国内精品久久久久精免费| 亚洲精品中文字幕在线视频| 亚洲精品色激情综合| 一本一本综合久久| 一级毛片女人18水好多| 亚洲乱码一区二区免费版| 给我免费播放毛片高清在线观看| 丁香六月欧美| 亚洲全国av大片| 国产黄a三级三级三级人| 欧美成狂野欧美在线观看| 国产精品一区二区免费欧美| 日本黄色视频三级网站网址| 亚洲18禁久久av| 欧美日韩一级在线毛片| 国产亚洲精品综合一区在线观看 | 人成视频在线观看免费观看| 国产精品久久电影中文字幕| 国产精品免费视频内射| 亚洲,欧美精品.| 最近在线观看免费完整版| 可以在线观看毛片的网站| 91字幕亚洲| 少妇裸体淫交视频免费看高清 | 亚洲精品av麻豆狂野| 嫩草影视91久久| 叶爱在线成人免费视频播放| 久久久久国产一级毛片高清牌| 又黄又爽又免费观看的视频| 99久久精品热视频| 好男人在线观看高清免费视频| 窝窝影院91人妻| 大型黄色视频在线免费观看| 欧美黑人精品巨大| 日韩大码丰满熟妇| 精品国产乱子伦一区二区三区| 日本免费一区二区三区高清不卡| 久久香蕉国产精品| 亚洲一区高清亚洲精品| 亚洲一区中文字幕在线| 国产欧美日韩一区二区三| 午夜精品在线福利| 黄色毛片三级朝国网站| 国产av在哪里看| 欧美zozozo另类| 制服人妻中文乱码| 久久久久久人人人人人| 99国产综合亚洲精品| 欧美中文综合在线视频| 亚洲成人免费电影在线观看| 亚洲av美国av| 亚洲人成伊人成综合网2020| 欧美日韩国产亚洲二区| 欧美黑人精品巨大| 午夜精品久久久久久毛片777| 国产熟女xx| 久久这里只有精品中国| 后天国语完整版免费观看| netflix在线观看网站| 久久精品国产亚洲av香蕉五月| 国产熟女午夜一区二区三区| 久久久久免费精品人妻一区二区| 国产av一区二区精品久久| 国产熟女午夜一区二区三区| 黑人操中国人逼视频| 两个人免费观看高清视频| av福利片在线| 日本熟妇午夜| 美女高潮喷水抽搐中文字幕| 老司机午夜十八禁免费视频| 老熟妇乱子伦视频在线观看| 国产区一区二久久| 国产在线精品亚洲第一网站| 亚洲av成人不卡在线观看播放网| 两个人视频免费观看高清| 日本在线视频免费播放| 午夜视频精品福利| 一夜夜www| xxx96com| av天堂在线播放| 麻豆av在线久日| 一级毛片高清免费大全| 2021天堂中文幕一二区在线观| 99久久无色码亚洲精品果冻| 国内精品久久久久精免费| 亚洲成人精品中文字幕电影| 国模一区二区三区四区视频 | 亚洲第一电影网av| 999久久久国产精品视频| 日日干狠狠操夜夜爽| 久久久精品大字幕| 亚洲欧美日韩无卡精品| 最近手机中文字幕大全| 人人妻人人看人人澡| 亚洲美女视频黄频| 国产美女午夜福利| 欧美最新免费一区二区三区| 国产一区二区激情短视频| 人妻久久中文字幕网| 欧美精品一区二区大全| 69av精品久久久久久| 免费人成在线观看视频色| 色噜噜av男人的天堂激情| av女优亚洲男人天堂| 久久鲁丝午夜福利片| 长腿黑丝高跟| 欧美成人精品欧美一级黄| 精品一区二区三区人妻视频| 欧美日韩乱码在线| 亚洲av成人精品一区久久| 日韩欧美精品v在线| 久久久色成人| 看免费成人av毛片| 哪里可以看免费的av片| 欧美丝袜亚洲另类| 婷婷精品国产亚洲av| 91久久精品国产一区二区三区| 亚洲第一电影网av| 国产极品天堂在线| 国模一区二区三区四区视频| 日韩成人伦理影院| 在线免费观看的www视频| 久久精品久久久久久久性| 成年女人看的毛片在线观看| 观看美女的网站| 欧美潮喷喷水| 最新中文字幕久久久久| 51国产日韩欧美| 欧美日韩综合久久久久久| 偷拍熟女少妇极品色| 中文字幕久久专区| 搡女人真爽免费视频火全软件| 少妇人妻精品综合一区二区 | 2021天堂中文幕一二区在线观| 禁无遮挡网站| 黄色一级大片看看| 亚洲精品国产成人久久av| 亚洲av中文字字幕乱码综合| 亚洲av.av天堂| 精品人妻一区二区三区麻豆| 日韩av不卡免费在线播放| 最后的刺客免费高清国语| 欧美激情国产日韩精品一区| 九色成人免费人妻av| 亚洲久久久久久中文字幕| 欧美又色又爽又黄视频| 毛片女人毛片| 网址你懂的国产日韩在线| 我要看日韩黄色一级片| 亚洲在线自拍视频| 最新中文字幕久久久久| 久久精品国产亚洲网站| 不卡一级毛片| 91久久精品国产一区二区成人| 国产成人午夜福利电影在线观看| 在线免费观看的www视频| 老司机福利观看| 午夜福利在线在线| 国产老妇伦熟女老妇高清| 国语自产精品视频在线第100页| 国产精品精品国产色婷婷| 大又大粗又爽又黄少妇毛片口| 久久久成人免费电影| 成人二区视频| 亚洲激情五月婷婷啪啪| 亚洲四区av| 精品不卡国产一区二区三区| 69人妻影院| 亚洲成av人片在线播放无| 天堂√8在线中文| 午夜久久久久精精品| 看非洲黑人一级黄片| 国产精品一二三区在线看| 久久99蜜桃精品久久| 国产亚洲av嫩草精品影院| 欧美最新免费一区二区三区| 91在线精品国自产拍蜜月| 午夜免费男女啪啪视频观看| 精品午夜福利在线看| 尤物成人国产欧美一区二区三区| av免费观看日本| 国产精品一二三区在线看| ponron亚洲| 在现免费观看毛片| 久久久久久久久久久免费av| 嫩草影院新地址| 久久久久免费精品人妻一区二区| 亚洲av成人av| 欧美3d第一页| 全区人妻精品视频| 国产精品久久视频播放| 欧美潮喷喷水| 黄片wwwwww| 亚洲国产精品国产精品| 亚洲自偷自拍三级| 日韩一区二区视频免费看| 国产精品爽爽va在线观看网站| 欧美+日韩+精品| 六月丁香七月| 欧美日韩一区二区视频在线观看视频在线 | 国产一级毛片在线| 蜜臀久久99精品久久宅男| 51国产日韩欧美| 久久精品人妻少妇| 可以在线观看的亚洲视频| 国产免费一级a男人的天堂| 久久九九热精品免费| 国产高清不卡午夜福利| 夜夜爽天天搞| 国产精品.久久久| 男女视频在线观看网站免费| 欧美激情久久久久久爽电影| 午夜亚洲福利在线播放| 国产午夜精品论理片| 最近视频中文字幕2019在线8| 亚洲人成网站在线观看播放| 色综合站精品国产| 赤兔流量卡办理| 久久精品综合一区二区三区| 一级毛片aaaaaa免费看小| 变态另类丝袜制服| 黄色欧美视频在线观看| 中文精品一卡2卡3卡4更新| 欧美一级a爱片免费观看看| 国产精品一区二区三区四区久久| 亚洲人成网站高清观看| 深夜精品福利| 中文字幕久久专区| 中文精品一卡2卡3卡4更新| 成人综合一区亚洲| 长腿黑丝高跟| 国产在线男女| 欧美性猛交黑人性爽| 日韩欧美在线乱码| 乱人视频在线观看| 麻豆成人av视频| av在线播放精品| .国产精品久久| 给我免费播放毛片高清在线观看| 寂寞人妻少妇视频99o| 成人av在线播放网站| 久久欧美精品欧美久久欧美| a级毛色黄片| 精品人妻视频免费看| 女的被弄到高潮叫床怎么办| 久久亚洲精品不卡| 99久久成人亚洲精品观看| 久久久久久久久大av| 尤物成人国产欧美一区二区三区| 色综合亚洲欧美另类图片| 国产在视频线在精品|