999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于DEA和Tobit模型的安徽省水資源利用效率測(cè)算與影響因素分析

2024-06-23 21:45:00李寶春
安徽農(nóng)業(yè)科學(xué) 2024年11期

摘要 利用超效率SBM-DEA模型、DEA-Malmquist 指數(shù)對(duì)2011—2020年安徽省水資源利用效率進(jìn)行測(cè)算,并采用Tobit模型分析其影響因素。結(jié)果表明,2011—2020年安徽省水資源利用效率未達(dá)到有效水平,利用效率平均值從0.931波動(dòng)下降為0.873,16個(gè)地市中僅有合肥、淮北、黃山的水資源利用效率達(dá)到DEA有效,其余13個(gè)城市的水資源利用效率均處于生產(chǎn)前沿面之下。安徽省水資源利用效率Malmquist 指數(shù)平均值為0.981,總體處于下降階段,技術(shù)效率變化、技術(shù)進(jìn)步變化均阻礙了效率的上升;合肥、淮北、蕪湖、銅陵、安慶和黃山的Malmquist 指數(shù)平均值處于上升階段,其余10個(gè)城市的Malmquist 指數(shù)值均處于下降階段。人均GDP、第一產(chǎn)業(yè)比重、農(nóng)業(yè)用水比重因素均對(duì)安徽省水資源利用效率的提升呈現(xiàn)出顯著的正向促進(jìn)作用;而人均用水量、人均污水排放量則均呈現(xiàn)出顯著負(fù)向作用。

關(guān)鍵詞 水資源利用效率;DEA模型;Malmquist 指數(shù);Tobit模型;安徽省

中圖分類號(hào) TV213? 文獻(xiàn)標(biāo)識(shí)碼 A? 文章編號(hào) 0517-6611(2024)11-0183-04

doi:10.3969/j.issn.0517-6611.2024.11.039

Estimation and Analysis of Influencing Factors of Water Resource Use Efficiency in Anhui Province Based on DEA and Tobit Models

LI Bao-chun

(Anhui and Huaihe River Institute of Hydraulic Research,Hefei,Anhui 230088)

Abstract The superefficiency SBM-DEA model and DEA-Malmquist index were used to measure the water resource use efficiency in Anhui Province from 2011 to 2020, and the Tobit model was used to analyze the influencing factors.The result showed that water resource use efficiency in Anhui Province has not reached the effective level, and the average value of water resource use efficiency fluctuates from 0.931 to 0.873,among the 16 cities, only Hefei, Huaibei and Huangshan achieve DEA efficiency in water use,and the remaining 13 cities, the water resource use efficiency is below the production frontier.The average Malmquist index of water resource use efficiency in Anhui Province is 0.981, which is generally in a declining stage,the change of technical efficiency and the change of technological progress have hindered the increase of water resources utilization efficiency in Anhui Province.The average Malmquist index of water resource use efficiency in Hefei, Huaibei, Wuhu, Tongling, Anqing and Huangshan are on the rise, showing a good growth trend, the remaining 10 cities are in the declining stage.The per capita GDP, the proportion of primary industry and the proportion of agricultural water use all have significant positive promoting effects on the improvement of water use efficiency in Anhui Province,both per capita water consumption and per capita sewage discharge have significant negative effects on water resources utilization efficiency.

Key words Water resource use efficiency;DEA model;Malmquist index;Tobit model;Anhui Province

作者簡(jiǎn)介 李寶春(1978—),男,安徽六安人,工程師,從事水利工程設(shè)計(jì)研究。

收稿日期 2023-11-03

伴隨著社會(huì)的發(fā)展,人類對(duì)水資源的需求量不斷增加,如何更加高效利用水資源已成為社會(huì)經(jīng)濟(jì)增長(zhǎng)和人們生活質(zhì)量協(xié)調(diào)可持續(xù)發(fā)展的重要內(nèi)容[1]。安徽省地處中國(guó)東部,屬暖溫帶與亞熱帶過(guò)渡地區(qū),水資源較為緊缺,且南北空間分布差異較大。近年來(lái),隨著安徽省城鎮(zhèn)化進(jìn)程的加快和工業(yè)化程度的快速提高,水資源需求不斷增長(zhǎng),再加上水利工程設(shè)施老化、水資源調(diào)蓄能力不足、灌溉用水效率低下等問(wèn)題凸顯,導(dǎo)致安徽省水資源利用效率偏低[2]。在此背景下,對(duì)安徽省水資源利用效率及其影響因素進(jìn)行分析研究,對(duì)于構(gòu)建節(jié)水型社會(huì)、推動(dòng)安徽省經(jīng)濟(jì)社會(huì)可持續(xù)發(fā)展具有重要意義。

水資源利用效率評(píng)價(jià)作為學(xué)界的研究熱點(diǎn),備受學(xué)者關(guān)注,目前在水資源利用效率的測(cè)算和影響因素領(lǐng)域已開(kāi)展了大量研究。鐘麗雯等[3]采用數(shù)據(jù)包絡(luò)分析方法(DEA)研究了廣西及其14個(gè)地級(jí)市2006—2017年水資源利用結(jié)構(gòu)、效率及其時(shí)空分布特征。孟鈺等[4]基于TOPSIS模糊物元-熵權(quán)評(píng)估模型,利用Morans I指數(shù)對(duì)黃河流域各省(區(qū))綜合水資源利用效率進(jìn)行了空間分異研究。孫才志等[5]采用SBM模型和Malmquist指數(shù)模型對(duì)我國(guó)31個(gè)省(區(qū)、市)2000—2014年水資源綠色效率進(jìn)行了測(cè)度,并分析了時(shí)空演變特征及其影響機(jī)理。汪杰等[6]應(yīng)用超效率DEA模型和Malmquist指數(shù)模型對(duì)湖南省環(huán)洞庭湖地區(qū)3市的2008—2017年農(nóng)業(yè)生產(chǎn)效率進(jìn)行研究,并采用Tobit模型分析其影響因素。雖然學(xué)者對(duì)水資源利用效率進(jìn)行了大量研究,但是對(duì)效率進(jìn)行解構(gòu)并對(duì)影響效率原因進(jìn)行深入分析的研究較為缺乏。因此,筆者采用DEA模型對(duì)安徽省水資源利用效率進(jìn)行測(cè)算,并采用Tobit模型對(duì)影響因素進(jìn)行分析,旨在為快速提升安徽省水資源利用效率和區(qū)域高質(zhì)量發(fā)展提供決策參考依據(jù)。

1 研究方法

1.1 超效率SBM-DEA模型

數(shù)據(jù)包絡(luò)分析法(DEA)是Charnes等于1987年提出的一種測(cè)算相對(duì)效率的方法,該方法在固定決策單元輸入或輸出不變的基礎(chǔ)上,根據(jù)決策單元偏離最優(yōu)前沿面的距離來(lái)測(cè)算相對(duì)效率[7]。該方法雖然可以避免人為主觀因素對(duì)計(jì)算結(jié)果的影響,但是易受到松弛變量影響。Tone于2001年提出了SBM-DEA模型,該模型很好地解決了非期望產(chǎn)出和投入與產(chǎn)出指標(biāo)等松弛變量導(dǎo)致的測(cè)算結(jié)果偏高的問(wèn)題[8]。該模型表達(dá)式為:

ρ=min1-1NNn=1Sxnxt′k′n1+1M+1Mm=1Symyt′k′m+Ii=1Sbibt′k′m

s.t.=xt′k′n=Tt=1Kk=1λtkxtkn+Sxn,n=1,2,…,N

yt′k′m=Tt=1Kk=1λtkytkm+Sym,m=1,2,…,M

bt′k′i=Tt=1Kk=1λtkbtki+Sbi,i=1,2,…,I

λtk≥0,Sxn≥0,Sym≥0,Sbi≥0,k=1,2,…,K

式中:ρ為評(píng)價(jià)單元的超效率值;N為評(píng)價(jià)單元的投入量;M、I分別為期望、非期望的產(chǎn)出量;Sxn、Sym、Sbi為松弛變量;xt′k′n、yt′k′m、bt′k′m為在第k′個(gè)生產(chǎn)單元的第t′期間的投入產(chǎn)出量;λtk表示權(quán)重。當(dāng)ρ<1時(shí),表明該評(píng)價(jià)單元DEA無(wú)效;當(dāng)ρ≥1時(shí),表明該評(píng)價(jià)單元DEA有效;ρ值越大則相對(duì)效率越高。

1.2 DEA-Malmquist 指數(shù)

DEA-Malmquist 指數(shù)最早由Malmquist 于1953年提出,后由Fare等于1994年將Malmquist 指數(shù)和DEA模型相結(jié)合,構(gòu)建可以計(jì)算評(píng)價(jià)單元的全要素生產(chǎn)率,以反映出評(píng)價(jià)單元生產(chǎn)效率在時(shí)間上的變化趨勢(shì),屬于一種動(dòng)態(tài)效率評(píng)價(jià)工具[9]。計(jì)算公式為:

CHEff=Dt+1(xt+1,yt+1)Dt(xt,yt)

CHTech=

Dt(xt+1,yt+1)Dt+1(xt+1,yt+1)×

Dt(xt,yt)Dt+1(xt,yt)

TFP=M(xt+1,yt+1,xt,yt)

=Dt(xt+1,yt+1)Dt+1(xt+1,yt+1)×

Dt(xt,yt)Dt+1(xt,yt)

×Dt+1(xt+1,yt+1)Dt(xt,yt)

式中:xt、yt、xt+1、yt+1表示t、t+1時(shí)期內(nèi)的投入與產(chǎn)出;Dt(xt,yt)代表效率水平用t時(shí)期內(nèi)的技術(shù)表示;Dt(xt+1,yt+1)代表t+1時(shí)期的效率水平用t時(shí)期的技術(shù)表示;Dt+1(xt ,yt)代表t時(shí)期的效率水平用t+1時(shí)期的技術(shù)表示;Dt+1(xt+1,yt+1)代表當(dāng)前效率水平用t+1時(shí)期的技術(shù)表示。TFP為t、t+1時(shí)期的生產(chǎn)率指數(shù),CHEff為技術(shù)效率變化指數(shù),CHTech為技術(shù)進(jìn)步變化指數(shù),CHPe為純技術(shù)效率變化,CHSe為規(guī)模效率。

1.3 Tobit模型

Tobit模型主要用于研究因變量含有截?cái)鄶?shù)據(jù)情況下的回歸問(wèn)題,是一種解決因變量觀察數(shù)據(jù)含有截?cái)鄶?shù)據(jù)的回歸模型方法。計(jì)算公式為:

y*=βxi+μi μi∶N(0,σ2)

yi=y* yi>0

yi=0 yi≤0

式中,yi為因變量,xi為因變量,β為系數(shù),誤差項(xiàng)μi獨(dú)立且服從正態(tài)分布。

2 實(shí)證分析

2.1 評(píng)價(jià)指標(biāo)選取及數(shù)據(jù)來(lái)源

該研究以安徽省16地市為研究區(qū)域,建立2011—2020年安徽省16個(gè)地市的面板數(shù)據(jù)。研究數(shù)據(jù)來(lái)源于2011—2020年《安徽省水資源公報(bào)》《安徽統(tǒng)計(jì)年鑒》。基于對(duì)已有文獻(xiàn)指標(biāo)的梳理,并考慮研究數(shù)據(jù)的統(tǒng)計(jì)口徑范圍,將用水總量、固定資產(chǎn)投資總額、從業(yè)人數(shù)作為投入指標(biāo),將市轄區(qū)GDP作為產(chǎn)出指標(biāo)。

2.2 SBM-DEA模型分析

采用MaxDEA8.0軟件對(duì)2011—2020年安徽省16個(gè)地市的水資源利用效率進(jìn)行測(cè)算,將收集的研究指標(biāo)數(shù)據(jù)輸入,結(jié)果見(jiàn)圖1、2。由圖1可知,2011—2020年安徽省16個(gè)地市中,僅有合肥、淮北、黃山3個(gè)城市的水資源利用效率平均值等于1,表明這3個(gè)城市的水資源利用水平較高,達(dá)到DEA有效;而其余13個(gè)城市的水資源利用效率均小于1,均位于生產(chǎn)前沿面之下,表示這13個(gè)城市存在水資源浪費(fèi)、水資源配置效率低下等問(wèn)題。16個(gè)地市中,池州(0.754)、宣城(0.722)的水資源利用效率平均值排名最后2位。由圖2可知,2011—2020年安徽省水資源利用效率平均值從0.931波動(dòng)下降為0.873,利用效率有所降低。具體分析發(fā)現(xiàn),2011—2015年安徽省水資源利用效率平均值均持續(xù)下降,由2011年的0.931降為2015年的0.871,之后又上升到2016年的0.906,此后又持續(xù)降低,由2017年的0.884降為2019年的0.854,而2020年又上升為0.873,呈上升趨勢(shì)。近年來(lái)安徽省經(jīng)濟(jì)和技術(shù)的發(fā)展,整體正向提升了水資源利用和管理能力。

2.3 DEA-Malmquist 指數(shù)分析

由表1可知,Malmquist 指數(shù)大于1的城市有合肥、淮北、蕪湖、銅陵、安慶和黃山,表明這6個(gè)城市的水資源利用效率在研究期內(nèi)均表現(xiàn)為上升,呈現(xiàn)出穩(wěn)步增長(zhǎng)態(tài)勢(shì);其余10個(gè)城市的Malmquist 指數(shù)值均小于1,表明這10個(gè)城市水資源利用效率在研究期呈現(xiàn)出下降態(tài)勢(shì);安徽省Malmquist 指數(shù)平均值為0.981,表明研究期內(nèi)安徽省水資源利用效率總體表現(xiàn)為下降態(tài)勢(shì),其中CHEff、CHTech平均值均小于1,表示安徽省水資源利用效率的提升同時(shí)受到技術(shù)效率變化、技術(shù)進(jìn)步變化的阻礙作用,安徽省水資源利用現(xiàn)有技術(shù)應(yīng)用水平有待提高。

進(jìn)一步分析可知,安徽省16個(gè)地市中,合肥市的Malmquist 指數(shù)值最大,為1.077,表明水資源利用效率處于上升階段;CHEff、CHTech分別為1.000和1.077,表明技術(shù)效率和技術(shù)進(jìn)步共同促進(jìn)了水資源利用效率的加快提升;CHPe和CHSe均為1.000,表明水資源利用管理水平和投入規(guī)模均較高。淮北市、安慶市和黃山市的Malmquist 指數(shù)值均大于1,且CHEff、CHTech、CHPe和CHSe均大于或等于1,表明技術(shù)效率和技術(shù)進(jìn)步均正向促進(jìn)了以上3個(gè)城市水資源利用效率的提升,且其水資源利用管理水平和投入規(guī)模均較高。蕪湖市的Malmquist 指數(shù)值為1.009,大于1,表明其水資源利用效率呈上升態(tài)勢(shì);CHEff、CHTech分別為0.996、1.013,表明技術(shù)進(jìn)步促進(jìn)了該市水資源利用效率的提升;CHPe和CHSe分別為1.001、0.995,可以看出存在投入規(guī)模不足的問(wèn)題。銅陵市的Malmquist 指數(shù)值為1.008,大于1,表明其水資源利用效率呈上升態(tài)勢(shì);CHEff、CHTech分別為0.983、1.025,表明該市水資源利用效率的提升是由于技術(shù)進(jìn)步與創(chuàng)新;CHPe和CHSe分別為0.998、0.985,因而存在水資源利用管理不善、投入規(guī)模不足等問(wèn)題。其余10個(gè)地市的Malmquist 指數(shù)值均小于1,且大部分地市的CHEff、CHSe均小于1,表明技術(shù)效率低下和投入規(guī)模不足是這10個(gè)地市Malmquist 指數(shù)降低的主要原因。

2.4 水資源利用效率影響因素分析

關(guān)于水資源利用效率影響因素的識(shí)別和篩選的研究較多,丁緒輝等[10]將二三產(chǎn)業(yè)規(guī)模、市場(chǎng)化水平、農(nóng)業(yè)種植結(jié)構(gòu)、工業(yè)廢水排放強(qiáng)度、用水強(qiáng)度等作為水資源利用效率的影響因素。武繼堯等[11]將水資源總量、農(nóng)業(yè)用水總量、工業(yè)用水比重、第三產(chǎn)業(yè)增加值、萬(wàn)元GDP用水量何交通基礎(chǔ)設(shè)施作為水資源利用效率的影響因素。應(yīng)卓暉等[12]將人均用水量、第一產(chǎn)業(yè)比重、第二產(chǎn)業(yè)比重、人均GDP、農(nóng)業(yè)用水量作為水資源利用效率的影響因素。何偉等[13]選取人均GDP、第二產(chǎn)業(yè)比重、第三產(chǎn)業(yè)從業(yè)人口比重、人均水資源量、萬(wàn)元GDP水耗、人均生活用水等作為水資源利用效率的影響因素。參考前人研究成果,并考慮相關(guān)數(shù)據(jù)的可比性、可獲得性,從經(jīng)濟(jì)發(fā)展水平、產(chǎn)業(yè)結(jié)構(gòu)、資源稟賦、水污染排放和農(nóng)業(yè)灌溉5個(gè)維度選取變量,其中,經(jīng)濟(jì)發(fā)展水平用人均GDP(元/人)進(jìn)行表征,產(chǎn)業(yè)結(jié)構(gòu)用第一產(chǎn)業(yè)比重(%)進(jìn)行表征,資源稟賦用人均用水量(m3/人)進(jìn)行表征,水污染排放采用人均污水排放量(t/人)表征,農(nóng)業(yè)灌溉用農(nóng)業(yè)用水比重(%)進(jìn)行表征。根據(jù)上述分析結(jié)果,借助Stata 16.0軟件進(jìn)行Tobit回歸分析,結(jié)果見(jiàn)表2。

由表2可知,人均GDP對(duì)安徽省水資源利用效率的提升具有正向促進(jìn)作用,且正向作用凸顯。人均GDP代表的是當(dāng)?shù)氐慕?jīng)濟(jì)發(fā)展水平,經(jīng)濟(jì)發(fā)展水平越高的地區(qū)就會(huì)積極引進(jìn)先進(jìn)的技術(shù)和管理機(jī)制,并逐漸形成區(qū)域聚集,當(dāng)發(fā)展至一定階段時(shí),技術(shù)創(chuàng)新水平的提升會(huì)促使低成本的資源配置,并促進(jìn)水資源的循環(huán)利用。第一產(chǎn)業(yè)比重、農(nóng)業(yè)用水比重對(duì)安徽省水資源利用效率的提升均表現(xiàn)出顯著的正向促進(jìn)作用,安徽省農(nóng)業(yè)較為發(fā)達(dá),水資源灌溉技術(shù)較為成熟,近年來(lái)噴灌、低壓管灌、微灌等高效節(jié)水灌溉技術(shù)也得到了推廣應(yīng)用,這些技術(shù)均對(duì)水資源利用效率具有積極的正向影響。人均用水量對(duì)安徽省水資源利用效率的影響較為顯著,且呈負(fù)向作用,這表明在一些水資源量相對(duì)豐富的地區(qū),居民缺乏節(jié)水意識(shí),浪費(fèi)水資源現(xiàn)象嚴(yán)重,而在水資源量相對(duì)缺乏的地區(qū),居民的節(jié)水意識(shí)較強(qiáng),再加上政府部門加大了水資源管理力度,完善了水權(quán)分配制度和節(jié)水法規(guī),一定程度上促進(jìn)了水資源的有效利用。人均污水排放量對(duì)安徽省水資源利用效率的影響較為顯著,且呈負(fù)向作用,從水資源投入方面看,可通過(guò)降低水污染物向水體排放,改善水環(huán)境質(zhì)量,增加水資源的經(jīng)濟(jì)產(chǎn)出,以提升水資源利用效率;從產(chǎn)出方面看,水污染物屬于一種非期望產(chǎn)出,通過(guò)實(shí)施降低水污染物排放政策,可有效提升水資源利用效率。

3 結(jié)論與討論

(1)2011—2020年安徽省水資源利用效率平均值從0.931波動(dòng)下降為0.873,利用效率有所降低。16個(gè)地市中,僅有合肥、淮北、黃山的水資源利用效率平均值等于1,達(dá)到DEA有效,其余13個(gè)城市的水資源利用效率均小于1,表示以上城市存在水資源浪費(fèi)、水資源配置效率低下等問(wèn)題。

(2)2011—2020年安徽省水資源利用效率Malmquist 指數(shù)平均值為0.981,總體處于下降階段,安徽省水資源利用效率的提升同時(shí)受到技術(shù)效率變化、技術(shù)進(jìn)步變化的阻礙作用,安徽省水資源利用現(xiàn)有技術(shù)應(yīng)用水平有待提高。16個(gè)地市中,合肥、淮北、蕪湖、銅陵、安慶和黃山的水資源利用效率Malmquist 指數(shù)平均值處于上升階段,其余10個(gè)城市的Malmquist 指數(shù)值均處于下降階段。

(3)人均GDP、第一產(chǎn)業(yè)比重、農(nóng)業(yè)用水比重因素均對(duì)安徽省水資源利用效率的提升表現(xiàn)出顯著的正向促進(jìn)作用;人均用水量、人均污水排放量均對(duì)安徽省水資源利用效率均表現(xiàn)出顯著負(fù)向影響。

(4)安徽省水資源利用效率具有很大的提升空間,應(yīng)從加大水資源利用的技術(shù)效率投入和完善管理政策兩方面入手,結(jié)合各地市具體情況不斷優(yōu)化水資源的投入產(chǎn)出,探索科學(xué)合理的資源配置方式,從而促進(jìn)安徽省水資源利用效率的提高。

參考文獻(xiàn)

[1] 張祥永,曹陽(yáng),吳文飛.我國(guó)省域水資源可持續(xù)性與生態(tài)系統(tǒng)服務(wù)耦合協(xié)調(diào)關(guān)系[J].人民黃河,2022,44(12):73-78,100.

[2] 任志安,卞曉桐,劉柏陽(yáng).合肥都市圈城鎮(zhèn)化質(zhì)量對(duì)巢湖流域水資源利用效率的影響研究[J].青島科技大學(xué)學(xué)報(bào)(社會(huì)科學(xué)版),2021,37(4):1-10.

[3] 鐘麗雯,張建兵,蔡蕓霜,等.廣西水資源利用效率及其時(shí)空格局[J].經(jīng)濟(jì)地理,2020,40(6):193-202.

[4] 孟鈺,杜瓊英,管新建,等.黃河流域水資源利用效率評(píng)估及空間分異研究[J].中國(guó)農(nóng)村水利水電,2020(10):12-16.

[5] 孫才志,馬奇飛,趙良仕.基于SBM-Malmquist生產(chǎn)率指數(shù)模型的中國(guó)水資源綠色效率變動(dòng)研究[J].資源科學(xué),2018,40(5):993-1005.

[6] 汪杰,李姣,周翠煙,等.基于DEA-Tobit模型的湖南省環(huán)洞庭湖地區(qū)農(nóng)業(yè)生產(chǎn)效率研究[J].天津農(nóng)業(yè)科學(xué),2019,25(12):48-55.

[7] 李世玉,劉笑冰.基于DEA模型的京津冀農(nóng)業(yè)水資源利用效率研究[J].天津農(nóng)業(yè)科學(xué),2021,27(9):69-74.

[8] 吉雪強(qiáng),尚杰.基于三階段SBM模型的中國(guó)農(nóng)業(yè)生態(tài)效率研究[J].中國(guó)農(nóng)業(yè)資源與區(qū)劃,2021,42(7):210-217.

[9] 張麗,鄭軍.基于DEA-Malmquist模型的山東省農(nóng)業(yè)生產(chǎn)效率測(cè)算研究[J].湖北農(nóng)業(yè)科學(xué),2019,58(21):219-223.

[10] 丁緒輝,賀菊花,王柳元.考慮非合意產(chǎn)出的省際水資源利用效率及驅(qū)動(dòng)因素研究:基于SE-SBM與Tobit模型的考察[J].中國(guó)人口·資源與環(huán)境,2018,28(1):157-164.

[11] 武繼堯,韓增林,吳慧歌.遼寧省水資源利用效率與空間溢出效應(yīng)研究[J].資源開(kāi)發(fā)與市場(chǎng),2017,33(4):417-422.

[12] 應(yīng)卓暉,趙衡,王富強(qiáng),等.基于DEA和Tobit模型的河南省水資源利用效率評(píng)價(jià)及影響因素[J].南水北調(diào)與水利科技,2021,19(2):255-262.

[13] 何偉,王語(yǔ)苓.黃河流域城市水資源利用效率測(cè)算及影響因素分析[J].環(huán)境科學(xué)學(xué)報(bào),2021,41(11):4760-4770.

主站蜘蛛池模板: 国产小视频免费| 狠狠干综合| 国产成人精品男人的天堂下载 | 精品国产中文一级毛片在线看| 日韩精品久久无码中文字幕色欲| 亚洲成a人片| 欧美成人一区午夜福利在线| 无码人妻免费| 国产毛片高清一级国语 | 精品国产成人a在线观看| 亚洲欧洲日本在线| 国产在线精品99一区不卡| 99热国产在线精品99| 亚洲人免费视频| 欧美另类视频一区二区三区| 天天色天天综合| 亚洲AⅤ永久无码精品毛片| 免费a级毛片视频| 内射人妻无码色AV天堂| 亚洲精品综合一二三区在线| 精品一区二区三区视频免费观看| 国产精品伦视频观看免费| 五月天在线网站| 亚洲系列无码专区偷窥无码| 亚洲国产欧美国产综合久久 | 国产黄色爱视频| 国产91精品久久| 精品国产成人av免费| 亚洲无码不卡网| 国产三级韩国三级理| aⅴ免费在线观看| 2021国产在线视频| 高清色本在线www| 国产成人精品亚洲77美色| 亚洲区一区| 在线亚洲天堂| 久久青草精品一区二区三区 | 成人在线观看一区| 国产18页| 午夜人性色福利无码视频在线观看| 国产福利一区在线| 波多野结衣久久高清免费| 中文一级毛片| 亚洲欧美不卡| 五月激情综合网| 99视频在线观看免费| 国产swag在线观看| 久久久久亚洲精品成人网| 另类欧美日韩| 国产嫖妓91东北老熟女久久一| 久久鸭综合久久国产| 蜜桃视频一区二区三区| 97se综合| 久久天天躁狠狠躁夜夜2020一| 日韩视频福利| 伊人91视频| 久久五月视频| 亚洲精品无码高潮喷水A| 亚洲人妖在线| 国产欧美中文字幕| 国产凹凸视频在线观看| 青青草国产一区二区三区| 国产午夜不卡| 国产一区二区三区日韩精品| 在线看片中文字幕| 麻豆AV网站免费进入| 国产男人的天堂| 欧美成人一级| 高清无码一本到东京热| 亚洲第一区欧美国产综合| 69视频国产| 爱色欧美亚洲综合图区| 麻豆精品在线视频| 欧美成人怡春院在线激情| 国产精品视频导航| 国产在线91在线电影| 中文字幕 91| 亚洲欧美精品一中文字幕| 全部无卡免费的毛片在线看| 国产网站黄| 国产精品第5页| 香蕉在线视频网站|