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

    山東半島藍(lán)色經(jīng)濟(jì)區(qū)經(jīng)濟(jì)—環(huán)境復(fù)合系統(tǒng)仿真與發(fā)展對(duì)策研究

    2017-04-01 16:32:21董會(huì)忠王格??
    軟科學(xué) 2017年2期
    關(guān)鍵詞:系統(tǒng)動(dòng)力學(xué)

    董會(huì)忠++王格??

    摘要:構(gòu)建了以山東半島藍(lán)色經(jīng)濟(jì)區(qū)為例的經(jīng)濟(jì)—環(huán)境復(fù)合系統(tǒng)發(fā)展系統(tǒng)動(dòng)力學(xué)模型,運(yùn)用VENSIM軟件對(duì)區(qū)域未來發(fā)展情景進(jìn)行定量化動(dòng)態(tài)研究與模擬。通過調(diào)節(jié)復(fù)合系統(tǒng)的決策變量,設(shè)定了4種不同的發(fā)展情景方案,預(yù)測(cè)在不同情景方案下山東半島藍(lán)色經(jīng)濟(jì)區(qū)的發(fā)展?fàn)顩r,得出山東半島藍(lán)色經(jīng)濟(jì)區(qū)經(jīng)濟(jì)—環(huán)境復(fù)合系統(tǒng)協(xié)調(diào)發(fā)展的最優(yōu)方案,為山東半島藍(lán)色經(jīng)濟(jì)區(qū)的規(guī)劃、治理和發(fā)展提供決策支持。

    關(guān)鍵詞:山東半島藍(lán)色經(jīng)濟(jì)區(qū);系統(tǒng)動(dòng)力學(xué);決策變量;情景方案

    DOI:10.13956/j.ss.1001-8409.2017.02.22

    中圖分類號(hào):F205;F127 文獻(xiàn)標(biāo)識(shí)碼:A 文章編號(hào):1001-8409(2017)02-0103-06

    Study on Simulation and Development Countermeasures of

    Economyenvironment Composite System in

    Shandong Peninsula Blue Economic Zone

    DONG Huizhong,WANG Ge

    (School of Business, Shandong University of Technology, Zibo 255012)

    Abstract:This paper constructed a system dynamics model of composite system of environment and economic development taking Shandong Peninsula Blue Economic Zone as an example. It used the VENSIM software to do quantitative study and dynamic simulation on the future development situation of the region. By adjusting the composite decision variables of the system, it set up 4 different kinds of development scenarios and predicted the development status of Shandong Peninsula Blue Economic Zone in different scenarios, then get the optimal solution of Shandong peninsula blue economic zone on complex system of economy and environment coordinated development to provide decision support for the planning, management and development of Shandong Peninsula blue Economic Zone.

    Key words:Shandong Peninsula Blue Economic Zone; system dynamics; decision variables; scenarios

    區(qū)域協(xié)調(diào)發(fā)展涉及到社會(huì)、經(jīng)濟(jì)、資源、環(huán)境等諸多因素[1],是多個(gè)子系統(tǒng)耦合演化的結(jié)果。許多地區(qū)在發(fā)展過程中僅僅關(guān)注經(jīng)濟(jì)發(fā)展速度,而忽視當(dāng)?shù)刭Y源與環(huán)境承載力,造成資源的嚴(yán)重短缺和生態(tài)環(huán)境的普遍惡化,區(qū)域發(fā)展處于失衡狀態(tài),難以持續(xù)發(fā)展。當(dāng)前,人們已普遍意識(shí)到以犧牲資源環(huán)境為代價(jià)的發(fā)展不可持續(xù),努力探求社會(huì)經(jīng)濟(jì)與生態(tài)環(huán)境協(xié)調(diào)發(fā)展的途徑已成為各國政府高度關(guān)注的焦點(diǎn)和眾多學(xué)者研究的熱點(diǎn)。作為一個(gè)復(fù)雜巨系統(tǒng),區(qū)域經(jīng)濟(jì)—環(huán)境復(fù)合系統(tǒng)具有多變量、非線性、動(dòng)態(tài)性和多重反饋的特性,系統(tǒng)動(dòng)力學(xué)是研究復(fù)雜巨系統(tǒng)的科學(xué)方法,能夠?qū)Σ煌榫胺桨赶碌膮^(qū)域發(fā)展態(tài)勢(shì)進(jìn)行有效研究,并對(duì)未來發(fā)展趨勢(shì)進(jìn)行科學(xué)預(yù)測(cè)。該方法在國外廣泛應(yīng)用于經(jīng)濟(jì)社會(huì)與生態(tài)環(huán)境研究領(lǐng)域[3-5],代表文獻(xiàn)有《增長(zhǎng)的限制》和《趨向全球平衡》(以WORD Ⅲ模型為基礎(chǔ))[2]。上世紀(jì)80年代初,系統(tǒng)動(dòng)力學(xué)方法被引入我國后普遍應(yīng)用于國民經(jīng)濟(jì)、區(qū)域合作、可持續(xù)發(fā)展、技術(shù)創(chuàng)新等領(lǐng)域[6-11]。本文采用系統(tǒng)動(dòng)力學(xué)仿真技術(shù)對(duì)山東半島藍(lán)色經(jīng)濟(jì)區(qū)環(huán)境經(jīng)濟(jì)復(fù)合系統(tǒng)進(jìn)行仿真模擬,并對(duì)系統(tǒng)狀態(tài)進(jìn)行評(píng)估,定量研究半島經(jīng)濟(jì)區(qū)經(jīng)濟(jì)—環(huán)境系統(tǒng)協(xié)調(diào)發(fā)展?fàn)顟B(tài),對(duì)各種情景的作用和效果進(jìn)行比較和分析,從而作出合理選擇,為山東半島藍(lán)色經(jīng)濟(jì)區(qū)發(fā)展規(guī)劃的制定提供科學(xué)依據(jù)。

    1系統(tǒng)動(dòng)力學(xué)模型構(gòu)建

    11研究區(qū)域邊界確定

    山東半島藍(lán)色經(jīng)濟(jì)區(qū)陸域面積64萬平方千米,范圍包括青島、東營、煙臺(tái)、濰坊、威海、日照6市及濱州的無棣、沾化2個(gè)沿??h所屬陸域,海域面積1595萬平方千米,包括山東全部海域。為研究方便,本文構(gòu)建的模型及相關(guān)數(shù)據(jù)分析界定在青島、東營、煙臺(tái)、濰坊、威海、日照及濱州市全部區(qū)域。

    12區(qū)域復(fù)合系統(tǒng)的結(jié)構(gòu)分析及子系統(tǒng)劃分

    根據(jù)經(jīng)濟(jì)環(huán)境協(xié)調(diào)發(fā)展理論[12],將山東半島藍(lán)色經(jīng)濟(jì)區(qū)經(jīng)濟(jì)環(huán)境復(fù)合系統(tǒng)劃分為經(jīng)濟(jì)、人口、資源、環(huán)境四個(gè)子系統(tǒng),各子系統(tǒng)內(nèi)部要素相互作用,相互關(guān)聯(lián),形成一個(gè)動(dòng)態(tài)復(fù)雜系統(tǒng),以系統(tǒng)各變量之間的相互關(guān)系為基本框架,建立山東半島藍(lán)色經(jīng)濟(jì)區(qū)復(fù)合系統(tǒng)模型,通過調(diào)整決策變量,對(duì)預(yù)設(shè)發(fā)展模式進(jìn)行仿真模擬,進(jìn)而找到山東半島藍(lán)色經(jīng)濟(jì)區(qū)發(fā)展的最佳模式。

    (1)經(jīng)濟(jì)子系統(tǒng)

    經(jīng)濟(jì)子系統(tǒng)把能夠反映經(jīng)濟(jì)發(fā)展水平的第一產(chǎn)業(yè)、工業(yè)、建筑業(yè)和第三產(chǎn)業(yè)產(chǎn)值作為流位變量(水平變量);把第一產(chǎn)業(yè)GDP增加值、第三產(chǎn)業(yè)GDP增加值、建筑業(yè)GDP增加值和工業(yè)GDP增加值作為流率變量(速率變量);把區(qū)域GDP、第一產(chǎn)業(yè)GDP增長(zhǎng)率、第三產(chǎn)業(yè)GDP增長(zhǎng)率、建筑業(yè)GDP增長(zhǎng)率和工業(yè)GDP增長(zhǎng)率等作為輔助變量。通過人均GDP、產(chǎn)業(yè)用水、工業(yè)COD排放量和工業(yè)氨氮排放量等輔助變量分別與人口、資源、環(huán)境子系統(tǒng)建立關(guān)聯(lián)?!叭a(chǎn)”比例是反映山東半島藍(lán)色經(jīng)濟(jì)區(qū)區(qū)域經(jīng)濟(jì)產(chǎn)業(yè)結(jié)構(gòu)合理性的指標(biāo),因此以第一產(chǎn)業(yè)GDP增長(zhǎng)率、工業(yè)GDP增長(zhǎng)率和第三產(chǎn)業(yè)GDP增長(zhǎng)率作為經(jīng)濟(jì)子系統(tǒng)的決策變量。通過以上指標(biāo)建立經(jīng)濟(jì)子系統(tǒng)模型。包括14個(gè)方程和4個(gè)表函數(shù),其中,主要方程有:

    第一產(chǎn)業(yè)GDP=INTEG(第一產(chǎn)業(yè)GDP增加值,70731,單位:億元)

    建筑業(yè)GDP=INTEG(建筑業(yè)GDP增加值,38204,,單位:億元)

    工業(yè)GDP=INTEG(工業(yè)GDP增加值,317624,單位:億元)

    第三產(chǎn)業(yè)GDP=INTEG(第三產(chǎn)業(yè)GDP增加值,21111,單位:億元)

    區(qū)域GDP=第一產(chǎn)業(yè)GDP+第三產(chǎn)業(yè)產(chǎn)業(yè)GDP+建筑業(yè)GDP+工業(yè)GDP(單位:億元)

    工業(yè)GDP增長(zhǎng)率=工業(yè)GDP增長(zhǎng)率 LOOKUP(Time)

    建筑業(yè)GDP增長(zhǎng)率=建筑業(yè)GDP增長(zhǎng)率 LOOKUP(Time)

    第一產(chǎn)業(yè)GDP增長(zhǎng)率=第一產(chǎn)業(yè)GDP增長(zhǎng)率 LOOKUP(Time)

    第三產(chǎn)業(yè)GDP增長(zhǎng)率=第三產(chǎn)業(yè)GDP增長(zhǎng)率 LOOKUP(Time)

    (2)資源子系統(tǒng)

    土地、水、能源及其他礦產(chǎn)等自然資源是人類社會(huì)發(fā)展的重要物質(zhì)基礎(chǔ),自然資源的合理開發(fā)和利用對(duì)于區(qū)域協(xié)調(diào)發(fā)展具有重要意義。在資源子系統(tǒng)中,把可用水資源量和耕地面積作為水平變量;國家基建占地、退耕還林還草占地、鄉(xiāng)村集體占地以及水資源可供量和工業(yè)用水量、農(nóng)業(yè)用水量、生活用水量、水資源循環(huán)利用率,再用水回用率等作為輔助變量。通過生活用水量以及人均耕地面積與人口子系統(tǒng)建立關(guān)聯(lián),以生活污水處理回用量作為資源子系統(tǒng)和環(huán)境子系統(tǒng)溝通的橋梁。通過引入耕地壓力指數(shù)模型[13]來考察山東半島藍(lán)色經(jīng)濟(jì)區(qū)耕地壓力情況。耕地壓力指數(shù)——最小人均耕地面積與實(shí)際人均耕地面積之比,其算法如下:

    K=SminS(1)

    式中K為耕地壓力指數(shù),Smin為最小人均耕地面積,即一定區(qū)域范圍內(nèi)保障食物需求的最小人均耕地面積,S為人均實(shí)際耕地面積。計(jì)算耕地壓力指數(shù)時(shí),在給出人均實(shí)際耕地面積的情況下,最小人均耕地面積是影響耕地壓力指數(shù)的關(guān)鍵變量,其算法如下:

    Smin=Gr×βP×Q×k(2)

    式中,Smin為山東半島藍(lán)色經(jīng)濟(jì)區(qū)最小人均耕地面積,Gr為山東半島藍(lán)色經(jīng)濟(jì)區(qū)區(qū)域人均糧食需求量,β為山東半島藍(lán)色經(jīng)濟(jì)區(qū)區(qū)域糧食自給率,P為山東半島藍(lán)色經(jīng)濟(jì)區(qū)區(qū)域糧食單產(chǎn),Q為山東半島藍(lán)色經(jīng)濟(jì)區(qū)區(qū)域糧食播種面積占總播種面積的比重,k為山東半島藍(lán)色經(jīng)濟(jì)區(qū)區(qū)域復(fù)種指數(shù)。

    當(dāng)K>1時(shí),表明耕地壓力明顯,容易出現(xiàn)糧食危機(jī),需要采取緊急必要措施降低耕地壓力指數(shù);當(dāng)K<1時(shí),表明耕地壓力不明顯,糧食生產(chǎn)處于安全范圍之內(nèi);K=1時(shí),最小人均耕地面積等于實(shí)際人均耕地面積,表明耕地壓力開始顯現(xiàn),應(yīng)當(dāng)密切關(guān)注耕地面積的變化,采取措施避免耕地壓力變大,確保糧食生產(chǎn)的安全。資源子系統(tǒng)方程有22個(gè),表函數(shù)6個(gè),主要方程有:

    水資源量=INTEG(生產(chǎn)量-消耗量,1101,單位:億立方米)

    生活用水量=農(nóng)村生活用水量+城鎮(zhèn)生活用水量(單位:億立方米)

    工業(yè)用水量=工業(yè)GDP×工業(yè)億元產(chǎn)值用水量(單位:億立方米)

    再用水回用量=生活污水處理量×再用水回用率(單位:億立方米)

    工業(yè)億元產(chǎn)值用水量=工業(yè)億元產(chǎn)值用水量 LOOKUP(Time)

    城鎮(zhèn)人均生活用水量=城鎮(zhèn)人均生活用水量 LOOKUP(Time)

    耕地面積=耕地增加量-耕地減少量(單位:公頃)

    耕地壓力指數(shù)=最小人均耕地面積/實(shí)際人均耕地面積

    (3)人口子系統(tǒng)

    人口子系統(tǒng)建模過程中選取總?cè)丝谧鳛榱魑蛔兞浚ㄋ?/p>

    平變量);以年自然增長(zhǎng)人口和年凈遷入人口為流率變量(速率變量);把農(nóng)村人口、城鎮(zhèn)人口、自然增長(zhǎng)率、人口凈遷入率作為輔助變量。人口子系統(tǒng)共有方程9個(gè),3個(gè)表函數(shù),主要方程有:

    總?cè)丝?INTEG(年自然增長(zhǎng)人口+年凈遷入人口,528628,單位:萬人)

    城鎮(zhèn)人口=總?cè)丝凇脸鞘谢剑▎挝唬喝f人)

    農(nóng)村人口=總?cè)丝?城鎮(zhèn)人口(單位:萬人)

    自然增長(zhǎng)率=自然增長(zhǎng)率 LOOKUP(Time)

    人口凈遷入率=人口凈遷入率 LOOKUP(Time)

    城鎮(zhèn)化水平=城鎮(zhèn)化水平 LOOKUP(Time)

    (4)環(huán)境子系統(tǒng)

    環(huán)境子系統(tǒng)建模過程中,選取COD存量和氨氮存量作為流位變量(水平變量);以COD產(chǎn)生量、COD削減量,氨氮增加量、氨氮削減量為流率變量(速率變量);以工業(yè)COD產(chǎn)生量、生活COD產(chǎn)生量、工業(yè)廢水COD削減量、工業(yè)氨氮產(chǎn)生量、工業(yè)廢水氨氮削減量、生活COD削減量、生活氨氮產(chǎn)生量、生活氨氮削減量、工業(yè)廢水處理量、工業(yè)廢水排放量、工業(yè)廢水治理費(fèi)用、工業(yè)億元產(chǎn)值廢水量等為輔助變量。選取COD和氨氮存量作為水環(huán)境污染程度的指示指標(biāo)。環(huán)境子系統(tǒng)共有方程22個(gè),表函數(shù)11個(gè),主要方程有:

    COD存量=INTEG(COD產(chǎn)生量-COD削減量,539793,單位:噸)

    氨氮存量=INTEG(氨氮產(chǎn)生量-氨氮削減量,53538,單位:噸)

    工業(yè)廢水產(chǎn)生量=工業(yè)GDP×工業(yè)億元產(chǎn)值廢水量(單位:噸)

    生活污水產(chǎn)生量=城鎮(zhèn)人口×城鎮(zhèn)居民人均生活污水產(chǎn)污系數(shù)(單位:萬噸)

    工業(yè)廢水治理費(fèi)用占工業(yè)GDP比重=工業(yè)廢水治理費(fèi)用占工業(yè)GDP比重 LOOKUP(Time)

    城鎮(zhèn)居民人均生活污水產(chǎn)污系數(shù)=城鎮(zhèn)居民人均生活污水產(chǎn)污系數(shù) LOOKUP(Time)

    13經(jīng)濟(jì)—環(huán)境復(fù)合系統(tǒng)模型構(gòu)建

    通過對(duì)山東半島藍(lán)色經(jīng)濟(jì)區(qū)經(jīng)濟(jì)、社會(huì)、資源、環(huán)境子系統(tǒng)的結(jié)構(gòu)分析,選取能夠合理描述各子系統(tǒng)的變量,建立復(fù)合系統(tǒng)動(dòng)力學(xué)模型(如圖1所示)。

    14模型有效性檢驗(yàn)

    系統(tǒng)動(dòng)力學(xué)通過Vensim軟件建立系統(tǒng)動(dòng)力學(xué)模型,對(duì)復(fù)合系統(tǒng)發(fā)展態(tài)勢(shì)進(jìn)行預(yù)測(cè),其中模型的有效性檢驗(yàn)是模型預(yù)測(cè)的關(guān)鍵環(huán)節(jié)。本模型模擬運(yùn)行時(shí)間為2013~2020年,以2013年為基準(zhǔn)年,仿真步長(zhǎng)為1年,模型有效性檢驗(yàn)采用的歷史統(tǒng)計(jì)數(shù)據(jù)主要來源于《山東省統(tǒng)計(jì)年鑒》以及各市統(tǒng)計(jì)年鑒。以2003~2012年歷史數(shù)據(jù)為檢驗(yàn)依據(jù),其中某些指標(biāo)如年凈遷入率、自然增長(zhǎng)率、城鎮(zhèn)居民人均COD產(chǎn)污系數(shù)等是根據(jù)2003~2012年所收集數(shù)據(jù)計(jì)算得出,某些指標(biāo)隨著時(shí)間的推移而變化,這些指標(biāo)采用表函數(shù)的形式給出,而變量間關(guān)系不明確的指標(biāo)也是采用表函數(shù)形式給出。本文采用相對(duì)誤差法對(duì)半島藍(lán)色經(jīng)濟(jì)區(qū)復(fù)合系統(tǒng)模型模擬的結(jié)果進(jìn)行檢驗(yàn),以總?cè)丝诤虲OD存量?jī)蓚€(gè)水平變量為主要例證,檢驗(yàn)結(jié)果如表1,通過表1可以看出,2003~2012年總?cè)丝诤虲OD存量的仿真值和歷史值的相對(duì)誤差介于-1%~51295%之間,擬合精確度較高,相對(duì)誤差檢驗(yàn)結(jié)果比較理想。

    2山東半島藍(lán)色經(jīng)濟(jì)區(qū)模型預(yù)測(cè)與分析

    21情景設(shè)定

    在對(duì)2003~2012年數(shù)據(jù)有效性檢驗(yàn)的基礎(chǔ)上,本文設(shè)定不同的決策變量來調(diào)控2013~2020年山東半島藍(lán)色經(jīng)濟(jì)區(qū)區(qū)域經(jīng)濟(jì)發(fā)展的模式。選取第一產(chǎn)業(yè)GDP增長(zhǎng)率、工業(yè)GDP增長(zhǎng)率、第三產(chǎn)業(yè)GDP增長(zhǎng)率作為調(diào)整經(jīng)濟(jì)結(jié)構(gòu)的決策變量,將工業(yè)廢水治理費(fèi)用占工業(yè)GDP的比重、再用水回用率以及農(nóng)業(yè)億元產(chǎn)值用水量作為調(diào)控環(huán)境和資源的決策變量。由于實(shí)際人均耕地面積與人口基數(shù)有很大關(guān)系,在耕地面積逐年減少、人口不斷增加的情況下,實(shí)際人均耕地面積呈現(xiàn)逐年減少的狀況,因而最小人均耕地面積成為影響耕地壓力指數(shù)的關(guān)鍵變量;在糧食自給率為1,糧食單產(chǎn)保持2012年的水平上,復(fù)種指數(shù)成為影響最小人均耕地面積的關(guān)鍵變量。根據(jù)山東省三次產(chǎn)業(yè)比重調(diào)整戰(zhàn)略,提高第三產(chǎn)業(yè)比重,降低第一、二產(chǎn)業(yè)比重的方向,設(shè)定出四種發(fā)展模式情景:

    情景Ⅰ:設(shè)定第一產(chǎn)業(yè)GDP增長(zhǎng)率降低25個(gè)百分點(diǎn),從2013年的11%下降到2020年的85%;工業(yè)GDP增長(zhǎng)率增加17個(gè)百分點(diǎn),從2013年53%增加到2020年的7%;第三產(chǎn)業(yè)GDP增長(zhǎng)率小幅度穩(wěn)步提高,增加29個(gè)百分點(diǎn),從2013年的127%增加到2020年的156%;工業(yè)廢水治理費(fèi)用占工業(yè)GDP的比重和復(fù)種指數(shù)以及農(nóng)業(yè)億元產(chǎn)值用水量維持在2012年的水平,再用水回用率從2013年到2020年提高6個(gè)百分點(diǎn)。

    情景Ⅱ:設(shè)定第一產(chǎn)業(yè)GDP增長(zhǎng)率2016年降低到情景Ⅰ2020年85%的水平,到2020年降低到62%;工業(yè)GDP增長(zhǎng)率略微降低,2013年到2020年降低08個(gè)百分點(diǎn);第三產(chǎn)業(yè)GDP增速與情景Ⅰ相比較明顯,2016年達(dá)到情景Ⅰ2020年156%的水平,到2020年第三產(chǎn)業(yè)GDP增長(zhǎng)率達(dá)到175%;工業(yè)廢水治理費(fèi)用占工業(yè)GDP的比重、再用水回用率、復(fù)種指數(shù)均穩(wěn)步提升;農(nóng)業(yè)億元產(chǎn)值用水量略微降低。

    情景Ⅲ:設(shè)定第一產(chǎn)業(yè)GDP增長(zhǎng)率降低幅度與情景Ⅰ相似;工業(yè)GDP增長(zhǎng)率明顯降低,2016年降低到情景Ⅱ2020年45%的水平,到2020年工業(yè)GDP增長(zhǎng)率降低到40%;第三產(chǎn)業(yè)GDP增速明顯提高,2016年達(dá)到164%,2020年達(dá)到185%;再用水回用率以及復(fù)種指數(shù)在情景Ⅱ的基礎(chǔ)上繼續(xù)提高,工業(yè)廢水治理費(fèi)用占工業(yè)GDP的比重略微提升,與情景Ⅱ相似,農(nóng)業(yè)億元產(chǎn)值用水量在情景Ⅱ的基礎(chǔ)上繼續(xù)降低。

    情景Ⅳ:第一產(chǎn)業(yè)GDP增長(zhǎng)率、工業(yè)GDP增長(zhǎng)率顯著降低,降幅較情景Ⅱ相比較大。其中,第一產(chǎn)業(yè)GDP增長(zhǎng)率到2020年降低到43%,工業(yè)GDP增長(zhǎng)率降低到35%;第三產(chǎn)業(yè)GDP增長(zhǎng)率顯著提高,到2020年第三產(chǎn)業(yè)GDP增長(zhǎng)率達(dá)到268%,工業(yè)廢水治理費(fèi)用占工業(yè)GDP比重顯著增加,再用水回用率、復(fù)種指數(shù)顯著提高,與情景Ⅲ相比提升幅度較大,農(nóng)業(yè)億元產(chǎn)值用水量明顯減少。

    22模擬結(jié)果分析

    情景模擬結(jié)果分析目的在于找出適合山東半島藍(lán)色經(jīng)濟(jì)區(qū)的較為理想的發(fā)展模式,因而指標(biāo)變量的選取至關(guān)重要,它能夠反映區(qū)域復(fù)合系統(tǒng)的基本情況。圖2至圖5為模型運(yùn)行得出的指標(biāo)變量變化趨勢(shì)圖,表3為模型運(yùn)行得出的各項(xiàng)指標(biāo)模擬值。

    從仿真模擬出的結(jié)果可以看出:

    情景Ⅰ發(fā)展模式是逐步降低第一產(chǎn)業(yè)GDP增長(zhǎng)率,略微增加工業(yè)GDP增長(zhǎng)率,同時(shí)第三產(chǎn)業(yè)GDP增長(zhǎng)率穩(wěn)步提高;工業(yè)廢水治理費(fèi)用占工業(yè)GDP比重、復(fù)種指數(shù)、農(nóng)業(yè)億元產(chǎn)值用水量維持在2012年的水平。從模擬結(jié)果可以看出,區(qū)域GDP增速從2016年開始增幅略高于目前發(fā)展模式;COD存量與目前發(fā)展模式相比無明顯差異;可用水資源量下滑趨勢(shì)明顯降低,低于目前發(fā)展模式下的水資源量,供水形勢(shì)嚴(yán)峻,2016年供水缺口為1375億立方米,2020年供水缺口達(dá)到2891億立方米;耕地壓力指數(shù)與現(xiàn)狀相比,無明顯變化,耕地壓力形勢(shì)嚴(yán)峻,容易產(chǎn)生糧食危機(jī)。

    情景Ⅱ發(fā)展模式是在穩(wěn)步降低第一、二產(chǎn)業(yè)GDP增長(zhǎng)率的前提下,在情景Ⅰ的基礎(chǔ)上提高第三產(chǎn)業(yè)GDP增長(zhǎng)率的速度;同時(shí)工業(yè)廢水治理費(fèi)用占工業(yè)GDP比重、復(fù)種指數(shù)穩(wěn)步提高,再用水回用率從2013年的40%提高到2016年的45%、再到2020年的48%,農(nóng)業(yè)億元產(chǎn)值用水量減少緩慢。從模擬結(jié)果可以看出,區(qū)域GDP與情景Ⅰ相比,無明顯差別,但增速略高于情景Ⅰ;COD存量增長(zhǎng)趨勢(shì)趨于平緩,明顯低于情景Ⅰ,到2020年COD存量為54741萬噸,明顯低于情景ⅠCOD存量2020年87036萬噸的水平;可用水資源量減少趨勢(shì)明顯低于情景Ⅰ,但供水缺口仍然存在,2016年供水缺口為105億立方米,2020年供水缺口為2211億立方米,耕地壓力指數(shù)呈現(xiàn)降低趨勢(shì)但依然高于1,耕地壓力依然明顯。

    情景Ⅲ發(fā)展模式是第一、二產(chǎn)業(yè)GDP增長(zhǎng)率降幅與情景Ⅰ相似,但低于情景Ⅱ,第三產(chǎn)業(yè)GDP增長(zhǎng)率較情景Ⅱ相比增幅較大,2016年第三產(chǎn)業(yè)GDP增長(zhǎng)率為164%,2020年為185%;同時(shí)工業(yè)廢水治理費(fèi)用比重提高,與情景Ⅱ相比,無明顯差別,再用水回用率2016年提高到48%,2020年達(dá)到50%,農(nóng)業(yè)億元產(chǎn)值用水量降低幅度高于情景Ⅱ。從模擬結(jié)果可以看出,區(qū)域GDP增速略高于情景Ⅰ和情景Ⅱ;COD存量增速明顯低于情景Ⅰ和現(xiàn)狀發(fā)展模式,但高于情景Ⅱ,可用水資源量降幅低于情景Ⅱ,耕地壓力指數(shù)明顯低于情景Ⅱ, 2020年耕地壓力指數(shù)為126,耕地壓力依然存在。

    情景Ⅳ發(fā)展模式是明顯降低第一、二產(chǎn)業(yè)增長(zhǎng)率,同時(shí)大幅提高第三產(chǎn)業(yè)GDP增長(zhǎng)率,到2016年,第三產(chǎn)業(yè)GDP增長(zhǎng)率為185%,2020年達(dá)到268%;工業(yè)廢水治理費(fèi)用比重明顯加大,再用水回用率、復(fù)種指數(shù)明顯提高,農(nóng)業(yè)億元產(chǎn)值用水量明顯降低。從模擬結(jié)果可以看出,區(qū)域GDP增長(zhǎng)幅度明顯高于情景Ⅰ、Ⅱ、Ⅲ,區(qū)域GDP規(guī)模2016年為417429億元,2020年達(dá)到942623億元,人均GDP明顯提高;COD存量大幅度降低,從2016年開始COD存量開始降低,到2020年,COD和氨氮治理明顯,水環(huán)境明顯改善, 2016年供水缺口為695億立方米,2020年供水缺口為1176億立方米,間接提高了水資源可用量;耕地壓力指數(shù)2020年為098,低于1,耕地壓力明顯改善。

    3結(jié)論與建議

    根據(jù)以上分析可以得出在情景Ⅰ中雖注重產(chǎn)業(yè)結(jié)構(gòu)調(diào)整,但產(chǎn)業(yè)結(jié)構(gòu)調(diào)整比例不明顯,COD和氨氮環(huán)境壓力形勢(shì)較現(xiàn)狀發(fā)展更加嚴(yán)峻,耕地壓力指數(shù)較大,耕地壓力明顯;情景Ⅱ調(diào)整了經(jīng)濟(jì)產(chǎn)業(yè)結(jié)構(gòu),第一產(chǎn)業(yè)增速降低,第二產(chǎn)業(yè)增速小幅增高,第三產(chǎn)業(yè)增速較情景Ⅰ穩(wěn)步提高,同時(shí)提高環(huán)境治理費(fèi)用,提高再用水回用率,環(huán)境壓力明顯減小,但可用水資源量的緊張程度依然明顯,供水缺口比情景Ⅰ有所減少,耕地壓力指數(shù)與情景Ⅰ相比雖有所下降,但依然高于1,耕地壓力明顯;情景Ⅲ第一、二產(chǎn)業(yè)增速逐步降低,第三產(chǎn)業(yè)增速有較大提高,環(huán)境治理費(fèi)用比重、再用水回用率、復(fù)種指數(shù)提高,農(nóng)業(yè)億元產(chǎn)值用水有較明顯的減少,環(huán)境壓力比情景Ⅰ明顯縮小,與情景Ⅱ相比有小幅增加,供水缺口低于情景Ⅰ、Ⅱ,水資源緊張程度依然存在;情景Ⅳ是著重加大產(chǎn)業(yè)結(jié)構(gòu)調(diào)整,明顯降低第一、二產(chǎn)業(yè)增長(zhǎng)速度,大幅提高第三產(chǎn)業(yè)增長(zhǎng)速度,環(huán)境治理投資比重增大,再用水回用率明顯提高,農(nóng)業(yè)億元產(chǎn)值用水量明顯較少,到2020年,區(qū)域經(jīng)濟(jì)產(chǎn)值達(dá)到942623億元;COD存量2016年為471253萬噸,到2020年COD存量基本消除,COD和氨氮治理明顯,水環(huán)境壓力明顯改善;供水缺口2016年為695億立方米,2020年為1176億立方米,有效緩解了水資源緊張程度。通過以上結(jié)論,推薦情景Ⅳ為山東半島藍(lán)色經(jīng)濟(jì)區(qū)發(fā)展方向的最優(yōu)情景。

    通過上述四種情景發(fā)展模式的仿真模擬結(jié)果,提出如下建議:(1)結(jié)合山東半島藍(lán)色經(jīng)濟(jì)區(qū)6市差異,進(jìn)一步推動(dòng)產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化升級(jí),降低第一、二產(chǎn)業(yè)的比重,大力提高第三產(chǎn)業(yè)比重,推動(dòng)宏觀產(chǎn)業(yè)結(jié)構(gòu)優(yōu)化,形成合理的產(chǎn)業(yè)結(jié)構(gòu)發(fā)展模式,以進(jìn)一步提高經(jīng)濟(jì)發(fā)展的質(zhì)量和效益;(2)以“保持經(jīng)濟(jì)平穩(wěn)增長(zhǎng)、水資源消耗低增長(zhǎng)”為發(fā)展模式,推進(jìn)節(jié)水減排計(jì)劃,合理分配水資源用水量在各產(chǎn)業(yè)結(jié)構(gòu)中的比例;(3)出臺(tái)相關(guān)的政策和提高排污標(biāo)準(zhǔn),制定相應(yīng)的評(píng)價(jià)體系,通過增加工業(yè)廢水治理費(fèi)用,大力發(fā)展廢水再回用水的循環(huán)利用率,從而達(dá)到減少水環(huán)境污染、增加水資源量的雙重效果。

    參考文獻(xiàn):

    [1]覃成林. 區(qū)域協(xié)調(diào)發(fā)展機(jī)制體系研究[J]. 經(jīng)濟(jì)學(xué)家, 2011, (4): 63-70.

    [2]王其藩. 系統(tǒng)動(dòng)力學(xué)[M].2009年修訂版. 上海:上海財(cái)經(jīng)大學(xué)出版社, 2009.

    [3]Forrester J W. The System Dynamics National Model: Macrobehavior from Microstructure Modeling Growth Strategy in a Biotechnology Startup Firm [J].System Dynamics Review, 1989, 7(2): 93-116.

    [4]S M Disney A,T Potter, B M, Cardner. The Impact of Vendor Managed Inventory on Transport Operation [J].Transportation Research E, 2003.

    [5]Marquez A C, Bianchi C, Gupta J N D. Operational and Financial Effectiveness of Ecollaboration Tools in Supply Chain Integration[J]. European Journal of Operational Research,2004.

    [6]張建慧,雷星暉,李金良. 基于系統(tǒng)動(dòng)力學(xué)城市低碳交通發(fā)展模式研究——以鄭州市為例[J]. 軟科學(xué), 2012, 26(4): 77-81.

    [7]陳國衛(wèi),金家善,耿俊豹. 系統(tǒng)動(dòng)力學(xué)應(yīng)用研究綜述[J]. 控制工程, 2012, 19 (6):921-928.

    [8]李維乾,解建倉,李建勛,等. 基于系統(tǒng)動(dòng)力學(xué)的閉環(huán)反饋水資源優(yōu)化配置研究[J]. 西北農(nóng)林科技大學(xué)學(xué)報(bào)(自然科學(xué)版), 2013, 41(11):209-216.

    [9]賈一偉. 基于系統(tǒng)動(dòng)力學(xué)的高??萍籍a(chǎn)業(yè)可持續(xù)發(fā)展研究[D]. 北京交通大學(xué), 2013.

    [10]譚術(shù)魁,張路,齊睿. 基于系統(tǒng)動(dòng)力學(xué)的區(qū)域耕地壓力指數(shù)研究[J]. 自然資源學(xué)報(bào), 2012, 27(5):757-765.

    [11]李華,蔡永立. 基于SD的生態(tài)安全指標(biāo)閾值的確定及應(yīng)用——以上海崇明島為例[J]. 生態(tài)學(xué)報(bào), 2010, 30(13): 3654-3664.

    [12]劉偉,陳超凡,桓漢青,等. 基于系統(tǒng)動(dòng)力學(xué)模型的低碳城市發(fā)展研究[J]. 環(huán)境污染與防治, 2014, 36(4): 86-91.

    [13]蔡運(yùn)龍,傅澤強(qiáng),戴爾阜. 區(qū)域最小人均耕地面積與耕地資源調(diào)控[J]. 地理學(xué)報(bào), 2002, 57(2): 127-134.

    (責(zé)任編輯:辜萍)

    猜你喜歡
    系統(tǒng)動(dòng)力學(xué)
    基于隨機(jī)函數(shù)Petri網(wǎng)的系統(tǒng)動(dòng)力學(xué)關(guān)聯(lián)分析模型
    大學(xué)生網(wǎng)游成癮的系統(tǒng)動(dòng)力學(xué)研究
    商(2016年32期)2016-11-24 17:43:22
    基于SD模型的文化產(chǎn)業(yè)發(fā)展科技需求仿真研究
    基于系統(tǒng)動(dòng)力學(xué)的山東省農(nóng)村消費(fèi)的研究
    商(2016年19期)2016-06-27 14:18:13
    基于系統(tǒng)動(dòng)力學(xué)的供應(yīng)鏈庫存控制系統(tǒng)
    商(2016年19期)2016-06-27 09:31:21
    基于系統(tǒng)動(dòng)力學(xué)供應(yīng)鏈的聯(lián)合庫存的研究
    商(2016年16期)2016-06-12 17:38:42
    基于系統(tǒng)動(dòng)力學(xué)的庫存管理研究
    商(2016年11期)2016-05-04 15:02:21
    產(chǎn)學(xué)研協(xié)同創(chuàng)新萌發(fā)階段研究
    商(2016年10期)2016-04-25 08:51:05
    基于系統(tǒng)動(dòng)力學(xué)的住宅能耗分析
    科技視界(2016年2期)2016-03-30 11:59:32
    聯(lián)盟組合成員多元化與企業(yè)創(chuàng)新績(jī)效之間的關(guān)系
    軟科學(xué)(2015年4期)2015-04-20 01:44:01
    欧美黄色片欧美黄色片| 久久国产精品大桥未久av| 免费观看av网站的网址| 亚洲国产av影院在线观看| 国产精品欧美亚洲77777| 99九九在线精品视频| 婷婷成人精品国产| 亚洲,欧美精品.| 亚洲欧美一区二区三区黑人 | 老熟女久久久| 亚洲av电影在线进入| 婷婷色麻豆天堂久久| 中文精品一卡2卡3卡4更新| 亚洲 欧美一区二区三区| 日韩一区二区三区影片| 久久久精品免费免费高清| 欧美xxⅹ黑人| 国产男人的电影天堂91| 男人爽女人下面视频在线观看| 18禁观看日本| 大香蕉久久网| 精品国产一区二区久久| 好男人视频免费观看在线| 如何舔出高潮| 女性生殖器流出的白浆| 久久久精品国产亚洲av高清涩受| 大片免费播放器 马上看| 日韩中文字幕欧美一区二区 | 乱人伦中国视频| 99re6热这里在线精品视频| 亚洲av欧美aⅴ国产| 18禁裸乳无遮挡动漫免费视频| 亚洲国产欧美网| 18禁动态无遮挡网站| 考比视频在线观看| 最新的欧美精品一区二区| 高清av免费在线| 一区二区三区乱码不卡18| 久久久久久久大尺度免费视频| 欧美黄色片欧美黄色片| 青草久久国产| 久久青草综合色| 国产精品一二三区在线看| 久久精品aⅴ一区二区三区四区 | 最近手机中文字幕大全| 丁香六月天网| 国产免费福利视频在线观看| 成人国语在线视频| 五月天丁香电影| 久久精品国产亚洲av天美| 久久毛片免费看一区二区三区| 国产福利在线免费观看视频| 亚洲激情五月婷婷啪啪| 国产熟女午夜一区二区三区| 在线观看美女被高潮喷水网站| 综合色丁香网| 国产亚洲一区二区精品| 99国产精品免费福利视频| 丝袜人妻中文字幕| av天堂久久9| av女优亚洲男人天堂| 在线观看人妻少妇| 亚洲成色77777| 久久久久精品久久久久真实原创| 91精品国产国语对白视频| 建设人人有责人人尽责人人享有的| 另类亚洲欧美激情| av.在线天堂| 亚洲成国产人片在线观看| 国产免费一区二区三区四区乱码| 午夜福利视频在线观看免费| 成人免费观看视频高清| 女人精品久久久久毛片| 在线 av 中文字幕| 国产成人午夜福利电影在线观看| 成人黄色视频免费在线看| 国产成人av激情在线播放| 热re99久久精品国产66热6| 亚洲国产精品999| 久久99精品国语久久久| 日日摸夜夜添夜夜爱| 在线亚洲精品国产二区图片欧美| 久久久久久久久久久免费av| 亚洲av成人精品一二三区| 国产成人a∨麻豆精品| 伊人久久大香线蕉亚洲五| 看免费成人av毛片| 亚洲色图 男人天堂 中文字幕| 久久久久久久久免费视频了| 麻豆精品久久久久久蜜桃| 18+在线观看网站| 日产精品乱码卡一卡2卡三| av一本久久久久| 女性被躁到高潮视频| 一个人免费看片子| 日韩av不卡免费在线播放| 免费av中文字幕在线| 久久久久国产精品人妻一区二区| 国产午夜精品一二区理论片| 制服诱惑二区| 黄色怎么调成土黄色| 久热久热在线精品观看| 美女脱内裤让男人舔精品视频| 桃花免费在线播放| 久久精品夜色国产| 伦精品一区二区三区| 少妇被粗大猛烈的视频| 欧美人与性动交α欧美软件| 热99久久久久精品小说推荐| 五月伊人婷婷丁香| 丰满少妇做爰视频| 久久久精品免费免费高清| 只有这里有精品99| 高清不卡的av网站| 91成人精品电影| 大香蕉久久网| 香蕉国产在线看| 精品第一国产精品| 侵犯人妻中文字幕一二三四区| 在线观看人妻少妇| 欧美日韩成人在线一区二区| 国产成人a∨麻豆精品| 涩涩av久久男人的天堂| 国产免费福利视频在线观看| 国产综合精华液| 国产一级毛片在线| 国产成人精品福利久久| 各种免费的搞黄视频| 国产精品熟女久久久久浪| 亚洲第一av免费看| 亚洲色图 男人天堂 中文字幕| av卡一久久| videos熟女内射| 亚洲av日韩在线播放| 亚洲av福利一区| 2021少妇久久久久久久久久久| 青春草视频在线免费观看| a级毛片在线看网站| 色视频在线一区二区三区| 久久久久久免费高清国产稀缺| 国产精品香港三级国产av潘金莲 | 久久国产精品大桥未久av| 可以免费在线观看a视频的电影网站 | 久久这里只有精品19| 久久 成人 亚洲| 97在线人人人人妻| 日韩av不卡免费在线播放| 色视频在线一区二区三区| 国产在线一区二区三区精| 99国产综合亚洲精品| 亚洲精品国产一区二区精华液| 九色亚洲精品在线播放| 啦啦啦在线免费观看视频4| 国产深夜福利视频在线观看| 边亲边吃奶的免费视频| 精品午夜福利在线看| 国产精品久久久久久精品电影小说| 宅男免费午夜| 午夜日本视频在线| 久久久久人妻精品一区果冻| 免费看av在线观看网站| 久久久久久久久久久免费av| 精品一区在线观看国产| 国产精品蜜桃在线观看| 大话2 男鬼变身卡| 色视频在线一区二区三区| 97在线人人人人妻| 一区二区三区精品91| 午夜免费鲁丝| 午夜福利乱码中文字幕| 菩萨蛮人人尽说江南好唐韦庄| av不卡在线播放| 国产午夜精品一二区理论片| 亚洲精品日韩在线中文字幕| 亚洲国产欧美在线一区| 亚洲国产精品成人久久小说| 久久久欧美国产精品| 久久精品久久久久久噜噜老黄| 99久久精品国产国产毛片| 国产av一区二区精品久久| 两性夫妻黄色片| 嫩草影院入口| kizo精华| 最近最新中文字幕免费大全7| 九色亚洲精品在线播放| 成人亚洲欧美一区二区av| av又黄又爽大尺度在线免费看| 少妇 在线观看| 在线天堂中文资源库| 亚洲av福利一区| 热re99久久精品国产66热6| 最近中文字幕2019免费版| 高清av免费在线| 亚洲欧美日韩另类电影网站| 最近中文字幕2019免费版| 激情五月婷婷亚洲| 亚洲图色成人| av在线观看视频网站免费| 久久国产亚洲av麻豆专区| 日本欧美视频一区| 国产精品偷伦视频观看了| 日日摸夜夜添夜夜爱| 亚洲国产看品久久| 视频区图区小说| 97人妻天天添夜夜摸| 人人妻人人爽人人添夜夜欢视频| 日韩av免费高清视频| 18禁动态无遮挡网站| 亚洲一码二码三码区别大吗| 王馨瑶露胸无遮挡在线观看| 精品人妻熟女毛片av久久网站| 又大又黄又爽视频免费| 精品少妇内射三级| 男人添女人高潮全过程视频| 久久精品久久久久久噜噜老黄| 下体分泌物呈黄色| 久久久久久久久久人人人人人人| 成人漫画全彩无遮挡| 亚洲av在线观看美女高潮| 亚洲欧美一区二区三区黑人 | 性高湖久久久久久久久免费观看| 777米奇影视久久| 啦啦啦啦在线视频资源| 边亲边吃奶的免费视频| 国语对白做爰xxxⅹ性视频网站| 两个人免费观看高清视频| 观看av在线不卡| 欧美 亚洲 国产 日韩一| 80岁老熟妇乱子伦牲交| 国产97色在线日韩免费| 蜜桃国产av成人99| 满18在线观看网站| 嫩草影院入口| 女人高潮潮喷娇喘18禁视频| 熟女电影av网| 欧美精品av麻豆av| 亚洲欧美色中文字幕在线| 国产精品一国产av| 九色亚洲精品在线播放| 中文字幕另类日韩欧美亚洲嫩草| 人人妻人人添人人爽欧美一区卜| 18+在线观看网站| 女人精品久久久久毛片| 成年av动漫网址| 伊人亚洲综合成人网| 啦啦啦视频在线资源免费观看| 一区福利在线观看| 又大又黄又爽视频免费| 丰满饥渴人妻一区二区三| 你懂的网址亚洲精品在线观看| 天美传媒精品一区二区| 国产日韩欧美在线精品| 日韩制服骚丝袜av| 香蕉精品网在线| 亚洲精品aⅴ在线观看| 毛片一级片免费看久久久久| 国产毛片在线视频| www.精华液| 新久久久久国产一级毛片| 日本av免费视频播放| 国产精品嫩草影院av在线观看| 久久精品国产亚洲av涩爱| 亚洲精品国产一区二区精华液| 韩国精品一区二区三区| 国产 一区精品| 久久久久网色| 国产精品av久久久久免费| 久久久久人妻精品一区果冻| 国产av一区二区精品久久| 女人精品久久久久毛片| 丝瓜视频免费看黄片| 日韩一本色道免费dvd| 久久午夜福利片| 久久亚洲国产成人精品v| 欧美精品国产亚洲| 国产精品麻豆人妻色哟哟久久| 欧美人与性动交α欧美精品济南到 | 国产亚洲av片在线观看秒播厂| 国产精品成人在线| 亚洲欧美一区二区三区国产| 边亲边吃奶的免费视频| 色婷婷久久久亚洲欧美| 国产av精品麻豆| 青青草视频在线视频观看| 巨乳人妻的诱惑在线观看| 国产麻豆69| 国产精品一二三区在线看| 男女下面插进去视频免费观看| 国产成人a∨麻豆精品| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产精品三级大全| 高清在线视频一区二区三区| 国产欧美日韩综合在线一区二区| 另类精品久久| 18在线观看网站| 国产老妇伦熟女老妇高清| 成人亚洲欧美一区二区av| 精品少妇久久久久久888优播| 成年女人在线观看亚洲视频| 日韩三级伦理在线观看| 亚洲少妇的诱惑av| 亚洲av成人精品一二三区| 婷婷色av中文字幕| av在线观看视频网站免费| 水蜜桃什么品种好| 纵有疾风起免费观看全集完整版| 国产老妇伦熟女老妇高清| 精品国产一区二区三区久久久樱花| 久久久久久久精品精品| 久久av网站| 久久久精品94久久精品| 日韩精品免费视频一区二区三区| 日韩中文字幕视频在线看片| 男女午夜视频在线观看| 日本wwww免费看| 国产精品蜜桃在线观看| 欧美日韩精品成人综合77777| 三级国产精品片| 边亲边吃奶的免费视频| av女优亚洲男人天堂| 狠狠婷婷综合久久久久久88av| 精品少妇内射三级| 成人亚洲精品一区在线观看| 亚洲成人手机| 精品少妇黑人巨大在线播放| 美女xxoo啪啪120秒动态图| 亚洲精品久久成人aⅴ小说| 色吧在线观看| 极品少妇高潮喷水抽搐| 国产97色在线日韩免费| 天天躁狠狠躁夜夜躁狠狠躁| 成人国语在线视频| 午夜影院在线不卡| 国产一区亚洲一区在线观看| 涩涩av久久男人的天堂| 又大又黄又爽视频免费| 熟妇人妻不卡中文字幕| 精品一区二区三区四区五区乱码 | 午夜激情久久久久久久| 精品人妻一区二区三区麻豆| 七月丁香在线播放| 精品人妻一区二区三区麻豆| 最近最新中文字幕免费大全7| 黑人巨大精品欧美一区二区蜜桃| 久久久久精品久久久久真实原创| 国产爽快片一区二区三区| 免费观看a级毛片全部| 日韩中文字幕欧美一区二区 | 黑丝袜美女国产一区| 美女福利国产在线| 五月开心婷婷网| 免费观看性生交大片5| 免费在线观看完整版高清| av天堂久久9| 久久久精品免费免费高清| 日韩av免费高清视频| 一级爰片在线观看| 777米奇影视久久| 午夜影院在线不卡| 免费久久久久久久精品成人欧美视频| 巨乳人妻的诱惑在线观看| 国产一区亚洲一区在线观看| 精品视频人人做人人爽| www.熟女人妻精品国产| 另类亚洲欧美激情| 美女xxoo啪啪120秒动态图| 欧美激情 高清一区二区三区| av天堂久久9| 亚洲精华国产精华液的使用体验| 久久综合国产亚洲精品| 成年人免费黄色播放视频| 最近中文字幕高清免费大全6| 不卡视频在线观看欧美| 亚洲精品av麻豆狂野| 国产精品二区激情视频| 亚洲精品久久午夜乱码| 国产一区二区三区综合在线观看| av国产精品久久久久影院| 亚洲久久久国产精品| 国产精品女同一区二区软件| 国产欧美亚洲国产| 国产免费一区二区三区四区乱码| 你懂的网址亚洲精品在线观看| 婷婷色综合www| 精品国产露脸久久av麻豆| 精品少妇内射三级| 肉色欧美久久久久久久蜜桃| 亚洲av免费高清在线观看| 丝袜美腿诱惑在线| 一本大道久久a久久精品| 精品国产一区二区三区久久久樱花| 日韩欧美精品免费久久| 国产精品av久久久久免费| 婷婷色综合www| 国产免费福利视频在线观看| 国产成人aa在线观看| 男的添女的下面高潮视频| 男人爽女人下面视频在线观看| 成人毛片a级毛片在线播放| 国产精品久久久久成人av| 免费女性裸体啪啪无遮挡网站| 2021少妇久久久久久久久久久| 亚洲av日韩在线播放| 丝袜脚勾引网站| 亚洲精品视频女| 亚洲五月色婷婷综合| 这个男人来自地球电影免费观看 | kizo精华| 亚洲 欧美一区二区三区| 女人久久www免费人成看片| 国产免费视频播放在线视频| 在线天堂最新版资源| 超碰成人久久| 交换朋友夫妻互换小说| 建设人人有责人人尽责人人享有的| 国产精品无大码| 欧美亚洲日本最大视频资源| 欧美日韩精品成人综合77777| 成年av动漫网址| 久久久久久伊人网av| 十八禁网站网址无遮挡| 久久国产亚洲av麻豆专区| 免费高清在线观看日韩| 五月天丁香电影| 久久久亚洲精品成人影院| 一二三四在线观看免费中文在| 亚洲一区二区三区欧美精品| 久热久热在线精品观看| 色婷婷av一区二区三区视频| 久久狼人影院| 国产亚洲最大av| 熟女av电影| 啦啦啦中文免费视频观看日本| 成人亚洲精品一区在线观看| 不卡av一区二区三区| 永久免费av网站大全| 少妇熟女欧美另类| 成人黄色视频免费在线看| 九九爱精品视频在线观看| 韩国精品一区二区三区| 久久久久精品久久久久真实原创| 性色av一级| av有码第一页| 看免费成人av毛片| 久久久久人妻精品一区果冻| 成人午夜精彩视频在线观看| 久久午夜福利片| 一区二区日韩欧美中文字幕| 国产97色在线日韩免费| 亚洲美女搞黄在线观看| 日本午夜av视频| 天堂俺去俺来也www色官网| 亚洲国产精品一区二区三区在线| 色播在线永久视频| 99国产精品免费福利视频| 成人影院久久| 少妇被粗大的猛进出69影院| 各种免费的搞黄视频| 十八禁高潮呻吟视频| 制服诱惑二区| 黄网站色视频无遮挡免费观看| 电影成人av| 五月天丁香电影| 肉色欧美久久久久久久蜜桃| 好男人视频免费观看在线| 成人免费观看视频高清| 色视频在线一区二区三区| 18+在线观看网站| 一级毛片电影观看| 日韩欧美一区视频在线观看| 国产精品一国产av| 久久精品久久久久久久性| 最近最新中文字幕大全免费视频 | h视频一区二区三区| 成人国产麻豆网| 另类精品久久| freevideosex欧美| 亚洲欧美一区二区三区国产| 纵有疾风起免费观看全集完整版| 国产精品久久久久成人av| 亚洲成色77777| 中文字幕av电影在线播放| 亚洲第一青青草原| 搡女人真爽免费视频火全软件| 欧美 日韩 精品 国产| 亚洲色图 男人天堂 中文字幕| 三级国产精品片| 电影成人av| a 毛片基地| 香蕉丝袜av| videossex国产| 国产在线视频一区二区| 国产日韩一区二区三区精品不卡| 欧美av亚洲av综合av国产av | 自拍欧美九色日韩亚洲蝌蚪91| 国产精品蜜桃在线观看| 制服诱惑二区| 亚洲精品在线美女| 人成视频在线观看免费观看| 成年人午夜在线观看视频| 男女免费视频国产| 久久 成人 亚洲| 久久鲁丝午夜福利片| 看免费av毛片| 国产成人午夜福利电影在线观看| 亚洲精品国产一区二区精华液| 成年动漫av网址| 男女边吃奶边做爰视频| 十八禁高潮呻吟视频| 菩萨蛮人人尽说江南好唐韦庄| 熟妇人妻不卡中文字幕| 人人妻人人澡人人看| 精品99又大又爽又粗少妇毛片| 亚洲美女黄色视频免费看| 国产一级毛片在线| av免费观看日本| 丰满饥渴人妻一区二区三| 青草久久国产| 三级国产精品片| 日韩av在线免费看完整版不卡| 一二三四在线观看免费中文在| a级毛片在线看网站| 777久久人妻少妇嫩草av网站| 久久精品国产亚洲av高清一级| 在线 av 中文字幕| 久久久久久久久久人人人人人人| 9191精品国产免费久久| 国产在线免费精品| 色视频在线一区二区三区| 五月开心婷婷网| 丝袜美足系列| 99国产综合亚洲精品| 99久久精品国产国产毛片| 亚洲欧美精品自产自拍| 亚洲综合精品二区| 久久国产精品大桥未久av| 午夜老司机福利剧场| 最近中文字幕2019免费版| 人妻人人澡人人爽人人| 肉色欧美久久久久久久蜜桃| 女性被躁到高潮视频| 亚洲中文av在线| 青青草视频在线视频观看| 高清视频免费观看一区二区| xxxhd国产人妻xxx| 18禁国产床啪视频网站| 搡老乐熟女国产| 亚洲欧美成人精品一区二区| 亚洲精品久久午夜乱码| 国产一区二区三区综合在线观看| 免费观看性生交大片5| 国产欧美亚洲国产| 卡戴珊不雅视频在线播放| 亚洲成人av在线免费| 黄色怎么调成土黄色| 亚洲在久久综合| 美女福利国产在线| 一级毛片黄色毛片免费观看视频| 日日爽夜夜爽网站| 亚洲成人av在线免费| 亚洲,一卡二卡三卡| 国产成人精品福利久久| 亚洲婷婷狠狠爱综合网| 免费女性裸体啪啪无遮挡网站| 丝袜脚勾引网站| 视频在线观看一区二区三区| av不卡在线播放| 电影成人av| 欧美成人午夜精品| 欧美日韩视频精品一区| 黑人巨大精品欧美一区二区蜜桃| 色吧在线观看| 午夜免费观看性视频| 黑人欧美特级aaaaaa片| 美女高潮到喷水免费观看| 一级毛片我不卡| 亚洲精品aⅴ在线观看| 黄网站色视频无遮挡免费观看| 少妇被粗大猛烈的视频| 国产亚洲av片在线观看秒播厂| 亚洲内射少妇av| 欧美黄色片欧美黄色片| 啦啦啦视频在线资源免费观看| 久久女婷五月综合色啪小说| 国产精品久久久久成人av| 久久久亚洲精品成人影院| 亚洲欧美日韩另类电影网站| 亚洲熟女精品中文字幕| 99久国产av精品国产电影| 日韩,欧美,国产一区二区三区| 青青草视频在线视频观看| 国产免费现黄频在线看| 王馨瑶露胸无遮挡在线观看| 制服丝袜香蕉在线| 欧美成人午夜精品| 视频区图区小说| 国产色婷婷99| 亚洲av综合色区一区| 亚洲av成人精品一二三区| 久久狼人影院| 久久久久人妻精品一区果冻| 在线天堂中文资源库| 亚洲欧美清纯卡通| 美女大奶头黄色视频| 精品少妇一区二区三区视频日本电影 | 国产精品偷伦视频观看了| 九九爱精品视频在线观看| 亚洲,欧美精品.| 久久韩国三级中文字幕| 亚洲av日韩在线播放| 建设人人有责人人尽责人人享有的| 1024视频免费在线观看| 亚洲欧美日韩另类电影网站| 欧美日韩一级在线毛片| 国产精品久久久av美女十八|