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

    小麥模型算法集成平臺構(gòu)建與算法比較

    2022-08-05 00:37:36陳先冠馮利平白慧卿王春雷余衛(wèi)東
    農(nóng)業(yè)機械學報 2022年6期
    關(guān)鍵詞:春化發(fā)育期實測值

    陳先冠 馮利平 白慧卿 王春雷 王 靖 余衛(wèi)東

    (1.中國農(nóng)業(yè)大學資源與環(huán)境學院,北京 100193;2.中國氣象局河南省農(nóng)業(yè)氣象保障與應(yīng)用技術(shù)重點實驗室,鄭州 450003)

    0 引言

    作物模型是利用系統(tǒng)分析和計算機技術(shù),綜合作物生理生態(tài)、農(nóng)業(yè)氣象、土壤和農(nóng)學等學科研究成果,將作物與其生態(tài)環(huán)境因子作為整體進行動態(tài)的定量化分析,能夠再現(xiàn)作物生長發(fā)育和產(chǎn)量形成的過程。作物模型技術(shù)自20世紀60年代開創(chuàng)以來,經(jīng)過多年發(fā)展,經(jīng)歷了從經(jīng)驗走向機理,從理論走向?qū)嵺`的過程[1-3]。目前已有一系列的模型產(chǎn)品,如荷蘭的WOFOST[4]、美國的DSSAT[5]和澳大利亞的APSIM[6]等。近年來,隨著國際AgMIP項目的開展,作物生長模型研究開始側(cè)重于現(xiàn)有模型的比較、改進和完善[7]。

    基于過程的作物模型可以綜合反映基因型、環(huán)境和管理措施對作物生長發(fā)育的影響,已經(jīng)廣泛用于氣候變化影響評估[8]、作物管理[9]、品種設(shè)計等多個方面[10-11],但模型模擬存在較大的不確定性[12-15]。由于建模目標以及對模型機理認知的不同,不同建模者構(gòu)建的模型不同,從而導(dǎo)致不同模型模擬作物生長發(fā)育和產(chǎn)量形成等過程差異較大[12-15]。氣候變化對作物產(chǎn)量的影響評估中,作物模型比氣候模式對產(chǎn)量變化的不確定性影響更大,模型結(jié)構(gòu)是氣候變化對產(chǎn)量影響不確定性的最主要來源[14-15]。因此,需要系統(tǒng)地比較不同模型的模擬性能并改進相應(yīng)算法,以降低模型結(jié)構(gòu)的不確定性。FLEISHER等[16]對9個馬鈴薯模型進行比較,評價了不同模型對氣候變化的響應(yīng),量化了模型的不確定性。RAN等[17]發(fā)現(xiàn)在充分灌溉和氣候年際變異高的情況下,光能驅(qū)動模型DSSAT比水分驅(qū)動模型AquaCrop更適合生物量和產(chǎn)量的模擬。PALOSUO等[18]比較了8個小麥模型在不同氣候條件下的模擬性能,表明沒有一個模型能夠完美地再現(xiàn)作物生長和產(chǎn)量形成的過程。前人已對作物模型做了大量比較,但由于模型的模擬結(jié)果是由多個過程綜合影響得到的,而只分析模擬結(jié)果(如生物量、產(chǎn)量和葉面積等)很難直接判斷模型的哪一部分需要改進,可見目前大部分的模型比較研究只闡明了模型的不確定性,但對模型改進的指導(dǎo)較為有限[19]。因此,為比較特定過程中不同算法的差異,需要采取控制試驗的方式控制模型的各個過程,在單一平臺進行算法比較。

    通過在單一平臺集成特定過程的不同算法,可以促進算法的比較、評估與改進,有利于降低模型結(jié)構(gòu)的不確定性。DONATELLI等[20]最先通過重用和組合已有算法構(gòu)建評估土壤溫度的通用框架。JIN等[19]進一步采用了“算法集成”的思想,比較了16種主要玉米模型的溫度脅迫和水分脅迫算法,為降低模型模擬的不確定性和提升產(chǎn)量的模擬精度提供了參考。本研究在模型模塊化的基礎(chǔ)上,在主要模塊中集成多種算法,旨在比較與改進現(xiàn)有的模擬算法、提升模型模擬性能以及提高模型開發(fā)效率及質(zhì)量。本文首先通過參考國內(nèi)外主流小麥模型CERES-Wheat[21]、APSIM-Wheat[6]、WOFOST[4]、SWAT[22]、WheatSM[23]的主要模塊算法,利用Python語言搭建標準平臺,將各模塊的不同算法整合到平臺中;再利用田間試驗數(shù)據(jù)對模型平臺進行檢驗;最后進行特定模塊的算法比較。

    1 平臺描述

    1.1 平臺框架

    模型平臺主要包括氣象、土壤水分、土壤養(yǎng)分、土壤溫度、發(fā)育期、生物量、器官建成、產(chǎn)量形成、田間管理等模塊。由于Python語言采用模塊化設(shè)計,方便添加及更新模塊,模塊之間傳遞變量簡單,本模型平臺使用Python語言將各個模塊有機地結(jié)合在一起,各模塊不僅可單獨運行,也可作為一個整體運行,小麥模型算法集成平臺(Wheat model algorithm integration platform,WMAIP)框架如圖1所示。

    圖1 小麥模型算法集成平臺框架Fig.1 Schematic of wheat model algorithm integration platform (WMAIP)

    1.2 發(fā)育期模塊

    發(fā)育期模塊包含小麥鐘模型法(Wheat clock,WC)與熱時法(Thermal time,TT)。

    1.2.1小麥鐘模型法

    小麥鐘模型[23]表達式為

    (1)

    其中

    (2)

    (3)

    式中M——發(fā)育期或發(fā)育階段的發(fā)育進程,完成時M為1

    t——時間

    k——基本發(fā)育系數(shù),由品種自身的遺傳特性決定,k越大,說明該品種發(fā)育速度快,是早熟品種,反之則發(fā)育遲緩,表現(xiàn)晚熟

    TE——溫度效應(yīng)因子,反映溫度對小麥發(fā)育的非線性影響程度

    p——溫度系數(shù),反映品種在某一發(fā)育階段內(nèi)對溫度的敏感性,p越大,說明在下限到最適溫度范圍內(nèi)時,增溫對發(fā)育速度的促進作用越大,反之則越小

    PE——光周期效應(yīng)因子,反映光照對小麥發(fā)育的非線性影響程度

    q——光周期反應(yīng)特征遺傳系數(shù),反映某品種在某一發(fā)育階段內(nèi)對光周期反應(yīng)的敏感性,q越大,說明該品種對日長的反應(yīng)越敏感,反之則越鈍感

    fEC——肥料、播種深度等可控栽培措施因子

    Tmean——日平均溫度,℃

    Tb——某一階段內(nèi)小麥生長發(fā)育的下限溫度,播種—出苗階段、出苗—拔節(jié)階段、拔節(jié)—開花階段、開花—成熟階段分別為1、3、3、9℃

    To——發(fā)育階段內(nèi)小麥生長的最適溫度,播種—開花階段為20℃,開花—成熟階段為22℃

    PL——每天的光照時長,包括曙暮光在內(nèi),h

    PLo——發(fā)育階段內(nèi)的最適光照時長,取18 h

    PLb——發(fā)育階段內(nèi)的臨界光照時長,取8 h

    本文對小麥鐘模型中春化作用的起始階段進行改進,由原來的出苗期改為播種期計算春化效應(yīng)。春化效應(yīng)因子VE為

    (4)

    式中VT——春化溫度,℃

    VTb——春化下限溫度,取-4℃

    VTo1——春化最適溫度下限,取3℃

    VTo2——春化最適溫度上限,取7℃

    VTm——春化上限溫度,取18℃

    1.2.2熱時法

    熱時法參考APSIM-Wheat發(fā)育期模塊[24],計算公式為

    TT,adj=∑(TT,Δmin(fD,fV))

    (5)

    其中

    (6)

    fD=1-0.002Rp(20-PL)2

    (7)

    fV=1-(0.005 454 5RV+0.000 3)(50-V)

    (8)

    式中TT,adj——光周期和春化作用訂正后的熱時,℃·d

    TT,Δ——每天的熱時,℃·d

    fD——光周期影響因子

    V——總春化作用

    fV——春化作用影響因子

    Tc——日平均冠層溫度,℃

    Rp——品種參數(shù),光周期敏感系數(shù)

    RV——品種參數(shù),春化作用敏感系數(shù)

    春化作用VΔ由日平均冠層溫度、日最高溫度Tmax和日最低溫度Tmin計算得來。

    (9)

    當日最高氣溫大于30℃并且總春化作用V小于10時,會發(fā)生反春化作用VΔ,d,計算公式為

    VΔ,d=min(0.5(Tmax-30),V)
    (Tmax>30℃且V<10)

    (10)

    其中

    V=∑(VΔ-VΔ,d)

    (11)

    1.3 生物量模塊

    生物量模塊包括二氧化碳同化率法(Carbon assimilation,CA)、光能利用效率法(Radiation use efficiency,RUE)以及群體光合作用法(Canopy photosynthesis,CP)3種算法。

    1.3.1二氧化碳同化率法

    二氧化碳同化率法采用WOFOST模型的干物質(zhì)生產(chǎn)算法[25]

    WΔ=Ce(A-Rm)

    (12)

    式中WΔ——日生物量生長速率,kg/(hm2·d)

    Ce——同化物轉(zhuǎn)化效率系數(shù),取0.7 kg/kg

    A——總同化速率,kg/(hm2·d)

    Rm——維持呼吸速率,kg/(hm2·d)

    冠層光合速率采用Gaussian 3點積分法,對葉片所在的時間和空間的瞬時同化速率進行積分,求得日總同化速率。相對冠層高度L處瞬時總同化速率是光飽和時葉片的CO2同化速率和所吸收的輻射量的函數(shù),以負指數(shù)形式表示為

    (13)

    式中AL——單位葉面積在相對高度L處的瞬時總同化速率,kg/(hm2·h)

    Am——光飽和時葉片的CO2同化速率,因白天溫度不同而異,kg/(hm2·h)

    EFF——單葉片光能利用效率,

    (kg/(hm2·h))/(J/(m2·s))

    Ia——吸收的總輻射量,J/(m2·s)

    維持呼吸量與日平均溫度有關(guān),維持呼吸速率計算公式為

    (14)

    Rm=Rm,TW

    (15)

    式中Rm,T——溫度T下維持呼吸系數(shù),kg/kg

    Q10——呼吸墑,取2

    Tr——參考溫度,取25℃

    Rm,Tr——參考溫度維持呼吸,抽穗前取0.02 kg/kg,抽穗后取0.01 kg/kg

    W——總干物質(zhì)量,kg/hm2

    1.3.2光能利用效率法

    光能利用效率法參考APSIM-Wheat模型[24]。在APSIM-Wheat中,基于光能利用效率(RUE)的概念,光能驅(qū)動干物質(zhì)生產(chǎn)的計算公式為

    WΔ=10IiRUEmin(fT,poto,fN,poto)fW,potofC

    (16)

    式中Ii——冠層截獲輻射量,MJ/(m2·d)

    RUE——光能利用效率,g/MJ

    fT,poto——溫度脅迫因子

    fN,poto——氮素脅迫因子

    fW,poto——水分脅迫因子

    fC——二氧化碳因子

    1.3.3群體光合作用法

    群體光合作用法參考WheatSM中的門司公式[23]

    (17)

    式中PGd——群體光合作用日總量,g/(m2·d)

    DL——日長,h

    Pmax——光飽和點下的最大凈光合作用速率,g/(m2·h)

    Pa——光響應(yīng)曲線的初始斜率,表示葉片在弱光條件下的光合能力,g/MJ

    I0——太陽總輻射量,MJ/(m2·d)

    Kec——群體消光系數(shù)

    σ——麥田反射率,取8%

    LAI——葉面積指數(shù)

    1.4 器官建成模塊

    1.4.1物質(zhì)分配

    小麥在生長過程中各個器官(根、莖、葉和穗)的增重與全株的增重的比率稱作器官的分配系數(shù),反映的是各器官的生長分配動態(tài)[23]

    Wi,d=WΔfi,d

    (18)

    式中Wi,d——器官i在時間d的分配量,kg/(hm2·d)

    fi,d——器官i在時間d的分配系數(shù)

    1.4.2葉面積

    葉面積指數(shù)采用比葉重法模擬[25],比葉重采用APSIM-Wheat模型方法計算[24]

    LAI=0.1WL,LSLA

    (19)

    其中

    (20)

    式中WL,L——活葉干物質(zhì)量,kg/hm2

    SLA——比葉重,m2/g

    SLA,max——最大比葉重,取0.027 m2/g

    SLA,min——最小比葉重,取0.022 m2/g

    1.4.3根系生長

    根系生長采用ASSENG等[26]提出的方法,計算公式為

    RLVΔ,j=PMLFRml,jRΔ,j

    (21)

    式中RLVΔ,j——土層j根長密度的日增加量,cm/(cm3·d)

    FRml,j——土層j中的土壤阻力因子

    PML——根生物量與根長的轉(zhuǎn)換系數(shù),取21 000 cm/g

    RΔ,j——土層j中根重密度增加量,g/(cm3·d)

    1.5 產(chǎn)量形成模塊

    產(chǎn)量形成模塊包括籽粒灌漿法(Grain filling,GF)、收獲指數(shù)法(Harvest index,HI)及生物量轉(zhuǎn)移法(Biomass remobilization,BR)。

    1.5.1籽粒灌漿法

    籽粒灌漿法參考APSIM-Wheat[24]。籽粒數(shù)由開花期的莖干物質(zhì)決定,有

    Ng=0.1RgWs

    (22)

    式中Ng——籽粒數(shù)

    Ws——開花期的莖干物質(zhì)量,kg/hm2

    Rg——每克莖干物質(zhì)的籽粒數(shù)

    籽粒需碳量計算公式為

    Dg=NgRfhg,TmeanfN,grain

    (23)

    其中

    (24)

    式中Dg——籽粒需碳量,g/(m2·d)

    Rf——潛在灌漿速率,開花-灌漿始期取0.001 g/d,灌漿始期-成熟取0.002 g/d

    hg,Tmean——日均溫Tmean對籽粒灌漿影響因子

    fN,grain——氮素脅迫對籽粒灌漿的影響因子

    最后,籽粒需碳量被最大籽粒質(zhì)量限制,有

    Dg=min(Dg,Dgm)

    (25)

    其中Dgm=NgSgm-Qmeal(Dgm≥0)

    (26)

    式中Dgm——籽粒最大需碳量,g/m2

    Qmeal——籽粒質(zhì)量,g/m2

    Sgm——最大籽粒質(zhì)量,品種參數(shù),g

    1.5.2收獲指數(shù)法

    收獲指數(shù)法參考SWAT模型[27],實際收獲指數(shù)計算公式為

    (27)

    式中HIact——實際收獲指數(shù)

    HImax——最大收獲指數(shù)

    HImin——最小收獲指數(shù)

    fW,ave——生長季平均水分脅迫系數(shù)

    S1、S2——水分脅迫影響參數(shù),分別取14.3、0.286

    1.5.3生物量轉(zhuǎn)移法

    生物量轉(zhuǎn)移法參考WheatSM模型[23]。在生長階段,小麥的生物量不斷積累并分配給小麥籽粒,最終形成產(chǎn)量,其基本模型為

    (28)

    式中Y——小麥產(chǎn)量,kg/hm2

    de——從播種至小麥抽穗期的時間,d

    dm——全生育期時間,d

    Ter1——抽穗前莖鞘存儲物向籽粒轉(zhuǎn)運的效率,kg/kg

    Ter2——抽穗后光合產(chǎn)物向籽粒轉(zhuǎn)運的效率,kg/kg

    WΔ,d——時間d時生物量生長速率,kg/(hm2·d)

    1.6 潛在蒸散量

    潛在蒸散量的計算公式為

    ETc=EToKc

    (29)

    式中ETc——潛在蒸散量,mm/d

    Kc——冬小麥不同發(fā)育階段的作物系數(shù)

    ETo——作物參考蒸散量,采用Penman-Monteith公式[28]計算,mm/d

    1.7 土壤水分模塊

    土壤水分模塊主要參考ALFASM模型[21,29],主要包括冠層截留、地表徑流、水分收入、水分入滲、根系吸水、水分再分配等過程。其中水分脅迫因子計算式為

    (30)

    其中

    (31)

    式中WRWU——根系吸水量,mm/d

    EP0——潛在蒸騰量,mm/d

    r1、r2、r3——方程參數(shù),分別取0.002 67、62、6.68

    θt,j——土層j的實際含水量,cm3/cm3

    θwp,j——土層j的凋萎含水量,cm3/cm3

    RLV,j——土層j的根長密度,cm/cm3

    ZΔ,j——土層j的厚度,cm

    N——根系達到的最深土層數(shù)

    1.8 土壤氮素模塊

    土壤氮素模塊主要參考CERES-Wheat模型[21,30],主要包括有機氮礦化、氮素固定、硝化、反硝化、尿素水解、氮吸收等過程。

    1.9 土壤溫度模塊

    土壤溫度采用土壤溫波方程[31]

    (32)

    其中

    ω=2π/τ

    (33)

    式中TZ,t——Z深度t時刻的土壤溫度,℃

    A0——地面溫波振幅,即最高溫度為Tmean+A0,最低溫度為Tmean-A0

    ω——溫度曲線的正弦角度

    τ——周期Z——觀測深度,m

    D——土溫衰減深度,m

    2 田間試驗與模型驗證

    2.1 田間試驗

    2.1.1研究區(qū)概況

    冬小麥分期播種試驗在河北省滄州市吳橋縣姚莊村中國農(nóng)業(yè)大學吳橋試驗站(37°37′N,116°23′E)進行。試驗地點位于海河平原黑龍港地區(qū)中部,海拔20 m,地下水位6~9 m,為暖溫帶季風氣候,多年平均年降水量542 mm,小麥季平均年降水量117 mm。試驗地土壤類型為黏壤土,0~200 cm土體的土壤參數(shù)見表1[32]。

    表1 試驗地土壤剖面物理參數(shù)Tab.1 Physical parameters of soil profiles at experimental field

    2.1.2數(shù)據(jù)來源

    逐日氣象資料來自吳橋縣氣象局,主要包括:日最高溫度(℃)、日最低溫度(℃)、日照時數(shù)(h)、日降水量 (mm)、平均風速(m/s)、空氣相對濕度(%)。

    冬小麥分期播種試驗于2017年10月開始至2019年6月結(jié)束,供試品種濟麥22(Jimai_22)為當?shù)刂髟云贩N。2017年10月—2018年6月小麥生長季降水量為347 mm,總輻射量為3 385.2 MJ/m2,日平均溫度為10.6℃,是典型濕潤年;2018年10月—2019年6月小麥生長季降水量為79 mm,總輻射量為3 576.3 MJ/m2,日平均溫度為10.7℃,是典型干旱年。其中2017—2018年分期播種試驗設(shè)置5個播期,每播期3個重復(fù),小區(qū)面積為3.6 m×9 m=32.4 m2;5個播期分別為10月6日(D1)、10月16日(D2)、10月23日(D3)、10月30日(D4)和11月6日(D5),播量分別為187.5、300、375、375、375 kg/hm2,行距15 cm,播前澆底墑水,春季澆拔節(jié)水,灌水定額75 mm。2018—2019年試驗增設(shè)1個播期處理(11月13日,D6),播量375 kg/hm2,其他處理同2017—2018年試驗。

    栽培方式參照小麥節(jié)水高產(chǎn)栽培技術(shù)[33],一次性施入底肥尿素225 kg/hm2、磷酸二銨300 kg/hm2、硫酸鉀225 kg/hm2,春季不再追肥。在小麥生長期間,于越冬期、返青期、起身期、拔節(jié)期、孕穗期、開花期、乳熟期、成熟期共8個時期在每個小區(qū)選取具有代表性的2段0.5 m小麥植株,在105℃下殺青30 min,然后在75℃下干燥至質(zhì)量不變,計算各時期小麥地上部生物量。收獲時每小區(qū)取1.8 m2冬小麥進行測產(chǎn),自然風干后脫粒測產(chǎn)求得產(chǎn)量,文中的產(chǎn)量均以含水率13%計。在上述關(guān)鍵發(fā)育期取土,從地表至地下200 cm處,每20 cm作為一層,采用土鉆法取土并用干燥法測定土壤含水量。

    利用WebPlotDigitizer軟件(http:∥arohatgi.info/WebPlotDigitizer/)提取文獻中吳橋試驗站2011—2014年期間利用濟麥22開展的播期耦合水分試驗數(shù)據(jù),主要包括關(guān)鍵發(fā)育期、地上部生物量和小麥產(chǎn)量。試驗采用裂區(qū)區(qū)組設(shè)計,以播期為主區(qū),水分為副區(qū)。4個播期處理為10月10日、10月20日、10月30日和11月10日。3個水分處理分別為春季不澆水、春季澆拔節(jié)水、春季澆拔節(jié)水和開花水,每次灌水定額75 mm[34]。

    2.2 模型構(gòu)成及參數(shù)率定

    2.2.1模擬模型組成

    基于模型平臺構(gòu)建了6種模擬模型,包括發(fā)育期模塊2種算法、生物量模塊3種算法及產(chǎn)量形成模塊3種算法,模擬模型的結(jié)構(gòu)組成如表2所示。其中,M1和M4用于發(fā)育期算法比較;M1、M2和M3控制了發(fā)育期模塊和產(chǎn)量形成模塊,而生物量模塊采用了不同的算法,用于生物量算法比較;M4、M5和M6控制了發(fā)育期模塊和生物量模塊,而產(chǎn)量形成模塊采用了不同的算法,用于產(chǎn)量形成算法比較。

    表2 不同模擬模型的結(jié)構(gòu)組成Tab.2 Structural composition of different simulation models

    2.2.2作物參數(shù)率定

    利用Python包sopt(https:∥github.com/Lyrichu/sopt)中的簡單遺傳算法(Simple genetic algorithm,SGA)進行參數(shù)自動率定,其中,種群數(shù)目參數(shù)population_size取30,進化代數(shù)參數(shù)generations取10,其余參數(shù)采用默認值?;?011—2012年、2012—2013年和2017—2018年3個年度的試驗數(shù)據(jù)進行調(diào)參,運用2013—2014年和2018—2019年兩個年度的試驗數(shù)據(jù)進行模型驗證。首先根據(jù)試驗觀測的發(fā)育期對發(fā)育期模塊進行參數(shù)率定;然后根據(jù)試驗觀測的干物質(zhì)量對生物量模塊進行參數(shù)率定;最后根據(jù)試驗觀測的產(chǎn)量對產(chǎn)量形成模塊進行參數(shù)率定(表3)。

    本文利用模型模擬值和實際觀測值的均方根誤差(Root mean squared error,RMSE)作為目標函數(shù)[35]。RMSE越小,模擬值和觀測值之間偏差越小,模擬值與觀測值一致性越好,模型模擬結(jié)果越準確可靠。

    2.3 模型集合及模擬效果評價

    采用6個模型模擬結(jié)果的平均值代表模型集合模擬結(jié)果。采用回歸系數(shù)(Regression coefficient,b0)、決定系數(shù) (Coefficient of determination,R2)、均方根誤差(RMSE)和歸一化均方根誤差(Normalized root mean squared error,NRMSE)來評價模型的模擬值與實測值的吻合程度。

    3 結(jié)果與分析

    3.1 模型平臺參數(shù)校準與驗證

    利用WMAIP平臺組成的6個模擬模型都能較好地模擬土壤貯水量和冬小麥的生物學指標。2 m深土層土壤貯水量的實測值和模擬值的差別較小,6個模擬模型的調(diào)參結(jié)果和驗證結(jié)果中均方根誤差(RMSE)分別在45.6~49.0 mm和55.1~59.4 mm之間,歸一化均方根誤差(NRMSE)分別在7.68%~8.26%和10.78%~11.63%之間(表4)。6個模型都能夠較好地模擬葉面積指數(shù)的動態(tài)變化,調(diào)參結(jié)果和驗證結(jié)果中RMSE分別在0.58~0.72 mm2/mm2和0.79~1.34 mm2/mm2之間,NRMSE分別在19.43%~22.28%和26.50%~37.75%之間(表4,圖2)。6個模型都能夠較好地模擬地上部生物量的動態(tài)變化,其中M1~M4的調(diào)參結(jié)果NRMSE分別為11.22%、12.34%、12.28%、11.28%,驗證結(jié)果中NRMSE分別為18.18%、18.72%、16.13%、17.48%(表4,圖2)。產(chǎn)量的實測值和模擬值的差別較小,6個模擬模型的調(diào)參結(jié)果和驗證結(jié)果中RMSE分別在946.1~1 219.8 kg/hm2和923.5~1 402.3 kg/hm2之間,NRMSE分別在12.27%~15.82%和12.48%~18.95%之間(表4)。模型集合有助于提高地上部生物量和產(chǎn)量的模擬精度,模型集合的調(diào)參結(jié)果和驗證結(jié)果中NRMSE都小于單個模擬模型。

    3.2 發(fā)育期算法比較

    冬小麥發(fā)育期模塊的調(diào)參結(jié)果顯示,模型模擬值與觀測值之間都具有較好的一致性。其中,小麥鐘模型播種-拔節(jié)期、播種-開花期和播種-成熟期的模擬值和實測值的RMSE分別為3.5、1.4、1.9 d,NRMSE分別為2.06%、0.69%、0.81%(圖3a);熱時法播種-拔節(jié)期、播種-開花期和播種-成熟期的模擬值和實測值的RMSE分別為2.6、2.6、1.7 d,NRMSE分別為1.50%、1.32%、0.75%(圖3c)。

    表3 冬小麥品種濟麥22的遺傳參數(shù)Tab.3 Genetic parameters for winter wheat variety Jimai_22

    表4 不同模型的調(diào)參與驗證結(jié)果Tab.4 Calibration and validation results for different simulation model

    圖2 葉面積指數(shù)和地上部生物量的實測值和模擬值比較Fig.2 Comparisons of observed and simulated leaf area index and above-ground biomass

    圖3 兩種發(fā)育期算法的冬小麥發(fā)育期調(diào)參和驗證結(jié)果Fig.3 Comparison between simulated and observed growth period of winter wheat with “wheat clock” model (CM) and thermal time (TT) approach for calibration and validation

    冬小麥發(fā)育期模塊的驗證結(jié)果表明,發(fā)育期模塊的2種算法都能夠很好地模擬冬小麥的發(fā)育期,熱時法模擬播種至拔節(jié)階段較優(yōu),小麥鐘模型模擬播種至開花階段和播種至成熟階段較優(yōu)。小麥鐘模型的播種-拔節(jié)期、播種-開花期和播種-成熟期模擬值和實測值的RMSE分別為6.5、1.9、1.3 d,NRMSE分別為4.00%、1.02%、0.56%(圖3b);熱時法的播種-拔節(jié)期、播種-開花期和播種-成熟期模擬值和實測值的RMSE分別為4.0、4.3、2.0 d,NRMSE分別為2.46%、2.24%、0.91%(圖3d)。

    3.3 生物量積累算法比較

    生物量模塊調(diào)參結(jié)果如圖4所示,群體光合作用法生物量模擬值和實測值的回歸系數(shù)b0和決定系數(shù)R2分別為0.96和0.94;光能利用效率法生物量模擬值和實測值模型的b0和R2分別為0.95和0.94;二氧化碳同化率法生物量模擬值和實測值的b0和R2分別為0.93和0.94。生物量模塊3種算法都能較好地模擬生物量動態(tài),模擬值與實測值有很好的線性關(guān)系,可以反映90%以上的生物量變異。

    生物量模塊的驗證結(jié)果中,群體光合作用法生物量模擬值和實測值的b0和R2分別為0.93和0.92;光能利用效率法生物量模擬值和實測值的b0和R2分別為0.80和0.92;二氧化碳同化率法生物量模擬值和實測值的b0和R2分別為0.82和0.94。3種生物量算法均為模擬小麥生物量的較佳模型。

    圖4 3種生物量算法的冬小麥生物量調(diào)參與驗證結(jié)果Fig.4 Comparison between observed and simulated above-ground biomass of winter wheat with different algorithm for calibration and validation

    3.4 產(chǎn)量形成算法比較

    圖5 3種產(chǎn)量算法的冬小麥產(chǎn)量調(diào)參和驗證結(jié)果Fig.5 Comparison between observed and simulated yield of winter wheat with different algorithms for calibration and validation

    產(chǎn)量模塊調(diào)參結(jié)果如圖5所示,生物量轉(zhuǎn)移法的產(chǎn)量模擬值和實測值的b0和R2分別為0.51和0.35;籽粒灌漿法的產(chǎn)量模擬值和實測值的b0和R2分別為0.70和0.31;收獲指數(shù)法的產(chǎn)量模擬值和實測值的b0和R2分別為0.63和0.44。3種產(chǎn)量算法都能夠捕捉冬小麥不同播期和水分下產(chǎn)量的變化趨勢,但生物量轉(zhuǎn)移法的模擬產(chǎn)量略低于實測產(chǎn)量,而籽粒灌漿法和收獲指數(shù)法的模擬產(chǎn)量略高于實測產(chǎn)量。

    產(chǎn)量形成模塊驗證結(jié)果中,3種算法都能夠很好地模擬冬小麥的產(chǎn)量,模擬和實測的NRMSE都在19%以下(表4)。其中,生物量轉(zhuǎn)移法產(chǎn)量模擬值和實測值的b0和R2分別為0.29和0.53;籽粒灌漿法的產(chǎn)量模擬值和實測值的b0和R2分別為0.51和0.44;收獲指數(shù)法的產(chǎn)量模擬值和實測值的b0和R2分別為0.31和0.38(圖5)。3種產(chǎn)量算法模擬的產(chǎn)量變化趨勢較為一致,但生物量轉(zhuǎn)移法效果略好。

    4 討論

    本研究集成多模塊的多種算法,組成6個模擬模型,采用5個小麥生長季播期水分的試驗數(shù)據(jù)對模型平臺進行測試,表明相同模塊不同算法都能夠有效模擬冬小麥的發(fā)育期、地上部生物量和產(chǎn)量,也能夠模擬葉面積指數(shù)的動態(tài)變化,反映WMAIP平臺可用于集成多種算法,比較特定過程不同算法的差異。WU等[36]發(fā)現(xiàn)小麥鐘模型在晚播條件下的模擬性能較適期播種差,表明該算法存在一定的不足。由于春化階段是在小麥出苗后開始,改進前的小麥鐘模型模擬超晚播小麥發(fā)育期可能會出現(xiàn)無法完成春化作用現(xiàn)象,因此本文針對小麥的發(fā)育特性,將春化作用模擬改進為從播種后開始計算[37-38],提升了小麥鐘模型發(fā)育期的模擬能力。熱時法和小麥鐘模型法都考慮了光周期敏感系數(shù)和溫度敏感系數(shù),都能較好地模擬小麥發(fā)育期。在評價土壤和作物指標的模擬效果時,劉建剛等[39]認為NRMSE小于等于15%時模擬效果很好,NRMSE大于15%小于等于30%則模擬效果較好,而當NRMSE大于30%則模擬效果較差。本研究中葉面積指數(shù)的調(diào)參結(jié)果中NRMSE在19.43%~22.28%之間,調(diào)參的效果較好,但驗證結(jié)果中NRMSE在26.50%~37.75%之間,主要由干旱條件下模擬效果較差導(dǎo)致,說明需要進一步提高干旱條件下葉面積指數(shù)的模擬精度。3種生物量算法都能較好地模擬冬小麥的生物量動態(tài),在小麥生長中前期,二氧化碳同化率和光能利用效率兩種算法模擬的生物量一致性較高,而群體光合作用法模擬的生物量較大;在小麥生長后期,群體光合作用和二氧化碳同化率模擬的生物量表現(xiàn)為增長緩慢,主要由于較低的光合生產(chǎn)與較高的維持呼吸導(dǎo)致,但光能利用效率法模擬的生物量增長較快,主要由于其未單獨考慮呼吸作用導(dǎo)致(圖2d)。2017—2018年小麥生長季,3種生物量算法差異較小,而2018—2019年小麥生長季,群體光合作用法模擬的生物量較光能利用效率法和二氧化碳同化率法高,表明群體光合作用法在高輻射條件下可模擬出較高的生物量,也反映出在典型氣候條件下,不同算法的差異較小,而在極端氣候條件下,不同算法的差異較大(圖2c、2d)。3種產(chǎn)量算法模擬效果相似,主要由于3種算法的原理相似,生物量轉(zhuǎn)移法和籽粒灌漿法都考慮花后物質(zhì)生產(chǎn)以及花前物質(zhì)轉(zhuǎn)移,而收獲指數(shù)法主要考慮同化產(chǎn)物在籽粒和營養(yǎng)器官上的分配比例及水分脅迫的影響。與PALOSUO等[18]的研究結(jié)果較為一致,在大多數(shù)情況下,很難找到無誤差或最佳的模型,而模型集合的效果優(yōu)于單一模型。

    另外,在特定過程的不同算法中,存在輸入輸出的差別,需要統(tǒng)一特定過程的輸入輸出才能對不同算法進行有效集成,同時今后還需增加試驗數(shù)據(jù)對土壤氮素、根系生長及土壤溫度等模塊進行驗證。

    5 結(jié)束語

    在借鑒國內(nèi)外作物模型的理論和方法的基礎(chǔ)上,對小麥模型進行模塊化設(shè)計,構(gòu)建了小麥模型算法集成平臺WMAIP。該平臺針對特定模塊集成多種算法,有助于現(xiàn)有算法的比較與改進,并且能進一步集成模擬,降低特定過程模擬的不確定性?;谀P推脚_組成了6個具有代表性的模擬模型,采用吳橋試驗站2011—2014年及2017—2019年播期水分試驗數(shù)據(jù)對模型進行參數(shù)校準與驗證,驗證結(jié)果表明,特定模塊不同算法的模擬值與實測值均吻合較好,模擬誤差在合理范圍之內(nèi),其中發(fā)育期、地上部生物量、產(chǎn)量、土壤貯水量模擬值和實測值的歸一化均方根誤差(NRMSE)分別在0.56%~4.00%、16.13%~18.72%、12.48%~18.95%、10.78%~11.63%之間。模型集合的模擬效果優(yōu)于單一模型。發(fā)育期模塊中熱時法模擬播種至拔節(jié)階段較優(yōu),小麥鐘模型模擬播種至開花階段和播種至成熟階段較優(yōu);生物量模塊中3種算法均為模擬小麥生物量的較佳模型,但在高輻射條件下,群體光合作用法模擬的生物量較高;產(chǎn)量模塊中3種算法模擬的產(chǎn)量變化趨勢較為一致,但生物量轉(zhuǎn)移法效果略好。以上結(jié)果表明,WMAIP平臺可用于集成特定模塊的多種算法,為算法比較、算法改進以及算法集成模擬提供了一種重要的工具。

    猜你喜歡
    春化發(fā)育期實測值
    春化過程中2個時鐘蛋白CCA1和LHY通過不同順式元件激活VIN3的轉(zhuǎn)錄
    蔬菜(2022年1期)2022-11-11 01:08:01
    光照對小麥的春化效應(yīng)及其機制初探
    ±800kV直流輸電工程合成電場夏季實測值與預(yù)測值比對分析
    常用高溫軸承鋼的高溫硬度實測值與計算值的對比分析
    哈爾濱軸承(2020年1期)2020-11-03 09:16:22
    設(shè)施葡萄果實發(fā)育期管理技術(shù)要點
    河北果樹(2020年4期)2020-01-09 16:06:08
    市售純牛奶和巴氏殺菌乳營養(yǎng)成分分析
    中國奶牛(2019年10期)2019-10-28 06:23:36
    一種基于實測值理論計算的導(dǎo)航臺電磁干擾分析方法
    電子制作(2018年23期)2018-12-26 01:01:22
    春化有效積溫對結(jié)球甘藍春化相關(guān)物質(zhì)的影響
    不同春化和授粉前處理對不結(jié)球白菜抽薹及種子產(chǎn)量的影響
    基于累積熱量單位的甜瓜幼苗發(fā)育期模擬
    麻豆一二三区av精品| 男女床上黄色一级片免费看| av在线播放免费不卡| 亚洲成人精品中文字幕电影| 村上凉子中文字幕在线| 久久人人精品亚洲av| 国产真人三级小视频在线观看| avwww免费| 国产精品久久久人人做人人爽| 极品教师在线免费播放| 97超级碰碰碰精品色视频在线观看| 欧美一级a爱片免费观看看 | 嫩草影院精品99| 免费看十八禁软件| 国产精品久久电影中文字幕| 手机成人av网站| 亚洲自拍偷在线| 90打野战视频偷拍视频| 欧美色欧美亚洲另类二区| 黄色毛片三级朝国网站| 999久久久国产精品视频| 久久久久久久久久黄片| 男人的好看免费观看在线视频 | 亚洲av成人av| 99精品久久久久人妻精品| 天堂影院成人在线观看| 久久久久久免费高清国产稀缺| 色综合亚洲欧美另类图片| av视频在线观看入口| 欧美3d第一页| 欧美激情久久久久久爽电影| 男人舔女人下体高潮全视频| 久久精品91无色码中文字幕| 91大片在线观看| 国产精品亚洲一级av第二区| 亚洲18禁久久av| 高清毛片免费观看视频网站| 午夜精品一区二区三区免费看| 一边摸一边做爽爽视频免费| 日本精品一区二区三区蜜桃| 亚洲人成电影免费在线| 久久精品国产综合久久久| 天天躁夜夜躁狠狠躁躁| 久久久水蜜桃国产精品网| 亚洲精品久久国产高清桃花| 国内精品久久久久精免费| 欧美成人午夜精品| av欧美777| 国产黄色小视频在线观看| 国产成人欧美在线观看| 成人国产综合亚洲| 国产精品久久久人人做人人爽| 日本 欧美在线| 国产v大片淫在线免费观看| 国产精品电影一区二区三区| 精品国内亚洲2022精品成人| 免费观看精品视频网站| 黄色毛片三级朝国网站| 欧美激情久久久久久爽电影| 久久精品91蜜桃| 久久香蕉精品热| 黑人巨大精品欧美一区二区mp4| 婷婷六月久久综合丁香| 欧美zozozo另类| 成人亚洲精品av一区二区| 天天一区二区日本电影三级| 啦啦啦免费观看视频1| 51午夜福利影视在线观看| 又爽又黄无遮挡网站| 搞女人的毛片| 国产一区在线观看成人免费| 国模一区二区三区四区视频 | 国产乱人伦免费视频| 亚洲片人在线观看| 高清毛片免费观看视频网站| 成人一区二区视频在线观看| 精品人妻1区二区| 成在线人永久免费视频| 成熟少妇高潮喷水视频| 国语自产精品视频在线第100页| 久久精品夜夜夜夜夜久久蜜豆 | 欧美乱码精品一区二区三区| 巨乳人妻的诱惑在线观看| 日本免费a在线| 又爽又黄无遮挡网站| 特级一级黄色大片| 亚洲成人中文字幕在线播放| 久久精品成人免费网站| 亚洲在线自拍视频| 久久精品91蜜桃| 国产精品香港三级国产av潘金莲| 舔av片在线| 亚洲欧美激情综合另类| 国产99久久九九免费精品| 每晚都被弄得嗷嗷叫到高潮| 亚洲国产欧美人成| 久久久精品国产亚洲av高清涩受| 亚洲一区二区三区不卡视频| 亚洲乱码一区二区免费版| 亚洲狠狠婷婷综合久久图片| www国产在线视频色| 欧美精品啪啪一区二区三区| 99热只有精品国产| 午夜精品在线福利| 午夜久久久久精精品| 国产午夜精品论理片| 亚洲自拍偷在线| 一个人免费在线观看电影 | 亚洲最大成人中文| 1024香蕉在线观看| 制服诱惑二区| 色综合婷婷激情| 波多野结衣高清无吗| 久久久久久大精品| 中出人妻视频一区二区| 夜夜看夜夜爽夜夜摸| 99精品欧美一区二区三区四区| 麻豆国产av国片精品| 精品高清国产在线一区| 精华霜和精华液先用哪个| 亚洲aⅴ乱码一区二区在线播放 | 九色成人免费人妻av| 一本一本综合久久| 一本综合久久免费| 久久天堂一区二区三区四区| 黄色片一级片一级黄色片| 国产亚洲精品综合一区在线观看 | 午夜久久久久精精品| 999精品在线视频| 哪里可以看免费的av片| 亚洲人成电影免费在线| 在线观看舔阴道视频| 无遮挡黄片免费观看| 黄频高清免费视频| 一本精品99久久精品77| 一进一出抽搐gif免费好疼| 亚洲人成77777在线视频| 国产69精品久久久久777片 | 国产精品一及| 久久精品成人免费网站| 热99re8久久精品国产| 亚洲一区二区三区不卡视频| 少妇粗大呻吟视频| 国产黄a三级三级三级人| 久久久久九九精品影院| √禁漫天堂资源中文www| 国产精品国产高清国产av| 日韩成人在线观看一区二区三区| 长腿黑丝高跟| 精品免费久久久久久久清纯| 校园春色视频在线观看| 国产一区二区三区视频了| 亚洲av成人av| e午夜精品久久久久久久| 99国产精品一区二区三区| 久久婷婷人人爽人人干人人爱| 成人手机av| 精品高清国产在线一区| 香蕉av资源在线| 别揉我奶头~嗯~啊~动态视频| 欧美日本视频| 黄片小视频在线播放| 久99久视频精品免费| 日韩国内少妇激情av| 国产欧美日韩一区二区三| 操出白浆在线播放| 亚洲专区国产一区二区| 久久精品影院6| 日本一区二区免费在线视频| 国产精品一区二区三区四区免费观看 | 国产av一区在线观看免费| 老司机靠b影院| 日韩欧美免费精品| 国产成人av激情在线播放| 国产精品影院久久| 色精品久久人妻99蜜桃| 国产精品野战在线观看| av欧美777| 亚洲一卡2卡3卡4卡5卡精品中文| 超碰成人久久| 亚洲18禁久久av| 成人国语在线视频| 麻豆久久精品国产亚洲av| 亚洲性夜色夜夜综合| 国产精品久久电影中文字幕| 99国产极品粉嫩在线观看| 老鸭窝网址在线观看| 免费观看精品视频网站| 最好的美女福利视频网| 日本 av在线| 成人高潮视频无遮挡免费网站| 亚洲专区中文字幕在线| 免费电影在线观看免费观看| 亚洲熟妇中文字幕五十中出| 99在线人妻在线中文字幕| av超薄肉色丝袜交足视频| 可以在线观看的亚洲视频| 日本在线视频免费播放| 一本大道久久a久久精品| 亚洲 国产 在线| 亚洲精品美女久久久久99蜜臀| 精品久久久久久久末码| 特大巨黑吊av在线直播| 韩国av一区二区三区四区| 伊人久久大香线蕉亚洲五| 99精品欧美一区二区三区四区| 在线观看午夜福利视频| 日本在线视频免费播放| 欧美日韩精品网址| 听说在线观看完整版免费高清| 日本在线视频免费播放| 日韩精品青青久久久久久| 日韩av在线大香蕉| 欧美最黄视频在线播放免费| 国产不卡一卡二| 国产熟女午夜一区二区三区| 成人国产综合亚洲| a级毛片a级免费在线| 国产精品永久免费网站| 成人三级黄色视频| 别揉我奶头~嗯~啊~动态视频| 激情在线观看视频在线高清| 天天添夜夜摸| 国产精品野战在线观看| 久久天堂一区二区三区四区| 成熟少妇高潮喷水视频| 色在线成人网| 非洲黑人性xxxx精品又粗又长| 国产亚洲av高清不卡| 19禁男女啪啪无遮挡网站| 国产区一区二久久| 精品不卡国产一区二区三区| 18禁国产床啪视频网站| 51午夜福利影视在线观看| 日韩欧美 国产精品| 久久99热这里只有精品18| 久久婷婷人人爽人人干人人爱| 麻豆av在线久日| 欧美最黄视频在线播放免费| 午夜福利在线在线| 欧美黑人欧美精品刺激| 一级作爱视频免费观看| 香蕉久久夜色| 男人舔女人下体高潮全视频| 国产欧美日韩一区二区精品| 国内精品久久久久精免费| www.自偷自拍.com| 国产精品,欧美在线| 精品欧美国产一区二区三| 黑人欧美特级aaaaaa片| 曰老女人黄片| 美女大奶头视频| 亚洲av第一区精品v没综合| 久久久国产成人精品二区| 精品国产乱子伦一区二区三区| 国产亚洲精品久久久久久毛片| 国产精品久久久人人做人人爽| 51午夜福利影视在线观看| 成人三级做爰电影| 亚洲午夜精品一区,二区,三区| 午夜福利18| 免费高清视频大片| 成年人黄色毛片网站| 欧美久久黑人一区二区| 精品少妇一区二区三区视频日本电影| 欧美黑人精品巨大| 精品国产乱码久久久久久男人| 欧美成狂野欧美在线观看| 亚洲国产中文字幕在线视频| 亚洲av片天天在线观看| 亚洲全国av大片| 日韩欧美三级三区| 老鸭窝网址在线观看| 色尼玛亚洲综合影院| 777久久人妻少妇嫩草av网站| 五月伊人婷婷丁香| а√天堂www在线а√下载| 久久草成人影院| 怎么达到女性高潮| 免费无遮挡裸体视频| 丁香欧美五月| 神马国产精品三级电影在线观看 | 夜夜夜夜夜久久久久| 窝窝影院91人妻| 国产爱豆传媒在线观看 | 国产精品电影一区二区三区| 国产三级在线视频| 亚洲欧美一区二区三区黑人| 法律面前人人平等表现在哪些方面| 久久国产精品影院| 男女床上黄色一级片免费看| 国产区一区二久久| 欧美 亚洲 国产 日韩一| 久久精品夜夜夜夜夜久久蜜豆 | 99热这里只有是精品50| www日本在线高清视频| 国产高清视频在线播放一区| 黄片小视频在线播放| 免费电影在线观看免费观看| 性色av乱码一区二区三区2| 欧美极品一区二区三区四区| 在线十欧美十亚洲十日本专区| 男女下面进入的视频免费午夜| 久久亚洲真实| xxx96com| 大型黄色视频在线免费观看| 1024香蕉在线观看| 麻豆成人午夜福利视频| 国产久久久一区二区三区| 国产精品国产高清国产av| 国产欧美日韩一区二区三| 久久久久国内视频| 色尼玛亚洲综合影院| www日本在线高清视频| 波多野结衣高清无吗| 麻豆国产97在线/欧美 | 成人18禁在线播放| 欧美绝顶高潮抽搐喷水| 色在线成人网| 天天添夜夜摸| 99久久精品国产亚洲精品| 国内精品久久久久精免费| 欧美日韩黄片免| 国产精品自产拍在线观看55亚洲| 久久久国产成人精品二区| 欧美丝袜亚洲另类 | 黑人巨大精品欧美一区二区mp4| 欧美性长视频在线观看| 国产人伦9x9x在线观看| 午夜激情av网站| 99re在线观看精品视频| 国产欧美日韩精品亚洲av| tocl精华| 高清在线国产一区| 亚洲精品一卡2卡三卡4卡5卡| 宅男免费午夜| 最近最新免费中文字幕在线| 国产精品 欧美亚洲| 熟女电影av网| 久久香蕉国产精品| 亚洲男人的天堂狠狠| 男人舔女人的私密视频| av在线天堂中文字幕| 黄色丝袜av网址大全| 日韩欧美精品v在线| 日韩av在线大香蕉| 99国产精品99久久久久| 久久人人精品亚洲av| 99国产精品一区二区蜜桃av| 国产又黄又爽又无遮挡在线| 天天躁夜夜躁狠狠躁躁| 99久久无色码亚洲精品果冻| 精品国内亚洲2022精品成人| 国产又色又爽无遮挡免费看| 久久精品成人免费网站| 色综合婷婷激情| 18禁国产床啪视频网站| 午夜免费激情av| 日韩大尺度精品在线看网址| 少妇粗大呻吟视频| 舔av片在线| 女人被狂操c到高潮| 岛国在线观看网站| 久久亚洲精品不卡| cao死你这个sao货| 99热这里只有是精品50| 亚洲av成人不卡在线观看播放网| 免费高清视频大片| 香蕉国产在线看| 成人18禁在线播放| 精品不卡国产一区二区三区| 国产精品,欧美在线| 久久久精品国产亚洲av高清涩受| 小说图片视频综合网站| 久久香蕉国产精品| 蜜桃久久精品国产亚洲av| 麻豆久久精品国产亚洲av| 亚洲精品美女久久久久99蜜臀| 精品福利观看| 美女高潮喷水抽搐中文字幕| 黄色片一级片一级黄色片| 欧美日韩瑟瑟在线播放| av国产免费在线观看| 桃红色精品国产亚洲av| 一本精品99久久精品77| 此物有八面人人有两片| 香蕉丝袜av| 久久久久久久久免费视频了| 少妇人妻一区二区三区视频| 久久亚洲真实| 国产精品 欧美亚洲| 国产精品亚洲一级av第二区| 午夜免费观看网址| 成人手机av| 午夜精品一区二区三区免费看| 国产乱人伦免费视频| 精品乱码久久久久久99久播| 精品一区二区三区四区五区乱码| 精品高清国产在线一区| 国产精品1区2区在线观看.| 日韩欧美一区二区三区在线观看| 国产精品精品国产色婷婷| 午夜福利成人在线免费观看| 真人做人爱边吃奶动态| 黑人巨大精品欧美一区二区mp4| 天天躁狠狠躁夜夜躁狠狠躁| 淫秽高清视频在线观看| 国产精品久久视频播放| 欧美另类亚洲清纯唯美| 99久久99久久久精品蜜桃| 不卡一级毛片| 亚洲激情在线av| 一个人免费在线观看电影 | 欧美日韩乱码在线| 一本一本综合久久| 在线国产一区二区在线| 亚洲中文字幕一区二区三区有码在线看 | 久久精品国产99精品国产亚洲性色| 精品久久久久久久久久久久久| 久久久久免费精品人妻一区二区| 悠悠久久av| 搞女人的毛片| 18禁观看日本| 一二三四社区在线视频社区8| 午夜福利免费观看在线| 两人在一起打扑克的视频| 啦啦啦观看免费观看视频高清| 全区人妻精品视频| 亚洲av日韩精品久久久久久密| 久久久久国产一级毛片高清牌| 曰老女人黄片| 国产成年人精品一区二区| 久99久视频精品免费| netflix在线观看网站| 久久99热这里只有精品18| 国产成人aa在线观看| 麻豆一二三区av精品| 少妇人妻一区二区三区视频| 欧美日韩亚洲综合一区二区三区_| 熟女少妇亚洲综合色aaa.| 久久精品91蜜桃| 精品久久久久久久毛片微露脸| www.精华液| 亚洲电影在线观看av| 久久久久久国产a免费观看| 在线观看日韩欧美| 中文字幕高清在线视频| 草草在线视频免费看| 97超级碰碰碰精品色视频在线观看| 不卡av一区二区三区| 亚洲真实伦在线观看| 老鸭窝网址在线观看| 黄色a级毛片大全视频| 最好的美女福利视频网| 这个男人来自地球电影免费观看| 男女那种视频在线观看| 日韩精品青青久久久久久| 国产高清有码在线观看视频 | 久久 成人 亚洲| 2021天堂中文幕一二区在线观| 人成视频在线观看免费观看| 成人一区二区视频在线观看| 国产一区二区在线观看日韩 | e午夜精品久久久久久久| 日韩精品中文字幕看吧| 亚洲激情在线av| 久久午夜亚洲精品久久| 日日夜夜操网爽| 两个人免费观看高清视频| 一本一本综合久久| 国产亚洲av高清不卡| 久久久久久大精品| 久久草成人影院| 全区人妻精品视频| 国产一区二区三区在线臀色熟女| 精品欧美一区二区三区在线| 午夜两性在线视频| 亚洲av日韩精品久久久久久密| 成人三级做爰电影| 夜夜看夜夜爽夜夜摸| 最近视频中文字幕2019在线8| 国产av一区在线观看免费| 一进一出抽搐gif免费好疼| 男女之事视频高清在线观看| 日本三级黄在线观看| 老汉色∧v一级毛片| 免费观看人在逋| 精品久久久久久久久久免费视频| 午夜免费成人在线视频| 国产成人精品久久二区二区91| 欧美久久黑人一区二区| 精品久久久久久,| 青草久久国产| 老司机深夜福利视频在线观看| 亚洲18禁久久av| 成人国语在线视频| 麻豆久久精品国产亚洲av| 51午夜福利影视在线观看| 久久国产精品人妻蜜桃| 午夜福利在线在线| 又大又爽又粗| 啦啦啦观看免费观看视频高清| 丁香欧美五月| 国产熟女午夜一区二区三区| 国产伦一二天堂av在线观看| 免费电影在线观看免费观看| 久久久久久亚洲精品国产蜜桃av| 亚洲电影在线观看av| 最近视频中文字幕2019在线8| 国内久久婷婷六月综合欲色啪| 女人高潮潮喷娇喘18禁视频| 亚洲精品久久国产高清桃花| 一本综合久久免费| 黄片小视频在线播放| 亚洲国产中文字幕在线视频| 久久久久久九九精品二区国产 | 天堂影院成人在线观看| 国产精品影院久久| 午夜免费激情av| 窝窝影院91人妻| 国产高清激情床上av| 中文字幕精品亚洲无线码一区| 亚洲午夜精品一区,二区,三区| 久久精品综合一区二区三区| 日本在线视频免费播放| 亚洲精品一区av在线观看| 日本一二三区视频观看| 日韩国内少妇激情av| 丰满人妻一区二区三区视频av | 国内精品久久久久精免费| 五月伊人婷婷丁香| 精品熟女少妇八av免费久了| 婷婷精品国产亚洲av| 午夜福利免费观看在线| 国产成人啪精品午夜网站| 久久精品人妻少妇| 中亚洲国语对白在线视频| 欧美日本视频| 亚洲精品国产精品久久久不卡| 日韩欧美国产在线观看| 日韩成人在线观看一区二区三区| 后天国语完整版免费观看| 午夜激情av网站| 91成年电影在线观看| 法律面前人人平等表现在哪些方面| 久久午夜综合久久蜜桃| 久久欧美精品欧美久久欧美| 母亲3免费完整高清在线观看| 波多野结衣巨乳人妻| 亚洲激情在线av| 97超级碰碰碰精品色视频在线观看| 男人舔奶头视频| 欧美成人一区二区免费高清观看 | 国产97色在线日韩免费| 99国产精品99久久久久| 国产精品香港三级国产av潘金莲| 热99re8久久精品国产| 国产视频一区二区在线看| 18禁美女被吸乳视频| 亚洲成av人片在线播放无| 亚洲精品美女久久av网站| 男女那种视频在线观看| 亚洲美女黄片视频| 久久久久久大精品| 一级毛片女人18水好多| 亚洲国产欧洲综合997久久,| 在线播放国产精品三级| 美女免费视频网站| 非洲黑人性xxxx精品又粗又长| 久久欧美精品欧美久久欧美| 免费在线观看黄色视频的| 女同久久另类99精品国产91| 欧美日韩亚洲综合一区二区三区_| 中文资源天堂在线| 欧美一区二区国产精品久久精品 | 亚洲国产欧美网| 久久久久久久久久黄片| 91麻豆av在线| 伊人久久大香线蕉亚洲五| 女警被强在线播放| 97碰自拍视频| 嫁个100分男人电影在线观看| 亚洲人成网站高清观看| 日韩 欧美 亚洲 中文字幕| 亚洲色图 男人天堂 中文字幕| 黄色成人免费大全| а√天堂www在线а√下载| 国产伦在线观看视频一区| 成在线人永久免费视频| 精品高清国产在线一区| 琪琪午夜伦伦电影理论片6080| 午夜日韩欧美国产| 在线观看美女被高潮喷水网站 | 男人的好看免费观看在线视频 | 69av精品久久久久久| 国产精品1区2区在线观看.| 丝袜人妻中文字幕| 国产一区二区三区视频了| 成人av在线播放网站| 国产av在哪里看| 久久久国产成人精品二区| 人人妻人人澡欧美一区二区| 国产欧美日韩一区二区精品| 妹子高潮喷水视频| 青草久久国产| 日本成人三级电影网站| 人成视频在线观看免费观看| 看黄色毛片网站| 久久久久久久久久黄片| 日本一本二区三区精品| 波多野结衣高清无吗| 国产成人aa在线观看| 国产精品98久久久久久宅男小说| 久久久国产精品麻豆| 精品久久久久久久人妻蜜臀av|