• <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角點檢測算法
    熟妇人妻不卡中文字幕| h视频一区二区三区| 久久久久久人妻| 日韩视频在线欧美| a级一级毛片免费在线观看| 亚洲精华国产精华液的使用体验| 国产精品国产av在线观看| 国产伦在线观看视频一区| 国产在线视频一区二区| h日本视频在线播放| 国产免费福利视频在线观看| 在线观看av片永久免费下载| 欧美精品一区二区免费开放| www.色视频.com| 丰满乱子伦码专区| av国产精品久久久久影院| 美女内射精品一级片tv| 亚洲国产色片| 国产亚洲一区二区精品| 久久久久久久久久久丰满| 大又大粗又爽又黄少妇毛片口| 日韩成人伦理影院| 一个人看的www免费观看视频| 午夜福利在线观看免费完整高清在| 91午夜精品亚洲一区二区三区| 免费观看a级毛片全部| 一本—道久久a久久精品蜜桃钙片| 亚洲国产欧美人成| 国产av国产精品国产| 中文字幕精品免费在线观看视频 | 国产一区有黄有色的免费视频| 久久毛片免费看一区二区三区| www.av在线官网国产| 久久人妻熟女aⅴ| 久久99蜜桃精品久久| 亚洲美女视频黄频| 国产精品免费大片| 最近中文字幕高清免费大全6| 国产亚洲精品久久久com| 建设人人有责人人尽责人人享有的 | 插阴视频在线观看视频| 久久精品夜色国产| 国产爽快片一区二区三区| 夜夜爽夜夜爽视频| 国产精品99久久久久久久久| 最近中文字幕2019免费版| 久久精品人妻少妇| 免费播放大片免费观看视频在线观看| 亚洲国产欧美在线一区| 超碰av人人做人人爽久久| 丰满少妇做爰视频| 色视频www国产| 涩涩av久久男人的天堂| kizo精华| 国产精品久久久久久久电影| 国产精品国产三级专区第一集| 日韩制服骚丝袜av| 精品一区在线观看国产| 国产成人a∨麻豆精品| 黄色视频在线播放观看不卡| 久久久久久久久久久丰满| kizo精华| 全区人妻精品视频| 精华霜和精华液先用哪个| av免费在线看不卡| 国产精品精品国产色婷婷| 欧美成人一区二区免费高清观看| 少妇 在线观看| 亚洲av在线观看美女高潮| 欧美变态另类bdsm刘玥| 亚洲av二区三区四区| 男女啪啪激烈高潮av片| 亚洲伊人久久精品综合| 噜噜噜噜噜久久久久久91| 天天躁日日操中文字幕| 久久人人爽人人片av| av黄色大香蕉| 亚洲国产精品999| 精品熟女少妇av免费看| 国产精品偷伦视频观看了| 各种免费的搞黄视频| 视频中文字幕在线观看| 国产老妇伦熟女老妇高清| 国产91av在线免费观看| 国产精品国产av在线观看| 精品国产露脸久久av麻豆| 亚洲美女黄色视频免费看| 精品午夜福利在线看| 偷拍熟女少妇极品色| 亚洲国产色片| 蜜桃久久精品国产亚洲av| 人人妻人人看人人澡| h日本视频在线播放| 亚洲av国产av综合av卡| 亚洲av男天堂| 久久精品久久精品一区二区三区| 麻豆成人午夜福利视频| 欧美高清性xxxxhd video| 亚洲人成网站在线观看播放| 久久久久久久久久人人人人人人| 在线观看免费日韩欧美大片 | 久久久久久久大尺度免费视频| 黄色欧美视频在线观看| 国产欧美亚洲国产| 22中文网久久字幕| 久久久午夜欧美精品| 2022亚洲国产成人精品| 永久免费av网站大全| 边亲边吃奶的免费视频| 国产真实伦视频高清在线观看| 精品一区二区免费观看| xxx大片免费视频| 日韩中文字幕视频在线看片 | 国产中年淑女户外野战色| 丝瓜视频免费看黄片| 亚洲丝袜综合中文字幕| 久久av网站| 日本免费在线观看一区| 国产真实伦视频高清在线观看| 日本猛色少妇xxxxx猛交久久| 青青草视频在线视频观看| 国产在视频线精品| 国产人妻一区二区三区在| 啦啦啦啦在线视频资源| 最近最新中文字幕大全电影3| 97超碰精品成人国产| 久久久久久久久久久免费av| 亚洲一区二区三区欧美精品| 七月丁香在线播放| 观看美女的网站| 中文乱码字字幕精品一区二区三区| 最近的中文字幕免费完整| 纯流量卡能插随身wifi吗| 亚洲图色成人| 一级爰片在线观看| 国产成人a区在线观看| 一级毛片aaaaaa免费看小| 人妻制服诱惑在线中文字幕| 亚洲电影在线观看av| 一区二区三区精品91| 日本vs欧美在线观看视频 | 国产欧美另类精品又又久久亚洲欧美| 一级毛片久久久久久久久女| 王馨瑶露胸无遮挡在线观看| 水蜜桃什么品种好| 晚上一个人看的免费电影| 亚洲丝袜综合中文字幕| 黄色一级大片看看| 国产女主播在线喷水免费视频网站| 看非洲黑人一级黄片| 成人国产av品久久久| 国产精品三级大全| 成人漫画全彩无遮挡| 亚洲综合精品二区| 男人狂女人下面高潮的视频| 成人美女网站在线观看视频| 亚洲中文av在线| 国产极品天堂在线| 国产伦精品一区二区三区视频9| 日韩成人av中文字幕在线观看| 亚洲人与动物交配视频| 黄色欧美视频在线观看| 你懂的网址亚洲精品在线观看| 九九久久精品国产亚洲av麻豆| 免费观看性生交大片5| 亚洲国产最新在线播放| 国产精品久久久久久精品电影小说 | 亚洲精华国产精华液的使用体验| 天堂中文最新版在线下载| 国产永久视频网站| 中国三级夫妇交换| 久久韩国三级中文字幕| 狂野欧美激情性bbbbbb| 亚洲精品日韩av片在线观看| 少妇的逼好多水| videossex国产| 欧美日韩亚洲高清精品| 国产av一区二区精品久久 | 美女国产视频在线观看| 美女福利国产在线 | 色5月婷婷丁香| 亚洲精品视频女| 精品久久国产蜜桃| 国产人妻一区二区三区在| 国产色婷婷99| 天堂中文最新版在线下载| 18禁裸乳无遮挡动漫免费视频| 精品久久久精品久久久| 国产日韩欧美亚洲二区| 波野结衣二区三区在线| 国内揄拍国产精品人妻在线| 久久 成人 亚洲| 如何舔出高潮| 精品酒店卫生间| 3wmmmm亚洲av在线观看| 日本爱情动作片www.在线观看| 三级国产精品片| 蜜桃在线观看..| 只有这里有精品99| 国产视频首页在线观看| 久久国内精品自在自线图片| av黄色大香蕉| 欧美日韩一区二区视频在线观看视频在线| 少妇被粗大猛烈的视频| 国产在线视频一区二区| 国产av码专区亚洲av| 一级毛片 在线播放| 少妇猛男粗大的猛烈进出视频| 免费看日本二区| 在线精品无人区一区二区三 | 一级av片app| 麻豆成人av视频| 日本午夜av视频| 天堂8中文在线网| a级毛色黄片| 国产精品不卡视频一区二区| 两个人的视频大全免费| 国产成人精品婷婷| 在现免费观看毛片| 免费人成在线观看视频色| 欧美最新免费一区二区三区| 在线天堂最新版资源| 日本vs欧美在线观看视频 | 噜噜噜噜噜久久久久久91| 日韩一区二区三区影片| 国产精品99久久久久久久久| 五月伊人婷婷丁香| 亚洲欧美日韩卡通动漫| 亚洲精品国产av蜜桃| 一个人看的www免费观看视频| 亚洲综合色惰| 18禁裸乳无遮挡免费网站照片| 国产色婷婷99| 在线观看av片永久免费下载| 国产精品爽爽va在线观看网站| 99热6这里只有精品| 中国美白少妇内射xxxbb| 天天躁夜夜躁狠狠久久av| 97在线视频观看| 亚洲欧美日韩另类电影网站 | 视频区图区小说| 91九色精品人成在线观看| 乱人伦中国视频| 精品国产乱码久久久久久男人| 精品欧美一区二区三区在线| 一边摸一边抽搐一进一出视频| 国产av一区二区精品久久| 一区二区三区精品91| 一区二区三区四区激情视频| 91字幕亚洲| 精品国产一区二区三区久久久樱花| 极品人妻少妇av视频| 一本大道久久a久久精品| 男人添女人高潮全过程视频| 又大又爽又粗| 91老司机精品| 自线自在国产av| 肉色欧美久久久久久久蜜桃| 日韩av免费高清视频| 国产视频一区二区在线看| 深夜精品福利| 国产国语露脸激情在线看| 亚洲精品国产av成人精品| 日韩精品免费视频一区二区三区| 男女边摸边吃奶| 久久久久久久精品精品| 国产在线一区二区三区精| 久久人妻福利社区极品人妻图片 | 亚洲精品久久成人aⅴ小说| 一级片免费观看大全| 人人澡人人妻人| 黑人巨大精品欧美一区二区蜜桃| 国产亚洲精品久久久久5区| √禁漫天堂资源中文www| 女人被躁到高潮嗷嗷叫费观| 精品人妻熟女毛片av久久网站| av福利片在线| 国产成人a∨麻豆精品| 涩涩av久久男人的天堂| 最近最新中文字幕大全免费视频 | 赤兔流量卡办理| 99精品久久久久人妻精品| 99久久综合免费| 国产深夜福利视频在线观看| 看十八女毛片水多多多| 亚洲av男天堂| 曰老女人黄片| 欧美老熟妇乱子伦牲交| 久9热在线精品视频| 精品少妇内射三级| 欧美精品高潮呻吟av久久| 国产女主播在线喷水免费视频网站| 两个人看的免费小视频| 999久久久国产精品视频| 黄频高清免费视频| 亚洲欧洲精品一区二区精品久久久| 亚洲精品中文字幕在线视频| 看十八女毛片水多多多| 美女扒开内裤让男人捅视频| 国产亚洲一区二区精品| 夫妻性生交免费视频一级片| 亚洲精品乱久久久久久| av天堂在线播放| 亚洲人成电影免费在线| 午夜视频精品福利| 免费高清在线观看视频在线观看| 国产视频一区二区在线看| 永久免费av网站大全| 又黄又粗又硬又大视频| 性色av乱码一区二区三区2| 国产成人a∨麻豆精品| 精品少妇黑人巨大在线播放| 欧美日韩视频高清一区二区三区二| 国产精品久久久人人做人人爽| 色综合欧美亚洲国产小说| 新久久久久国产一级毛片| 国产极品粉嫩免费观看在线| 欧美日韩亚洲高清精品| 亚洲人成电影免费在线| 国产爽快片一区二区三区| 后天国语完整版免费观看| 久久人人爽人人片av| 伊人亚洲综合成人网| 热99国产精品久久久久久7| 日韩熟女老妇一区二区性免费视频| 黄片播放在线免费| 十八禁网站网址无遮挡| 人人妻人人添人人爽欧美一区卜| 一二三四社区在线视频社区8| 亚洲精品国产色婷婷电影| 青春草视频在线免费观看| 国产一区亚洲一区在线观看| 久久国产精品男人的天堂亚洲| 十分钟在线观看高清视频www| 韩国精品一区二区三区| 亚洲人成电影免费在线| 好男人电影高清在线观看| 99精品久久久久人妻精品| 超碰97精品在线观看| 超色免费av| 国产视频一区二区在线看| 色精品久久人妻99蜜桃| 久热爱精品视频在线9| av网站免费在线观看视频| 亚洲精品日本国产第一区| 51午夜福利影视在线观看| 夜夜骑夜夜射夜夜干| 亚洲人成77777在线视频| 精品高清国产在线一区| 色网站视频免费| 王馨瑶露胸无遮挡在线观看| 久久中文字幕一级| 亚洲自偷自拍图片 自拍| 最近中文字幕2019免费版| 精品卡一卡二卡四卡免费| 91九色精品人成在线观看| 精品久久久久久久毛片微露脸 | 国产老妇伦熟女老妇高清| 亚洲七黄色美女视频| 亚洲精品国产av成人精品| 精品国产一区二区三区久久久樱花| 好男人视频免费观看在线| 久久女婷五月综合色啪小说| 国语对白做爰xxxⅹ性视频网站| 国产av精品麻豆| 人妻人人澡人人爽人人| 大话2 男鬼变身卡| 亚洲国产精品成人久久小说| 中国美女看黄片| 又紧又爽又黄一区二区| 亚洲五月婷婷丁香| 欧美成人精品欧美一级黄| 久久九九热精品免费| 激情视频va一区二区三区| 免费高清在线观看视频在线观看| 蜜桃国产av成人99| 成年人黄色毛片网站| 丰满饥渴人妻一区二区三| 国产色视频综合| 无限看片的www在线观看| 精品视频人人做人人爽| 国产精品 欧美亚洲| 久久久久久久大尺度免费视频| 国产精品久久久久久人妻精品电影 | 亚洲国产av影院在线观看| 人人妻,人人澡人人爽秒播 | 老汉色∧v一级毛片| 欧美成狂野欧美在线观看| 在线观看www视频免费| 欧美日韩视频精品一区| 老司机靠b影院| 国产熟女午夜一区二区三区| 伊人久久大香线蕉亚洲五| 国产麻豆69| 极品少妇高潮喷水抽搐| 国产男人的电影天堂91| 国产av一区二区精品久久| av天堂久久9| 欧美黑人欧美精品刺激| 免费一级毛片在线播放高清视频 | www.999成人在线观看| 久久精品久久精品一区二区三区| 国产亚洲欧美在线一区二区| 99九九在线精品视频| 在线天堂中文资源库| 国产片内射在线| 新久久久久国产一级毛片| 亚洲色图 男人天堂 中文字幕| 亚洲午夜精品一区,二区,三区| 国产成人一区二区在线| 亚洲成国产人片在线观看| 侵犯人妻中文字幕一二三四区| 国产三级黄色录像| 国产亚洲欧美精品永久| www.自偷自拍.com| 欧美日韩成人在线一区二区| 免费在线观看视频国产中文字幕亚洲 | 亚洲精品国产区一区二| 亚洲中文av在线| 欧美国产精品一级二级三级| 精品亚洲成国产av| 操出白浆在线播放| 亚洲欧美激情在线| 99国产精品一区二区蜜桃av | 久久人人97超碰香蕉20202| 1024视频免费在线观看| 久久精品亚洲熟妇少妇任你| 午夜福利在线免费观看网站| 高潮久久久久久久久久久不卡| 青春草视频在线免费观看| 操美女的视频在线观看| 中文字幕亚洲精品专区| 这个男人来自地球电影免费观看| 在线 av 中文字幕| 日本91视频免费播放| 欧美人与性动交α欧美精品济南到| 亚洲精品一卡2卡三卡4卡5卡 | 国产一区二区三区综合在线观看| 久久久久精品国产欧美久久久 | 亚洲五月色婷婷综合| 最近最新中文字幕大全免费视频 | 狠狠婷婷综合久久久久久88av| 美女福利国产在线| 国产一区二区三区综合在线观看| 中文字幕最新亚洲高清| 午夜福利乱码中文字幕| 在线av久久热| 亚洲少妇的诱惑av| 18禁黄网站禁片午夜丰满| 美女视频免费永久观看网站| 亚洲人成电影观看| 黄色视频在线播放观看不卡| 亚洲欧美精品自产自拍| 王馨瑶露胸无遮挡在线观看| 欧美激情极品国产一区二区三区| 两个人看的免费小视频| 国产又色又爽无遮挡免| 午夜久久久在线观看| 欧美在线黄色| 首页视频小说图片口味搜索 | av片东京热男人的天堂| 妹子高潮喷水视频| 日韩av在线免费看完整版不卡| 在线 av 中文字幕| 伊人亚洲综合成人网| 777久久人妻少妇嫩草av网站| 色网站视频免费| 亚洲成人手机| 成人免费观看视频高清| 久久99热这里只频精品6学生| 亚洲国产精品999| 伦理电影免费视频| 午夜免费观看性视频| 成人三级做爰电影| 亚洲第一av免费看| 少妇的丰满在线观看| 精品人妻1区二区| 菩萨蛮人人尽说江南好唐韦庄| 成年av动漫网址| 69精品国产乱码久久久| 蜜桃国产av成人99| 成年人黄色毛片网站| 国产亚洲一区二区精品| 汤姆久久久久久久影院中文字幕| 久热爱精品视频在线9| 亚洲国产欧美一区二区综合| 曰老女人黄片| 伊人亚洲综合成人网| 欧美日韩av久久| 男的添女的下面高潮视频| 国产精品久久久久久精品电影小说| 亚洲欧美一区二区三区黑人| 夜夜骑夜夜射夜夜干| 国产精品av久久久久免费| 一级毛片我不卡| 叶爱在线成人免费视频播放| 精品国产国语对白av| 成在线人永久免费视频| 亚洲国产欧美网| 国产高清不卡午夜福利| 中文字幕高清在线视频| 女人久久www免费人成看片| 久久国产亚洲av麻豆专区| 国产三级黄色录像| 丝袜美足系列| 天天添夜夜摸| 日韩免费高清中文字幕av| 精品卡一卡二卡四卡免费| 成人免费观看视频高清| 国产成人免费无遮挡视频| 成年人免费黄色播放视频| 日日夜夜操网爽| 大香蕉久久网| 美女福利国产在线| 老司机影院毛片| 一级毛片我不卡| 男女之事视频高清在线观看 | 亚洲欧美色中文字幕在线| 中文字幕色久视频| 久久国产精品男人的天堂亚洲| 大码成人一级视频| 在线观看一区二区三区激情| 久久久久视频综合| 少妇精品久久久久久久| 91麻豆av在线| 国产主播在线观看一区二区 | 久久毛片免费看一区二区三区| av一本久久久久| 高清视频免费观看一区二区| 国产99久久九九免费精品| 成人免费观看视频高清| netflix在线观看网站| 中文字幕人妻丝袜制服| 免费日韩欧美在线观看| 欧美在线一区亚洲| 又紧又爽又黄一区二区| 成人手机av| 亚洲国产成人一精品久久久| 欧美精品一区二区免费开放| 精品一区二区三区四区五区乱码 | av又黄又爽大尺度在线免费看| 波多野结衣av一区二区av| 中文字幕色久视频| 国产精品香港三级国产av潘金莲 | 久久精品aⅴ一区二区三区四区| 日本猛色少妇xxxxx猛交久久| 国产欧美日韩一区二区三 | 国产一区亚洲一区在线观看| 男女无遮挡免费网站观看| 高清av免费在线| 国产午夜精品一二区理论片| 久久热在线av| 九草在线视频观看| 99九九在线精品视频| 性色av乱码一区二区三区2| 美女国产高潮福利片在线看| 伊人亚洲综合成人网| 亚洲中文字幕日韩| 69精品国产乱码久久久| 久久精品国产a三级三级三级| 久久久久国产一级毛片高清牌| 精品一区二区三卡| 18禁国产床啪视频网站| 黄色a级毛片大全视频| 亚洲av日韩精品久久久久久密 | 国产免费又黄又爽又色| 久久久精品94久久精品| 国产成人一区二区在线| 久久 成人 亚洲| 性高湖久久久久久久久免费观看| 王馨瑶露胸无遮挡在线观看| 51午夜福利影视在线观看| 侵犯人妻中文字幕一二三四区| 亚洲五月婷婷丁香| 久久久久久人人人人人| 操美女的视频在线观看| 日韩 欧美 亚洲 中文字幕| 国产人伦9x9x在线观看| 99国产精品免费福利视频| 国产成人免费无遮挡视频| 香蕉丝袜av| 欧美日韩av久久| 国产色视频综合| 久久精品国产a三级三级三级| 国产精品人妻久久久影院| 成年动漫av网址| 久久久久久久久免费视频了| 看十八女毛片水多多多| 久久精品亚洲熟妇少妇任你| 亚洲av国产av综合av卡| 国产爽快片一区二区三区| 婷婷丁香在线五月| 国产91精品成人一区二区三区 | 91国产中文字幕| 中文字幕制服av| 日本vs欧美在线观看视频| 久久国产精品男人的天堂亚洲| 久久ye,这里只有精品| 视频区欧美日本亚洲| 成年人午夜在线观看视频| 欧美亚洲 丝袜 人妻 在线| 在线亚洲精品国产二区图片欧美| 首页视频小说图片口味搜索 | 黄片播放在线免费| 国产高清不卡午夜福利| 午夜激情av网站| 最黄视频免费看| 免费少妇av软件| 国产成人av激情在线播放| 美女脱内裤让男人舔精品视频| 久热爱精品视频在线9| 女警被强在线播放| 免费少妇av软件| 久久久久久久久免费视频了| 99国产综合亚洲精品|