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

    平泉油松建筑材林立地類(lèi)型劃分及立地質(zhì)量評(píng)價(jià)*

    2021-12-01 01:50:40張曉文于青軍羅桂生吳丹妮賈忠奎
    林業(yè)科學(xué) 2021年9期
    關(guān)鍵詞:平泉油松人工林

    張曉文 于青軍 羅桂生 賈 茜 吳丹妮 賈忠奎

    (1.北京林業(yè)大學(xué)林學(xué)院 國(guó)家林業(yè)和草原局油松工程技術(shù)研究中心 北京 100083;2.河北省平泉市國(guó)有黃土梁子林場(chǎng) 平泉 067506)

    當(dāng)前,我國(guó)林業(yè)處于一個(gè)重要的轉(zhuǎn)折期,正以木材生產(chǎn)為主向以生態(tài)建設(shè)為轉(zhuǎn)變,以采伐天然林為主向以采伐人工林為主轉(zhuǎn)變。根據(jù)森林起源,我國(guó)人工林面積比重較低,質(zhì)量低下,而且全面實(shí)施天然林保護(hù)政策后,國(guó)內(nèi)將停止7 000多萬(wàn)hm2國(guó)有天然林商業(yè)性采伐,每年會(huì)減少5 000萬(wàn)m3木材供給。在這樣的背景下,我國(guó)木材資源供給遠(yuǎn)遠(yuǎn)不能滿足日益增長(zhǎng)的木材資源需求,木材進(jìn)口量激增,木材對(duì)外依存度持續(xù)增長(zhǎng),因北,發(fā)展更多的建筑材林對(duì)解決木材供需矛盾、保障國(guó)家生產(chǎn)安全具有重要意義。建筑材林是速生豐產(chǎn)林的一種,不同于一般用材林,其具有定向、速生、豐產(chǎn)、優(yōu)質(zhì)、高效和穩(wěn)定6大特性,以為建筑、家具和裝修等行業(yè)提供原料或大徑材為經(jīng)營(yíng)目標(biāo)。

    我國(guó)建筑材林研究起步較晚,20世紀(jì)60年代開(kāi)始對(duì)速生豐產(chǎn)人工林木材性質(zhì)進(jìn)行過(guò)一些探討,并對(duì)建筑材優(yōu)良種源進(jìn)行了評(píng)估。21世紀(jì)以來(lái),我國(guó)學(xué)者對(duì)建筑材的研究主要集中在木材性質(zhì)方面,如童方平等(2003)探究間伐對(duì)濕地松(Pinuselliottii)木材材性的影響效應(yīng),結(jié)果發(fā)現(xiàn)間伐對(duì)木材干縮率和抗彎強(qiáng)度影響顯著;歐陽(yáng)林男等(2019)分析6個(gè)桉樹(shù)(Eucalyptusspp.)大徑材樹(shù)種木材力學(xué)性質(zhì)差異,結(jié)果表明大花序桉(E.cloeziana)木材綜合性能最優(yōu),其次為細(xì)葉桉(E.tereticornis)、赤桉(E.camaldulensis)、尾葉桉(E.urophylla)、粗皮桉(E.pellita)和巨桉(E.grandis)。

    在國(guó)外,學(xué)者對(duì)各類(lèi)建筑材樹(shù)種木材性質(zhì)進(jìn)行了大量研究。Bhat等(2010)分析化學(xué)改性對(duì)馬占相思(Acaciamangium)木材物理力學(xué)性質(zhì)的影響發(fā)現(xiàn),丙酸酐和琥珀酸酐均能增強(qiáng)木材的耐腐性和機(jī)械性能。Bal等(2013)探究熱處理對(duì)尾巨桉幼齡木和成熟木一些力學(xué)性能的影響表明,在相同處理下,幼齡木的力學(xué)性能下降程度大于成熟木。Kaygin等(2016)研究發(fā)現(xiàn),海拔對(duì)歐洲赤松(Pinussylvestris)木材性質(zhì)影響顯著。Kiaei等(2018)研究得出,鵝耳櫪葉槭(Acercarpinifolium)被抑制木與優(yōu)勢(shì)木在物理性能方面無(wú)顯著差異,但在力學(xué)性能方面差異顯著。

    目前,我國(guó)已對(duì)馬尾松(Pinusmassoniana)(程曦依等, 2017; 丁貴杰等, 2000; 李改云等, 2010)、火炬松(Pinustaeda)(蔣開(kāi)彬等, 2018; 徐有明等, 2007)、杉木(Cunninghamialanceolata)(程琳等, 2017; 尹擁君等, 2001)、濕地松(羅真付等, 2012; 趙奮成等, 2018)、落葉松(Larixspp.)(馬成武等, 2018; 朱玉杰等, 2016)等一系列樹(shù)種的建筑材林培育技術(shù)進(jìn)行了研究,油松(Pinustabulaeformis)作為我國(guó)北方地區(qū)的主要造林樹(shù)種之一,對(duì)其進(jìn)行建筑材林方面的研究鮮見(jiàn)報(bào)道。油松為松科針葉常綠喬木,樹(shù)皮灰褐色,喜光,耐干旱瘠薄,且材質(zhì)堅(jiān)韌細(xì)致,耐腐耐久用,抗壓力強(qiáng),可供建筑、車(chē)船和器具等使用,是我國(guó)特有樹(shù)種(鄭萬(wàn)鈞, 1983)。油松開(kāi)始旺盛生長(zhǎng)較早,而衰退慢,這種生長(zhǎng)特性既可培育中、小徑材,也可培育大徑材(鄭均寶, 1998)。因此,開(kāi)展油松建筑材林高效培育技術(shù)研究具有重大意義。

    立地質(zhì)量評(píng)價(jià)即對(duì)某一立地單元在采用某種利用方式時(shí)的最大自然生產(chǎn)潛力進(jìn)行評(píng)價(jià)(詹昭寧, 1989),其是研究、掌握森林生長(zhǎng)環(huán)境對(duì)森林生產(chǎn)力影響的一種重要手段,對(duì)林地經(jīng)營(yíng)和決策至關(guān)重要(Louwetal., 2002),是實(shí)現(xiàn)科學(xué)造林和經(jīng)營(yíng)森林的關(guān)鍵。立地質(zhì)量評(píng)價(jià)方法通常可分為生物因子法和地理因子法兩大類(lèi),其中生物因子法中的地位級(jí)法和立地指數(shù)法最常用,立地指數(shù)法多用于評(píng)價(jià)人工林(吳恒等, 2015)。我國(guó)學(xué)者對(duì)森林立地質(zhì)量評(píng)價(jià)開(kāi)展過(guò)一些研究,部分學(xué)者采用逐步回歸分析(張超等, 2015)、人工神經(jīng)網(wǎng)絡(luò)模型(鞏垠熙等, 2013)、建立數(shù)學(xué)模型(朱光玉等, 2010)和遙感(李明澤等, 2017; 王永昌等, 2007)等方法對(duì)我國(guó)主要森林立地質(zhì)量進(jìn)行了評(píng)價(jià)。但隨著國(guó)家對(duì)林業(yè)需求的轉(zhuǎn)變,森林立地質(zhì)量評(píng)價(jià)要求更加精細(xì),需要進(jìn)一步量化到二級(jí)林種上,建筑材林作為人工用材林的二級(jí)林種,目前對(duì)其的研究多集中于高效培育技術(shù)或木材基本物理力學(xué)性質(zhì)方面(徐有明等, 2007; 周海賓等, 2012),未將立地質(zhì)量評(píng)價(jià)與木材材性結(jié)合起來(lái)。鑒于此,本研究對(duì)河北省平泉市油松人工林進(jìn)行立地類(lèi)型劃分和立地質(zhì)量評(píng)價(jià),并將立地質(zhì)量分級(jí),以期篩選出適合培育油松建筑材林的林地。

    1 研究區(qū)概況

    河北平泉(118°21′—119°15′E, 40°24′—40°40′N(xiāo))是河北、遼寧、內(nèi)蒙古三省(自治區(qū))交界地帶,位于燕山山脈末端,平均海拔660 m。該市地處我國(guó)半濕潤(rùn)地區(qū),年均降水量540 mm,屬溫帶季風(fēng)性氣候,四季分明,春季風(fēng)大干旱、回暖較快,夏季炎熱多雨,冬季寒冷干燥,年均氣溫7.3 ℃。7月平均氣溫22~32 ℃,降水占全年降水量的68%~74%; 年日照時(shí)數(shù)2 000~2 900 h,年日照率65%,全年無(wú)霜期120~130天。土壤類(lèi)型以褐土為主,森林類(lèi)型為天然次生林,主要樹(shù)種有油松、落葉松、刺槐(Robiniapseudoacacia)等,灌木層主要有胡枝子(Lespedezabicolor)、錦雞兒(Caraganasinica)、繡線菊(Spiraeasalicifolia)、山杏(Armeniacasibirica)等,常見(jiàn)草本有菊葉委陵菜(Potentillatanacetifolia)、細(xì)葉薹草(Carexduriusculasubsp.stenophylloides)和玉竹(Polygonatumodoratum)等。

    2 研究方法

    2.1 樣地布設(shè)與調(diào)查

    在平泉市馬家南溝、婁南溝、蘇達(dá)營(yíng)子等地布設(shè)92塊樣地,使其能夠代表平泉油松人工林的生長(zhǎng)狀況。為了保證調(diào)查數(shù)據(jù)的代表性和精確性,并綜合考慮密度、林齡等因素,調(diào)查樣地包含幼齡林、中齡林、近熟林和成熟林,且郁閉度均在0.5以上。同時(shí)要求布設(shè)樣地的造林技術(shù)和經(jīng)營(yíng)措施基本一致,未發(fā)生過(guò)蟲(chóng)害、火災(zāi)等自然災(zāi)害。

    在每塊樣地內(nèi)進(jìn)行每木檢尺,實(shí)測(cè)樹(shù)高、胸徑、活枝下高、東西冠幅、南北冠幅,并記錄海拔、坡度、坡向、土層厚度、腐殖質(zhì)厚度、枯落物厚度和土壤類(lèi)型等因子。采用調(diào)查林場(chǎng)小班與生長(zhǎng)錐相結(jié)合的方法確定樣地林齡。每塊樣地內(nèi)選取胸徑最大且最高的5株樹(shù)木取樹(shù)高平均值(段劼等, 2009)。測(cè)量過(guò)程中,每株優(yōu)勢(shì)木最少測(cè)高2次,2次誤差不超過(guò)0.1 m時(shí)取其平均值作為優(yōu)勢(shì)木高。

    2.2 樣地?cái)?shù)據(jù)整理

    將92塊樣地調(diào)查數(shù)據(jù)匯總成表,按齡階(段劼等, 2009)分別統(tǒng)計(jì)出平均林齡和優(yōu)勢(shì)木平均高,并計(jì)算各齡階的樹(shù)高標(biāo)準(zhǔn)差:

    (1)

    式中:Si為第i齡階樹(shù)高標(biāo)準(zhǔn)差;Hij為第i齡階中第j株優(yōu)勢(shì)木高(j=1,2,3,…,ni);ni為第i齡階中優(yōu)勢(shì)木株數(shù)。

    各齡階以平均優(yōu)勢(shì)木高為準(zhǔn),以3倍樹(shù)高標(biāo)準(zhǔn)差為范圍,剔除各齡階內(nèi)數(shù)據(jù)異常的樣地或優(yōu)勢(shì)木,重新統(tǒng)計(jì)樣地?cái)?shù)、優(yōu)勢(shì)木平均高。經(jīng)檢驗(yàn),在92塊樣地中選取的優(yōu)勢(shì)木均符合條件,滿足立地指數(shù)表編制要求。編表數(shù)據(jù)按齡階整理后,優(yōu)勢(shì)木平均高6~14 m,具體結(jié)果見(jiàn)表1。

    2.3 立地類(lèi)型劃分方法

    2.3.1 立地因子分級(jí) 采用主導(dǎo)因子分析法(馬明東, 2006)篩選主導(dǎo)因子,其中類(lèi)目范圍劃分以其對(duì)油松生長(zhǎng)影響是否顯著、實(shí)踐中是否便于調(diào)查和應(yīng)用為依據(jù)(唐誠(chéng)等, 2018)。實(shí)地調(diào)查時(shí)發(fā)現(xiàn),平泉油松人工林主要分布在海拔720~820 m,且土壤類(lèi)型以褐土為主,故排除海拔和土壤類(lèi)型對(duì)林地的影響;油松人工林分布坡位差異不顯著,故排除坡位對(duì)林地的影響。將坡向、坡度、土層厚度、腐殖質(zhì)厚度、枯落物厚度等立地因子進(jìn)行分級(jí)并賦值(表2)。

    2.3.2 主導(dǎo)因子篩選及立地類(lèi)型劃分 采用Forstat軟件,將表2中立地因子與優(yōu)勢(shì)木平均高進(jìn)行相關(guān)性分析,結(jié)果如表3所示。

    根據(jù)Pr>F判斷因子顯著性,當(dāng)Pr>F<0.01時(shí)極顯著,0.01≤Pr>F<0.05時(shí)顯著,Pr>F≥0.05時(shí)不顯著。由表3可知,坡向和土層厚度對(duì)油松生長(zhǎng)影響顯著,故選擇坡向和土層厚度作為主導(dǎo)因子,據(jù)此將平泉油松人工林立地類(lèi)型劃分為陽(yáng)坡薄土、陽(yáng)坡中土、陽(yáng)坡厚土、陰坡薄土、陰坡中土和陰坡厚土,土層厚度劃分依據(jù)見(jiàn)表2。

    表3 立地因子與優(yōu)勢(shì)木平均高的相關(guān)性Tab.3 The correlation between site factors and dominant tree mean height

    2.4 立地質(zhì)量評(píng)價(jià)方法

    2.4.1 基準(zhǔn)年齡的確定 基準(zhǔn)年齡(A0)指林分生長(zhǎng)趨于穩(wěn)定時(shí)的年齡,是立地質(zhì)量評(píng)價(jià)的重要指標(biāo)(陳永富, 2010)。確定基準(zhǔn)年齡應(yīng)考慮以下幾方面: 1) 基準(zhǔn)年齡時(shí),林分樹(shù)高生長(zhǎng)已趨于穩(wěn)定,且能反映林分生長(zhǎng)的立地差異(吳中能等, 2010); 2) 確定基準(zhǔn)年齡前應(yīng)先對(duì)樹(shù)種的生長(zhǎng)規(guī)律進(jìn)行分析,還要考慮樹(shù)種的培育目的。段劼等(2009)研究北京地區(qū)油松基準(zhǔn)年齡時(shí),將其定為20年。本研究在確定基準(zhǔn)年齡時(shí),首先采用隨機(jī)抽樣方法,在不同立地類(lèi)型林分中各選取1株解析木,共計(jì)42株,分析樹(shù)高0~35年的生長(zhǎng)過(guò)程(圖1),然后對(duì)全部解析木求樹(shù)高連年生長(zhǎng)量和平均生長(zhǎng)量的算數(shù)平均值,再繪制其隨年齡的變化曲線(圖2)。

    由圖1可知, 20年后油松樹(shù)高生長(zhǎng)逐漸減緩,趨于穩(wěn)定;由圖2可知,油松樹(shù)高平均生長(zhǎng)量在20年左右達(dá)到峰值,其后逐漸減緩,趨于穩(wěn)定,而連年生長(zhǎng)量在10、20年達(dá)到峰值,其后下降。所以,本研究將平泉油松建筑材林的基準(zhǔn)年齡定為20年。

    圖1 各立地類(lèi)型樹(shù)高生長(zhǎng)曲線Fig.1 Tree height growth curves of various site types

    圖2 解析木樹(shù)高生長(zhǎng)曲線Fig.2 Analytic tree height growth curves

    2.4.2 指數(shù)級(jí)距的確定 指數(shù)級(jí)距指在基準(zhǔn)年齡時(shí)優(yōu)勢(shì)木平均高變幅與指數(shù)級(jí)個(gè)數(shù)的比值。同一森林植被地帶同一樹(shù)種的指數(shù)級(jí)距應(yīng)基本一致,一般來(lái)說(shuō),集約經(jīng)營(yíng)水平較高、以速生豐產(chǎn)為主的地區(qū),級(jí)距可以適當(dāng)小些,管理水平較低、立地條件較差、以防護(hù)或其他效能為主的地區(qū),級(jí)距可以大些(鄭聰慧等, 2013)。本研究用以下2種方法確定指數(shù)級(jí)距。

    1) 重點(diǎn)分析油松建筑材林20年(基準(zhǔn)年齡)時(shí)林分優(yōu)勢(shì)木平均高的分布情況,各樣地優(yōu)勢(shì)木平均高變幅為3.2~9.6 m,絕對(duì)變幅ΔH為6.4 m,一般指數(shù)級(jí)個(gè)數(shù)不超過(guò)10個(gè),采用式(2)計(jì)算指數(shù)級(jí)距:

    S=ΔH/k。

    (2)

    根據(jù)計(jì)算結(jié)果,確定指數(shù)級(jí)距(S)為1 m,指數(shù)級(jí)個(gè)數(shù)(k)為7,分別為4、5、6、7、8、9、10 m。

    2) 根據(jù)唐瑞德(1986)給出的立地指數(shù)級(jí)距計(jì)算公式S=1.5+1/2lnΔA/A0,(ΔA=5,為立地指數(shù)表的齡階;A0=20),求得指數(shù)級(jí)距S=0.8 m,依據(jù)指數(shù)級(jí)距最小為1 m的原則,確定指數(shù)級(jí)距為1 m。

    以上2種方法確定的指數(shù)級(jí)距均為1 m,因此本研究將平泉油松建筑材林的指數(shù)級(jí)距確定為1 m。

    2.4.3 導(dǎo)向曲線的擬合 編制立地指數(shù)表時(shí),導(dǎo)向曲線的選擇至關(guān)重要,其與立地指數(shù)表的精度密切相關(guān)。選擇導(dǎo)向曲線時(shí)既要考慮樹(shù)高生長(zhǎng)規(guī)律,又要保證擬合效果。本研究根據(jù)全部樣地優(yōu)勢(shì)木解析資料,利用各公式對(duì)齡階優(yōu)勢(shì)木平均高數(shù)據(jù)進(jìn)行曲線擬合,常用的導(dǎo)向曲線擬合公式及其參數(shù)見(jiàn)表4。

    表4 導(dǎo)向曲線擬合公式及其參數(shù)①Tab.4 Guiding curve fitting formula and its parameters

    由表4可知,線性、對(duì)數(shù)曲線、拋物線、冪函數(shù)、等比級(jí)數(shù)曲線和指數(shù)曲線的R2均超過(guò)0.8,其中R2最大的是冪函數(shù)曲線(0.855),這是因?yàn)閮绾瘮?shù)只有2個(gè)參數(shù),冪函數(shù)圖像在林齡為0時(shí)樹(shù)高為0,且冪函數(shù)曲線趨勢(shì)更符合油松生長(zhǎng)規(guī)律。所以本研究選擇冪函數(shù)曲線,導(dǎo)向曲線方程為:

    H=0.524A0.815。

    (3)

    2.4.4 立地指數(shù)表的編制 本研究采用相對(duì)優(yōu)勢(shì)高法編制立地指數(shù)表,具體步驟如下:

    1) 計(jì)算各齡階樹(shù)高理論值Hi,將各齡階階值分別代入導(dǎo)向曲線,即可得到樹(shù)高理論值;

    2) 計(jì)算基準(zhǔn)年齡時(shí)的樹(shù)高理論值H0;

    3) 計(jì)算各齡階調(diào)整系數(shù)Ki:

    Ki=Hi/H0×100%;

    (4)

    4) 計(jì)算各齡階各指數(shù)級(jí)樹(shù)高值Hik:

    Hik=Ki×H0k。

    (5)

    式中:H0k為標(biāo)準(zhǔn)年齡時(shí)第k指數(shù)級(jí)的樹(shù)高值。

    由式(3)和(4)計(jì)算得出各齡階樹(shù)高理論值及調(diào)整系數(shù)(表5),由式(5)和表5整理得到油松人工林立地指數(shù)表(表6)。

    表5 各齡階調(diào)整系數(shù)及樹(shù)高理論值Tab.5 Adjust coeffident and tree height in theory of different age classes

    表6 平泉油松人工林立地指數(shù)表(A0=20年)Tab.6 Site index table of P. tabulaeformis plantations in Pingquan (A0=20 a) m

    2.4.5 立地指數(shù)表的檢驗(yàn) 立地指數(shù)表編制完成必須進(jìn)行檢驗(yàn),在保證精度和適用性后才能應(yīng)用于實(shí)際生產(chǎn)中。本研究采用卡方檢驗(yàn),首先利用現(xiàn)有42株優(yōu)勢(shì)木解析木樹(shù)高和年齡,在立地指數(shù)表(表6)中查到相對(duì)應(yīng)的立地指數(shù),然后將優(yōu)勢(shì)木解析木林齡分別代入導(dǎo)向曲線得到各優(yōu)勢(shì)木解析木樹(shù)高理論值Ht,與優(yōu)勢(shì)木解析木樹(shù)高實(shí)際值Hp進(jìn)行χ2檢驗(yàn),公式為:

    (6)

    2.5 油松木材物理力學(xué)性質(zhì)測(cè)定

    立地質(zhì)量對(duì)林木生長(zhǎng)影響很大,且其影響將伴隨林木一生,因此,在人工林培育中立地選擇十分重要。實(shí)際應(yīng)用中,木材體積干縮率、抗彎強(qiáng)度、抗彎彈性模量和順紋抗壓強(qiáng)度是優(yōu)良建筑材選擇的重要參數(shù)(羅彬等, 2008; 周海賓等, 2012)。為了探究不同立地指數(shù)下木材物理力學(xué)性質(zhì)是否存在顯著差異,本研究選取立地條件為差、中、優(yōu)3種類(lèi)型樣地(立地條件劃分標(biāo)準(zhǔn)見(jiàn)表7),每種類(lèi)型樣地3塊,在樣地中伐取1株優(yōu)勢(shì)木作為樣木,對(duì)其進(jìn)行物理力學(xué)性質(zhì)測(cè)試。

    表7 立地條件劃分標(biāo)準(zhǔn)Tab.7 Site condition division standard

    3 結(jié)果與分析

    3.1 油松人工林立地類(lèi)型劃分

    根據(jù)篩選的主導(dǎo)因子,將平泉油松人工林立地劃分為陽(yáng)坡厚土、陽(yáng)坡中土、陽(yáng)坡薄土、陰坡厚土、陰坡中土和陰坡薄土6種類(lèi)型?;贏rcGIS地理信息系統(tǒng),結(jié)合立地類(lèi)型劃分結(jié)果,繪制平泉油松人工林立地類(lèi)型分布圖(圖3)。從圖3可以看出,平泉油松人工林主要分布在薄土層和中土層,厚土層分布較少,其中薄土層面積14 239.15 hm2,占油松人工林總面積的27.9%,中土層面積33 355.98 hm2,占油松人工林總面積的65.5%,厚土層面積3 304.67 hm2,占油松人工林總面積的6.6%。由圖3還可以看出,平泉油松人工林陰坡分布較多,陽(yáng)坡分布較少,其中陰坡面積38 976.01 hm2,占油松人工林總面積的76.6%,陽(yáng)坡面積11 920.04 hm2,占油松人工林總面積的23.4%。

    圖3 平泉油松人工林立地類(lèi)型分布Fig.3 Site type distribution of P. tabulaeformis plantations in Pingquan

    3.2 油松人工林立地質(zhì)量評(píng)價(jià)

    3.2.1 擬合優(yōu)勢(shì)木平均高和林分平均高曲線 為了將平泉油松人工林立地指數(shù)表與油松二類(lèi)調(diào)查數(shù)據(jù)相結(jié)合,需要擬合平泉油松人工林優(yōu)勢(shì)木平均高和林分平均高曲線,以探究油松優(yōu)勢(shì)木平均高與林分平均高的關(guān)系。本研究基于平泉油松人工林優(yōu)勢(shì)木平均高和林分平均高調(diào)查數(shù)據(jù),以林分平均高為自變量(x)、優(yōu)勢(shì)木平均高為因變量(y),采用直線型、二次曲線型、三次曲線型、對(duì)數(shù)曲線型、冪函數(shù)型、S型曲線、復(fù)合函數(shù)型、指數(shù)曲線型和邏輯斯特曲線型9種模型擬合平泉油松人工林優(yōu)勢(shì)木平均高與林分平均高的關(guān)系。各曲線模型的相關(guān)系數(shù)和參數(shù)見(jiàn)表8。

    表8 優(yōu)勢(shì)木平均高和林分平均高曲線擬合參數(shù)Tab.8 Fitting parameters of dominant tree mean height and stand mean height

    通過(guò)R2比較,選擇二次曲線型,得到平泉油松人工林優(yōu)勢(shì)木平均高-林分平均高擬合曲線,公式為:

    y=0.512 26+1.563 74x-0.034 67x2。

    (7)

    3.2.2 立地指數(shù)等級(jí)分布圖的繪制 依據(jù)式(7)和表6,結(jié)合平泉油松二類(lèi)調(diào)查數(shù)據(jù),利用ArcGIS軟件,繪制立地指數(shù)等級(jí)分布圖(圖4)。統(tǒng)計(jì)分析發(fā)現(xiàn),立地指數(shù)為9、10的地塊占地面積14 277.665 hm2,占油松人工林總面積的28.1%,其中黃土梁子鎮(zhèn)和柳溪滿族鄉(xiāng)占地面積13 645.96 hm2,說(shuō)明這兩個(gè)地區(qū)立地質(zhì)量較好。立地指數(shù)為6、7、8的地塊占地面積33 648.51 hm2,占油松人工林總面積的66.1%,而立地指數(shù)為4、5的僅有2 973.63 hm2,僅占油松人工林總面積的5.8%。這說(shuō)明河北平泉地區(qū)立地質(zhì)量狀況總體較好,適合培育油松建筑材林。

    3.3 油松木材物理力學(xué)性質(zhì)分析

    分別計(jì)算各立地質(zhì)量樣地各木材物理力學(xué)指標(biāo)測(cè)定值的最小值、最大值、均值和標(biāo)準(zhǔn)差,結(jié)果見(jiàn)表9。

    圖4 平泉油松立地指數(shù)等級(jí)分布Fig.4 Site index distribution of P. tabulaeformis in Pingquan

    表9 各立地質(zhì)量樣地木材物理力學(xué)性質(zhì)①Tab.9 Physical and mechanical properties of wood in various site quality samples

    體積干縮率反映木材在干燥過(guò)程中的收縮程度,與木材中射線細(xì)胞數(shù)量和分布有關(guān)(Mohammadetal., 2012)。由表9可知,各樣地體積干縮率均值為0.097%~0.125%,極差為0.028%,立地質(zhì)量?jī)?yōu)的樣地體積干縮率最高,立地質(zhì)量中的樣地體積干縮率次之,立地質(zhì)量差的樣地體積干縮率最低,且兩兩之間在0.05水平上差異顯著,說(shuō)明木材干縮性能隨立地質(zhì)量提升而增強(qiáng)。

    抗彎彈性模量和抗彎強(qiáng)度代表木材的抗彎性能(Ruelleetal., 2007)。由表9可知,各樣地抗彎彈性模量均值為10.27~13.79 GPa,極差為3.52 GPa,立地質(zhì)量?jī)?yōu)的樣地抗彎彈性模量最高,立地質(zhì)量中的樣地抗彎彈性模量次之,立地質(zhì)量差的樣地抗彎彈性模量最低,立地質(zhì)量差、立地質(zhì)量中的樣地與立地質(zhì)量?jī)?yōu)的樣地在0.05水平上均差異顯著,說(shuō)明抗彎彈性模量隨著立地質(zhì)量提升而提高,在一定范圍內(nèi)差異顯著。各樣地抗彎強(qiáng)度均值為98.61~119.94 MPa,極差為21.33 MPa,立地質(zhì)量?jī)?yōu)的樣地抗彎強(qiáng)度最高,立地質(zhì)量中的樣地抗彎強(qiáng)度次之,立地質(zhì)量差的樣地抗彎強(qiáng)度最低,立地質(zhì)量差和立地質(zhì)量中的樣地均與立地質(zhì)量?jī)?yōu)的樣地在0.05水平上差異顯著,說(shuō)明木材抗彎強(qiáng)度隨立地質(zhì)量提升而增強(qiáng)。

    各樣地順紋抗壓強(qiáng)度均值為56.623~78.612 MPa,極差為21.989 MPa,立地質(zhì)量?jī)?yōu)的樣地順紋抗壓強(qiáng)度最高,立地質(zhì)量中的樣地順紋抗壓強(qiáng)度次之,立地質(zhì)量差的樣地順紋抗壓強(qiáng)度最低,立地質(zhì)量差和立地質(zhì)量中的樣地均與立地質(zhì)量?jī)?yōu)的樣地在0.05水平上差異顯著,而立地質(zhì)量差與立地質(zhì)量中的樣地之間差異不顯著。

    綜上所述,木材物理力學(xué)性質(zhì)各項(xiàng)指標(biāo)均隨著立地質(zhì)量提升而提高,說(shuō)明立地質(zhì)量好的林地(立地指數(shù)為9、10)更適合培育大、中徑材,而立地質(zhì)量中的林地(立地指數(shù)為6、7、8)適合培育中、小徑材,立地質(zhì)量差的林地(立地指數(shù)為4、5)在木材材性方面不夠穩(wěn)定,不適合培育建筑材,可用作生態(tài)公益林。

    3.4 油松建筑材培育分布圖的繪制

    基于ArcGIS地理信息系統(tǒng),從圖4中篩選出立地指數(shù)9、10的林地為適宜培育大、中徑材的林地, 立地指數(shù)6、7、8的林地為適宜培育中、小徑材的林地, 立地指數(shù)4、5的林地為不適宜培育建筑材的林地,繪制出油松建筑材培育分布區(qū)(圖5)。經(jīng)篩選統(tǒng)計(jì),可用來(lái)培育大、中徑材的地塊占地面積13 645.96 hm2,占油松人工林總面積的28.1%;可用來(lái)培育中、小徑材的地塊占地面積33 648.51 hm2,占油松人工林總面積的66.1%;生態(tài)公益林占地面積2 973.63 hm2,僅占油松人工林總面積的5.8%。這說(shuō)明,河北平泉地區(qū)適合用來(lái)培育建筑材。可用于培育大、中徑材的地塊中,以黃土梁子鎮(zhèn)分布最多,共計(jì)3 697.61 hm2,其次是柳溪滿族鄉(xiāng)占地面積1 897.39 hm2。黃土梁子鎮(zhèn)和柳溪滿族鄉(xiāng)所占比例明顯高于其他地區(qū),兩地所占比例為油松人工林總面積的42.9%??捎糜谂嘤?、小徑材的地塊中,黃土梁子鎮(zhèn)占地面積5 570.86 hm2,其次是柳溪滿族鄉(xiāng)占地面積4 770.64 hm2,黃土梁子鎮(zhèn)和柳溪滿族鄉(xiāng)所占比例為油松人工林總面積的30.7%。因此若要在河北平泉地區(qū)培育油松建筑材林,應(yīng)優(yōu)先考慮黃土梁子鎮(zhèn)和柳溪滿族鄉(xiāng)。

    圖5 平泉油松建筑材林適宜培育分布區(qū)Fig.5 P. tabulaeformis construction timber cultivation area in Pingquan

    4 討論

    本研究對(duì)油松建筑材林進(jìn)行立地類(lèi)型劃分和立地質(zhì)量評(píng)價(jià),基于ArcGIS地理信息系統(tǒng)繪制平泉油松人工林立地類(lèi)型分布圖、立地指數(shù)等級(jí)分布圖以及建筑材林適宜培育區(qū)圖,可以直觀反映出油松人工林立地類(lèi)型分布、立地指數(shù)分布以及建筑材林適宜培育區(qū)分布情況,能為林業(yè)實(shí)踐提供指導(dǎo)。

    采用主導(dǎo)因子分析法篩選出影響油松人工林優(yōu)勢(shì)木平均高的主導(dǎo)因子為坡向和土層厚度,據(jù)此將油松人工林立地類(lèi)型劃分為6種。將坡向劃分為陽(yáng)坡和陰坡,并沒(méi)有細(xì)化至半陰坡和半陽(yáng)坡,以后研究中應(yīng)加入該部分探索。本研究中,海拔作為立地類(lèi)型劃分的關(guān)鍵因子并沒(méi)有被篩選出來(lái),以后研究中應(yīng)注重海拔對(duì)立地質(zhì)量影響方面的分析。

    段劼等(2009)選擇雙曲線為導(dǎo)向曲線,采用相對(duì)優(yōu)勢(shì)高法對(duì)北京低山地區(qū)油松進(jìn)行立地質(zhì)量評(píng)價(jià),并編制立地指數(shù)表。本研究同樣采用相對(duì)優(yōu)勢(shì)高法,在10種曲線中,經(jīng)對(duì)比分析選擇冪函數(shù)式為導(dǎo)向曲線,并進(jìn)行立地指數(shù)表編制,經(jīng)卡方檢驗(yàn)后合格。前人在進(jìn)行立地質(zhì)量評(píng)價(jià)時(shí),多止步于立地指數(shù)表的編制,本研究利用ArcGIS地理信息系統(tǒng),將立地指數(shù)表與油松小班調(diào)查數(shù)據(jù)結(jié)合,得到包含立地指數(shù)信息的油松小班調(diào)查屬性表,由該表能夠快捷、方便查出某一小班的立地類(lèi)型和立地質(zhì)量。這一結(jié)果在實(shí)際應(yīng)用中能夠大幅度提高工作效率,具有很強(qiáng)的現(xiàn)實(shí)意義。

    木材化學(xué)性質(zhì)也是影響建筑材性能的重要指標(biāo),楊昇等(2019)研究熱處理對(duì)進(jìn)口輻射松(Pinusradiata)木材化學(xué)性質(zhì)的影響,結(jié)果發(fā)現(xiàn), 熱處理會(huì)導(dǎo)致輻射松木材細(xì)胞壁中木質(zhì)素分子結(jié)構(gòu)發(fā)生變化,且隨著熱處理溫度升高,木質(zhì)素結(jié)構(gòu)改變程度逐漸增加。秦理哲等(2019)研究季銨銅防腐劑對(duì)馬尾松木材化學(xué)性質(zhì)的影響,結(jié)果表明,防腐處理會(huì)引起馬尾松木材酸堿性質(zhì)、化學(xué)元素組成和官能團(tuán)強(qiáng)度的變化,影響其化學(xué)性質(zhì)。本研究?jī)H從物理力學(xué)性質(zhì)方面對(duì)油松木材進(jìn)行研究,沒(méi)有考慮化學(xué)性質(zhì)以及解剖性能,以后研究中可以加入化學(xué)性質(zhì)和解剖性能的探索。

    5 結(jié)論

    1)河北平泉油松人工林主要分布在薄土層和中土層,厚土層分布較少,其中薄土層面積占油松人工林總面積的27.9%,中土層面積占65.5%,厚土層面積僅占6.6%;陰坡分布較多,陽(yáng)坡分布較少,其中陰坡面積占油松人工林總面積的76.6%,陽(yáng)坡占23.4%。由此可知,平泉地區(qū)的立地特點(diǎn)適合油松林生長(zhǎng)。

    2)平泉地區(qū)油松人工林立地指數(shù)為9、10的地塊面積占油松人工林面積的28.1%,立地指數(shù)為6、7、8的地塊面積占66.1%,立地指數(shù)為4、5的地塊面積占5.8%,說(shuō)明河北平泉地區(qū)立地質(zhì)量狀況較好,適合培育油松建筑材林。

    3)木材物理力學(xué)性質(zhì)隨立地質(zhì)量提升而提高,說(shuō)明立地質(zhì)量好的林地更適合培育優(yōu)質(zhì)的大、中徑級(jí)建筑材林,立地質(zhì)量中的林地適合培育中、小徑級(jí)建筑材林,立地質(zhì)量差的林地不適合培育建筑材林。黃土梁子鎮(zhèn)和柳溪滿族鄉(xiāng)可作為平泉地區(qū)培育油松建筑材林的主要地點(diǎn)。

    猜你喜歡
    平泉油松人工林
    平泉市萬(wàn)壽菊帶農(nóng)增收致富
    平泉市人大常委會(huì)多措并舉做實(shí)聯(lián)動(dòng)監(jiān)督“后半篇文章”
    公民與法治(2022年8期)2022-02-04 07:13:05
    油松容器育苗技術(shù)
    平泉市農(nóng)村信用合作聯(lián)社
    北方油松育苗技術(shù)
    我國(guó)人工林生態(tài)輪伐期探討
    人工林生產(chǎn)力提升對(duì)策探析
    山東林業(yè)科技(2016年5期)2016-07-05 00:43:04
    平泉發(fā)改局:打好立縣富民“三張牌”
    公民與法治(2016年7期)2016-05-17 04:11:26
    濕地松人工林生長(zhǎng)規(guī)律
    x7x7x7水蜜桃| 精品乱码久久久久久99久播| 国产精品 国内视频| 两人在一起打扑克的视频| 欧美最黄视频在线播放免费| 国产精品99久久99久久久不卡| 999久久久国产精品视频| 神马国产精品三级电影在线观看 | 亚洲 国产 在线| 日本一区二区免费在线视频| 成在线人永久免费视频| 欧美日韩乱码在线| 在线观看免费日韩欧美大片| 免费在线观看黄色视频的| 一卡2卡三卡四卡精品乱码亚洲| 久久人妻av系列| 搡老熟女国产l中国老女人| 色哟哟哟哟哟哟| 女同久久另类99精品国产91| 在线视频色国产色| 久久亚洲精品不卡| 宅男免费午夜| 俄罗斯特黄特色一大片| 欧美成人性av电影在线观看| 国产精品久久久久久人妻精品电影| 成人手机av| 亚洲成a人片在线一区二区| 午夜久久久久精精品| 一卡2卡三卡四卡精品乱码亚洲| 一本精品99久久精品77| 91老司机精品| 9191精品国产免费久久| 国产精品一区二区精品视频观看| 欧美三级亚洲精品| 婷婷精品国产亚洲av在线| avwww免费| 久久精品影院6| 香蕉丝袜av| 18禁裸乳无遮挡免费网站照片 | 亚洲中文日韩欧美视频| 亚洲av五月六月丁香网| 久久天堂一区二区三区四区| 国产av又大| 国产精品一区二区免费欧美| 国产爱豆传媒在线观看 | 午夜精品久久久久久毛片777| 中文在线观看免费www的网站 | 国产欧美日韩一区二区三| 免费人成视频x8x8入口观看| 高清在线国产一区| 国产精品1区2区在线观看.| 亚洲成人精品中文字幕电影| 色播亚洲综合网| 精品无人区乱码1区二区| 亚洲在线自拍视频| 亚洲成人久久性| 十八禁网站免费在线| 高清毛片免费观看视频网站| 啦啦啦韩国在线观看视频| 国产精品亚洲一级av第二区| 美女高潮喷水抽搐中文字幕| 真人一进一出gif抽搐免费| 欧美久久黑人一区二区| 国产爱豆传媒在线观看 | 国产真实乱freesex| 免费观看精品视频网站| 日本五十路高清| 俄罗斯特黄特色一大片| 狂野欧美激情性xxxx| 欧美乱色亚洲激情| 美女大奶头视频| 国产伦人伦偷精品视频| 黄色丝袜av网址大全| 久热爱精品视频在线9| 一区二区日韩欧美中文字幕| 午夜福利在线观看吧| 久久久久精品国产欧美久久久| 国语自产精品视频在线第100页| 男人舔女人下体高潮全视频| 久久午夜亚洲精品久久| 亚洲国产欧美网| av中文乱码字幕在线| 国产精品精品国产色婷婷| 桃色一区二区三区在线观看| 白带黄色成豆腐渣| 国产蜜桃级精品一区二区三区| 亚洲精华国产精华精| 久久伊人香网站| 亚洲专区国产一区二区| 国产三级黄色录像| 99久久综合精品五月天人人| 午夜两性在线视频| 欧美激情极品国产一区二区三区| 成人国语在线视频| 久久久久国产精品人妻aⅴ院| 午夜福利欧美成人| 桃红色精品国产亚洲av| 午夜日韩欧美国产| avwww免费| 这个男人来自地球电影免费观看| 久久国产亚洲av麻豆专区| 自线自在国产av| 成在线人永久免费视频| www.熟女人妻精品国产| 国产99白浆流出| 国产av不卡久久| 精品熟女少妇八av免费久了| 熟妇人妻久久中文字幕3abv| 免费在线观看完整版高清| 久9热在线精品视频| 亚洲中文日韩欧美视频| 国产成人系列免费观看| 亚洲黑人精品在线| 啦啦啦观看免费观看视频高清| 18禁观看日本| 精品国产乱码久久久久久男人| 国产v大片淫在线免费观看| 亚洲色图av天堂| 欧美午夜高清在线| 热99re8久久精品国产| 欧美乱码精品一区二区三区| 婷婷六月久久综合丁香| 午夜精品久久久久久毛片777| 中文字幕高清在线视频| 国产又色又爽无遮挡免费看| 国产熟女午夜一区二区三区| 亚洲国产精品合色在线| 欧美黑人欧美精品刺激| 久久久久久亚洲精品国产蜜桃av| 亚洲精品一卡2卡三卡4卡5卡| 亚洲成人国产一区在线观看| 两人在一起打扑克的视频| 美女国产高潮福利片在线看| 国产av不卡久久| xxxwww97欧美| 亚洲五月色婷婷综合| 亚洲男人天堂网一区| 成人三级做爰电影| 成人特级黄色片久久久久久久| 19禁男女啪啪无遮挡网站| 美女 人体艺术 gogo| 好看av亚洲va欧美ⅴa在| 国产亚洲精品久久久久久毛片| 久久婷婷人人爽人人干人人爱| 性色av乱码一区二区三区2| 中文字幕人妻熟女乱码| 亚洲国产精品合色在线| 亚洲精品粉嫩美女一区| 性欧美人与动物交配| 国产精品98久久久久久宅男小说| 久热这里只有精品99| 在线观看免费日韩欧美大片| 亚洲自偷自拍图片 自拍| 一区二区日韩欧美中文字幕| 午夜老司机福利片| 黄色成人免费大全| 女生性感内裤真人,穿戴方法视频| 国产亚洲欧美98| 婷婷精品国产亚洲av| 亚洲久久久国产精品| 午夜激情福利司机影院| 在线观看66精品国产| 久久久水蜜桃国产精品网| 免费在线观看成人毛片| 在线观看免费视频日本深夜| 亚洲性夜色夜夜综合| 国产亚洲欧美精品永久| 久久精品91蜜桃| 淫秽高清视频在线观看| 欧美另类亚洲清纯唯美| 999久久久国产精品视频| 欧美激情高清一区二区三区| 精品国产乱子伦一区二区三区| 最好的美女福利视频网| 国产亚洲av嫩草精品影院| 成人欧美大片| 草草在线视频免费看| 久久午夜亚洲精品久久| 国产成+人综合+亚洲专区| 国产欧美日韩一区二区三| 国产成人一区二区三区免费视频网站| 国产精品国产高清国产av| 日本精品一区二区三区蜜桃| 欧美一级毛片孕妇| 好男人电影高清在线观看| 最近最新中文字幕大全免费视频| 亚洲在线自拍视频| 国内揄拍国产精品人妻在线 | 在线观看66精品国产| 国产一区二区激情短视频| 国产v大片淫在线免费观看| 久久精品91蜜桃| 亚洲九九香蕉| 99精品在免费线老司机午夜| 少妇被粗大的猛进出69影院| 成年免费大片在线观看| 国产高清视频在线播放一区| 又大又爽又粗| 亚洲自偷自拍图片 自拍| 99热只有精品国产| 国产亚洲欧美精品永久| 韩国精品一区二区三区| 男人的好看免费观看在线视频 | 最好的美女福利视频网| 久9热在线精品视频| 少妇熟女aⅴ在线视频| 欧美黄色片欧美黄色片| 中出人妻视频一区二区| 国产精品久久久久久亚洲av鲁大| 男人的好看免费观看在线视频 | 美女免费视频网站| 99热6这里只有精品| 欧美精品亚洲一区二区| 成人免费观看视频高清| 亚洲在线自拍视频| 男女做爰动态图高潮gif福利片| 一级作爱视频免费观看| 国产精品亚洲一级av第二区| 久久久久久久精品吃奶| 美国免费a级毛片| 宅男免费午夜| 欧美中文综合在线视频| 波多野结衣高清无吗| 亚洲精华国产精华精| 亚洲av成人一区二区三| 在线天堂中文资源库| 亚洲欧美精品综合久久99| 亚洲精品国产区一区二| 麻豆av在线久日| 91九色精品人成在线观看| 久久天堂一区二区三区四区| 人妻久久中文字幕网| 午夜视频精品福利| 国产成人影院久久av| 成人永久免费在线观看视频| 国产精品九九99| www.精华液| 此物有八面人人有两片| 男人操女人黄网站| 久久国产精品影院| 日日干狠狠操夜夜爽| 一级毛片女人18水好多| 国产欧美日韩一区二区三| 国产亚洲欧美98| 身体一侧抽搐| 欧美成人一区二区免费高清观看 | 亚洲精品av麻豆狂野| 国产不卡一卡二| 18禁黄网站禁片午夜丰满| 又紧又爽又黄一区二区| 两个人视频免费观看高清| 成人欧美大片| 波多野结衣高清作品| 特大巨黑吊av在线直播 | 精品欧美一区二区三区在线| 波多野结衣高清无吗| 在线观看免费日韩欧美大片| 一级黄色大片毛片| 免费搜索国产男女视频| 久久精品亚洲精品国产色婷小说| 欧美日韩中文字幕国产精品一区二区三区| 在线播放国产精品三级| 亚洲片人在线观看| 精品无人区乱码1区二区| 亚洲三区欧美一区| 中国美女看黄片| 日本熟妇午夜| av在线天堂中文字幕| 淫妇啪啪啪对白视频| 男女下面进入的视频免费午夜 | 少妇粗大呻吟视频| 欧美成狂野欧美在线观看| 亚洲国产精品合色在线| 日本免费a在线| 丁香六月欧美| 十八禁网站免费在线| 国内久久婷婷六月综合欲色啪| 最好的美女福利视频网| 在线观看日韩欧美| 亚洲中文字幕一区二区三区有码在线看 | 亚洲成人精品中文字幕电影| 99精品久久久久人妻精品| 少妇 在线观看| 精品国产乱码久久久久久男人| 欧美一级a爱片免费观看看 | 亚洲自拍偷在线| 国产av不卡久久| 久久性视频一级片| 看片在线看免费视频| 国产不卡一卡二| 大型黄色视频在线免费观看| 国产精品九九99| 制服人妻中文乱码| 欧美午夜高清在线| 级片在线观看| 国产又色又爽无遮挡免费看| 欧美成人免费av一区二区三区| 婷婷精品国产亚洲av在线| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美日韩黄片免| 免费看十八禁软件| 麻豆成人av在线观看| 大型av网站在线播放| 黑人操中国人逼视频| 99在线人妻在线中文字幕| 精品久久久久久久久久久久久 | 精品无人区乱码1区二区| 亚洲精品在线美女| 俺也久久电影网| 欧美日韩精品网址| 色播亚洲综合网| 老司机深夜福利视频在线观看| 免费高清视频大片| xxx96com| 亚洲成av片中文字幕在线观看| 黄色丝袜av网址大全| 国产又色又爽无遮挡免费看| 欧美另类亚洲清纯唯美| 99re在线观看精品视频| 成人亚洲精品av一区二区| 久久久国产欧美日韩av| 国内精品久久久久久久电影| 婷婷精品国产亚洲av在线| 国产极品粉嫩免费观看在线| 亚洲 国产 在线| 欧美人与性动交α欧美精品济南到| 国产视频一区二区在线看| 亚洲国产欧美一区二区综合| 在线十欧美十亚洲十日本专区| 国产精品av久久久久免费| 久久国产亚洲av麻豆专区| 国产av不卡久久| 50天的宝宝边吃奶边哭怎么回事| 亚洲成av片中文字幕在线观看| 久久午夜综合久久蜜桃| 欧美黄色片欧美黄色片| 高潮久久久久久久久久久不卡| 亚洲激情在线av| 日本三级黄在线观看| 99re在线观看精品视频| www.精华液| av在线播放免费不卡| 亚洲第一欧美日韩一区二区三区| 美女大奶头视频| 在线av久久热| 成人三级黄色视频| 国产v大片淫在线免费观看| 村上凉子中文字幕在线| 首页视频小说图片口味搜索| 久久国产乱子伦精品免费另类| 亚洲成a人片在线一区二区| 亚洲精品在线观看二区| 久久久久久亚洲精品国产蜜桃av| 亚洲国产精品999在线| xxxwww97欧美| 国产在线观看jvid| 久久久久国产一级毛片高清牌| 2021天堂中文幕一二区在线观 | 国产精品久久久久久精品电影 | 欧美日韩中文字幕国产精品一区二区三区| 热re99久久国产66热| 欧美激情极品国产一区二区三区| 亚洲中文日韩欧美视频| 久久中文字幕人妻熟女| 亚洲va日本ⅴa欧美va伊人久久| 男女之事视频高清在线观看| 欧美激情久久久久久爽电影| 亚洲精华国产精华精| 一二三四在线观看免费中文在| 亚洲精华国产精华精| 变态另类成人亚洲欧美熟女| 大香蕉久久成人网| 成人av一区二区三区在线看| 可以在线观看的亚洲视频| 999久久久国产精品视频| 国产国语露脸激情在线看| 久久精品人妻少妇| 97超级碰碰碰精品色视频在线观看| 婷婷精品国产亚洲av在线| 日韩一卡2卡3卡4卡2021年| 韩国精品一区二区三区| 欧美黄色淫秽网站| 亚洲第一av免费看| 日本三级黄在线观看| 色综合亚洲欧美另类图片| 中文在线观看免费www的网站 | 真人一进一出gif抽搐免费| 身体一侧抽搐| 美女 人体艺术 gogo| 午夜福利在线在线| 国产精品亚洲av一区麻豆| 在线视频色国产色| 欧美中文综合在线视频| 久久天堂一区二区三区四区| 人妻丰满熟妇av一区二区三区| 变态另类丝袜制服| 波多野结衣巨乳人妻| videosex国产| 国产伦一二天堂av在线观看| 狠狠狠狠99中文字幕| 老熟妇乱子伦视频在线观看| 久久久久九九精品影院| 亚洲精品中文字幕在线视频| 午夜老司机福利片| 国产成人系列免费观看| 日本在线视频免费播放| 在线免费观看的www视频| 黄色a级毛片大全视频| АⅤ资源中文在线天堂| 精品高清国产在线一区| 久久人妻福利社区极品人妻图片| www日本在线高清视频| 母亲3免费完整高清在线观看| 国产精品99久久99久久久不卡| 在线观看午夜福利视频| 免费在线观看影片大全网站| 亚洲人成伊人成综合网2020| 啦啦啦免费观看视频1| 午夜成年电影在线免费观看| 亚洲av美国av| 日日摸夜夜添夜夜添小说| 亚洲性夜色夜夜综合| 亚洲va日本ⅴa欧美va伊人久久| 免费人成视频x8x8入口观看| 成人特级黄色片久久久久久久| 最近在线观看免费完整版| 18禁观看日本| 成年女人毛片免费观看观看9| 精华霜和精华液先用哪个| www.精华液| 亚洲,欧美精品.| 亚洲精品av麻豆狂野| 日韩视频一区二区在线观看| 国产精品国产高清国产av| 视频在线观看一区二区三区| 成人av一区二区三区在线看| 在线av久久热| 丰满的人妻完整版| 国内精品久久久久精免费| 久久久久久久午夜电影| 中文字幕精品免费在线观看视频| 麻豆av在线久日| 国产一区二区三区视频了| 好男人在线观看高清免费视频 | 亚洲国产精品久久男人天堂| 99热这里只有精品一区 | 国产亚洲精品久久久久久毛片| 最近在线观看免费完整版| 日韩欧美国产一区二区入口| 成人亚洲精品av一区二区| 宅男免费午夜| 热re99久久国产66热| 搡老妇女老女人老熟妇| 在线观看66精品国产| 好男人在线观看高清免费视频 | 久久伊人香网站| 欧美黑人欧美精品刺激| 91在线观看av| 亚洲精华国产精华精| 身体一侧抽搐| 十分钟在线观看高清视频www| av欧美777| 成人18禁高潮啪啪吃奶动态图| 亚洲精品美女久久av网站| 精品卡一卡二卡四卡免费| 满18在线观看网站| 黄网站色视频无遮挡免费观看| 女警被强在线播放| 午夜福利欧美成人| 亚洲欧美日韩高清在线视频| 精品国产乱码久久久久久男人| 精品福利观看| 中文字幕人妻熟女乱码| 丁香欧美五月| 黄片播放在线免费| 叶爱在线成人免费视频播放| 一级黄色大片毛片| 99久久国产精品久久久| 无遮挡黄片免费观看| 欧洲精品卡2卡3卡4卡5卡区| 变态另类丝袜制服| 色婷婷久久久亚洲欧美| 午夜福利视频1000在线观看| 女人高潮潮喷娇喘18禁视频| 日本三级黄在线观看| 欧美日本视频| 亚洲五月天丁香| 日韩欧美 国产精品| 国产单亲对白刺激| 欧美黄色片欧美黄色片| 好看av亚洲va欧美ⅴa在| 亚洲精华国产精华精| 一进一出抽搐gif免费好疼| 桃色一区二区三区在线观看| 国产欧美日韩一区二区三| 国产1区2区3区精品| 亚洲av熟女| 亚洲 欧美 日韩 在线 免费| 成年免费大片在线观看| 亚洲成人国产一区在线观看| 91在线观看av| 亚洲精品久久国产高清桃花| 天天躁夜夜躁狠狠躁躁| 久久精品国产99精品国产亚洲性色| 午夜免费观看网址| netflix在线观看网站| av超薄肉色丝袜交足视频| 欧美三级亚洲精品| 成人国语在线视频| 亚洲男人天堂网一区| 国产亚洲精品一区二区www| 久久精品成人免费网站| av天堂在线播放| 黄色丝袜av网址大全| 久久久久久大精品| 亚洲专区国产一区二区| 午夜老司机福利片| 国产免费av片在线观看野外av| 日日夜夜操网爽| 精品免费久久久久久久清纯| 欧美日韩黄片免| 欧美中文日本在线观看视频| 亚洲精品色激情综合| 韩国av一区二区三区四区| 久久久久国内视频| 神马国产精品三级电影在线观看 | 欧美成人午夜精品| 精品一区二区三区四区五区乱码| 草草在线视频免费看| 99国产精品99久久久久| 亚洲国产精品成人综合色| 少妇粗大呻吟视频| 亚洲国产中文字幕在线视频| 美女高潮到喷水免费观看| 曰老女人黄片| 国产精品美女特级片免费视频播放器 | 色婷婷久久久亚洲欧美| 欧美人与性动交α欧美精品济南到| 18禁黄网站禁片免费观看直播| 人妻丰满熟妇av一区二区三区| 18禁美女被吸乳视频| 久久婷婷成人综合色麻豆| 精品午夜福利视频在线观看一区| 变态另类丝袜制服| 午夜激情av网站| 国产精品美女特级片免费视频播放器 | 曰老女人黄片| 18禁黄网站禁片午夜丰满| 男女下面进入的视频免费午夜 | 宅男免费午夜| 亚洲国产日韩欧美精品在线观看 | 成人特级黄色片久久久久久久| 欧美黄色淫秽网站| www.999成人在线观看| 国产熟女xx| 最近最新免费中文字幕在线| 国产成人欧美| 成熟少妇高潮喷水视频| 啦啦啦 在线观看视频| 欧美日韩精品网址| 很黄的视频免费| 国产精品久久久久久精品电影 | 亚洲aⅴ乱码一区二区在线播放 | 国产三级黄色录像| 国产乱人伦免费视频| 国产国语露脸激情在线看| 他把我摸到了高潮在线观看| 热re99久久国产66热| 国产成人影院久久av| 无遮挡黄片免费观看| 久久青草综合色| 69av精品久久久久久| 999久久久国产精品视频| 国产片内射在线| 亚洲性夜色夜夜综合| 妹子高潮喷水视频| 婷婷丁香在线五月| 成人免费观看视频高清| 亚洲熟女毛片儿| 女人被狂操c到高潮| 久久久水蜜桃国产精品网| 白带黄色成豆腐渣| 在线看三级毛片| 中文字幕人妻丝袜一区二区| 欧美绝顶高潮抽搐喷水| 可以在线观看的亚洲视频| 久久婷婷人人爽人人干人人爱| 搡老熟女国产l中国老女人| 午夜福利视频1000在线观看| 美女午夜性视频免费| 热re99久久国产66热| 18禁裸乳无遮挡免费网站照片 | 人人妻人人澡人人看| 久久精品aⅴ一区二区三区四区| 中文资源天堂在线| 欧美色视频一区免费| 国产一区在线观看成人免费| 看黄色毛片网站| 午夜免费激情av| 亚洲国产精品sss在线观看| 脱女人内裤的视频| 中亚洲国语对白在线视频| 97碰自拍视频| 黄片播放在线免费| 亚洲黑人精品在线| 欧美又色又爽又黄视频| 国产成人av激情在线播放| 18禁美女被吸乳视频| www.自偷自拍.com| 最新美女视频免费是黄的| 亚洲欧美一区二区三区黑人| 亚洲 国产 在线| 国产片内射在线| 狠狠狠狠99中文字幕| 一区二区三区高清视频在线|