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

    基于Penman-Monteith模型的低丘紅壤區(qū)稻田蒸散模擬?

    2020-04-22 15:13:34文建川景元書韓麗娟
    中國農(nóng)業(yè)氣象 2020年4期
    關(guān)鍵詞:冠層實(shí)測(cè)值表達(dá)式

    文建川,景元書**,韓麗娟

    (1.南京信息工程大學(xué)氣象災(zāi)害預(yù)報(bào)預(yù)警與評(píng)估協(xié)同創(chuàng)新中心/應(yīng)用氣象學(xué)院,南京 210044;2.國家氣象中心,北京 100081)

    生態(tài)系統(tǒng)中植被發(fā)育以及系統(tǒng)生產(chǎn)力都受到蒸散的直接影響,同時(shí)它不僅是土壤-植被-大氣連續(xù)體中水分傳輸中的重要環(huán)節(jié),也是生態(tài)系統(tǒng)中能量平衡的主要組成部分[1]。據(jù)統(tǒng)計(jì),以蒸散的形式消耗農(nóng)業(yè)用水占到90%[2-3],蒸散的準(zhǔn)確估算有利于作物水分生產(chǎn)力提升以及改善農(nóng)田或流域灌溉管理[4-5],因此,蒸散觀測(cè)和模擬一直受到生態(tài)氣象領(lǐng)域?qū)W者的廣泛關(guān)注。

    近幾十年來,蒸散估算模型眾多,在該研究領(lǐng)域較為有名的是以能量平衡和微氣候?qū)W方法為原理的Penman-Monteith模型(簡稱PM模型),該模型模擬精度主要取決于冠層阻力這一關(guān)鍵性參數(shù)[7]。國內(nèi)外研究者基于冠層阻力公式,對(duì)其與PM模型結(jié)合后模擬蒸散量開展了不少研究。Irmak 采用非線性回歸方法設(shè)計(jì)了7個(gè)形式相同,氣象因子數(shù)量不同導(dǎo)致復(fù)雜程度不一樣的模型,之后利用該模型對(duì)大豆蒸散量進(jìn)行模擬,取得了較好的結(jié)果[8-9]。吳林等[10]將Irmak 公式應(yīng)用到PM模型中,可以很好地估算黑河綠洲區(qū)玉米半小時(shí)尺度上的蒸散量。趙華等[11]采用Jarvis 公式計(jì)算了南京地區(qū)2012年和2013年水稻拔節(jié)-成熟期的冠層阻力,并將其耦合到PM模型中模擬水稻該時(shí)期的蒸散,兩年的擬合度都超過0.95。Li 等[12]采用PM_Jarvis模型時(shí),基于玉米葉面積指數(shù)大小將整個(gè)生育期劃分為兩個(gè)階段,分段擬合模型中的經(jīng)驗(yàn)系數(shù),對(duì)比后發(fā)現(xiàn)以葉面積指數(shù)為0.5 分段擬合的系數(shù)進(jìn)行蒸散模擬的效果最佳。Srivastava等[13]比較了Monteith、Katerji-Perrier、Todorovic和Jarvis 四種冠層阻力公式應(yīng)用在PM模型中的模擬效果,結(jié)果表明Jarvis和Todorovic 都較好地估算了半濕潤地區(qū)玉米的蒸散。Jarvis和Irmak 公式中包含了較多經(jīng)驗(yàn)系數(shù),有學(xué)者分別將其它研究的經(jīng)驗(yàn)系數(shù)值和通過非線性擬合得到的系數(shù)值代入到模型中估算蒸散量,比較后發(fā)現(xiàn)借用其它研究的經(jīng)驗(yàn)系數(shù)得到的模擬值與實(shí)測(cè)值相關(guān)性較差[14]。

    因?yàn)椴煌貐^(qū)氣候差異大,這兩種模型同屬為經(jīng)驗(yàn)?zāi)P停琂arvis和Irmak 公式應(yīng)用在PM模型中估算蒸散量,能否取得較高的模擬精度,是否適用于低丘紅壤區(qū)稻田蒸散研究,有待進(jìn)一步探索。因此,需要針對(duì)試驗(yàn)區(qū)及作物,挑選出最為合適的冠層阻力公式,得到模擬精度高的蒸散模型。另外,Jarvis模型中各脅迫函數(shù)存在多種表達(dá)式,導(dǎo)致有多種組合形式,較多學(xué)者在使用該模型時(shí),基本上都是選取其中一種組合形式進(jìn)行模擬,但是不同的表達(dá)式所涉及的參數(shù)數(shù)量不同,可能導(dǎo)致最后的模擬精度有差別。趙華[15]采用了Jarvis模型中的幾種組合形式,對(duì)冬小麥和水稻蒸散量進(jìn)行研究,但只選取了各生育期中典型晴天數(shù)據(jù),未對(duì)作物全生育期進(jìn)行驗(yàn)證。因此,本研究以波文比儀監(jiān)測(cè)該地區(qū)稻田的蒸散量作為實(shí)測(cè)值,對(duì)比分析PM_Jarvis1-8(脅迫函數(shù)的8種組合形式)與PM_Irmak 共9個(gè)模型的蒸散模擬效果,旨在選出更適用于低丘紅壤區(qū)稻田的實(shí)際蒸散模型,準(zhǔn)確估算蒸散量,以利于制定更合理的灌溉計(jì)劃和提高稻田的水分利用效率。

    1 材料與方法

    1.1 試驗(yàn)站點(diǎn)概況

    試驗(yàn)站點(diǎn)位于江西省鷹潭市中國科學(xué)院紅壤生態(tài)試驗(yàn)站(116°55'E,28°15'N),該站地處武夷山至潘陽湖平原的過渡地帶鷹潭盆地,屬亞熱帶濕潤季風(fēng)氣候,年平均氣溫17.6℃,降水量1788.8mm,日照時(shí)數(shù)1739.4h。試驗(yàn)觀測(cè)時(shí)間為2014年和2015年晚稻生長季(7-10月),水稻移栽后,依照當(dāng)?shù)卮筇锷a(chǎn)進(jìn)行田間管理。觀察期內(nèi)水稻生育期見表1。

    表1 試驗(yàn)區(qū)2014年和2015年水稻生育期Table1 Growth stage of rice in experimental area in 2014 and 2015

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

    波文比儀架設(shè)在稻田中部,該系統(tǒng)由NR-Lite凈輻射傳感器(Kipp&Zonen,Netherland)、兩層空氣溫濕度傳感器( HMP155A)、CR1000 數(shù)據(jù)采集器、風(fēng)速儀( 010C-1)組成,自動(dòng)采集1.5m 高度差之間的溫濕度差,計(jì)算得到波文比值。土壤熱通量由埋深為 5cm的土壤熱通量板(Hukseflux,HFP01,Netherland)測(cè)得。氣象數(shù)據(jù)由布設(shè)在農(nóng)田中部的自動(dòng)氣象站(HOBO U30,ONSET,USA)監(jiān)測(cè),觀測(cè)內(nèi)容包括相對(duì)濕度、氣溫、2m 風(fēng)速、氣壓等小氣候數(shù)據(jù),觀測(cè)頻率為20min。

    1.3 數(shù)據(jù)處理方法

    波文比測(cè)定水熱通量是以地表能量平衡方程為依據(jù),即

    式中,Rn 為凈輻射通量(W·m-2),λET 為潛熱通量(W·m-2),H 為感熱通量(W·m-2),G 為土壤熱通量(W·m-2)。其中,Rn和G可以實(shí)測(cè)得到,λET和H 通過計(jì)算得到。

    波文比(β)為感熱通量和潛熱通量的比值,根據(jù)莫寧-奧布霍夫理論,假定水汽和熱量的湍流擴(kuò)散系數(shù)Kh與Kw相等,則有

    式中,γ 為干濕表常數(shù),取值為0.067kPa·℃;ΔT、Δe 分別為波文比系統(tǒng)的兩層觀測(cè)高度的溫度差和水汽壓差。

    波文比儀測(cè)定出兩層高度的溫度差和水汽壓差、凈輻射通量以及土壤熱通量板測(cè)定的土壤熱通量數(shù)據(jù),結(jié)合式(1)和式(2),可得潛熱通量λET 為

    式中,λ 為水的汽化潛熱,取值為2.45MJ·kg-1,ET 為蒸散量(mm)。

    因?yàn)槿粘?、日落以及降水天氣,可能造成誤差,所以根據(jù)Perez 等[16]提供的方法確定無效β 值,將其剔除并進(jìn)行插補(bǔ)。計(jì)算式為

    式中,δ1和δ2分別為溫、濕度傳感器的測(cè)量精度。

    1.4 模型介紹

    1.4.1 冠層阻力計(jì)算模型

    冠層阻力rc,采用兩種模型進(jìn)行計(jì)算。

    (1)PM_Jarvis模型

    Jarvis 認(rèn)為冠層阻力與環(huán)境因子具有協(xié)同關(guān)系,而且不同環(huán)境因子對(duì)冠層阻力的影響?yīng)毩ⅲ虼?,提出了冠層阻力與環(huán)境因子之間的關(guān)系式,即

    式中,rmin為最小冠層阻力,根據(jù)前人[17]研究結(jié)果,水稻在拔節(jié)期、孕穗期、抽穗期、乳熟期的rmin分別取為40、45、50、60s·m-1,LAIe為有效葉面積指數(shù),F(xiàn)(Rn)、F(D)、F(T)分別為太陽輻射脅迫函數(shù)、飽和水汽壓差脅迫函數(shù)、溫度脅迫函數(shù),因以上3種脅迫函數(shù)均有多種表達(dá)式,導(dǎo)致存在不同組合形式,不同的表達(dá)形式有可能影響模擬精度,因此,每個(gè)脅迫函數(shù)均選取其中兩種表達(dá)式。

    有效葉面積指數(shù)LAIe為[18]

    太陽輻射脅迫函數(shù)分別采用線性關(guān)系式[19]和指數(shù)關(guān)系式[12],即

    飽和水汽壓差脅迫函數(shù)分別采用雙曲線關(guān)系式[19]和線性關(guān)系式[12],即

    空氣溫度脅迫函數(shù)分別采用以下兩種關(guān)系式[12,19],即

    式中,Rn為凈輻射通量(W·m-2),D 為飽和水汽壓差(kPa),T 為空氣溫度(℃),a1、a2、a3、a4、b1、b2、b3、c1、c2、c3、c4為模型經(jīng)驗(yàn)系數(shù)。

    基于上述各脅迫函數(shù)的響應(yīng)表達(dá)式,得到8種不同的Jarvis模型組合形式(表2)。

    (2)PM_Irmak模型

    Irmak 等采用廣義非線性回歸方法,構(gòu)建了關(guān)于氣象因子、葉面積指數(shù)與空氣動(dòng)力學(xué)阻力的冠層阻力模型,即

    式中,rc為冠層阻力(s·m-1),Rn 為凈輻射(W·m-2),T 為空氣溫度(℃),RH 為相對(duì)濕度,U 為風(fēng)速(m·s-1),LAI 為葉面積指數(shù),ra為空氣動(dòng)力學(xué)阻力(s·m-1),按Rana 等提出的公式計(jì)算[20],即

    表2 PM_Jarvis 冠層阻力模型中3種脅迫函數(shù)表達(dá)式的不同組合Table2 Different combinations of three stress function expressions of the PM_Jarvis canopy resistance model

    式中,z 為參考高度(m),取值為1.5m;h 為冠層高度(m),以實(shí)測(cè)株高代替;k 為Von Karman常數(shù),取值為0.41;u 為風(fēng)速(m·s-1),z0=0.13h,d=0.63h。

    1.4.2 蒸散量計(jì)算模型

    Monteith 在Penman 基礎(chǔ)上,引入了空氣動(dòng)力學(xué)阻力和冠層阻力等參數(shù),提出以能量平衡和水汽擴(kuò)散理論為基礎(chǔ)的適用于作物蒸騰計(jì)算的阻力模式,即Penman-Monteith模型。

    式中,λET 為潛熱通量(W·m-2);Rn 為凈輻射通量(W·m-2);G 為土壤熱通量(W·m-2);Δ 為飽和水汽壓曲線斜率(kPa·K-1);ρ 為空氣密度,取值為 1.29kg·m-3;Cp 為空氣定壓比熱,取值為1004J·kg-1·K-1;D 為飽和水汽壓差(kPa);γ 為干濕表常數(shù),取值為0.067kPa·K-1;rc為冠層阻力(s·m-1);ra為空氣動(dòng)力學(xué)阻力(s·m-1)。

    1.5 模型精度分析

    以波文比儀測(cè)定的蒸散量ET 作為參照,分別把不同冠層阻力模型應(yīng)用到Penman-Monteith 公式中,再將模擬得到的ET 與實(shí)測(cè)值比較。模型的模擬精度評(píng)價(jià)指標(biāo)包括均方根誤差(RMSE)、平均絕對(duì)誤差(MAE)、納什效率系數(shù)(NSE)和一致性指數(shù)(IA)、確定系數(shù)(R2)。

    式中,Pi為模擬蒸散量(mm·h-1),Oi為實(shí)測(cè)蒸散量(mm·h-1),為實(shí)測(cè)蒸散量的平均值(mm·h-1),為模擬蒸散量的平均值(mm·h-1),n 為個(gè)數(shù)。

    RMSE和MAE 反映模擬值與實(shí)測(cè)值之間的差異,數(shù)值越接近0,表明誤差越?。籒SE 越接近1,表明模型可信度越高;IA 反映模擬值與實(shí)測(cè)值的符合程度,數(shù)值越接近1,說明模型的預(yù)測(cè)效果越好。R2越接近1,表明擬合程度越好。

    2 結(jié)果與分析

    2.1 PM_Jarvis和PM_Irmak 冠層阻力模型模擬效果

    基于2014年波文比儀監(jiān)測(cè)的水稻全生育期實(shí)測(cè)蒸散量與同階段的氣象數(shù)據(jù),將表2中組合好的8種Jarvis模型及Irmak模型分別與Penman-Monteith模型結(jié)合,簡稱PM_Jarvis1-8模型(脅迫函數(shù)8種組合形式)、PM_Irmak模型,利用1stOpt 軟件進(jìn)行非線性擬合(http: //www.7d-soft.com/en/),通過全局優(yōu)化算法,得到各模型的經(jīng)驗(yàn)系數(shù)值(表3)。統(tǒng)計(jì)這9種蒸散模型的誤差及相關(guān)性統(tǒng)計(jì)參量,結(jié)果見表4。從表4中可知,不同組合形式的PM_Jarvis模型精度指標(biāo)都存在一定的差異,其中PM_Jarvis6模型精度最低,其RMSE、MAE 最高,NSE、IA 及R2最低。PM_Irmak模型的RMSE、MAE 在9種模型中最低,NSE、IA 及R2較其余模型更接近1,表明該模型模擬精度最高。就PM_Jarvis模型而言,8種組合形式模型的模擬精度存在差別,將模型進(jìn)行兩兩對(duì)比,結(jié)果表明,對(duì)于太陽輻射脅迫函數(shù),F(xiàn)1(Rn)的表達(dá)式在應(yīng)用于該模型中較F2(Rn)效果更佳;對(duì)于飽和水汽壓差脅迫函數(shù)和空氣溫度脅迫函數(shù),F(xiàn)1(D)、F1(T)的表達(dá)式分別優(yōu)于F2(D)和F2(T)。

    表3 PM_Jarvis(式7-12)和PM_Irmak(式13)模型中各參數(shù)的擬合值Table3 Simulation value of parameters in PM_Jarvis(Equation7-12) and PM_Irmak models(Equation13)

    表4 PM_Jarvis和PM_Irmak 模擬2014年逐時(shí)蒸散的統(tǒng)計(jì)參量Table4 Statistical parameters of hourly evapotranspiration simulated by PM _Jarvis and PM_Irmak models

    2.2 PM_Jarvis和PM_Irmak 模擬結(jié)果驗(yàn)證

    利用2015年7月23日-10月31日的氣象資料和水稻生理數(shù)據(jù),代入9種模型計(jì)算蒸散量,并與波文比儀監(jiān)測(cè)的該時(shí)段蒸散量比較,進(jìn)一步驗(yàn)證模型。從圖1看出,PM_Irmak模型模擬的蒸散量與實(shí)測(cè)蒸散量基本沿1:1 線分布,且散點(diǎn)集中,表明模擬值與實(shí)測(cè)值較接近;其余8種PM_Jarvis模型散點(diǎn)大都位于1:1 線下方,總體上表現(xiàn)出偏低趨勢(shì)。為了更全面地體現(xiàn)模型模擬效果,在全生育期及各生育階段進(jìn)行精度分析,結(jié)果見表5和圖2。由表5可以看出,PM_Irmak模型在估算水稻全生育期蒸散量中,各評(píng)價(jià)指標(biāo)整體均優(yōu)于其余模型,同樣,就PM_Jarvis模型而言,脅迫函數(shù)表達(dá)式不同導(dǎo)致8種PM_Jarvis模型精度不一,與2014年結(jié)果一致。

    圖1 PM-Jarvis和PM-Irmak模型模擬逐時(shí)蒸散量與實(shí)測(cè)蒸散量對(duì)比Fig.1 Comparison of hourly evapotranspiration( ET) between measured and simulated by PM-Jarvis and PM-Irmak models

    表5 PM-Jarvis和PM-Irmak模型模擬的稻田全生育期蒸散量精度驗(yàn)證(2015年)Table5 Accuracy of evapotranspiration during the whole growth stage of rice paddy simulated by PM-Jarvis and PM-Irmak models (2015)

    將水稻的全生育期分為4個(gè)階段,由圖2可見,8個(gè)PM_Jarvis模型的4個(gè)評(píng)價(jià)指標(biāo)在返青-分蘗期與后 3個(gè)生育期有明顯差距,在該生育階段,PM_Irmak模型的均方根誤差(RMSE)和平均絕對(duì)誤差(MAE)顯著低于8個(gè)PM_Jarvis模型,納什效率系數(shù)(NSE)和一致性指數(shù)(IA)明顯較高;在后3個(gè)生育時(shí)期,PM_Irmak模型的RMSE和MAE雖然不是最低,NSE和IA 也不是最高,但與其余模型差距較小。所以認(rèn)為,PM_Irmak模型總體上對(duì)水稻整個(gè)生長階段蒸散量的模擬效果較好。

    圖2 PM-Jarvis和PM-Irmak模型模擬的稻田不同生育期蒸散量精度驗(yàn)證Fig.2 Accuracy of evapotranspiration in the different growth stages of rice paddy simulated by PM-Jarvis and PM-Irmak models

    2.3 PM-Jarvis和PM-Irmak模型模擬蒸散量的日動(dòng)態(tài)變化

    在水稻4個(gè)生育時(shí)期中,每一階段均選擇連續(xù)9d 用來觀察9個(gè)模型估算的蒸散與實(shí)測(cè)蒸散的日變化,結(jié)果見圖3(下頁)。由圖可以看出,在返青-分蘗期(7月23-31日),PM_Jarvis1-8模型呈現(xiàn)出低估現(xiàn)象,并且低估程度嚴(yán)重,蒸散估算值不及實(shí)際蒸散量的50%;而PM_Irmak模型在該生育時(shí)期,除了在個(gè)別時(shí)間點(diǎn)(7月27日)蒸散模擬值稍偏高,其余時(shí)段均與實(shí)測(cè)值較接近,模擬效果明顯優(yōu)于PM_Jarvis1-8模型。在后3個(gè)生育時(shí)期,各模型模擬的日變化動(dòng)態(tài)趨勢(shì)同步,數(shù)值上存在差異,但差異不大,均與實(shí)測(cè)蒸散的日動(dòng)態(tài)曲線基本吻合??偟膩碚f,PM_Irmak模型在水稻整個(gè)生長階段都能較好地估算蒸散量,所以更適用于紅壤地區(qū)稻田蒸散量模擬。

    3 結(jié)論與討論

    3.1 結(jié)論

    本研究利用不同冠層阻力公式的PM模型估算蒸散量,與波文比儀監(jiān)測(cè)的實(shí)測(cè)蒸散量對(duì)比,PM_Jarvis 1-8和PM_Irmak 共9種模型中,PM_Irmak模型的均方根誤差、平均絕對(duì)誤差比8種PM_Jarvis模型的均偏小,確定系數(shù)、納什效率系數(shù)和一致性指數(shù)相比更接近 1,整體模擬效果最好。8種PM_Jarvis模型在水稻返青分蘗期的初始階段存在明顯低估現(xiàn)象,而PM_Irmak模型在水稻所有生育階段的模擬值均與蒸散日變化動(dòng)態(tài)值接近。

    3.2 討論

    圖3 PM_Jarvis和PM_Irmak模型模擬值與實(shí)測(cè)值的蒸散日變化動(dòng)態(tài)Fig.3 Daily variation of observed and simulated evapotranspiration by PM_Jarvis and PM_Irmak models

    因Jarvis 公式脅迫函數(shù)表達(dá)式不止一種,本研究每個(gè)函數(shù)各選取兩種,組成了8種組合形式,發(fā)現(xiàn)PM_Jarvis1-8 模擬效果存在不同程度的差異,其中PM_Jarvis6模型(3個(gè)脅迫函數(shù)分別采用F2(Rn)、F2(D)、 F2(T)的表達(dá)式)精度最低。雖然這8種模型為同一形式,但脅迫函數(shù)的表達(dá)式不同,對(duì)最后蒸散量的模擬產(chǎn)生了較大影響。趙華[15]在進(jìn)行水稻和冬小麥蒸散量研究時(shí),也將不同組合形式的Jarvis 公式應(yīng)用到PM模型中,發(fā)現(xiàn)不同模型的均方根誤差和擬合度相差較大。8種PM_Jarvis模型,在返青期均呈現(xiàn)出明顯低于實(shí)測(cè)蒸散的現(xiàn)象。Li 等[21]在利用PM_Jarvis模型對(duì)葡萄園蒸散進(jìn)行模擬研究時(shí),同樣指出在葡萄生長前期階段模擬蒸散也明顯低于實(shí)測(cè)值,隨著葡萄逐漸生長發(fā)育,二者的差異在減小,擬合效果才得以提高。因?yàn)镴arvis 把每個(gè)環(huán)境變量對(duì)作物氣孔產(chǎn)生的影響假定為一種相互獨(dú)立的行為,并認(rèn)為作物所有葉片的氣孔開閉都同時(shí)具有一致性,基于這樣一種假設(shè)提出了該模型[22],但實(shí)際上環(huán)境變量之間存在著互相影響的關(guān)系,并且較為復(fù)雜,所以這有可能導(dǎo)致了該模型在作物部分生育階段蒸散模擬效果較差。

    PM_Irmak模型在返青-分蘗期各評(píng)價(jià)指標(biāo)明顯優(yōu)于8種PM_Jarvis模型,后3個(gè)生育時(shí)期該評(píng)價(jià)指標(biāo)雖不是最優(yōu),但9種模型之間差異較小,所以最終 PM_Irmak模型仍在水稻全生育期中較 8種PM_Jarvis模型有更好的模擬效果。Li 等在研究干旱地區(qū)玉米蒸散中發(fā)現(xiàn),PM_Irmak模型在玉米葉面積指數(shù)小于2的時(shí)期模擬蒸散量的誤差要小于葉面積指數(shù)大于2的時(shí)期,雖然在LAI 大于2時(shí)期中模型的精度略低于該研究中的其它冠層阻力模型,但是就整個(gè)生育期而言,PM_Irmak模型的均方根誤差和平均相對(duì)誤差都明顯要小[23],這與本研究結(jié)果一致。PM_Irmak模型在水稻4個(gè)生育階段估算的蒸散量與實(shí)測(cè)值的日變化動(dòng)態(tài)均較為吻合且數(shù)值接近,精度總體來說較高,說明該模型適用于本試驗(yàn)區(qū),能夠更好地模擬蒸散量。

    猜你喜歡
    冠層實(shí)測(cè)值表達(dá)式
    基于低空遙感的果樹冠層信息提取方法研究
    基于激光雷達(dá)的樹形靶標(biāo)冠層葉面積探測(cè)模型研究
    ±800kV直流輸電工程合成電場(chǎng)夏季實(shí)測(cè)值與預(yù)測(cè)值比對(duì)分析
    安徽省淮南森林冠層輻射傳輸過程的特征
    常用高溫軸承鋼的高溫硬度實(shí)測(cè)值與計(jì)算值的對(duì)比分析
    哈爾濱軸承(2020年1期)2020-11-03 09:16:22
    一個(gè)混合核Hilbert型積分不等式及其算子范數(shù)表達(dá)式
    表達(dá)式轉(zhuǎn)換及求值探析
    市售純牛奶和巴氏殺菌乳營養(yǎng)成分分析
    中國奶牛(2019年10期)2019-10-28 06:23:36
    淺析C語言運(yùn)算符及表達(dá)式的教學(xué)誤區(qū)
    一種基于實(shí)測(cè)值理論計(jì)算的導(dǎo)航臺(tái)電磁干擾分析方法
    電子制作(2018年23期)2018-12-26 01:01:22
    久久精品国产亚洲网站| 午夜福利视频精品| 国产一区二区三区av在线| 日韩一本色道免费dvd| 亚洲国产日韩一区二区| 亚洲无线观看免费| 老司机亚洲免费影院| 啦啦啦在线观看免费高清www| 你懂的网址亚洲精品在线观看| 色94色欧美一区二区| 久久久久久久久久久免费av| 黄色日韩在线| 日韩精品免费视频一区二区三区 | 亚洲av免费高清在线观看| 免费观看无遮挡的男女| 亚洲欧美成人精品一区二区| 一区二区三区免费毛片| 亚洲av.av天堂| 日韩制服骚丝袜av| 国产探花极品一区二区| 91成人精品电影| 曰老女人黄片| 黄色毛片三级朝国网站 | 青春草国产在线视频| 久久6这里有精品| av福利片在线| 亚洲av欧美aⅴ国产| 国产av码专区亚洲av| 乱码一卡2卡4卡精品| 水蜜桃什么品种好| 一区二区三区四区激情视频| 成人国产av品久久久| 国产成人午夜福利电影在线观看| 亚洲精品国产av蜜桃| 麻豆成人午夜福利视频| 黄色怎么调成土黄色| 丰满人妻一区二区三区视频av| 日日啪夜夜撸| 久久久久久久久久人人人人人人| 午夜免费男女啪啪视频观看| a级片在线免费高清观看视频| a级片在线免费高清观看视频| 亚洲国产精品国产精品| 国产美女午夜福利| 99久国产av精品国产电影| 亚洲色图综合在线观看| 激情五月婷婷亚洲| av天堂久久9| 一级二级三级毛片免费看| 久久精品久久久久久噜噜老黄| 国产成人精品久久久久久| 久久亚洲国产成人精品v| 亚洲欧美成人精品一区二区| 欧美激情国产日韩精品一区| 中文资源天堂在线| 极品教师在线视频| 国产黄色视频一区二区在线观看| 热re99久久国产66热| 成人漫画全彩无遮挡| 国产精品欧美亚洲77777| 久久ye,这里只有精品| 观看美女的网站| 大片免费播放器 马上看| 国产精品一区二区在线不卡| 国产一区二区三区av在线| 中文字幕人妻丝袜制服| 久久久久久久久久人人人人人人| 成年美女黄网站色视频大全免费 | 亚洲精品一二三| 国产精品不卡视频一区二区| 丝袜喷水一区| 日韩 亚洲 欧美在线| 人人妻人人添人人爽欧美一区卜| 国产成人a∨麻豆精品| 一本久久精品| 亚洲四区av| 久久午夜综合久久蜜桃| 女人精品久久久久毛片| 国产亚洲5aaaaa淫片| 亚洲电影在线观看av| 少妇猛男粗大的猛烈进出视频| 亚洲av福利一区| 日韩精品有码人妻一区| 久久99一区二区三区| 欧美日韩一区二区视频在线观看视频在线| kizo精华| 欧美精品高潮呻吟av久久| 久久国产精品男人的天堂亚洲 | 不卡视频在线观看欧美| 亚洲丝袜综合中文字幕| 成人美女网站在线观看视频| 国产伦精品一区二区三区四那| 午夜福利,免费看| 久久久久国产网址| 久久国产精品男人的天堂亚洲 | 久久99蜜桃精品久久| 日韩av在线免费看完整版不卡| 久久99蜜桃精品久久| 国产一区二区在线观看日韩| 亚洲国产av新网站| 少妇被粗大猛烈的视频| 两个人的视频大全免费| 黄色日韩在线| 国产精品99久久久久久久久| 日韩电影二区| 涩涩av久久男人的天堂| 女的被弄到高潮叫床怎么办| 少妇熟女欧美另类| 男人爽女人下面视频在线观看| 精品久久久精品久久久| 男男h啪啪无遮挡| 日本av手机在线免费观看| av一本久久久久| 在线观看三级黄色| 日日撸夜夜添| 99热国产这里只有精品6| 如日韩欧美国产精品一区二区三区 | 一级,二级,三级黄色视频| 中文字幕人妻丝袜制服| 色网站视频免费| 亚洲成人av在线免费| 熟女电影av网| av.在线天堂| 这个男人来自地球电影免费观看 | 久久影院123| √禁漫天堂资源中文www| 嫩草影院新地址| 夜夜骑夜夜射夜夜干| 韩国高清视频一区二区三区| 一级黄片播放器| 国产精品久久久久久精品电影小说| 男人狂女人下面高潮的视频| 亚洲精品一区蜜桃| 成人午夜精彩视频在线观看| 另类亚洲欧美激情| 国产成人aa在线观看| 视频区图区小说| 亚洲精品乱码久久久久久按摩| 精品人妻偷拍中文字幕| 日本黄大片高清| 亚洲欧美日韩东京热| 不卡视频在线观看欧美| 国产欧美亚洲国产| www.av在线官网国产| 黄色配什么色好看| 性高湖久久久久久久久免费观看| 大陆偷拍与自拍| 亚洲成人av在线免费| 男女啪啪激烈高潮av片| av免费观看日本| 亚洲精品日韩av片在线观看| 免费高清在线观看视频在线观看| 久久热精品热| 男女啪啪激烈高潮av片| 日韩av在线免费看完整版不卡| 一级,二级,三级黄色视频| 日韩av不卡免费在线播放| 天天躁夜夜躁狠狠久久av| 国产成人精品久久久久久| 男女免费视频国产| 寂寞人妻少妇视频99o| 日日撸夜夜添| 卡戴珊不雅视频在线播放| 亚洲欧美一区二区三区国产| av在线app专区| 欧美老熟妇乱子伦牲交| 精品人妻熟女毛片av久久网站| 中国美白少妇内射xxxbb| 人人妻人人爽人人添夜夜欢视频 | 99视频精品全部免费 在线| 在线看a的网站| 日日啪夜夜撸| 久久国产乱子免费精品| 亚洲精品视频女| 日本av免费视频播放| 观看av在线不卡| 午夜激情久久久久久久| 色视频在线一区二区三区| 在现免费观看毛片| 亚洲精品乱码久久久v下载方式| 欧美+日韩+精品| 肉色欧美久久久久久久蜜桃| 国产成人精品久久久久久| 永久网站在线| 免费观看av网站的网址| 久久热精品热| 如日韩欧美国产精品一区二区三区 | 国产欧美日韩综合在线一区二区 | 男女边摸边吃奶| 国产成人91sexporn| 亚洲精品国产成人久久av| 啦啦啦视频在线资源免费观看| 最后的刺客免费高清国语| 久久综合国产亚洲精品| 国产成人精品婷婷| 久久久精品94久久精品| av天堂久久9| 人人妻人人添人人爽欧美一区卜| 综合色丁香网| 18禁在线无遮挡免费观看视频| 人妻人人澡人人爽人人| 国产精品熟女久久久久浪| 熟女人妻精品中文字幕| 久久精品国产亚洲网站| 国产精品蜜桃在线观看| 黑人巨大精品欧美一区二区蜜桃 | 3wmmmm亚洲av在线观看| 日日撸夜夜添| 九色成人免费人妻av| 99九九线精品视频在线观看视频| 亚洲精品久久午夜乱码| 欧美日韩国产mv在线观看视频| 9色porny在线观看| 午夜老司机福利剧场| 人妻少妇偷人精品九色| 国产乱人偷精品视频| 蜜桃在线观看..| 免费大片黄手机在线观看| 人人妻人人澡人人爽人人夜夜| 免费av不卡在线播放| 欧美精品一区二区大全| 五月玫瑰六月丁香| 国产成人a∨麻豆精品| 晚上一个人看的免费电影| 欧美高清成人免费视频www| 黑丝袜美女国产一区| 久久av网站| 熟女电影av网| 日产精品乱码卡一卡2卡三| 国产真实伦视频高清在线观看| 亚洲欧洲日产国产| 欧美精品亚洲一区二区| 国产日韩欧美在线精品| 波野结衣二区三区在线| 精华霜和精华液先用哪个| 各种免费的搞黄视频| 亚洲欧洲精品一区二区精品久久久 | 亚洲美女视频黄频| 九色成人免费人妻av| 黄色怎么调成土黄色| 日韩免费高清中文字幕av| 免费黄色在线免费观看| 中文字幕亚洲精品专区| 国产色婷婷99| 国产又色又爽无遮挡免| 色94色欧美一区二区| 男的添女的下面高潮视频| 日本色播在线视频| 黑人高潮一二区| 亚洲情色 制服丝袜| 最后的刺客免费高清国语| 99热6这里只有精品| 国产深夜福利视频在线观看| h视频一区二区三区| 在线精品无人区一区二区三| 黑丝袜美女国产一区| 午夜日本视频在线| 免费播放大片免费观看视频在线观看| 我要看日韩黄色一级片| 美女cb高潮喷水在线观看| 熟妇人妻不卡中文字幕| 看免费成人av毛片| 国产免费一级a男人的天堂| 黄片无遮挡物在线观看| 成人漫画全彩无遮挡| 秋霞在线观看毛片| 少妇的逼水好多| 人人妻人人看人人澡| 国产高清不卡午夜福利| 日产精品乱码卡一卡2卡三| 亚洲一级一片aⅴ在线观看| 韩国高清视频一区二区三区| 一级av片app| 久久久久久久久久成人| 国产免费一级a男人的天堂| 3wmmmm亚洲av在线观看| 99九九线精品视频在线观看视频| 国产一级毛片在线| 久久久国产精品麻豆| 少妇被粗大猛烈的视频| 久久 成人 亚洲| 国产欧美日韩一区二区三区在线 | 亚洲va在线va天堂va国产| 51国产日韩欧美| 久久狼人影院| 三级国产精品欧美在线观看| 国产一区亚洲一区在线观看| av一本久久久久| 日韩欧美精品免费久久| 晚上一个人看的免费电影| a级一级毛片免费在线观看| 十八禁网站网址无遮挡 | 欧美精品国产亚洲| 啦啦啦在线观看免费高清www| 黄色毛片三级朝国网站 | 国产欧美亚洲国产| 国产真实伦视频高清在线观看| 国产欧美另类精品又又久久亚洲欧美| 久久精品久久精品一区二区三区| 国产精品99久久99久久久不卡 | 久久久久精品久久久久真实原创| av一本久久久久| 我要看黄色一级片免费的| 国产欧美日韩一区二区三区在线 | 久久久a久久爽久久v久久| 80岁老熟妇乱子伦牲交| 成人亚洲欧美一区二区av| 国产欧美日韩综合在线一区二区 | 国产日韩一区二区三区精品不卡 | 大陆偷拍与自拍| 爱豆传媒免费全集在线观看| 色婷婷久久久亚洲欧美| 成人美女网站在线观看视频| 日韩欧美 国产精品| 少妇人妻久久综合中文| 亚洲内射少妇av| 一本大道久久a久久精品| 国产深夜福利视频在线观看| 国产成人精品福利久久| 男人舔奶头视频| 多毛熟女@视频| a级毛色黄片| av有码第一页| av免费在线看不卡| 精品一区二区三区视频在线| 国产亚洲av片在线观看秒播厂| 国产精品99久久久久久久久| 国产永久视频网站| 女人精品久久久久毛片| 人人澡人人妻人| 精品99又大又爽又粗少妇毛片| av天堂中文字幕网| 丰满人妻一区二区三区视频av| 国产精品一区www在线观看| 欧美高清成人免费视频www| 亚洲精品aⅴ在线观看| 亚洲国产欧美在线一区| 少妇 在线观看| 欧美+日韩+精品| 精品国产一区二区久久| 男人和女人高潮做爰伦理| 免费看不卡的av| 欧美亚洲 丝袜 人妻 在线| 日韩三级伦理在线观看| 91精品国产国语对白视频| 亚洲国产欧美日韩在线播放 | 成人国产麻豆网| 亚洲色图综合在线观看| 热re99久久精品国产66热6| kizo精华| 亚洲av国产av综合av卡| 免费久久久久久久精品成人欧美视频 | 国产欧美另类精品又又久久亚洲欧美| 国产片特级美女逼逼视频| 亚洲精华国产精华液的使用体验| 亚洲精品日韩av片在线观看| 国产精品一区二区在线不卡| 欧美精品亚洲一区二区| 久热这里只有精品99| 欧美精品人与动牲交sv欧美| 夫妻午夜视频| 男人狂女人下面高潮的视频| 国产一区二区三区综合在线观看 | 美女视频免费永久观看网站| 美女中出高潮动态图| 亚洲欧美一区二区三区黑人 | 亚洲国产欧美在线一区| 春色校园在线视频观看| 国产欧美日韩一区二区三区在线 | 波野结衣二区三区在线| 久久久国产一区二区| 国产精品国产三级国产专区5o| 国产伦在线观看视频一区| 亚洲av国产av综合av卡| 夫妻午夜视频| 国产黄色视频一区二区在线观看| 久热这里只有精品99| 亚洲激情五月婷婷啪啪| 久久ye,这里只有精品| 如何舔出高潮| 亚洲一级一片aⅴ在线观看| 天堂8中文在线网| 精品亚洲成国产av| 亚洲精品久久久久久婷婷小说| 久久6这里有精品| 黄色欧美视频在线观看| 天堂中文最新版在线下载| 亚洲精品乱久久久久久| 国产精品熟女久久久久浪| 国产日韩欧美在线精品| av在线观看视频网站免费| 久久国产亚洲av麻豆专区| 夜夜爽夜夜爽视频| 欧美日韩在线观看h| 王馨瑶露胸无遮挡在线观看| 亚洲第一av免费看| 在线观看av片永久免费下载| 国产精品一区二区在线不卡| 久久精品国产亚洲av天美| 欧美精品人与动牲交sv欧美| av在线播放精品| 日韩欧美一区视频在线观看 | 日日摸夜夜添夜夜爱| 一级av片app| 丁香六月天网| 欧美日韩综合久久久久久| 自拍偷自拍亚洲精品老妇| 亚洲欧洲精品一区二区精品久久久 | 性色avwww在线观看| 久久这里有精品视频免费| 在线观看免费日韩欧美大片 | av国产精品久久久久影院| 亚洲精品一区蜜桃| 三级国产精品欧美在线观看| 一边亲一边摸免费视频| 高清欧美精品videossex| 亚洲精华国产精华液的使用体验| 亚洲综合色惰| 亚洲四区av| 亚洲无线观看免费| 看非洲黑人一级黄片| 午夜老司机福利剧场| 夜夜爽夜夜爽视频| 亚洲av国产av综合av卡| 日本爱情动作片www.在线观看| 99久久人妻综合| 在线看a的网站| 搡老乐熟女国产| 国产精品久久久久久久久免| 老司机影院毛片| 日韩av免费高清视频| 建设人人有责人人尽责人人享有的| 国产精品99久久久久久久久| 婷婷色综合大香蕉| 一区二区三区免费毛片| 青春草视频在线免费观看| av在线观看视频网站免费| 91精品国产国语对白视频| 国产高清三级在线| 国产黄色免费在线视频| 日韩成人av中文字幕在线观看| 亚洲av在线观看美女高潮| 永久免费av网站大全| 国产精品99久久久久久久久| 青青草视频在线视频观看| 国产黄色视频一区二区在线观看| 久久这里有精品视频免费| 久久精品国产自在天天线| 国产在视频线精品| 一个人看视频在线观看www免费| 一本色道久久久久久精品综合| 成人二区视频| 边亲边吃奶的免费视频| 国产免费视频播放在线视频| 国产淫语在线视频| 久久精品国产自在天天线| 欧美精品亚洲一区二区| 国产有黄有色有爽视频| 亚洲色图综合在线观看| 午夜影院在线不卡| 少妇人妻一区二区三区视频| 老司机亚洲免费影院| 性高湖久久久久久久久免费观看| 九色成人免费人妻av| 国产一级毛片在线| 国内少妇人妻偷人精品xxx网站| 丰满少妇做爰视频| 日韩,欧美,国产一区二区三区| 黑人猛操日本美女一级片| 久久久久久久大尺度免费视频| 日韩欧美 国产精品| 丰满少妇做爰视频| 国产在线免费精品| 国产黄片视频在线免费观看| 久久精品熟女亚洲av麻豆精品| 亚洲精品色激情综合| 伊人久久国产一区二区| 精品国产乱码久久久久久小说| 韩国高清视频一区二区三区| 欧美精品人与动牲交sv欧美| 久久97久久精品| 在线观看免费视频网站a站| 视频区图区小说| a级毛色黄片| 成人毛片a级毛片在线播放| 日韩欧美精品免费久久| 免费大片黄手机在线观看| 人人澡人人妻人| 久久久欧美国产精品| 少妇人妻久久综合中文| 国产一级毛片在线| 久久精品国产亚洲av天美| 日本黄色日本黄色录像| 伦理电影免费视频| 亚洲三级黄色毛片| 一级毛片 在线播放| 久久久久久久大尺度免费视频| 国产黄色免费在线视频| 下体分泌物呈黄色| 国产有黄有色有爽视频| 日本av手机在线免费观看| 一级毛片aaaaaa免费看小| 国产极品粉嫩免费观看在线 | 汤姆久久久久久久影院中文字幕| 亚洲精品国产av成人精品| 成年人免费黄色播放视频 | 成年人午夜在线观看视频| 国产白丝娇喘喷水9色精品| 久久影院123| 在线观看一区二区三区激情| 少妇丰满av| 一级av片app| 日韩av在线免费看完整版不卡| 久久国内精品自在自线图片| 国产成人freesex在线| 大片电影免费在线观看免费| 亚洲av不卡在线观看| 国产黄色免费在线视频| 美女cb高潮喷水在线观看| 国产在线免费精品| 免费观看性生交大片5| 人妻一区二区av| 成人漫画全彩无遮挡| 国产精品人妻久久久久久| 一级a做视频免费观看| 精品一区二区三卡| 中文精品一卡2卡3卡4更新| 国产乱人偷精品视频| 亚洲熟女精品中文字幕| 国产美女午夜福利| 中文天堂在线官网| 亚洲精品日本国产第一区| 成年美女黄网站色视频大全免费 | 精品一区在线观看国产| 国产在线男女| 欧美少妇被猛烈插入视频| 精品一品国产午夜福利视频| 中文精品一卡2卡3卡4更新| 国产高清三级在线| 亚洲精品乱久久久久久| 新久久久久国产一级毛片| 26uuu在线亚洲综合色| 久久99热6这里只有精品| 伊人亚洲综合成人网| 大片免费播放器 马上看| 国产成人精品福利久久| 久久青草综合色| 黄色一级大片看看| 在线观看免费高清a一片| 黄色一级大片看看| 黄色毛片三级朝国网站 | 精品久久国产蜜桃| 新久久久久国产一级毛片| 看非洲黑人一级黄片| 91精品国产九色| 如日韩欧美国产精品一区二区三区 | .国产精品久久| 99久久综合免费| 久久久久久久久久久免费av| 啦啦啦中文免费视频观看日本| 伦精品一区二区三区| 亚洲av电影在线观看一区二区三区| 欧美变态另类bdsm刘玥| 亚洲精品色激情综合| 国产一区有黄有色的免费视频| 国产一区二区在线观看av| 精品99又大又爽又粗少妇毛片| 99热国产这里只有精品6| 国产精品蜜桃在线观看| av播播在线观看一区| 99热这里只有是精品在线观看| 亚洲精品自拍成人| 国产成人aa在线观看| 午夜福利网站1000一区二区三区| 美女cb高潮喷水在线观看| 久久人妻熟女aⅴ| 国内精品宾馆在线| 亚洲国产精品一区二区三区在线| 国产精品久久久久成人av| 天堂中文最新版在线下载| 五月天丁香电影| 日韩一本色道免费dvd| 韩国av在线不卡| 免费高清在线观看视频在线观看| 韩国av在线不卡| 欧美亚洲 丝袜 人妻 在线| 日本黄色片子视频| av福利片在线| 精品人妻一区二区三区麻豆| 99热网站在线观看| 在线观看人妻少妇| 99久久中文字幕三级久久日本| 少妇精品久久久久久久| 久久人人爽人人爽人人片va| 少妇被粗大的猛进出69影院 | 免费看日本二区| 久久精品国产自在天天线| 国产欧美日韩一区二区三区在线 | 亚洲av国产av综合av卡| 日韩一区二区三区影片| 日本av免费视频播放| 国产欧美日韩一区二区三区在线 | 国产精品秋霞免费鲁丝片| 久久精品国产鲁丝片午夜精品| 亚洲av.av天堂| 久久久久久久亚洲中文字幕| 各种免费的搞黄视频| 成人综合一区亚洲| 一级毛片电影观看| 亚洲无线观看免费| av国产精品久久久久影院| 大香蕉久久网| 一区二区三区乱码不卡18| 国产精品一区二区性色av| 狂野欧美激情性bbbbbb|