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

    單位涌水量的換算方法對比與討論

    2024-03-08 10:49:12周立濤
    地下水 2024年1期
    關(guān)鍵詞:井徑曲線擬合富水

    周立濤,李 旭

    (安徽理工大學(xué) 地球與環(huán)境學(xué)院,安徽 淮南 232001)

    0 前言

    礦區(qū)的含水層富水性的等級(jí)評價(jià)可按單位涌水量來進(jìn)行判定,依據(jù)《礦區(qū)水文地質(zhì)工程地質(zhì)勘探規(guī)范》[1]、《煤炭礦井防治水設(shè)計(jì)規(guī)范》[2]中規(guī)定標(biāo)準(zhǔn)的鉆孔單位涌水量是以孔徑91 mm、抽水水位降深10 m為準(zhǔn)。但是在滿足“一孔多用”及研究區(qū)特定水文地質(zhì)條件下,多數(shù)鉆孔的孔徑并非91 mm、降深也很難精確到10 m,致使單位涌水量很難進(jìn)行對比分析,更難以評價(jià)礦區(qū)的含水層富水性。因此,需要對單位涌水量進(jìn)行標(biāo)準(zhǔn)化的換算。

    奚硯濤等人指出根據(jù)抽水試驗(yàn)規(guī)程對抽水試驗(yàn)數(shù)據(jù)進(jìn)行單位涌水量換算具有必要性、科學(xué)性和實(shí)踐性[3]。薛建坤等人對岱莊礦區(qū)單位涌水量進(jìn)行誤差統(tǒng)計(jì)得出如果降深值太小或太大,會(huì)導(dǎo)致計(jì)算的 q 值偏差很多,從而導(dǎo)致人們對含水層富水性產(chǎn)生錯(cuò)誤的認(rèn)識(shí)[4]。陳陽基于VB6.0軟件對單位涌水量計(jì)算方法進(jìn)行探討得出采用VB6.0軟件可以有效數(shù)據(jù)的處理效率和參數(shù)的準(zhǔn)確性[5]。孫尚云等[6]和莊保輝等[7]提出采用 Excel函數(shù)和圖表功能處理換算過程中的 Q-S 曲線擬合問題,快速、合理地解決了曲線類型判斷及單位涌水量 q 的計(jì)算。劉埔等人提出Q-S曲線的MATLAB擬合單位涌水量換算方法,通過比較擬合結(jié)果的多個(gè)誤差評價(jià)參數(shù),較好的解決了曲線 擬合的精確度問題,使得Q-S曲線擬合效果易于分析對比,回歸方程的確定更加準(zhǔn)確[8]。

    目前,針對單位涌水量換算問題,總體思路是對降深進(jìn)行修正和井徑換算,可以采用VB6.0、MATLAB和Excel軟件進(jìn)行輔助換算。然而,在實(shí)際操作計(jì)算中,為簡化計(jì)算通常將R孔和R91進(jìn)行等效處理,即R91的數(shù)值用R孔的數(shù)值代替。實(shí)際上R孔和R91是存在差異的,若等效處理,這必然會(huì)影響含水層富水性評價(jià)的精度。為此,本研究基于前人研究的3次抽水試驗(yàn)實(shí)測數(shù)據(jù),首先采用迭代方法進(jìn)行井徑的換算,并分別計(jì)算R91和R孔,之后采用MATLAB軟件進(jìn)行Q-S曲線的擬合得出單位涌水量并進(jìn)行含水層富水性評價(jià)。此方法避免了修正好的降深發(fā)生改變的問題,并且提高了計(jì)算精度。此外,將計(jì)算的結(jié)果與傳統(tǒng)的方法進(jìn)行對比,分析影響半徑變化對計(jì)算結(jié)果誤差的影響。

    1 單位涌水量的換算方法

    1.1 傳統(tǒng)的計(jì)算方法

    1.1.1 降深修正

    根據(jù)抽水試驗(yàn)獲得Q和S的三組數(shù)據(jù),利用最小二乘法或圖解法確定Q-S曲線的類型,有直線型、拋物線型、冪函數(shù)型或?qū)?shù)型。目前可以借助VB6.0、Excel和MATLAB軟件進(jìn)行Q-S曲線的擬合,得到Q-S的數(shù)學(xué)關(guān)系式。將S=10 m帶入到獲得的Q-S曲線,計(jì)算得到實(shí)際孔徑時(shí)降深為10 m時(shí)的涌水量。

    1.1.2 井徑修正

    降深修正后,由于井徑的差異,因此要進(jìn)一步對井徑進(jìn)行修正。地下水向承壓水井的穩(wěn)定運(yùn)動(dòng)理論滿足承壓含水層完整井穩(wěn)定流時(shí)Dupuit公式[10],如式(1):

    (1)

    式中:S井中水位降深,m;Q為抽水井流量,m3/d;M為含水層厚度,m;K為滲透系數(shù),m/d;r為井的半徑,m;R為影響半徑,m。

    針對同一含水層不同流量抽水時(shí)的K值,在計(jì)算公式一致的情況下,理論上差別不大。一般情況下采用集哈爾特經(jīng)驗(yàn)方程[11],如下式(2):

    (2)

    根據(jù)抽水試驗(yàn)已知水位降深、抽水井流量、含水層厚度和井的半徑,采用迭代方法計(jì)算滲透系數(shù)、降深修正后井孔的影響半徑和孔徑為91 mm時(shí)的影響半徑。具體步驟如下:首先初步確定K0值,代入式(2)可以求出R1,再將R1代入式(1)得到K1,之后再將得到的K1值代入式(2)得到R2,將得到的R2代入式(1)……直到K和R不再發(fā)生變化即可結(jié)束迭代,此時(shí)可以確定K和R的數(shù)值。將R孔和R91進(jìn)行等效處理為R,應(yīng)用下式(3)進(jìn)行流量的換算便可得出單位涌水量。

    (3)

    式中:Q91、R91、r91分別為孔徑為91 mm時(shí)的涌水量、影響半徑和鉆孔半徑;Q孔、R孔、r孔分別為經(jīng)過降深修正后井孔的穩(wěn)定抽水量、影響半徑和鉆井半徑。

    1.2 改進(jìn)的計(jì)算方法

    1.2.1 井徑與影響半徑修正

    為了能夠分別獲得R孔和R91,需要對抽水試驗(yàn)3次降深的實(shí)測數(shù)據(jù)分別進(jìn)行迭代換算出井徑為91 mm時(shí)的Q-S數(shù)據(jù)。迭代計(jì)算的方法如下:

    (1)第一步,根據(jù)抽水試驗(yàn)1次降深獲得的Q1和S1,另外初步假定含水層的滲透系數(shù)為K0,代入式(2)可以求出R孔-0,再將R孔-0代入式(1)得到K1,之后再將得到的K1值代入式(2)得到R孔-1,將得到的R孔-1代入式(1)……直到K和R孔不再發(fā)生變化即可結(jié)束迭代。此時(shí)可以確定K和R孔的數(shù)值。

    (2)第二步,如果井徑修正為91 mm,由于含水層滲透系數(shù)不變,但井徑變小,致使出水能力變?nèi)酢T谙嗤髁肯?井徑降深便會(huì)增加。降深增加致使降落漏斗變大,從而影響半徑增大。此時(shí)已知滲透系數(shù)、抽水井流量、含水層厚度和井徑,采用迭代法求取孔徑為91 mm時(shí)的涌水量的影響半徑和孔徑為91 mm時(shí)的降深。首先初步假定S91-0,代入式(2)得到R91-0,再將R91-0代入式(1)得到S91-1,之后再將得到的S91-1值代入式(2)得到R91-1,將得到的R91-1代入式(1)……直到S91和R91不再發(fā)生變化即可結(jié)束迭代。得到S91和R91的數(shù)值。

    (3)第三步,應(yīng)用式(3)進(jìn)行流量的換算。將上述得到的值代入式(3)即可得到井徑修正后的流量。此時(shí)便將實(shí)際鉆孔的Q-S數(shù)值換算成井徑為91 mm時(shí)的Q-S數(shù)據(jù)。

    上述步驟是對1次降深得到的Q-S數(shù)據(jù)進(jìn)行換算,所以需要進(jìn)行3次上述操作,將3次抽水試驗(yàn)實(shí)際Q-S數(shù)據(jù)換算成井徑為91 mm時(shí)的數(shù)據(jù)。

    1.2.2 降深修正

    根據(jù)上述換算方法可得出井徑為91 mm時(shí)3次降深的抽水試驗(yàn)數(shù)據(jù),利用最小二乘法或圖解法確定Q-S曲線的類型。本研究采用MATLAB軟件進(jìn)行Q-S曲線的擬合,獲取得到擬合效果最好的數(shù)學(xué)關(guān)系式,之后進(jìn)行單位涌水量的計(jì)算和含水層富水性的評價(jià)。

    2 實(shí)例分析

    本研究采用李超峰等[9]在亭南礦區(qū)進(jìn)行抽水試驗(yàn)的數(shù)據(jù),如下表1展開單位涌水量換算的實(shí)例分析。

    表1 亭南礦區(qū)2-1水文孔抽水試驗(yàn)數(shù)據(jù)表

    2.1 傳統(tǒng)計(jì)算方法

    2.1.1 Q-S曲線擬合

    本研究采用MATLAB軟件進(jìn)行Q-S曲線擬合,擬合步驟是首先將三組抽水試驗(yàn)數(shù)據(jù)輸入到MATLAB的工作區(qū),之后調(diào)用cftool曲線擬合工具箱,在Curve Fitting Tool窗口中把對應(yīng)的 x(Q) 、y(S)代入,使用Fitting進(jìn)行曲線擬合。最終生成擬合曲線和參數(shù)結(jié)果。通過對上述參數(shù)的綜合分析可以確定最佳擬合的曲線,進(jìn)而確定曲線的方程,計(jì)算降深為10 m時(shí)鉆孔的流量。

    采用MATLAB軟件分別進(jìn)行直線型、拋物線型、對數(shù)型和冪函數(shù)型四種曲線形態(tài)進(jìn)行擬合得出圖1。

    圖1 傳統(tǒng)計(jì)算方法的MATLAB擬合曲線圖

    根據(jù)MATLAB軟件獲取得到的四種曲線類型擬合效果參數(shù)結(jié)果統(tǒng)計(jì)得到表2,四種曲線類型擬合方程參數(shù)見表3??梢?拋物線型和冪函數(shù)型的擬合效果更好,因此在進(jìn)行單位涌水量的換算時(shí)使用拋物線型和冪函數(shù)型進(jìn)行計(jì)算。

    表2 傳統(tǒng)計(jì)算方法的曲線擬合效果表

    表3 傳統(tǒng)計(jì)算方法的曲線類型擬合參數(shù)表

    2.1.2 井徑換算

    迭代方法的計(jì)算步驟詳見1.1.2,此時(shí)我們將R孔≈R91采用Excel軟件可以節(jié)省計(jì)算時(shí)間,其計(jì)算精度也可以保證。采用Excel進(jìn)行的迭代計(jì)算結(jié)果如表4。

    表4 傳統(tǒng)計(jì)算方法的井徑換算計(jì)算表

    根據(jù)表4的數(shù)據(jù)可知,拋物線型和冪函數(shù)型得出的單位涌水量的數(shù)值在保留小數(shù)點(diǎn)后三位的情況下是一致的,為0.148。

    根據(jù)《礦區(qū)水文地質(zhì)工程地質(zhì)勘探規(guī)范》[1]、《煤炭礦井防治水設(shè)計(jì)規(guī)范》[2]的規(guī)定,按單位涌水量q劃分含水層富水性標(biāo)準(zhǔn)見表5。結(jié)合表5的含水層富水性劃分表可知,該區(qū)屬于中等富水性。

    表5 含水層富水性劃分表

    2.2 改進(jìn)的計(jì)算方法

    2.2.1 井徑換算

    根據(jù)抽水試驗(yàn)的三組數(shù)據(jù),采用1.2.1的步驟進(jìn)行迭代計(jì)算,此時(shí)將R孔和R91的數(shù)值分開計(jì)算得到表6。

    表6 更正計(jì)算方法的井徑換算計(jì)算表

    2.2.2 降深換算

    同樣采用2.1.1的方法進(jìn)行Q-S曲線的擬合,分別進(jìn)行直線型、拋物線型、對數(shù)型和冪函數(shù)型四種曲線形態(tài)進(jìn)行擬合得出圖2。

    圖2 更正計(jì)算方法的MATLAB擬合曲線圖

    根據(jù)MATLAB軟件獲取得到的四種曲線類型擬合效果參數(shù)結(jié)果統(tǒng)計(jì)成表7,四種曲線類型擬合方程參數(shù)見表8。

    表7 更正計(jì)算方法的曲線擬合效果表

    表8 更正計(jì)算方法的曲線類型擬合參數(shù)表

    根據(jù)表8可知,拋物線型的擬合效果相對更好,即Q=-0.008 987S2+11.35S-0.021 85。將S=10 m代入式中可知Q=112.579 5 m3/d。換算成單位涌水量即q=0.130 3 L/(s·m)。根據(jù)表5可知,該區(qū)屬于中等富水性。

    2.3 誤差分析

    通過系列的計(jì)算,我們得出傳統(tǒng)的換算方法計(jì)算出的單位涌水量為0.148 L/(s·m),而改進(jìn)過后的方法計(jì)算出的單位涌水量為0.130 3 L/(s·m) ,其相對誤差為12%。結(jié)果分析表明兩種計(jì)算方法存在一定的差異,因此有必要在單位涌水量換算過程中考慮影響半徑差異帶來的誤差。

    3 結(jié)語

    本文利用抽水實(shí)驗(yàn)數(shù)據(jù),采用兩種不同的方法對含水層單位涌水量進(jìn)行計(jì)算,分析影響半徑變化對計(jì)算結(jié)果誤差的影響,得出以下結(jié)論。

    (1)單位涌水量是評價(jià)含水層富水性的重要指標(biāo),改進(jìn)的單位涌水量換算方能夠分別計(jì)算出R91和R孔,這對單位涌水量的準(zhǔn)確計(jì)算提供了一定的理論支撐。

    (2)本研究采用Excel軟件展開迭代計(jì)算,使用MATLAB軟件進(jìn)行Q-S曲線的擬合得出單位涌水量的結(jié)果。該方法可以合理地找到最佳曲線類型,并計(jì)算結(jié)果的精度能夠得到保證。

    (3)在討論分析傳統(tǒng)計(jì)算方法和更正后的方法時(shí),兩者計(jì)算結(jié)果對比發(fā)現(xiàn)存在一定的差異,因此有必要在單位涌水量換算過程中考慮影響半徑差異帶來的誤差。

    猜你喜歡
    井徑曲線擬合富水
    富水砂卵石地層錨索施工工藝工法
    超級(jí)高密度電法在新疆某礦區(qū)富水性評價(jià)應(yīng)用研究
    曲線擬合的方法
    基于曲線擬合的投棄式剖面儀電感量算法
    電子測試(2017年12期)2017-12-18 06:35:46
    大井徑方位密度測井儀
    Matlab曲線擬合工具箱在地基沉降預(yù)測模型中的應(yīng)用
    Matlab曲線擬合法在地基沉降預(yù)測中的應(yīng)用
    掘進(jìn)巷道構(gòu)造富水性電法探測綜合應(yīng)用研究
    基于最小二乘橢圓擬合的井徑測量面偏角修正方法
    淺埋偏壓富水隧道掘進(jìn)支護(hù)工藝分析
    河南科技(2015年4期)2015-02-27 14:21:05
    久久精品国产自在天天线| 精品无人区乱码1区二区| 久久精品国产99精品国产亚洲性色| 成人综合一区亚洲| 麻豆久久精品国产亚洲av| 久久精品国产亚洲av香蕉五月| 床上黄色一级片| 色吧在线观看| 99久久精品一区二区三区| 尤物成人国产欧美一区二区三区| 亚洲天堂国产精品一区在线| 黄色日韩在线| 99热网站在线观看| 免费观看的影片在线观看| 亚洲国产精品久久男人天堂| av福利片在线观看| 亚洲一级一片aⅴ在线观看| 好男人视频免费观看在线| 国产精品伦人一区二区| 欧美激情国产日韩精品一区| 91av网一区二区| 嫩草影院精品99| 亚洲人与动物交配视频| kizo精华| av免费观看日本| 偷拍熟女少妇极品色| 午夜激情欧美在线| 日日摸夜夜添夜夜添av毛片| 日产精品乱码卡一卡2卡三| 男女边吃奶边做爰视频| 午夜激情福利司机影院| 成人综合一区亚洲| 寂寞人妻少妇视频99o| 我的女老师完整版在线观看| 久久九九热精品免费| 高清在线视频一区二区三区 | 国产av一区在线观看免费| 亚洲国产精品合色在线| 久久精品久久久久久噜噜老黄 | 一级av片app| 亚洲av中文字字幕乱码综合| 国内精品一区二区在线观看| 美女大奶头视频| 男女做爰动态图高潮gif福利片| 丰满乱子伦码专区| 人妻系列 视频| 久久久久久久久中文| 国内少妇人妻偷人精品xxx网站| 五月玫瑰六月丁香| 中国美女看黄片| 欧洲精品卡2卡3卡4卡5卡区| 插逼视频在线观看| 亚洲内射少妇av| 亚洲最大成人手机在线| 日韩一区二区视频免费看| av在线蜜桃| 国产精品一区二区三区四区久久| 最近2019中文字幕mv第一页| 亚洲国产日韩欧美精品在线观看| 18禁黄网站禁片免费观看直播| 日本免费a在线| 亚洲成a人片在线一区二区| 欧美高清成人免费视频www| 国产乱人视频| 偷拍熟女少妇极品色| 久久99热6这里只有精品| 丝袜喷水一区| 直男gayav资源| www日本黄色视频网| 国产亚洲精品av在线| 啦啦啦观看免费观看视频高清| 日韩国内少妇激情av| 久久久久久九九精品二区国产| 久久人人精品亚洲av| 九色成人免费人妻av| 激情 狠狠 欧美| 91久久精品国产一区二区成人| 久久久色成人| 99久久成人亚洲精品观看| 午夜免费男女啪啪视频观看| av福利片在线观看| 日本五十路高清| 亚洲自偷自拍三级| 国产精品国产三级国产av玫瑰| 亚洲av第一区精品v没综合| 一级毛片aaaaaa免费看小| 欧美成人一区二区免费高清观看| 少妇裸体淫交视频免费看高清| 一个人观看的视频www高清免费观看| 插阴视频在线观看视频| 能在线免费观看的黄片| 少妇熟女欧美另类| 一级av片app| 色5月婷婷丁香| 成人毛片a级毛片在线播放| 久久精品影院6| av女优亚洲男人天堂| 啦啦啦观看免费观看视频高清| 成人毛片60女人毛片免费| 天堂√8在线中文| 精品欧美国产一区二区三| 99热这里只有是精品50| 久久精品人妻少妇| 观看美女的网站| 麻豆av噜噜一区二区三区| 亚洲欧美日韩东京热| 国产高清激情床上av| 男女做爰动态图高潮gif福利片| 六月丁香七月| 午夜精品一区二区三区免费看| 色哟哟哟哟哟哟| 色哟哟·www| 一进一出抽搐动态| 日韩,欧美,国产一区二区三区 | 国产蜜桃级精品一区二区三区| 亚洲七黄色美女视频| 欧美成人精品欧美一级黄| 黄色日韩在线| 精品久久国产蜜桃| 99久国产av精品| 国产精品一区二区三区四区免费观看| 日韩中字成人| 久久精品国产亚洲av香蕉五月| 男人舔奶头视频| 只有这里有精品99| 国产一区二区三区av在线 | 一卡2卡三卡四卡精品乱码亚洲| 国产高清三级在线| 看片在线看免费视频| 白带黄色成豆腐渣| 午夜福利高清视频| 男女下面进入的视频免费午夜| 日韩欧美三级三区| 欧美一区二区国产精品久久精品| 69av精品久久久久久| 狂野欧美白嫩少妇大欣赏| 天堂影院成人在线观看| 亚洲性久久影院| 人体艺术视频欧美日本| 国产精品爽爽va在线观看网站| 日本黄色片子视频| 亚洲人成网站在线播放欧美日韩| 成人特级av手机在线观看| 日韩 亚洲 欧美在线| 精品久久久噜噜| 99久久久亚洲精品蜜臀av| 91久久精品国产一区二区三区| 婷婷色av中文字幕| 精品熟女少妇av免费看| 日本免费a在线| 欧洲精品卡2卡3卡4卡5卡区| 又粗又爽又猛毛片免费看| 综合色丁香网| 天堂网av新在线| 99国产极品粉嫩在线观看| 69人妻影院| 直男gayav资源| 国产女主播在线喷水免费视频网站 | 在线国产一区二区在线| 热99在线观看视频| 在现免费观看毛片| 欧美激情久久久久久爽电影| 国产精品三级大全| 亚洲图色成人| 国产av在哪里看| 嫩草影院入口| 人人妻人人看人人澡| 91av网一区二区| 最近手机中文字幕大全| 久久久久久久久久成人| 床上黄色一级片| 嘟嘟电影网在线观看| 全区人妻精品视频| 欧美性感艳星| 欧美一级a爱片免费观看看| 亚洲精品日韩在线中文字幕 | 色哟哟哟哟哟哟| 我要搜黄色片| 免费观看精品视频网站| 国产伦精品一区二区三区视频9| 亚洲最大成人手机在线| www.av在线官网国产| 亚洲精品色激情综合| 成人av在线播放网站| 91狼人影院| 26uuu在线亚洲综合色| 国产精华一区二区三区| 成人高潮视频无遮挡免费网站| 日韩欧美三级三区| 国产91av在线免费观看| 91午夜精品亚洲一区二区三区| 国产高清有码在线观看视频| 亚洲av成人精品一区久久| 男女啪啪激烈高潮av片| 一本久久中文字幕| 亚洲最大成人中文| 亚洲国产欧洲综合997久久,| 欧美另类亚洲清纯唯美| 只有这里有精品99| 禁无遮挡网站| 蜜桃久久精品国产亚洲av| 插逼视频在线观看| 日本免费a在线| 国产在线精品亚洲第一网站| 91aial.com中文字幕在线观看| 一级毛片我不卡| 在线国产一区二区在线| 亚洲国产欧洲综合997久久,| 直男gayav资源| 波多野结衣高清无吗| 美女内射精品一级片tv| 小蜜桃在线观看免费完整版高清| 99热只有精品国产| 成人二区视频| 你懂的网址亚洲精品在线观看 | 久久人人爽人人片av| 欧美日韩综合久久久久久| 中国国产av一级| 一本精品99久久精品77| av国产免费在线观看| 国产淫片久久久久久久久| 国产毛片a区久久久久| 变态另类丝袜制服| av在线老鸭窝| 久久99精品国语久久久| 精品一区二区三区视频在线| 九九在线视频观看精品| 欧美在线一区亚洲| 色5月婷婷丁香| 中文精品一卡2卡3卡4更新| 男人狂女人下面高潮的视频| 日日啪夜夜撸| 老司机影院成人| 亚洲最大成人中文| 一级二级三级毛片免费看| 色尼玛亚洲综合影院| 人人妻人人澡欧美一区二区| 国产真实乱freesex| 国内精品一区二区在线观看| 国产精品女同一区二区软件| 午夜视频国产福利| 亚洲18禁久久av| 91精品国产九色| 久久热精品热| 麻豆国产av国片精品| 伦理电影大哥的女人| 亚洲欧洲日产国产| 啦啦啦观看免费观看视频高清| 最新中文字幕久久久久| 美女高潮的动态| 1000部很黄的大片| 国产精品久久视频播放| 岛国毛片在线播放| 综合色av麻豆| a级一级毛片免费在线观看| 亚洲综合色惰| 亚洲av电影不卡..在线观看| 久久精品人妻少妇| 好男人在线观看高清免费视频| 欧美不卡视频在线免费观看| 超碰av人人做人人爽久久| 国产日韩欧美在线精品| 99热网站在线观看| 精品99又大又爽又粗少妇毛片| 91午夜精品亚洲一区二区三区| 久久久久九九精品影院| 少妇高潮的动态图| 狠狠狠狠99中文字幕| 国产av麻豆久久久久久久| 最新中文字幕久久久久| 久久精品国产99精品国产亚洲性色| av福利片在线观看| 亚洲国产欧洲综合997久久,| 人妻少妇偷人精品九色| 午夜精品国产一区二区电影 | 天堂√8在线中文| 精品久久久久久成人av| 欧美日韩精品成人综合77777| 男人和女人高潮做爰伦理| 久久久久国产网址| 黄色视频,在线免费观看| 啦啦啦观看免费观看视频高清| avwww免费| 久久久久性生活片| 久久久午夜欧美精品| 高清毛片免费观看视频网站| 国产精品一区二区在线观看99 | 91av网一区二区| 简卡轻食公司| 成年女人看的毛片在线观看| 欧美成人一区二区免费高清观看| 99久久九九国产精品国产免费| 国产精品日韩av在线免费观看| 观看免费一级毛片| 婷婷六月久久综合丁香| 性欧美人与动物交配| 人妻夜夜爽99麻豆av| 18禁裸乳无遮挡免费网站照片| 不卡视频在线观看欧美| 亚洲最大成人av| 亚洲人成网站高清观看| av免费在线看不卡| 亚洲国产精品久久男人天堂| 色噜噜av男人的天堂激情| 99热全是精品| 日本免费一区二区三区高清不卡| 免费观看的影片在线观看| 伦理电影大哥的女人| 日韩成人av中文字幕在线观看| 国产成人精品一,二区 | 国产在线男女| 亚洲欧洲日产国产| 亚洲丝袜综合中文字幕| 99在线人妻在线中文字幕| 国产在视频线在精品| 国产精品无大码| 国产人妻一区二区三区在| 国产一级毛片七仙女欲春2| 国产国拍精品亚洲av在线观看| 美女高潮的动态| ponron亚洲| 我要搜黄色片| 国产色婷婷99| 日日撸夜夜添| 中文字幕人妻熟人妻熟丝袜美| 看非洲黑人一级黄片| 麻豆av噜噜一区二区三区| 神马国产精品三级电影在线观看| 国产成人aa在线观看| 你懂的网址亚洲精品在线观看 | 国产在视频线在精品| 日韩中字成人| 中出人妻视频一区二区| 久久热精品热| 精品无人区乱码1区二区| 国产女主播在线喷水免费视频网站 | 搡女人真爽免费视频火全软件| 男人狂女人下面高潮的视频| 亚洲精品日韩在线中文字幕 | 黄片无遮挡物在线观看| 国产精品一区www在线观看| 一级黄色大片毛片| 国产精品久久久久久av不卡| 波野结衣二区三区在线| 欧美在线一区亚洲| 色5月婷婷丁香| 听说在线观看完整版免费高清| 亚洲三级黄色毛片| 亚洲国产欧美在线一区| 久久久久久久久久久丰满| 色噜噜av男人的天堂激情| 青春草国产在线视频 | 一本久久精品| 欧美日韩一区二区视频在线观看视频在线 | 国产精品一二三区在线看| 日韩av不卡免费在线播放| 亚洲精品国产av成人精品| 一个人看视频在线观看www免费| 欧美性猛交黑人性爽| 大又大粗又爽又黄少妇毛片口| www.色视频.com| 国产久久久一区二区三区| 亚洲国产精品sss在线观看| 国产中年淑女户外野战色| 国产精品一区二区性色av| 国产成人影院久久av| av福利片在线观看| 久久6这里有精品| 成人永久免费在线观看视频| 偷拍熟女少妇极品色| 亚洲欧美成人精品一区二区| 丰满乱子伦码专区| 午夜视频国产福利| 我要看日韩黄色一级片| 成人永久免费在线观看视频| 免费大片18禁| 欧美人与善性xxx| 欧美性猛交╳xxx乱大交人| 可以在线观看毛片的网站| 深夜精品福利| 亚洲在线观看片| 久久久久久久久久黄片| 国产久久久一区二区三区| 26uuu在线亚洲综合色| 女的被弄到高潮叫床怎么办| 91aial.com中文字幕在线观看| 久久人人爽人人片av| 国产精品久久久久久亚洲av鲁大| 美女高潮的动态| 村上凉子中文字幕在线| 日本黄色片子视频| 久久精品国产自在天天线| 中出人妻视频一区二区| 国内揄拍国产精品人妻在线| 乱码一卡2卡4卡精品| 一区福利在线观看| 国产精品蜜桃在线观看 | 精品一区二区三区视频在线| 天天躁夜夜躁狠狠久久av| 国产精品久久久久久精品电影小说 | 悠悠久久av| 亚洲国产色片| 久久热精品热| 中文字幕精品亚洲无线码一区| 亚洲激情五月婷婷啪啪| 免费看av在线观看网站| 色噜噜av男人的天堂激情| 晚上一个人看的免费电影| 18禁黄网站禁片免费观看直播| 欧美极品一区二区三区四区| 精品国产三级普通话版| 国产av一区在线观看免费| 99热只有精品国产| 啦啦啦观看免费观看视频高清| 一边亲一边摸免费视频| 97超碰精品成人国产| 99九九线精品视频在线观看视频| 欧美人与善性xxx| 欧美在线一区亚洲| .国产精品久久| 免费人成视频x8x8入口观看| 免费观看在线日韩| 男人的好看免费观看在线视频| 国产高清不卡午夜福利| 男人的好看免费观看在线视频| 国产精品一区二区性色av| 久久久色成人| 久久久久久九九精品二区国产| 久久精品国产亚洲网站| 免费看a级黄色片| 成年免费大片在线观看| 亚洲欧美成人综合另类久久久 | 成人亚洲欧美一区二区av| 能在线免费看毛片的网站| 夜夜看夜夜爽夜夜摸| 成人性生交大片免费视频hd| 国产极品天堂在线| 免费看日本二区| 婷婷六月久久综合丁香| 国产 一区 欧美 日韩| 国产高清视频在线观看网站| 亚洲美女搞黄在线观看| 久久久久久九九精品二区国产| 99国产精品一区二区蜜桃av| 国产 一区 欧美 日韩| 久久99精品国语久久久| 一边摸一边抽搐一进一小说| 国产高清不卡午夜福利| 99久国产av精品国产电影| 日韩一本色道免费dvd| 国产成人aa在线观看| 亚洲av中文av极速乱| 少妇熟女欧美另类| 久久精品国产自在天天线| 在线播放国产精品三级| 国产成年人精品一区二区| 欧美一级a爱片免费观看看| 淫秽高清视频在线观看| av在线亚洲专区| 国产黄片视频在线免费观看| 成人特级黄色片久久久久久久| av女优亚洲男人天堂| 在线观看av片永久免费下载| 91精品国产九色| 久久99热这里只有精品18| 老女人水多毛片| 亚洲国产欧洲综合997久久,| 91精品一卡2卡3卡4卡| 毛片女人毛片| 成人永久免费在线观看视频| 自拍偷自拍亚洲精品老妇| 91aial.com中文字幕在线观看| 男人的好看免费观看在线视频| 国产日韩欧美在线精品| 国产一区二区在线观看日韩| 热99re8久久精品国产| 亚洲av二区三区四区| 成人漫画全彩无遮挡| 国产伦精品一区二区三区四那| 国产精品国产三级国产av玫瑰| 欧美xxxx性猛交bbbb| 网址你懂的国产日韩在线| 美女黄网站色视频| 国产一级毛片在线| 男女下面进入的视频免费午夜| 国产精品麻豆人妻色哟哟久久 | 两个人的视频大全免费| 午夜激情欧美在线| 免费观看a级毛片全部| 乱码一卡2卡4卡精品| 久久精品影院6| 在线国产一区二区在线| 91久久精品国产一区二区成人| 亚洲国产高清在线一区二区三| 亚洲欧美中文字幕日韩二区| 日本与韩国留学比较| 久久久久久久久久久丰满| 五月玫瑰六月丁香| 国产伦在线观看视频一区| 日韩欧美精品免费久久| 国产精品综合久久久久久久免费| 色视频www国产| 26uuu在线亚洲综合色| 又黄又爽又刺激的免费视频.| 最近视频中文字幕2019在线8| 国产免费男女视频| 男女视频在线观看网站免费| 久久热精品热| 成人亚洲精品av一区二区| 长腿黑丝高跟| 91精品一卡2卡3卡4卡| 亚洲高清免费不卡视频| 永久网站在线| 欧美色欧美亚洲另类二区| 精品国产三级普通话版| 一级毛片aaaaaa免费看小| 国产精品女同一区二区软件| av.在线天堂| 狂野欧美白嫩少妇大欣赏| 97人妻精品一区二区三区麻豆| 亚洲天堂国产精品一区在线| 久久久久久国产a免费观看| 成人永久免费在线观看视频| 99久久精品国产国产毛片| 天天躁夜夜躁狠狠久久av| 啦啦啦啦在线视频资源| 深夜精品福利| 欧美日本视频| 亚洲高清免费不卡视频| 神马国产精品三级电影在线观看| 一卡2卡三卡四卡精品乱码亚洲| 99在线人妻在线中文字幕| 一个人看视频在线观看www免费| 97热精品久久久久久| 欧美+亚洲+日韩+国产| 欧美成人免费av一区二区三区| 我的女老师完整版在线观看| 国产午夜精品论理片| 一本久久精品| 欧美日本亚洲视频在线播放| 日本一本二区三区精品| 国产精品女同一区二区软件| 国产一级毛片七仙女欲春2| 亚洲成a人片在线一区二区| 最新中文字幕久久久久| 国产午夜精品一二区理论片| 亚洲aⅴ乱码一区二区在线播放| 久久精品91蜜桃| 免费看日本二区| 1024手机看黄色片| 久久久午夜欧美精品| 国产91av在线免费观看| 亚洲av成人av| 尾随美女入室| 午夜精品在线福利| 免费看a级黄色片| 久久久久久久亚洲中文字幕| 国产视频首页在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 精品一区二区三区视频在线| 少妇被粗大猛烈的视频| 可以在线观看毛片的网站| 天天躁日日操中文字幕| 91久久精品电影网| 中出人妻视频一区二区| 亚洲欧洲日产国产| 黄色配什么色好看| 成人性生交大片免费视频hd| 黄片wwwwww| 成人综合一区亚洲| 天天躁日日操中文字幕| 男的添女的下面高潮视频| 国产精品久久电影中文字幕| 成人特级av手机在线观看| 变态另类成人亚洲欧美熟女| 亚洲一区高清亚洲精品| 五月玫瑰六月丁香| 国产精品.久久久| 熟女人妻精品中文字幕| 国产三级在线视频| 十八禁国产超污无遮挡网站| 黄色视频,在线免费观看| 国产精品久久久久久久电影| 欧美成人精品欧美一级黄| 九色成人免费人妻av| 变态另类成人亚洲欧美熟女| 美女 人体艺术 gogo| 丰满的人妻完整版| 欧美高清成人免费视频www| 中国美女看黄片| 精品人妻视频免费看| 啦啦啦韩国在线观看视频| 久久韩国三级中文字幕| 亚洲av中文av极速乱| 久久久久久久亚洲中文字幕| 一个人免费在线观看电影| 两个人的视频大全免费| 人妻久久中文字幕网| 97在线视频观看| 欧美日韩一区二区视频在线观看视频在线 | 国产大屁股一区二区在线视频| 亚洲精品国产成人久久av| 亚洲人成网站在线播放欧美日韩| 久久久国产成人精品二区| 波多野结衣高清作品| 人人妻人人澡人人爽人人夜夜 | 成人漫画全彩无遮挡| 亚洲中文字幕一区二区三区有码在线看| 如何舔出高潮| 亚洲内射少妇av| 中文亚洲av片在线观看爽| 成人鲁丝片一二三区免费| 变态另类丝袜制服| 在线观看66精品国产| 国产精品一二三区在线看|