• <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角點檢測算法
    亚洲激情在线av| 精品免费久久久久久久清纯| 久久国产亚洲av麻豆专区| 满18在线观看网站| 午夜福利高清视频| 少妇熟女aⅴ在线视频| 一级毛片精品| 欧美一级毛片孕妇| 国产欧美日韩综合在线一区二区| 久久久久久久久久久久大奶| 中文字幕另类日韩欧美亚洲嫩草| 午夜福利免费观看在线| 亚洲国产看品久久| 午夜视频精品福利| 午夜激情av网站| 又黄又粗又硬又大视频| 免费一级毛片在线播放高清视频 | 9191精品国产免费久久| 亚洲第一电影网av| 国内精品久久久久久久电影| 麻豆av在线久日| 久久精品国产亚洲av高清一级| 在线天堂中文资源库| 日本 av在线| 搡老岳熟女国产| 国产精品综合久久久久久久免费 | 久久久久久免费高清国产稀缺| 麻豆成人av在线观看| 韩国av一区二区三区四区| 午夜久久久久精精品| 性少妇av在线| 啪啪无遮挡十八禁网站| 日本a在线网址| 岛国视频午夜一区免费看| 国产成人系列免费观看| 麻豆久久精品国产亚洲av| 精品高清国产在线一区| 免费观看精品视频网站| 91大片在线观看| 国产aⅴ精品一区二区三区波| 精品日产1卡2卡| 成人国语在线视频| 婷婷丁香在线五月| 久久久久国产一级毛片高清牌| 婷婷丁香在线五月| www.精华液| 天堂影院成人在线观看| 波多野结衣一区麻豆| av视频在线观看入口| 亚洲七黄色美女视频| 免费女性裸体啪啪无遮挡网站| 19禁男女啪啪无遮挡网站| 久久久水蜜桃国产精品网| 国产一区二区在线av高清观看| videosex国产| 亚洲精品国产精品久久久不卡| 日本精品一区二区三区蜜桃| 丁香欧美五月| 涩涩av久久男人的天堂| 亚洲精品一区av在线观看| 精品人妻1区二区| 久久久国产欧美日韩av| 亚洲国产精品sss在线观看| 久久人人97超碰香蕉20202| 亚洲精品国产区一区二| 激情视频va一区二区三区| 精品少妇一区二区三区视频日本电影| 麻豆成人av在线观看| 女同久久另类99精品国产91| 国语自产精品视频在线第100页| 侵犯人妻中文字幕一二三四区| 国产精品一区二区精品视频观看| 亚洲欧美日韩另类电影网站| 身体一侧抽搐| 久久欧美精品欧美久久欧美| 午夜免费成人在线视频| 人人妻人人澡欧美一区二区 | 亚洲精品国产区一区二| 亚洲 欧美一区二区三区| 欧美久久黑人一区二区| 亚洲欧美日韩高清在线视频| 成人三级做爰电影| 久久精品国产99精品国产亚洲性色 | 欧美 亚洲 国产 日韩一| 久久 成人 亚洲| 国产aⅴ精品一区二区三区波| xxx96com| 国产精品电影一区二区三区| 久久国产精品人妻蜜桃| 日本在线视频免费播放| www国产在线视频色| 欧美一级毛片孕妇| 欧美一级a爱片免费观看看 | 久久国产精品影院| 一本大道久久a久久精品| 色婷婷久久久亚洲欧美| 欧美国产精品va在线观看不卡| 91精品国产国语对白视频| 精品国内亚洲2022精品成人| 免费久久久久久久精品成人欧美视频| 女人精品久久久久毛片| 最新在线观看一区二区三区| 日本vs欧美在线观看视频| 性少妇av在线| 99久久综合精品五月天人人| 欧美黄色淫秽网站| 成年女人毛片免费观看观看9| 18禁黄网站禁片午夜丰满| 亚洲成av片中文字幕在线观看| 天天躁夜夜躁狠狠躁躁| 国产精品久久久久久亚洲av鲁大| 国产不卡一卡二| 在线永久观看黄色视频| 亚洲片人在线观看| 国产精品二区激情视频| 亚洲av日韩精品久久久久久密| 日韩av在线大香蕉| 亚洲五月天丁香| 亚洲成av人片免费观看| 国产极品粉嫩免费观看在线| 免费不卡黄色视频| 精品日产1卡2卡| 久久影院123| 免费在线观看日本一区| 真人做人爱边吃奶动态| 亚洲一区二区三区不卡视频| 一边摸一边抽搐一进一出视频| 精品国产亚洲在线| 午夜福利18| 男女床上黄色一级片免费看| av有码第一页| 免费观看人在逋| 一区二区三区精品91| 国产成人欧美在线观看| 久久久久精品国产欧美久久久| 国产成人av激情在线播放| 美国免费a级毛片| 亚洲 欧美 日韩 在线 免费| 99国产极品粉嫩在线观看| 一边摸一边抽搐一进一出视频| 999久久久精品免费观看国产| 好男人电影高清在线观看| 97人妻天天添夜夜摸| 动漫黄色视频在线观看| 亚洲国产日韩欧美精品在线观看 | 免费人成视频x8x8入口观看| 国产精华一区二区三区| 99久久国产精品久久久| 老司机午夜福利在线观看视频| 91老司机精品| 亚洲国产日韩欧美精品在线观看 | 90打野战视频偷拍视频| 欧美日韩乱码在线| 日韩av在线大香蕉| 香蕉国产在线看| 十分钟在线观看高清视频www| 法律面前人人平等表现在哪些方面| 热re99久久国产66热| 色婷婷久久久亚洲欧美| 欧美亚洲日本最大视频资源| 久久国产乱子伦精品免费另类| 两个人免费观看高清视频| 国产精品爽爽va在线观看网站 | 曰老女人黄片| 午夜精品国产一区二区电影| 亚洲国产看品久久| 亚洲精品美女久久av网站| 18禁美女被吸乳视频| 亚洲人成电影观看| 夜夜夜夜夜久久久久| 国产成人啪精品午夜网站| 色播在线永久视频| 精品少妇一区二区三区视频日本电影| 欧美+亚洲+日韩+国产| 在线观看免费视频网站a站| 美女高潮到喷水免费观看| 国产成人精品无人区| 可以在线观看的亚洲视频| 中出人妻视频一区二区| 亚洲欧美日韩高清在线视频| 少妇熟女aⅴ在线视频| 国产视频一区二区在线看| 人人妻人人澡欧美一区二区 | 极品人妻少妇av视频| 国产精品亚洲av一区麻豆| www.精华液| 在线国产一区二区在线| 久热爱精品视频在线9| 亚洲欧美日韩无卡精品| 美女高潮喷水抽搐中文字幕| 99国产综合亚洲精品| 如日韩欧美国产精品一区二区三区| 法律面前人人平等表现在哪些方面| 看黄色毛片网站| 亚洲人成网站在线播放欧美日韩| 嫩草影视91久久| 午夜免费观看网址| 欧美激情高清一区二区三区| 亚洲成av人片免费观看| 日本在线视频免费播放| 真人做人爱边吃奶动态| 精品不卡国产一区二区三区| 国产av在哪里看| 日韩欧美三级三区| 久久久久久久精品吃奶| 极品人妻少妇av视频| 黄片播放在线免费| 在线观看66精品国产| 最好的美女福利视频网| 男人舔女人下体高潮全视频| 国产精品九九99| 1024视频免费在线观看| 亚洲av电影不卡..在线观看| 亚洲精品av麻豆狂野| 在线观看免费视频日本深夜| 男女午夜视频在线观看| www.999成人在线观看| 两个人视频免费观看高清| 午夜老司机福利片| 午夜久久久久精精品| 国产精品免费视频内射| 天天躁狠狠躁夜夜躁狠狠躁| 欧美激情极品国产一区二区三区| 亚洲精品久久国产高清桃花| 精品欧美一区二区三区在线| 丝袜美足系列| 国产一级毛片七仙女欲春2 | 日本欧美视频一区| 精品乱码久久久久久99久播| 久久人妻熟女aⅴ| 一进一出抽搐动态| 50天的宝宝边吃奶边哭怎么回事| 一级a爱视频在线免费观看| 国产xxxxx性猛交| 久久狼人影院| 天天一区二区日本电影三级 | √禁漫天堂资源中文www| 久久精品亚洲熟妇少妇任你| 亚洲美女黄片视频| 亚洲色图av天堂| 好男人电影高清在线观看| 天堂动漫精品| 日韩大码丰满熟妇| 久久久久国产一级毛片高清牌| 免费一级毛片在线播放高清视频 | 欧美在线一区亚洲| 精品无人区乱码1区二区| 久久人人爽av亚洲精品天堂| 久久国产精品影院| 国产成人系列免费观看| 国产亚洲欧美在线一区二区| 日本 欧美在线| 欧美日韩亚洲国产一区二区在线观看| 午夜日韩欧美国产| 一区二区三区激情视频| 久久国产精品男人的天堂亚洲| 欧美日韩亚洲综合一区二区三区_| 一区在线观看完整版| 丰满人妻熟妇乱又伦精品不卡| 中文字幕精品免费在线观看视频| 好男人电影高清在线观看| 黄片小视频在线播放| 一边摸一边抽搐一进一小说| 可以免费在线观看a视频的电影网站| 怎么达到女性高潮| 亚洲 欧美 日韩 在线 免费| 一区二区日韩欧美中文字幕| 搡老妇女老女人老熟妇| 老司机午夜十八禁免费视频| 我的亚洲天堂| 不卡av一区二区三区| 12—13女人毛片做爰片一| 欧美午夜高清在线| 高清黄色对白视频在线免费看| 一级作爱视频免费观看| 欧美大码av| 久热爱精品视频在线9| 久热这里只有精品99| 久久精品国产亚洲av高清一级| 久久这里只有精品19| 99精品在免费线老司机午夜| 欧美日韩黄片免| 日韩欧美三级三区| 国产色视频综合| 夜夜躁狠狠躁天天躁| 一区二区三区激情视频| 午夜免费观看网址| 久久亚洲真实| 中文字幕色久视频| 亚洲自偷自拍图片 自拍| 波多野结衣高清无吗| 十八禁网站免费在线| 高清毛片免费观看视频网站| 欧美国产日韩亚洲一区| 精品福利观看| 亚洲精品久久成人aⅴ小说| 国产亚洲精品久久久久5区| 国产精品久久电影中文字幕| 免费在线观看视频国产中文字幕亚洲| 自线自在国产av| 欧美激情极品国产一区二区三区| tocl精华| 高潮久久久久久久久久久不卡| 国产精品亚洲一级av第二区| av视频免费观看在线观看| 亚洲国产欧美一区二区综合| 婷婷丁香在线五月| 成人永久免费在线观看视频| 免费看a级黄色片| 免费观看人在逋| 欧美色欧美亚洲另类二区 | 亚洲男人天堂网一区| 一本久久中文字幕| 国产成人啪精品午夜网站| 国产精品98久久久久久宅男小说| 少妇被粗大的猛进出69影院| 久久亚洲真实| 正在播放国产对白刺激| 国产精品亚洲一级av第二区| 他把我摸到了高潮在线观看| a级毛片在线看网站| 午夜福利成人在线免费观看| 手机成人av网站| avwww免费| 母亲3免费完整高清在线观看| 国产精品 欧美亚洲| 777久久人妻少妇嫩草av网站| 日本欧美视频一区| 亚洲九九香蕉| 性欧美人与动物交配| netflix在线观看网站| 国产三级在线视频| 久久精品国产亚洲av高清一级| 老司机福利观看| 岛国在线观看网站| 99re在线观看精品视频| 亚洲aⅴ乱码一区二区在线播放 | 欧美+亚洲+日韩+国产| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲久久久国产精品| 三级毛片av免费| 少妇熟女aⅴ在线视频| 久久精品国产亚洲av香蕉五月| 亚洲五月婷婷丁香| 中文字幕av电影在线播放| 精品电影一区二区在线| 琪琪午夜伦伦电影理论片6080| 国产极品粉嫩免费观看在线| av电影中文网址| 国产不卡一卡二| 给我免费播放毛片高清在线观看| 黄色毛片三级朝国网站| 午夜精品在线福利| 人妻久久中文字幕网| 欧美乱码精品一区二区三区| 日韩精品免费视频一区二区三区| 一区二区三区精品91| 日韩国内少妇激情av| 亚洲,欧美精品.| 亚洲成a人片在线一区二区| 久热爱精品视频在线9| 99精品欧美一区二区三区四区| 久久久国产成人精品二区| 亚洲精品在线美女| а√天堂www在线а√下载| 精品国内亚洲2022精品成人| 国产真人三级小视频在线观看| 好男人在线观看高清免费视频 | 国产精品98久久久久久宅男小说| videosex国产| 久久久久精品国产欧美久久久| 成人三级黄色视频| 无人区码免费观看不卡| 免费无遮挡裸体视频| 中文字幕人妻丝袜一区二区| 午夜久久久在线观看| 97人妻精品一区二区三区麻豆 | av天堂久久9| 国产精品免费一区二区三区在线| 一级毛片女人18水好多| 久久久久久大精品| 又黄又粗又硬又大视频| av网站免费在线观看视频| 亚洲av五月六月丁香网| 亚洲人成电影观看| 每晚都被弄得嗷嗷叫到高潮| 在线观看日韩欧美| 日本五十路高清| 女人被狂操c到高潮| 搡老岳熟女国产| 校园春色视频在线观看| 亚洲激情在线av| 级片在线观看| 天堂√8在线中文| 国产1区2区3区精品| 国产精品爽爽va在线观看网站 | 一级黄色大片毛片| 如日韩欧美国产精品一区二区三区| 久久婷婷成人综合色麻豆| 99香蕉大伊视频| 亚洲 国产 在线| 97超级碰碰碰精品色视频在线观看| 亚洲少妇的诱惑av| 熟女少妇亚洲综合色aaa.| 久久精品国产亚洲av香蕉五月| 亚洲精品国产色婷婷电影| 亚洲视频免费观看视频| 欧美最黄视频在线播放免费| 大陆偷拍与自拍| 亚洲欧美精品综合一区二区三区| 免费高清视频大片| 在线国产一区二区在线| 亚洲七黄色美女视频| 国产精品99久久99久久久不卡| 丝袜美腿诱惑在线| 九色国产91popny在线| 国产精华一区二区三区| www国产在线视频色| 欧美成人午夜精品| 色综合站精品国产| 亚洲在线自拍视频| 久久人人精品亚洲av| 色综合站精品国产| 国产精品一区二区在线不卡| 亚洲国产毛片av蜜桃av| 久久香蕉激情| av免费在线观看网站| 天堂√8在线中文| 多毛熟女@视频| 黑人巨大精品欧美一区二区蜜桃| 国产精华一区二区三区| 97人妻天天添夜夜摸| 少妇裸体淫交视频免费看高清 | 午夜免费成人在线视频| 麻豆国产av国片精品| 国产视频一区二区在线看| 免费搜索国产男女视频| 欧美乱色亚洲激情| 淫秽高清视频在线观看| 看免费av毛片| 午夜福利,免费看| 欧美国产日韩亚洲一区| 啦啦啦免费观看视频1| 日本撒尿小便嘘嘘汇集6| 日本免费一区二区三区高清不卡 | 免费看a级黄色片| 午夜老司机福利片| 搡老岳熟女国产| 国产片内射在线| 久久久水蜜桃国产精品网| 欧美日本视频| 亚洲国产精品成人综合色| 亚洲中文字幕日韩| 女警被强在线播放| 亚洲国产看品久久| 我的亚洲天堂| 国内毛片毛片毛片毛片毛片| 国产精品99久久99久久久不卡| 成年人黄色毛片网站| 91字幕亚洲| 久久草成人影院| 亚洲激情在线av| 国产精品野战在线观看| 亚洲欧美一区二区三区黑人| 亚洲 欧美 日韩 在线 免费| 老司机午夜十八禁免费视频| 国产高清激情床上av| 欧美一级a爱片免费观看看 | 久久久国产精品麻豆| 天天一区二区日本电影三级 | 久久欧美精品欧美久久欧美| 欧美日韩亚洲国产一区二区在线观看| 日本精品一区二区三区蜜桃| 美国免费a级毛片| 一级,二级,三级黄色视频| 每晚都被弄得嗷嗷叫到高潮| 老司机午夜福利在线观看视频| 精品国产美女av久久久久小说| 国内精品久久久久精免费| 久久中文字幕人妻熟女| 黑人巨大精品欧美一区二区mp4| 99精品久久久久人妻精品| 久久国产乱子伦精品免费另类| 在线播放国产精品三级| 国产精品九九99| 在线十欧美十亚洲十日本专区| 高清毛片免费观看视频网站| 99国产极品粉嫩在线观看| 国产激情久久老熟女| 亚洲av成人av| 亚洲va日本ⅴa欧美va伊人久久| 一本久久中文字幕| 精品人妻1区二区| svipshipincom国产片| 19禁男女啪啪无遮挡网站| 欧美性长视频在线观看| 国产精华一区二区三区| 亚洲伊人色综图| 日本三级黄在线观看| 亚洲精品粉嫩美女一区| 国产成+人综合+亚洲专区| 最近最新免费中文字幕在线| 久久人人97超碰香蕉20202| 亚洲avbb在线观看| 亚洲精品美女久久av网站| 日韩国内少妇激情av| 给我免费播放毛片高清在线观看| 色在线成人网| 91大片在线观看| 黑丝袜美女国产一区| 一级a爱视频在线免费观看| 可以在线观看的亚洲视频| 国产又色又爽无遮挡免费看| 99精品在免费线老司机午夜| 国产成人精品久久二区二区91| 999精品在线视频| 成年女人毛片免费观看观看9| 在线免费观看的www视频| 精品久久久精品久久久| 亚洲五月婷婷丁香| 亚洲欧美日韩无卡精品| 夜夜躁狠狠躁天天躁| 黄频高清免费视频| 色哟哟哟哟哟哟| 中出人妻视频一区二区| 天堂√8在线中文| www日本在线高清视频| 两性夫妻黄色片| 亚洲国产毛片av蜜桃av| 亚洲欧美日韩高清在线视频| 国产高清有码在线观看视频 | 一个人免费在线观看的高清视频| av免费在线观看网站| 国产麻豆成人av免费视频| 一区二区日韩欧美中文字幕| 丝袜美足系列| 丝袜人妻中文字幕| 91成人精品电影| 国产在线精品亚洲第一网站| 欧美日韩亚洲综合一区二区三区_| 亚洲熟女毛片儿| 亚洲一区高清亚洲精品| 亚洲成av片中文字幕在线观看| 超碰成人久久| 在线观看日韩欧美| 精品人妻在线不人妻| 亚洲中文日韩欧美视频| 又大又爽又粗| 无遮挡黄片免费观看| 黄色a级毛片大全视频| 精品免费久久久久久久清纯| 女人被躁到高潮嗷嗷叫费观| 法律面前人人平等表现在哪些方面| 色综合站精品国产| 在线观看免费视频日本深夜| 十分钟在线观看高清视频www| 国产成年人精品一区二区| 亚洲va日本ⅴa欧美va伊人久久| 一区在线观看完整版| 日韩精品免费视频一区二区三区| 男男h啪啪无遮挡| 日韩大码丰满熟妇| 村上凉子中文字幕在线| 高潮久久久久久久久久久不卡| 亚洲av熟女| 999久久久国产精品视频| 亚洲第一电影网av| 免费看美女性在线毛片视频| 久久婷婷成人综合色麻豆| 亚洲精品国产一区二区精华液| 国产一区二区在线av高清观看| 亚洲欧美精品综合一区二区三区| av片东京热男人的天堂| 亚洲久久久国产精品| 曰老女人黄片| aaaaa片日本免费| 99精品欧美一区二区三区四区| 岛国在线观看网站| 美女 人体艺术 gogo| 久久中文看片网| 亚洲最大成人中文| 人人妻人人澡人人看| 亚洲视频免费观看视频| 99国产精品一区二区三区| 久99久视频精品免费| 国产精品久久久久久精品电影 | 国产精品久久久av美女十八| 亚洲精品av麻豆狂野| 这个男人来自地球电影免费观看| 亚洲国产毛片av蜜桃av| 免费一级毛片在线播放高清视频 | 国产精品久久久久久亚洲av鲁大| 精品第一国产精品| 男女下面插进去视频免费观看| 日韩欧美国产一区二区入口| 亚洲情色 制服丝袜| cao死你这个sao货| 日本五十路高清| 午夜久久久久精精品| 国产亚洲欧美在线一区二区| 欧美乱色亚洲激情| 99riav亚洲国产免费| 啦啦啦观看免费观看视频高清 | 亚洲人成电影观看| 亚洲国产欧美日韩在线播放| 国产区一区二久久| 日韩 欧美 亚洲 中文字幕| 久久人人爽av亚洲精品天堂| 成人亚洲精品一区在线观看| 熟妇人妻久久中文字幕3abv| 久久九九热精品免费| 免费看a级黄色片| 久久人妻av系列| 国产成人精品在线电影| 亚洲精品一卡2卡三卡4卡5卡|