• <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天堂久久9| 性少妇av在线| 国产不卡av网站在线观看| 亚洲欧美一区二区三区久久| 久久99热这里只频精品6学生| av有码第一页| 十八禁网站网址无遮挡| 亚洲少妇的诱惑av| 精品高清国产在线一区| 国产人伦9x9x在线观看| av一本久久久久| 久久国产精品人妻蜜桃| 伊人亚洲综合成人网| 欧美黄色片欧美黄色片| 丁香六月天网| 不卡av一区二区三区| 久久久久久久精品精品| 国产成人啪精品午夜网站| 日本撒尿小便嘘嘘汇集6| 精品国产乱码久久久久久男人| 精品乱码久久久久久99久播| 精品人妻1区二区| 亚洲成人免费电影在线观看| 亚洲欧美清纯卡通| 日本av手机在线免费观看| 亚洲国产欧美在线一区| 国产成人一区二区三区免费视频网站| 国产熟女午夜一区二区三区| 久久人妻福利社区极品人妻图片| 久久久欧美国产精品| 久久精品aⅴ一区二区三区四区| 在线亚洲精品国产二区图片欧美| 国产精品熟女久久久久浪| 水蜜桃什么品种好| 91精品伊人久久大香线蕉| 制服人妻中文乱码| 亚洲成国产人片在线观看| 菩萨蛮人人尽说江南好唐韦庄| 日本一区二区免费在线视频| 51午夜福利影视在线观看| 日韩三级视频一区二区三区| 一本一本久久a久久精品综合妖精| 国产免费av片在线观看野外av| 久久久国产精品麻豆| 欧美av亚洲av综合av国产av| 亚洲精品一区蜜桃| 精品少妇一区二区三区视频日本电影| 真人做人爱边吃奶动态| 国产片内射在线| 精品第一国产精品| 成人国语在线视频| tocl精华| 男女免费视频国产| 亚洲av欧美aⅴ国产| 久久人妻熟女aⅴ| 国产精品.久久久| 五月开心婷婷网| 国产免费现黄频在线看| 国产精品一区二区在线观看99| 亚洲av电影在线进入| 亚洲视频免费观看视频| 久久精品熟女亚洲av麻豆精品| 亚洲精品久久成人aⅴ小说| 国产日韩欧美在线精品| 十八禁高潮呻吟视频| 美女高潮喷水抽搐中文字幕| 男女之事视频高清在线观看| av福利片在线| 亚洲欧洲精品一区二区精品久久久| 看免费av毛片| 亚洲精品一卡2卡三卡4卡5卡 | 欧美精品啪啪一区二区三区 | 精品一区二区三卡| 午夜福利一区二区在线看| 一区福利在线观看| 亚洲国产看品久久| 男女床上黄色一级片免费看| 女人久久www免费人成看片| www.999成人在线观看| 午夜福利在线免费观看网站| 叶爱在线成人免费视频播放| 国产有黄有色有爽视频| 69av精品久久久久久 | a级毛片在线看网站| 一区二区三区精品91| 中文字幕人妻熟女乱码| 中文欧美无线码| 精品福利永久在线观看| 麻豆av在线久日| 亚洲av日韩在线播放| 男女边摸边吃奶| av超薄肉色丝袜交足视频| 美女高潮喷水抽搐中文字幕| 黄片小视频在线播放| 成人免费观看视频高清| 亚洲精品日韩在线中文字幕| 日韩一卡2卡3卡4卡2021年| 精品欧美一区二区三区在线| 亚洲免费av在线视频| 男女国产视频网站| 精品国产乱子伦一区二区三区 | 日韩精品免费视频一区二区三区| 国产精品国产av在线观看| 一二三四社区在线视频社区8| 男女床上黄色一级片免费看| 久久久欧美国产精品| a级片在线免费高清观看视频| 中文字幕人妻丝袜一区二区| 国产免费av片在线观看野外av| 91麻豆精品激情在线观看国产 | 香蕉国产在线看| 91麻豆av在线| 亚洲av片天天在线观看| 精品福利观看| 动漫黄色视频在线观看| 日韩欧美国产一区二区入口| 欧美变态另类bdsm刘玥| 黄色视频不卡| www.av在线官网国产| 国产三级黄色录像| videosex国产| 亚洲av日韩在线播放| 国产成人精品久久二区二区91| 中文字幕另类日韩欧美亚洲嫩草| 91av网站免费观看| 人人妻人人澡人人爽人人夜夜| 丝袜美腿诱惑在线| 女性生殖器流出的白浆| 亚洲九九香蕉| 国产精品久久久久久人妻精品电影 | 日韩制服骚丝袜av| 久久人人爽av亚洲精品天堂| 久久国产精品影院| 91字幕亚洲| 日韩精品免费视频一区二区三区| 男女床上黄色一级片免费看| 亚洲av国产av综合av卡| 午夜激情久久久久久久| 老司机影院成人| 免费少妇av软件| 男女之事视频高清在线观看| 三级毛片av免费| 青春草亚洲视频在线观看| 亚洲专区国产一区二区| 青春草视频在线免费观看| 精品国产一区二区三区久久久樱花| 视频区图区小说| 9色porny在线观看| 美女国产高潮福利片在线看| 男人舔女人的私密视频| 99国产极品粉嫩在线观看| 久久中文看片网| 国产精品av久久久久免费| 欧美精品一区二区免费开放| 日本猛色少妇xxxxx猛交久久| 天天躁日日躁夜夜躁夜夜| 成人国产一区最新在线观看| 久久久久精品人妻al黑| 91国产中文字幕| 男女高潮啪啪啪动态图| av免费在线观看网站| 久久亚洲精品不卡| av网站免费在线观看视频| 各种免费的搞黄视频| 人妻人人澡人人爽人人| 国产野战对白在线观看| 成在线人永久免费视频| 午夜老司机福利片| 久久国产亚洲av麻豆专区| 亚洲国产欧美日韩在线播放| 人人妻人人爽人人添夜夜欢视频| 亚洲一区二区三区欧美精品| 高清在线国产一区| 国产亚洲一区二区精品| 亚洲第一青青草原| 两性夫妻黄色片| 美女主播在线视频| 亚洲精品国产av成人精品| 91成人精品电影| videosex国产| 国产精品欧美亚洲77777| 纯流量卡能插随身wifi吗| 满18在线观看网站| 亚洲中文日韩欧美视频| 国产成人系列免费观看| 亚洲熟女精品中文字幕| 国产精品免费视频内射| 午夜福利免费观看在线| svipshipincom国产片| 午夜视频精品福利| 色94色欧美一区二区| 成人手机av| 久久国产精品人妻蜜桃| 亚洲欧美日韩另类电影网站| 亚洲视频免费观看视频| 久久久久久久久久久久大奶| 久久久国产一区二区| 啦啦啦视频在线资源免费观看| 国产亚洲一区二区精品| 久久久久精品人妻al黑| 亚洲人成电影观看| 首页视频小说图片口味搜索| 国产成人免费无遮挡视频| 一级黄色大片毛片| 国产在线一区二区三区精| 久久精品国产亚洲av香蕉五月 | 国产精品影院久久| 久久天堂一区二区三区四区| 亚洲美女黄色视频免费看| 啦啦啦视频在线资源免费观看| 亚洲少妇的诱惑av| 99国产精品一区二区蜜桃av | 电影成人av| 淫妇啪啪啪对白视频 | 亚洲av日韩在线播放| 亚洲欧洲精品一区二区精品久久久| 97人妻天天添夜夜摸| 搡老岳熟女国产| 国产欧美亚洲国产| 两性夫妻黄色片| 国产日韩一区二区三区精品不卡| 老熟妇乱子伦视频在线观看 | 日韩欧美免费精品| 老司机午夜十八禁免费视频| 一本综合久久免费| av线在线观看网站| 日韩三级视频一区二区三区| 男女下面插进去视频免费观看| av片东京热男人的天堂| 2018国产大陆天天弄谢| 女性被躁到高潮视频| 国产亚洲精品第一综合不卡| 国产免费现黄频在线看| 欧美成狂野欧美在线观看| 各种免费的搞黄视频| 亚洲欧美日韩高清在线视频 | 99九九在线精品视频| 黑人巨大精品欧美一区二区mp4| 新久久久久国产一级毛片| 久久久久精品人妻al黑| 久久99一区二区三区| e午夜精品久久久久久久| 精品一区二区三区av网在线观看 | 国产男人的电影天堂91| 亚洲色图综合在线观看| av不卡在线播放| 老熟女久久久| 1024香蕉在线观看| 国内毛片毛片毛片毛片毛片| 国产成人av教育| 美女扒开内裤让男人捅视频| 亚洲 欧美一区二区三区| 欧美人与性动交α欧美软件| 妹子高潮喷水视频| 丝袜喷水一区| 老汉色av国产亚洲站长工具| 国产一区二区三区在线臀色熟女 | 久久这里只有精品19| 精品国产国语对白av| 91麻豆精品激情在线观看国产 | 美女视频免费永久观看网站| 男女边摸边吃奶| 国产日韩欧美亚洲二区| 啦啦啦中文免费视频观看日本| 国产av一区二区精品久久| 亚洲 欧美一区二区三区| 交换朋友夫妻互换小说| 夜夜骑夜夜射夜夜干| 啪啪无遮挡十八禁网站| 夫妻午夜视频| av视频免费观看在线观看| 亚洲精品在线美女| 亚洲欧美日韩另类电影网站| 成年人午夜在线观看视频| 亚洲一区中文字幕在线| 少妇猛男粗大的猛烈进出视频| 在线观看免费午夜福利视频| 一二三四社区在线视频社区8| 国产日韩欧美亚洲二区| 午夜免费成人在线视频| 老熟女久久久| 在线亚洲精品国产二区图片欧美| 巨乳人妻的诱惑在线观看| 亚洲av日韩精品久久久久久密| 美女午夜性视频免费| 国产精品一区二区在线观看99| 视频区图区小说| 青青草视频在线视频观看| 国产一区二区 视频在线| 老汉色∧v一级毛片| 国产在线视频一区二区| 亚洲伊人久久精品综合| 在线观看舔阴道视频| 国产成人系列免费观看| 亚洲一区二区三区欧美精品| 欧美少妇被猛烈插入视频| 天堂俺去俺来也www色官网| www.999成人在线观看| 欧美日韩一级在线毛片| 一本大道久久a久久精品| 日本欧美视频一区| a在线观看视频网站| 男男h啪啪无遮挡| 天天影视国产精品| 1024视频免费在线观看| 国产片内射在线| 亚洲精品中文字幕在线视频| 性高湖久久久久久久久免费观看| 亚洲免费av在线视频| 国产成人精品久久二区二区免费| 国产日韩欧美亚洲二区| 交换朋友夫妻互换小说| 精品国内亚洲2022精品成人 | 中文字幕另类日韩欧美亚洲嫩草| 日本猛色少妇xxxxx猛交久久| 黄网站色视频无遮挡免费观看| 亚洲精品美女久久久久99蜜臀| 日韩一卡2卡3卡4卡2021年| 欧美人与性动交α欧美软件| 欧美老熟妇乱子伦牲交| 精品熟女少妇八av免费久了| 夜夜骑夜夜射夜夜干| 久久天堂一区二区三区四区| 国产欧美日韩综合在线一区二区| 欧美激情久久久久久爽电影 | 国产在线观看jvid| 又大又爽又粗| 男女高潮啪啪啪动态图| 一级黄色大片毛片| 人人妻人人爽人人添夜夜欢视频| 久久久国产成人免费| 午夜福利视频精品| 国产日韩欧美在线精品| 久久精品成人免费网站| 精品少妇久久久久久888优播| 人人妻,人人澡人人爽秒播| 成人免费观看视频高清| 国产精品影院久久| 性高湖久久久久久久久免费观看| 亚洲免费av在线视频| 又黄又粗又硬又大视频| 国产在线一区二区三区精| 亚洲九九香蕉| 亚洲综合色网址| 国产日韩欧美视频二区| 大香蕉久久网| 美女扒开内裤让男人捅视频| 亚洲精品国产区一区二| 亚洲av电影在线进入| 久久热在线av| 国产精品久久久av美女十八| 国产黄色免费在线视频| 大片免费播放器 马上看| 日韩精品免费视频一区二区三区| 久9热在线精品视频| 99国产精品一区二区三区| 午夜福利,免费看| 久久久久久亚洲精品国产蜜桃av| 97精品久久久久久久久久精品| 欧美日本中文国产一区发布| 久久久久视频综合| 超碰97精品在线观看| 国产亚洲精品久久久久5区| 人人妻人人添人人爽欧美一区卜| 嫩草影视91久久| 中国美女看黄片| 夜夜骑夜夜射夜夜干| 久久人妻福利社区极品人妻图片| 国产亚洲精品一区二区www | 水蜜桃什么品种好| 国产黄色免费在线视频| 超碰97精品在线观看| 欧美日韩视频精品一区| 成人av一区二区三区在线看 | 丝袜喷水一区| 大片电影免费在线观看免费| 久热这里只有精品99| 亚洲国产中文字幕在线视频| 久久久欧美国产精品| 欧美另类一区| 免费在线观看完整版高清| 制服人妻中文乱码| 99国产精品免费福利视频| 黄片播放在线免费| 亚洲精品久久成人aⅴ小说| svipshipincom国产片| 国产精品国产av在线观看| 男女高潮啪啪啪动态图| 99re6热这里在线精品视频| 黄色a级毛片大全视频| 国产精品国产av在线观看| 久久影院123| 国产精品国产av在线观看| 国产日韩一区二区三区精品不卡| 精品少妇黑人巨大在线播放| 日日摸夜夜添夜夜添小说| 久久久久精品人妻al黑| 美女扒开内裤让男人捅视频| av在线老鸭窝| 亚洲人成电影免费在线| 精品少妇一区二区三区视频日本电影| 亚洲国产精品一区二区三区在线| 天天躁日日躁夜夜躁夜夜| 欧美日韩中文字幕国产精品一区二区三区 | 两个人看的免费小视频| 精品国产乱码久久久久久小说| 久久热在线av| 日韩电影二区| 两性夫妻黄色片| 色播在线永久视频| 国产欧美日韩一区二区三 | 欧美精品一区二区免费开放| 欧美激情久久久久久爽电影 | 国产亚洲一区二区精品| 国产男女超爽视频在线观看| 国产高清国产精品国产三级| 亚洲专区中文字幕在线| 免费女性裸体啪啪无遮挡网站| 久久人人爽av亚洲精品天堂| 大码成人一级视频| 一边摸一边做爽爽视频免费| 国产在线免费精品| 精品一区在线观看国产| 老熟女久久久| 免费看十八禁软件| 大片免费播放器 马上看| 热99久久久久精品小说推荐| 少妇人妻久久综合中文| 激情视频va一区二区三区| 男女下面插进去视频免费观看| 一级毛片女人18水好多| 美女午夜性视频免费| tocl精华| 中文字幕精品免费在线观看视频| h视频一区二区三区| 免费观看人在逋| 法律面前人人平等表现在哪些方面 | 亚洲 欧美一区二区三区| 99久久综合免费| 高潮久久久久久久久久久不卡| 国产精品一二三区在线看| 精品人妻熟女毛片av久久网站| 成人三级做爰电影| 久久精品国产a三级三级三级| 国产av又大| 色老头精品视频在线观看| 国产精品.久久久| 国产精品免费大片| 精品亚洲成国产av| 色精品久久人妻99蜜桃| 精品第一国产精品| 成人av一区二区三区在线看 | 亚洲av电影在线观看一区二区三区| 黄网站色视频无遮挡免费观看| 十八禁网站免费在线| 色精品久久人妻99蜜桃| 啦啦啦在线免费观看视频4| videos熟女内射| 久久久久久久久久久久大奶| 丁香六月欧美| 成年动漫av网址| 亚洲av国产av综合av卡| 国产av国产精品国产| 极品少妇高潮喷水抽搐| 永久免费av网站大全| 久久人人爽av亚洲精品天堂| 天天躁夜夜躁狠狠躁躁| 黄频高清免费视频| 大片电影免费在线观看免费| 国产成人欧美| 国产一区二区激情短视频 | 国产精品国产三级国产专区5o| 欧美日本中文国产一区发布| 中文字幕色久视频| 女人被躁到高潮嗷嗷叫费观| 91av网站免费观看| 欧美国产精品一级二级三级| 久久久久久人人人人人| 国产亚洲精品久久久久5区| 9色porny在线观看| 涩涩av久久男人的天堂| 这个男人来自地球电影免费观看| av国产精品久久久久影院| 欧美精品高潮呻吟av久久| 亚洲一卡2卡3卡4卡5卡精品中文| 精品少妇内射三级| 人妻 亚洲 视频| 99精品久久久久人妻精品| av免费在线观看网站| 99精国产麻豆久久婷婷| 两个人看的免费小视频| 啦啦啦免费观看视频1| 欧美人与性动交α欧美精品济南到| 亚洲欧美清纯卡通| 99re6热这里在线精品视频| √禁漫天堂资源中文www| 久久精品国产a三级三级三级| 亚洲一区二区三区欧美精品| 天天躁日日躁夜夜躁夜夜| 少妇被粗大的猛进出69影院| 国产97色在线日韩免费| 亚洲精品第二区| 欧美日韩亚洲综合一区二区三区_| 一级片'在线观看视频| 国产欧美日韩一区二区三区在线| 中国美女看黄片| 国产成+人综合+亚洲专区| 老汉色av国产亚洲站长工具| 国产激情久久老熟女| 久久久精品区二区三区| 51午夜福利影视在线观看| 正在播放国产对白刺激| 女人高潮潮喷娇喘18禁视频| 久久性视频一级片| 在线十欧美十亚洲十日本专区| 国产一区二区激情短视频 | 12—13女人毛片做爰片一| 中文字幕另类日韩欧美亚洲嫩草| 成人18禁高潮啪啪吃奶动态图| 国产精品免费视频内射| 亚洲成av片中文字幕在线观看| 国产av又大| 真人做人爱边吃奶动态| 热99久久久久精品小说推荐| 电影成人av| 人人妻人人添人人爽欧美一区卜| 亚洲国产看品久久| 99热国产这里只有精品6| 欧美精品亚洲一区二区| 久久女婷五月综合色啪小说| 亚洲av电影在线观看一区二区三区| 巨乳人妻的诱惑在线观看| 一区二区av电影网| 国产精品免费视频内射| 国内毛片毛片毛片毛片毛片| 中文字幕最新亚洲高清| 伊人久久大香线蕉亚洲五| 免费看十八禁软件| 91九色精品人成在线观看| 少妇被粗大的猛进出69影院| 亚洲av男天堂| 激情视频va一区二区三区| 交换朋友夫妻互换小说| 99久久综合免费| 欧美黄色片欧美黄色片| 国产精品秋霞免费鲁丝片| 国产亚洲一区二区精品| 满18在线观看网站| 中文字幕高清在线视频| 一区福利在线观看| 日本撒尿小便嘘嘘汇集6| 黄网站色视频无遮挡免费观看| 亚洲中文字幕日韩| 国产免费福利视频在线观看| 国产av国产精品国产| 国产一卡二卡三卡精品| 夫妻午夜视频| 亚洲视频免费观看视频| 色综合欧美亚洲国产小说| 男人添女人高潮全过程视频| av又黄又爽大尺度在线免费看| 日本五十路高清| 久久香蕉激情| 亚洲欧美激情在线| 人成视频在线观看免费观看| 一级毛片女人18水好多| 国产亚洲精品第一综合不卡| 美女主播在线视频| 精品少妇一区二区三区视频日本电影| 91av网站免费观看| 久久国产亚洲av麻豆专区| 宅男免费午夜| 法律面前人人平等表现在哪些方面 | 18禁观看日本| 香蕉国产在线看| 少妇粗大呻吟视频| 嫁个100分男人电影在线观看| 欧美成人午夜精品| 亚洲精品久久久久久婷婷小说| 午夜激情av网站| 免费女性裸体啪啪无遮挡网站| 日日夜夜操网爽| 亚洲成av片中文字幕在线观看| 国产精品一区二区在线不卡| 亚洲第一av免费看| 成人国语在线视频| 纯流量卡能插随身wifi吗| 制服诱惑二区| 性色av一级| 亚洲美女黄色视频免费看| 亚洲欧美日韩另类电影网站| 久久天躁狠狠躁夜夜2o2o| 免费日韩欧美在线观看| 国产亚洲精品久久久久5区| 9191精品国产免费久久| 欧美激情极品国产一区二区三区| 2018国产大陆天天弄谢| 热99久久久久精品小说推荐| 欧美另类一区| 成年动漫av网址| 高清av免费在线| 美女午夜性视频免费| 大码成人一级视频| 欧美97在线视频| 国产激情久久老熟女| 一区福利在线观看| 成人影院久久| 人人妻人人爽人人添夜夜欢视频| 国产1区2区3区精品| 欧美人与性动交α欧美精品济南到| 建设人人有责人人尽责人人享有的| a在线观看视频网站|