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

    任意形狀工作面沉陷預計計算方法

    2014-08-02 03:57:44王臨清
    金屬礦山 2014年5期
    關鍵詞:角點梯形形狀

    許 冬 王臨清 吳 侃

    (1.中國礦業(yè)大學環(huán)境與測繪學院,江蘇 徐州 221116;2.環(huán)境保護部環(huán)境工程評估中心,北京 100012)

    任意形狀工作面沉陷預計計算方法

    許 冬1王臨清2吳 侃1

    (1.中國礦業(yè)大學環(huán)境與測繪學院,江蘇 徐州 221116;2.環(huán)境保護部環(huán)境工程評估中心,北京 100012)

    為彌補現(xiàn)有任意形狀工作面沉陷預計計算方法存在的不足,更加準確地預計煤炭開采引起的地表移動和變形,基于概率積分法理論,將工作面按照角點劃分為若干梯形,再利用變步長Simpson二重積分法在各梯形區(qū)域按預計公式進行積分計算,最后將各梯形單元產(chǎn)生的影響進行疊加處理,實現(xiàn)了任意形狀工作面的沉陷預計計算。另外,通過研究拐點偏移距對工作面位置的影響,探討了任意形狀工作面實際計算邊界角點坐標的解算方法。在此基礎上設計、開發(fā)了開采沉陷預計程序,實現(xiàn)了地表移動變形預計、預計結果分析及可視化輸出。以某礦區(qū)實際開采工作面為例,對計算方法及預計程序進行了驗證。結果表明:程序設計合理,預計結果準確可靠,符合實際情況,能夠滿足大范圍沉陷預計的需要,可為“三下”采煤、土地復墾規(guī)劃提供必要的技術支持。

    任意形狀工作面 開采沉陷 概率積分法 移動變形 變步長Simpson方法 拐點偏移距

    煤炭資源的開采在我國高速發(fā)展的國民經(jīng)濟中發(fā)揮了重要的作用,同時也給礦區(qū)生態(tài)環(huán)境造成了比較嚴重的影響[1-2]。礦區(qū)開采沉陷預計是礦山開采沉陷學科的核心內(nèi)容之一[3],目前針對任意形狀工作面沉陷預計計算的問題,國內(nèi)學者從不同角度出發(fā),提出了各具特色的解決辦法[4-12]。文獻[4-6]利用數(shù)學手段將面積分轉(zhuǎn)化為線積分進行非矩形工作面的沉陷預計,雖然實現(xiàn)了任意煤層傾角任意形狀工作面的地表移動與變形計算,但存在數(shù)據(jù)準備繁瑣,不能快速實現(xiàn)多工作面預計計算等問題。文獻[7-8]提出將不規(guī)則工作面按照角點劃分成若干三角形,用多重積分的高斯方法求解地表移動和變形的方法,雖然預計結果精度較高,但存在積分上下限確定困難,計算速度慢等問題。文獻[9-10]提出采用將非矩形工作面劃分為若干個近似小矩形工作面,最后將預計結果進行疊加處理的方法進行預計計算。雖然方法操作簡單、便于理解,但存在預計結果帶有近似性、不能處理凹多邊形工作面、下沉等值線在開采邊界處呈現(xiàn)“鋸齒狀”等問題。

    筆者借鑒前人研究成果,應用概率積分法的基本原理,提出了任意形狀工作面沉陷預計計算的具體算法,并開發(fā)了一個能夠?qū)θ我庑螤疃喙ぷ髅孢M行沉陷預計與分析的程序,以期為評價“三下”采煤的可行性提供依據(jù)。

    1 理論模型及關鍵算法

    1.1 理論模型

    概率積分法因其所用的移動和變形預計公式中含有概率積分而得名[2],是以隨機介質(zhì)理論為基礎,以正態(tài)分布函數(shù)為影響函數(shù),用積分式準確合理描述下沉盆地的方法。基于概率積分法的全盆地移動與變形預計計算公式如下[11]。

    下沉:

    (1)

    水平移動:

    (2)

    (3)

    式中,W(x,y)為地表任意點(坐標為x、y)的下沉值,mm;Ux(x,y)為地表任意點走向方向的水平移動值,mm;Uy(x,y)為地表任意點傾斜方向的水平移動值,mm;Wcm、Ucm為地表充分采動的最大下沉值和最大水平移動值,mm;r為走向主斷面上采空區(qū)邊界,m;x,y為計算點相對坐標(考慮拐點偏移距);D為開采煤層區(qū)域;η,ξ分別為x,y方向的變形值,mm。

    1.2 梯形劃分與積分計算

    由式(1)~式(3)可以看出,按照全盆地移動與變形公式進行預計計算主要是在區(qū)域D上進行概率積分運算。對于這類復雜積分問題的求解,只有給出被積函數(shù)和積分上、下限才可進行數(shù)值積分計算。積分限與開采區(qū)域D有關,當開采區(qū)域是非矩形工作面時很難直接確定其積分限[8]。本研究嘗試采用將工作面按照角點劃分為若干梯形的方法來解決概率積分計算困難的問題,具體方法如下。

    設開采區(qū)域D為任意凸多邊形工作面(如圖1(a)所示),從A~G各頂點向X軸做垂線,按照A—B—C—E—F—G—A的角點順序(逆時針方向)將開采區(qū)域劃分為AA′B′B、BB′C′C、CC′EE′、EE′FF′、FF′GG′、GG′AA′共6個梯形單元,對于每個梯形單元可以唯一確定其積分的上下限,如AA′B′B單元的積分限如式(4)所示:

    (4)

    圖1 積分區(qū)域劃分Fig.1 Integral area division

    每個梯形單元的積分上、下限確定后,分別采用高效的變步長Simpson二重積分法進行數(shù)值積分計算,計算完成后將各梯形單元產(chǎn)生的影響進行疊加處理來計算整個采空區(qū)造成的地表移動變形。由于設定工作面角點按照逆時針順序(A—B—C—E—F—G—A)排列,根據(jù)正負判斷法則(如果x1x2,則S11′2′2前符號為正),開采區(qū)域D的移動變形值可用下式表示:

    一處狹窄的巷子入口,巷子里面黑黢黢的,黑的像是要把我吞沒。入口處一盞昏暗的路燈,仿佛隨時要被夜風吹滅。踩踏板的腳頓了一下。但我終究還是向這吃人般的黑暗騎去。

    (5)

    當工作面為凹多邊形時,仍可按凸多邊形的處理方法進行處理,凹多邊形工作面的梯形劃分如圖1(b)所示。因此,無論工作面是凸多邊形還是凹多邊形,均可使用梯形劃分方法,確定積分限進行概率積分計算,求解地表移動變形,使預計計算不再受開采工作面布置形式的影響。

    1.3 工作面計算邊界求取

    拐點偏移距是開采沉陷預計過程中的一個重要參數(shù),它的大小直接影響到地表下沉盆地的形狀、范圍以及移動變形值的大小[13-15]。考慮拐點偏移距后,工作面的開采邊界將會發(fā)生變化,求取任意形狀工作面實際開采邊界的主要步驟是:①根據(jù)工作面角點坐標計算工作面的重心;②以重心為交點分別做X、Y軸的平行線,將工作面劃分為4個子區(qū)域(Ⅰ~Ⅳ);③先判斷角點屬于哪一個子區(qū)域,再根據(jù)其所屬區(qū)域的不同選擇相應的計算方法,如圖2所示。圖2中S1~S4分別為走向左、右邊界拐點偏移距和上、下山邊界拐點偏移距。若在采空區(qū)一側(cè),取正值,在煤柱一側(cè),取負值[2]。

    圖2 分區(qū)域計算方法Fig.2 Calculation method of partition

    由于煤層存在一定的傾角,使計算邊界與考慮拐點偏移距后的實際開采邊界不完全一致,兩者在煤層傾斜方向存在一個偏移量,該偏移量d可按式(6)計算:

    (6)

    式中,h為下邊界采深,m;S4為下山邊界拐點偏移距,m;α為煤層傾角,(°);θ為開采影響傳播角,(°)。

    當偏移量d為正值時,從實際下山邊界向下山方

    2 程序?qū)崿F(xiàn)流程

    為了能夠更加方便、準確地進行開采沉陷預計計算,應用軟件工程思想并結合工程實際需求,以“數(shù)據(jù)”為核心,按照數(shù)據(jù)準備、處理和分析的原則對任意形狀工作面沉陷預計程序進行設計與開發(fā),主要實現(xiàn)地表移動變形預計、預計結果分析及可視化輸出等功能。程序?qū)崿F(xiàn)的基本流程是:①構建預計文件。從采掘工程平面圖上獲得工作面尺寸信息、開采參數(shù)、沉陷預計參數(shù),并圈定預計范圍、設定格網(wǎng)間距。②沉陷預計計算。根據(jù)預計文件中參數(shù)進行坐標轉(zhuǎn)換、求取工作面計算邊界,將計算邊界按角點劃分為若干梯形,分別對每一個梯形進行二重積分將積分結果進行疊加處理,得到下沉、水平移動等預計結果,進而繪制移動變形等值線、云圖等圖件;其中,對于多個工作面的沉陷預計采用將每個工作面產(chǎn)生的影響進行疊加的方法得到最終的地表移動變形值[16-17]。③沉陷結果分析。根據(jù)得到的預計結果可計算煤炭開采導致的沉陷體積和沉陷面積、制作開采影響后的礦區(qū)地形圖、評價區(qū)域內(nèi)建筑物的損壞等級、分析沉陷范圍內(nèi)的積水情況[18-19]。程序?qū)崿F(xiàn)的流程如圖3所示。

    圖3 程序?qū)崿F(xiàn)流程Fig.3 The program implementation process

    3 實例應用

    3.1 實驗區(qū)概況

    為測試預計結果的準確性,選取某礦區(qū)實際開采工作面進行任意形狀、多工作面的沉陷預計計算,工作面如圖4所示。該測試算例既包含凹多邊形工作面(如工作面1、2)又有凸多邊形工作面(如工作面3),能夠滿足程序測試的要求。

    圖4 工作面布置Fig.4 Layout of working face

    結合礦區(qū)實測資料,通過綜合分析,確定了沉陷預計時所使用的預計參數(shù),見表1。

    3.2 沉陷預計結果

    按照構建預計文件、沉陷預計計算、預計結果分析的流程對實驗區(qū)的3個工作面進行處理,生成下沉、傾斜、曲率、水平移動、水平變形等值線及下沉云圖。由于篇幅有限,只給出部分結果圖件,其中圖5為下沉等值線圖,圖6為下沉云圖,預計數(shù)值結果見表2。實例應用表明,通過采用基于梯形劃分的任意形狀工作面沉陷預計計算方法可使預計結果與實際情況更加符合。

    4 結 論

    (1)本研究采用將工作面按照角點劃分為若干梯形單元,用預計計算公式計算每一個梯形單元造成的地表移動變形,最后再將各小單元的影響進行疊加來計算整個采空區(qū)造成的地表移動變形的方法實現(xiàn)了任意形狀工作面的沉陷預計計算,解決了概率積分求解困難的問題,使預計計算不再受開采工作面布置形式的影響,為任意形狀工作面開采沉陷預計計算提供了新的解決思路。

    (2)通過研究拐點偏移距對工作面實際位置的影響,探討了求取任意形狀工作面實際計算邊界角點坐標的方法,使預計結果更符合地表移動與變形的實際情況。實現(xiàn)了多工作面沉陷預計計算,使預計結果能夠滿足礦區(qū)大范圍開采沉陷預計的需要。

    表1 開采沉陷預計參數(shù)Table 1 Prediction parameters of mining subsidence

    圖5 下沉等值線 (單位:m)Fig.5 Subsidence equivalent curve

    圖6 下沉云圖(單位:mm)Fig.6 Subsidence nephogram表2 部分預計數(shù)值結果

    Table 2 Part of prediction values mm

    (3)通過研究相關沉陷預計理論、算法,設計、開發(fā)了沉陷預計程序,實現(xiàn)了地表移動變形預計、預計結果分析及可視化輸出。通過實例應用表明:沉陷預計結果準確、可靠,程序具有良好的實用性與可擴展性,能夠簡單方便地進行開采沉陷預計,為“三下”采煤方案優(yōu)化設計、開展礦區(qū)土地復墾、保護礦區(qū)生態(tài)環(huán)境提供科學依據(jù)。

    [1] 呂偉才,楊世清.開采沉陷對礦區(qū)生態(tài)環(huán)境的影響評價[J].江蘇環(huán)境科技,2006,19(S1):57-59. Lu Weicai,Yang Shiqing.The effect evaluation on the mining area ecological environment by the mining subsidence[J].Jiangsu Environmental Science and Techonlogy,2006,19(S1):57-59.

    [2] 何國清,楊 倫,凌賡娣,等.礦山開采沉陷學[M].徐州:中國礦業(yè)大學出版社,1994. He Guoqing,Yang Lun,Ling Gengdi,et al.Mining Subsidence Theory[M].Xuzhou:China University of Mining and Technology Press,1994.

    [3] 吳 侃,黃珍珍,韓奎峰.礦區(qū)土地破壞預測[J].山東煤炭科技,2006(3):48-49. Wu Kan,Huang Zhenzhen,Han Kuifeng.Prediction of the destruction on mining land[J].Shandong Coal Science and Technology,2006(3):48-49.

    [4] 張華興.地表移動計算程序的可視化[J].煤炭科學技術,1998(11):29-30. Zhang Huaxing.The visualization of earth′s surface mobile computing program[J].Coal Science and Technology,1998(11):29-30.

    [5] 張華興.地表移動與變形的等價轉(zhuǎn)換線積分計算[J].煤礦開采,1991(1):14-19. Zhang Huaxing.The movement and deformation of ground equivalent conversion line integral calculation[J].Coal Mining Technology,1991(1):14-19.

    [6] 張華興,焦傳武.地表移動精確預計法——線積分運算[J].礦山測量,1988(12):21-29. Zhang Huaxing,Jiao Chuanwu.Precise method of surface movement expected:line integral operation[J].Mine Surveying,1988(12):21-29.

    [7] 康建榮.采動覆巖動態(tài)移動破壞規(guī)律及開采沉陷預計系統(tǒng)(MSPS)研究[D].北京:中國礦業(yè)大學,1999. Kang Jianrong.The Rule of Overburden Dynamic Movement and Damage and Mining Subsidence Prediction System[D].Beijing:China University of Mining and Technology,1999.

    [8] 康建榮,王金莊,溫澤民.任意形多工作面多線段開采沉陷預計系統(tǒng)(MSPS)[J].礦山測量,2000,1(3):24-27. Kang Jianrong,Wang Jinzhuang,Wen Zemin.Arbitrarily shaped multi-face multi-line of mining subsidence prediction system (MSPS)[J].Mine Surveying,2000,1(3):24-27.

    [9] 吳 侃,周 鳴.開采沉陷預測預報系統(tǒng)[M].徐州:中國礦業(yè)大學出版社,1999. Wu Kan,Zhou Ming.Mining Subsidence Prediction System[M].Xuzhou:China University of Mining and Technology Press,1999.

    [10] 吳 侃,葛家新,王玲丁,等.開采沉陷預計一體化方法[M].徐州:中國礦業(yè)大學出版社,1998. Wu Kan,Ge Jiaxin,Wang Lingding,et al.Integration Method of Mining Subsidence Prediction[M].Xuzhou:China University of Mining and Technology Press,1998.

    [11] 國家煤炭工業(yè)局.建筑物、水體、鐵路及主要井巷煤柱留設與壓煤開采規(guī)程[S].北京:煤炭工業(yè)出版社,2004. State Bureau of Coal Industry.Mining Regulations of Buildings,Water,Railway and Main Roadway Pillar Design and Pressure Coal[S].Beijing:Coal Industry Press,2004.

    [12] 顧 偉,譚志祥,鄧喀中.基于Active X Automation 的三維可視開采沉陷預計[J].金屬礦山,2013(3):119-122. Gu Wei,Tan Zhixiang,Deng Kazhong.3D visualization prediction of mining subsidence system based on the Active X Automation[J].Metal Mine,2013(3):119-122.

    [13] 陳 巖,奚硯濤.基于Surfer的金屬礦山開采沉陷預測和三維可視化[J].金屬礦山,2012(5):149-151. Chen Yan,Xi Yantao.Metal mine mining subsidence prediction and 3D visualization based on surfer[J].Metal Mine,2012(5):149-151.

    [14] 周萬茂,張華興,何瑞華.任意形狀工作面拐點移動距求取方法[J].煤礦開采,2000(12):13-16. Zhou Wanmao,Zhang Huaxing,He Ruihua.A method for determining the displacement distance of the inflection point of the mining faces with arbitrary shapes[J].Coal Mining Technology,2000(12):13-16.

    [15] 張合兵,王世東,李春意.煤礦區(qū)土地復墾中土地損毀預測方法研究[J].金屬礦山,2013(9):142-147. Zhang Hebing,Wang Shidong,Li Chunyi.Research on land destruction prediction methods of land reclamation in coal mining area[J].Metal Mine,2013(9):142-147.

    [16] 朱劉娟,陳俊杰,鄒友峰.任意形狀工作面開采地表移動變形預計的算法實現(xiàn)[J].遼寧工程技術大學學報:自然科學版,2005,24(3):337-339. Zhu Liujuan,Chen Junjie,Zou Youfeng.Algorithm of ground movement and deformation prediction in mining of random shape face[J].Journal of Liaoning Technical University:Natural Science Edition,2005,24(3):337-339.

    [17] 李 明,李琰慶,王紅梅.任意形狀工作面開采沉陷預計系統(tǒng)開發(fā)[J].礦業(yè)安全與環(huán)保,2008,35(5):18-21. Li Ming,Li Yanqing,Wang Hongmei.Development of predicting system for mining subsidence of arbitrarily shaped face[J].Mining Safety & Environmental Protection,2008,35(5):18-21.

    [18] 劉 虎,吳 侃,陳冉麗.基于AutoCAD的開采沉陷預計分析軟件研究[J].礦業(yè)安全與環(huán)保,2009,36(2):22-23. Liu Hu,Wu Kan,Chen Ranli.Research on mining subsidence prediction and analysis system based on AutoCAD[J].Mining Safety & Environmental Protection,2009,36(2):22-23.

    [19] 石秀偉,李 晶,趙換新,等.基于ArcGIS的煤礦開采沉陷預測與可視化分析[J].金屬礦山,2012(9):103-106. Shi Xiuwei,Li Jing,Zhao Huanxin,et al.The prediction and visualization of coal mining subsidence based on ArcGIS[J].Metal Mine,2012(9):103-106.

    (責任編輯 徐志宏)

    Mining Subsidence Prediction Calculation Methods of Random Shape Working Face

    Xu Dong1Wang Linqing2Wu Kan1

    (1.ChinaUniversityofMiningandTechnology,SchoolofEnvironmentalScienceandSpatialInformaticsXuzhou221116,China;2.AppraisalCenterforEnvironment&EngineeringMinistryofEnvironmentalProtection,Beijing100012,China)

    In order to make up for the shortage in mining subsidence prediction calculation methods of existing random shape working face,and more accurately predict the value of surface movement and deformation caused by coal mining,the working face is divided into several trapezoids according to corners based on the theory of probability integral method.The variable step-size Simpson double integral method is adopted to conduct the integration in every trapezoidal region according to the prediction formula,and deal with the effects which were produced by every trapezoidal unit to realize the subsidence prediction calculation of random shape working face through superposition method.In addition,through the study of inflection point displacement effect on the working face position,the solution method for actually calculating angular coordinate at boundary of the random shape working face is discussed.On this basis,the mining subsidence prediction program is designed and developed,realizing prediction of surface movement and distortion,analysis of prediction results and visual output.The actual mine working faces in a mining area was taken as a case to verify the calculation method and prediction program.The results showed that the program design is reasonable,and the prediction result is accurate and reliable,and accords with the fact.It can meet the needs for a wide range of subsidence prediction and provide the necessary technical support for “under three” coal mining and land reclamation planning.

    Random shape working face,Mining subsidence,Probability integral method,Movement and deformation,Variable step-size Simpson-method,Inflection point displacement

    2014-01-02

    “十二五”國家科技支撐計劃項目(編號:2012BAC04B03)。

    許 冬(1990—),男,碩士研究生。

    TD325+.2

    A

    1001-1250(2014)-05-055-05

    猜你喜歡
    角點梯形形狀
    玩轉(zhuǎn)梯形
    幼兒100(2023年37期)2023-10-23 11:38:58
    挖藕 假如悲傷有形狀……
    梯形達人
    一類變延遲中立型微分方程梯形方法的漸近估計
    你的形狀
    基于FAST角點檢測算法上對Y型與X型角點的檢測
    梯形
    啟蒙(3-7歲)(2017年6期)2017-11-27 09:34:55
    基于邊緣的角點分類和描述算法
    電子科技(2016年12期)2016-12-26 02:25:49
    看到的是什么形狀
    基于圓環(huán)模板的改進Harris角點檢測算法
    成人漫画全彩无遮挡| 国产精品99久久久久久久久| 国产精品久久久久久精品古装| 国产成人aa在线观看| 久久韩国三级中文字幕| 三级男女做爰猛烈吃奶摸视频| 国产白丝娇喘喷水9色精品| kizo精华| 波野结衣二区三区在线| 色视频www国产| 亚洲国产av新网站| 各种免费的搞黄视频| 女人十人毛片免费观看3o分钟| 婷婷色综合大香蕉| 午夜精品国产一区二区电影 | 日韩电影二区| 欧美日韩在线观看h| 美女高潮的动态| 亚洲精品日本国产第一区| av一本久久久久| 97在线视频观看| 亚洲精品国产色婷婷电影| 免费观看性生交大片5| 久久国内精品自在自线图片| 免费观看无遮挡的男女| 国产伦精品一区二区三区四那| 国产精品伦人一区二区| 韩国高清视频一区二区三区| 插阴视频在线观看视频| 亚洲精品国产av成人精品| av女优亚洲男人天堂| 又爽又黄无遮挡网站| 国产久久久一区二区三区| 成人亚洲精品av一区二区| 亚洲不卡免费看| 国产精品麻豆人妻色哟哟久久| 中文字幕久久专区| 一级毛片aaaaaa免费看小| 欧美激情国产日韩精品一区| 亚洲成人av在线免费| 中国美白少妇内射xxxbb| 日韩国内少妇激情av| 日韩欧美 国产精品| 久久99蜜桃精品久久| 九草在线视频观看| 精品久久久噜噜| 久久精品国产亚洲av涩爱| 亚洲三级黄色毛片| 国产 一区 欧美 日韩| 成人午夜精彩视频在线观看| 免费大片黄手机在线观看| 永久网站在线| 国产亚洲5aaaaa淫片| 久久久久久久大尺度免费视频| 日韩成人av中文字幕在线观看| 男人和女人高潮做爰伦理| 欧美成人精品欧美一级黄| 熟妇人妻不卡中文字幕| 伊人久久国产一区二区| 国产成人a∨麻豆精品| 亚洲国产高清在线一区二区三| 中文资源天堂在线| 成年女人在线观看亚洲视频 | 亚洲av日韩在线播放| 色吧在线观看| 老司机影院成人| av在线亚洲专区| 交换朋友夫妻互换小说| 久久久精品94久久精品| 久久精品久久久久久久性| 色5月婷婷丁香| 男女啪啪激烈高潮av片| 日日撸夜夜添| 高清日韩中文字幕在线| 一级片'在线观看视频| 黑人高潮一二区| 欧美日韩视频高清一区二区三区二| 久久久久性生活片| 久久久久国产网址| 亚洲av.av天堂| 欧美成人a在线观看| 国产美女午夜福利| 99热这里只有是精品在线观看| 亚洲精品第二区| 久久久久久久精品精品| 91在线精品国自产拍蜜月| 国产精品久久久久久精品电影小说 | 亚洲精品成人久久久久久| 亚洲精品久久午夜乱码| 高清毛片免费看| 99精国产麻豆久久婷婷| 色吧在线观看| 成年人午夜在线观看视频| 一级毛片aaaaaa免费看小| 日韩av在线免费看完整版不卡| 亚洲精品国产成人久久av| 韩国av在线不卡| 亚洲精华国产精华液的使用体验| 久久精品久久久久久噜噜老黄| 天天躁日日操中文字幕| 国产成人精品一,二区| 国内精品美女久久久久久| 亚洲内射少妇av| 欧美xxxx性猛交bbbb| 免费看不卡的av| 街头女战士在线观看网站| 嫩草影院新地址| 欧美高清性xxxxhd video| 国产有黄有色有爽视频| 国产黄a三级三级三级人| 最近2019中文字幕mv第一页| 麻豆成人av视频| 国产一区二区在线观看日韩| 免费观看的影片在线观看| 亚洲三级黄色毛片| 免费看a级黄色片| 国产黄片视频在线免费观看| 涩涩av久久男人的天堂| 午夜福利高清视频| av在线播放精品| 别揉我奶头 嗯啊视频| 九九在线视频观看精品| 欧美少妇被猛烈插入视频| 在线免费十八禁| 国产又色又爽无遮挡免| 日韩精品有码人妻一区| 亚洲精品国产av成人精品| 成年女人在线观看亚洲视频 | 亚洲图色成人| av.在线天堂| 2018国产大陆天天弄谢| 亚洲精品456在线播放app| 亚洲欧美日韩东京热| 精品久久久久久电影网| 麻豆成人午夜福利视频| 国产精品人妻久久久久久| 国产精品国产三级专区第一集| 黄色怎么调成土黄色| 国产精品国产三级国产专区5o| 国产高清三级在线| 日韩成人伦理影院| 欧美人与善性xxx| h日本视频在线播放| 边亲边吃奶的免费视频| 亚洲怡红院男人天堂| 禁无遮挡网站| 赤兔流量卡办理| 欧美日韩在线观看h| 美女内射精品一级片tv| 美女被艹到高潮喷水动态| 日本与韩国留学比较| 一二三四中文在线观看免费高清| 别揉我奶头 嗯啊视频| 中国美白少妇内射xxxbb| 亚洲经典国产精华液单| 只有这里有精品99| 精品久久久久久久久亚洲| 亚洲自偷自拍三级| 啦啦啦中文免费视频观看日本| 少妇熟女欧美另类| 国产色婷婷99| 欧美最新免费一区二区三区| 亚洲精品亚洲一区二区| 国内揄拍国产精品人妻在线| 国产精品久久久久久精品电影| 国产男女内射视频| 亚洲va在线va天堂va国产| 精品99又大又爽又粗少妇毛片| 人人妻人人爽人人添夜夜欢视频 | 国产精品熟女久久久久浪| 亚洲av中文字字幕乱码综合| 九九久久精品国产亚洲av麻豆| 亚洲av免费高清在线观看| 亚洲精品一二三| eeuss影院久久| 男插女下体视频免费在线播放| 少妇的逼水好多| 国产一区二区在线观看日韩| 亚洲人与动物交配视频| 成人免费观看视频高清| 永久网站在线| 久久国产乱子免费精品| 国产精品一区二区性色av| 亚洲精品乱码久久久久久按摩| 色吧在线观看| 亚洲av一区综合| 在现免费观看毛片| 免费黄网站久久成人精品| 色5月婷婷丁香| 一区二区三区精品91| av一本久久久久| 最近2019中文字幕mv第一页| 国产精品不卡视频一区二区| 亚洲欧洲国产日韩| 建设人人有责人人尽责人人享有的 | 国产色婷婷99| 日韩不卡一区二区三区视频在线| 日韩 亚洲 欧美在线| 国产成人精品福利久久| 国产av不卡久久| 欧美3d第一页| 色网站视频免费| 国国产精品蜜臀av免费| 极品少妇高潮喷水抽搐| 久久国产乱子免费精品| 国产黄色视频一区二区在线观看| 美女脱内裤让男人舔精品视频| 搡老乐熟女国产| 九九爱精品视频在线观看| 边亲边吃奶的免费视频| 成人无遮挡网站| 18+在线观看网站| 99久久精品热视频| 少妇 在线观看| h日本视频在线播放| 久久久久久国产a免费观看| 中文字幕亚洲精品专区| 亚洲一区二区三区欧美精品 | 我的老师免费观看完整版| 久久综合国产亚洲精品| 又黄又爽又刺激的免费视频.| 国产亚洲5aaaaa淫片| 亚洲国产最新在线播放| 看黄色毛片网站| 亚洲国产精品专区欧美| 极品教师在线视频| 午夜福利高清视频| 精品久久久精品久久久| 亚洲精品一二三| 日日撸夜夜添| 成年av动漫网址| 91午夜精品亚洲一区二区三区| 国精品久久久久久国模美| 国产探花在线观看一区二区| 免费av毛片视频| 韩国高清视频一区二区三区| 午夜福利高清视频| 亚洲最大成人手机在线| 亚洲色图av天堂| 国产亚洲午夜精品一区二区久久 | 一级a做视频免费观看| 亚洲精品国产av蜜桃| 中文字幕av成人在线电影| 午夜福利视频精品| 日日摸夜夜添夜夜添av毛片| 久久久国产一区二区| 久久久久久久国产电影| 2022亚洲国产成人精品| 日本欧美国产在线视频| 国产毛片a区久久久久| 夫妻午夜视频| 国产日韩欧美在线精品| 精品国产三级普通话版| 国产精品蜜桃在线观看| freevideosex欧美| 大香蕉97超碰在线| 最近中文字幕2019免费版| 欧美一区二区亚洲| 亚洲人成网站在线观看播放| 一级片'在线观看视频| 欧美极品一区二区三区四区| 国产午夜精品久久久久久一区二区三区| 99久国产av精品国产电影| 国产成人freesex在线| 亚洲欧洲日产国产| 夜夜爽夜夜爽视频| 最近最新中文字幕免费大全7| 欧美3d第一页| 欧美亚洲 丝袜 人妻 在线| 亚洲精品国产av成人精品| 久久久久国产网址| 观看美女的网站| 日韩视频在线欧美| 最近最新中文字幕免费大全7| 又爽又黄无遮挡网站| 男的添女的下面高潮视频| 欧美xxxx黑人xx丫x性爽| 久久精品国产亚洲av涩爱| 亚洲国产欧美人成| 青春草视频在线免费观看| 少妇丰满av| 麻豆乱淫一区二区| 日韩人妻高清精品专区| 男人和女人高潮做爰伦理| 一个人看的www免费观看视频| 成人漫画全彩无遮挡| 美女被艹到高潮喷水动态| 最近最新中文字幕免费大全7| 欧美 日韩 精品 国产| 男的添女的下面高潮视频| 精品久久久精品久久久| 五月开心婷婷网| 青青草视频在线视频观看| 秋霞伦理黄片| 日本猛色少妇xxxxx猛交久久| 国产日韩欧美亚洲二区| 菩萨蛮人人尽说江南好唐韦庄| 国产一级毛片在线| 国产午夜精品久久久久久一区二区三区| 国产 精品1| 国产老妇女一区| 最近中文字幕高清免费大全6| 欧美高清成人免费视频www| 高清av免费在线| 国产大屁股一区二区在线视频| 欧美日韩精品成人综合77777| 人妻系列 视频| 午夜福利网站1000一区二区三区| 国产老妇伦熟女老妇高清| 插逼视频在线观看| 大片免费播放器 马上看| 不卡视频在线观看欧美| 在线观看美女被高潮喷水网站| 另类亚洲欧美激情| 精品一区在线观看国产| 国产在线男女| 久久久久久久午夜电影| 久热这里只有精品99| 亚洲av免费在线观看| 国产乱人偷精品视频| 日韩精品有码人妻一区| 欧美极品一区二区三区四区| 欧美日韩精品成人综合77777| 青春草国产在线视频| 日韩一本色道免费dvd| 国产精品嫩草影院av在线观看| 一个人看的www免费观看视频| 国产高清有码在线观看视频| 亚洲性久久影院| 综合色av麻豆| 国产精品国产av在线观看| 2018国产大陆天天弄谢| 日韩一区二区视频免费看| 在线观看三级黄色| 大又大粗又爽又黄少妇毛片口| 亚洲av男天堂| 久久久a久久爽久久v久久| 国产精品久久久久久精品电影| 成年av动漫网址| 少妇被粗大猛烈的视频| 免费黄网站久久成人精品| 亚洲精品成人久久久久久| 久久精品久久久久久久性| 日产精品乱码卡一卡2卡三| 狂野欧美白嫩少妇大欣赏| 交换朋友夫妻互换小说| 亚洲四区av| 日本免费在线观看一区| 一本一本综合久久| 亚洲精品aⅴ在线观看| 亚洲最大成人手机在线| 欧美区成人在线视频| 欧美zozozo另类| 日本熟妇午夜| 91精品伊人久久大香线蕉| 99久久九九国产精品国产免费| 91在线精品国自产拍蜜月| 日韩国内少妇激情av| 午夜亚洲福利在线播放| 亚洲在久久综合| 少妇裸体淫交视频免费看高清| 亚洲精华国产精华液的使用体验| 亚州av有码| 欧美 日韩 精品 国产| 亚洲最大成人手机在线| 波多野结衣巨乳人妻| 午夜免费观看性视频| av国产免费在线观看| 精品人妻熟女av久视频| 成人高潮视频无遮挡免费网站| 日本wwww免费看| 热re99久久精品国产66热6| 成人毛片a级毛片在线播放| 午夜福利视频精品| 亚洲精品一区蜜桃| 国产午夜精品久久久久久一区二区三区| 能在线免费看毛片的网站| 国产精品国产av在线观看| 日本与韩国留学比较| 熟女人妻精品中文字幕| 波野结衣二区三区在线| 男女国产视频网站| 久久久久久久国产电影| 久久热精品热| 中文乱码字字幕精品一区二区三区| 日日摸夜夜添夜夜添av毛片| 亚洲av免费在线观看| 女人被狂操c到高潮| 成人免费观看视频高清| 97在线视频观看| 亚洲美女搞黄在线观看| 交换朋友夫妻互换小说| 欧美zozozo另类| 97超碰精品成人国产| 免费电影在线观看免费观看| 国产黄色视频一区二区在线观看| 麻豆久久精品国产亚洲av| 九色成人免费人妻av| 蜜桃久久精品国产亚洲av| 男的添女的下面高潮视频| 欧美日韩国产mv在线观看视频 | 狠狠精品人妻久久久久久综合| 天天躁夜夜躁狠狠久久av| 亚洲国产精品成人综合色| 欧美潮喷喷水| av免费观看日本| 婷婷色综合www| h日本视频在线播放| 噜噜噜噜噜久久久久久91| 久久久久久久亚洲中文字幕| 精品人妻视频免费看| 精品一区二区三区视频在线| 免费看光身美女| 久久影院123| 国产精品麻豆人妻色哟哟久久| 男女下面进入的视频免费午夜| 国产 精品1| 久久久久国产网址| 欧美性感艳星| 精品久久久噜噜| 亚洲av.av天堂| 亚洲不卡免费看| 亚洲欧美日韩无卡精品| 简卡轻食公司| 色网站视频免费| 永久免费av网站大全| 亚洲最大成人中文| 亚洲成人中文字幕在线播放| 午夜福利视频精品| 最近最新中文字幕大全电影3| 国内少妇人妻偷人精品xxx网站| 狂野欧美激情性xxxx在线观看| 国产亚洲一区二区精品| 国产探花极品一区二区| 亚洲欧洲国产日韩| 亚洲av免费高清在线观看| 国产在线男女| 亚洲精品,欧美精品| 搡老乐熟女国产| 免费黄色在线免费观看| 免费人成在线观看视频色| 交换朋友夫妻互换小说| 国产成人福利小说| 丝袜脚勾引网站| 人妻 亚洲 视频| av网站免费在线观看视频| av又黄又爽大尺度在线免费看| 成人无遮挡网站| 日本三级黄在线观看| av福利片在线观看| 一区二区av电影网| 少妇猛男粗大的猛烈进出视频 | 亚洲丝袜综合中文字幕| 交换朋友夫妻互换小说| 精品久久国产蜜桃| 亚洲国产精品成人久久小说| 中文字幕久久专区| 男女啪啪激烈高潮av片| 日韩强制内射视频| 午夜激情福利司机影院| 少妇的逼水好多| 只有这里有精品99| 亚洲精品影视一区二区三区av| 久久精品国产亚洲av天美| 国产男女超爽视频在线观看| www.色视频.com| 午夜福利视频1000在线观看| 久久久久久久亚洲中文字幕| 91狼人影院| 亚洲av中文av极速乱| 久久久久九九精品影院| 青春草视频在线免费观看| 国产有黄有色有爽视频| 亚洲国产精品专区欧美| 18禁在线无遮挡免费观看视频| 久久精品国产a三级三级三级| 免费av不卡在线播放| 亚洲欧美成人综合另类久久久| 欧美日韩视频高清一区二区三区二| 久久久久性生活片| 伦精品一区二区三区| 26uuu在线亚洲综合色| 六月丁香七月| 最近最新中文字幕大全电影3| 国产精品久久久久久久电影| 97超视频在线观看视频| 精品一区二区三区视频在线| 熟女av电影| 国内揄拍国产精品人妻在线| 日韩欧美精品免费久久| 国产成人a区在线观看| 中文字幕av成人在线电影| 青春草亚洲视频在线观看| 免费看日本二区| 99久久九九国产精品国产免费| 七月丁香在线播放| 男人爽女人下面视频在线观看| 日韩成人伦理影院| 欧美bdsm另类| 高清欧美精品videossex| 成人午夜精彩视频在线观看| 欧美日韩视频精品一区| 一区二区三区免费毛片| 免费人成在线观看视频色| 亚洲欧美日韩无卡精品| 最后的刺客免费高清国语| 国产探花在线观看一区二区| 免费看a级黄色片| 一区二区三区免费毛片| 国产成人一区二区在线| 国产成人精品婷婷| 国产欧美日韩一区二区三区在线 | 肉色欧美久久久久久久蜜桃 | 男人添女人高潮全过程视频| 亚洲精品色激情综合| 伊人久久国产一区二区| 国产亚洲5aaaaa淫片| 国产中年淑女户外野战色| 国产永久视频网站| 超碰av人人做人人爽久久| 国产69精品久久久久777片| 日韩国内少妇激情av| 又大又黄又爽视频免费| 亚洲精品影视一区二区三区av| 成人亚洲精品av一区二区| 三级经典国产精品| 色吧在线观看| 亚洲国产色片| 午夜免费男女啪啪视频观看| 亚洲激情五月婷婷啪啪| 国产精品女同一区二区软件| 国产成人freesex在线| 视频区图区小说| 插阴视频在线观看视频| 男人添女人高潮全过程视频| 日韩一本色道免费dvd| 国产精品女同一区二区软件| 18禁动态无遮挡网站| 欧美日韩综合久久久久久| 免费不卡的大黄色大毛片视频在线观看| 五月天丁香电影| 亚洲精品一二三| 亚洲熟女精品中文字幕| 少妇猛男粗大的猛烈进出视频 | 你懂的网址亚洲精品在线观看| 国产国拍精品亚洲av在线观看| 国产黄色视频一区二区在线观看| 久久亚洲国产成人精品v| 精品人妻视频免费看| 欧美成人a在线观看| 老女人水多毛片| 欧美精品人与动牲交sv欧美| av一本久久久久| 亚洲精品中文字幕在线视频 | 国产 一区 欧美 日韩| 边亲边吃奶的免费视频| 久久ye,这里只有精品| 香蕉精品网在线| 精品国产露脸久久av麻豆| 麻豆国产97在线/欧美| 免费看av在线观看网站| 三级经典国产精品| 啦啦啦在线观看免费高清www| 最近2019中文字幕mv第一页| 免费看光身美女| 亚州av有码| 国产精品一二三区在线看| 三级经典国产精品| 少妇的逼好多水| 国产v大片淫在线免费观看| 九九在线视频观看精品| 久久精品综合一区二区三区| 国产一级毛片在线| 老司机影院毛片| 久久久久久国产a免费观看| 晚上一个人看的免费电影| 免费观看在线日韩| 欧美性猛交╳xxx乱大交人| 久久精品国产亚洲av天美| 午夜福利网站1000一区二区三区| 黄色怎么调成土黄色| 国产精品麻豆人妻色哟哟久久| 国产伦精品一区二区三区视频9| 伊人久久国产一区二区| 日韩三级伦理在线观看| 一级毛片aaaaaa免费看小| 亚洲最大成人av| 99re6热这里在线精品视频| 欧美zozozo另类| 黄色一级大片看看| 丰满人妻一区二区三区视频av| 亚洲图色成人| 亚洲av男天堂| 久久精品国产亚洲av天美| 26uuu在线亚洲综合色| 欧美成人午夜免费资源| 97人妻精品一区二区三区麻豆| 91精品国产九色| 一区二区三区精品91| 国产精品国产三级国产专区5o| 汤姆久久久久久久影院中文字幕| 亚洲欧美中文字幕日韩二区| 男女啪啪激烈高潮av片| 亚洲四区av| 天堂中文最新版在线下载 | 久久人人爽人人片av| 中文天堂在线官网| 人妻少妇偷人精品九色| 国产熟女欧美一区二区| 国产成人a∨麻豆精品| 久久久久久久久久人人人人人人| 国产成人a∨麻豆精品| 色5月婷婷丁香| 激情 狠狠 欧美| 欧美日本视频|