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

    基于ANSYS的復(fù)雜地形條件下電法數(shù)據(jù)地形改正研究

    2023-08-23 12:44:00尚建閣李晨暉蘇麗娟邢尚鑫
    礦產(chǎn)與地質(zhì) 2023年3期
    關(guān)鍵詞:坡度電阻率數(shù)值

    尚建閣,李晨暉,蘇麗娟,李 冰,謝 崇,邢尚鑫

    (1.河南省有色金屬礦產(chǎn)探測工程技術(shù)研究中心,河南 鄭州 450016;2.河南省有色金屬地質(zhì)礦產(chǎn)局第二地質(zhì)大隊(duì),河南 鄭州 450016;3.河南省自然資源科技創(chuàng)新中心(礦山生態(tài)環(huán)境保護(hù)修復(fù)研究),河南 鄭州 450016;4.河南有色地質(zhì)礦產(chǎn)集團(tuán)有限公司,河南 鄭州 450016;5.河南省有色金屬地質(zhì)勘查總院,河南 鄭州 450052;6.河南省有色金屬深部找礦勘查技術(shù)研究重點(diǎn)實(shí)驗(yàn)室,河南 鄭州 450052)

    0 引言

    隨著全球能源礦產(chǎn)的需求量不斷增加,找礦工作的程度也越來越深入。直流激發(fā)極化法是一種常規(guī)的物探方法,也是目前應(yīng)用最為廣泛的物探方法之一,通過該方法可以快速的獲取地球內(nèi)部介質(zhì)的極化率、電阻率等物理參數(shù),但地球表面并不是理論上的一個(gè)光滑的平面,而是一個(gè)地形高低不平、連綿起伏的復(fù)雜曲面,特別是在山區(qū),地形變化更加復(fù)雜,正是由于這些復(fù)雜地形的存在,在使用直流激發(fā)極化法開展找礦過程中,發(fā)現(xiàn)的異常會(huì)產(chǎn)生畸變,可能導(dǎo)致淹沒地質(zhì)體的有用信息,呈現(xiàn)假異常,使研究者不能做出正確的推斷解釋,降低了物探綜合參數(shù)利用水平[1-3]。因此,復(fù)雜地形條件下電阻率地形改正方法的研究,成為物探工作者急需解決的關(guān)鍵問題之一。ANSYS軟件是一款大型的融電磁學(xué)、結(jié)構(gòu)力學(xué)、流體力學(xué)、仿聲學(xué)等于一體的大型通用有限元分析軟件[4]。ANSYS具有豐富的單元庫,提供了對(duì)各種物理場量的分析功能[5-6],為地球物理工作者在進(jìn)行數(shù)值模擬分析時(shí)提供了強(qiáng)大的數(shù)學(xué)分析功能。前人運(yùn)用ANSYS開展了一系列的地球物理電法數(shù)據(jù)的數(shù)值模擬,對(duì)典型的地電模型全空間介質(zhì)影響進(jìn)行模擬,驗(yàn)證了ANSYS開展電法數(shù)據(jù)模擬的正確性,均取得了較為理想的理論效果[7-9],但是對(duì)復(fù)雜地形條件下的模擬應(yīng)用較少。

    文章通過介紹運(yùn)用ANSYS開展電法地形改正數(shù)值模擬,總結(jié)經(jīng)驗(yàn)并運(yùn)用到實(shí)際生產(chǎn)找礦工作中,取得了良好應(yīng)用效果,為廣大地球物理工作者在復(fù)雜地形條件下開展直流激發(fā)極化法數(shù)據(jù)解釋推斷提供技術(shù)參考。

    1 方法原理

    地形影響問題是從20世紀(jì)50年代伴隨著直流電法在礦產(chǎn)資源勘查等方面的大規(guī)模應(yīng)用開始受到重視的[10]。復(fù)雜的地形對(duì)直流激發(fā)極化法中電阻率參數(shù)的影響是致命的,因此,正確消除起伏地形對(duì)電阻率參數(shù)的影響,對(duì)后期地球物理成果解釋推斷的真實(shí)性起到了至關(guān)重要的作用[11]。

    研究地形對(duì)視電阻率的影響,目前主要使用的有三種方法:物理模擬法、解析法和數(shù)值模擬法。所謂直流電法的地形改正,就是通過上述三種方法將在起伏地形上測量得到的電阻率數(shù)據(jù)轉(zhuǎn)換為平地形影響下數(shù)據(jù)的過程,從而消除不同地形所引起的誤差。由于物理模擬法和解析法不適用于場源和物性分布復(fù)雜的正演情況,因此,數(shù)值模擬法則成為電阻率法正演最主要的方法[12]。隨著計(jì)算機(jī)技術(shù)的發(fā)展和物探解釋方法的進(jìn)步,數(shù)值模擬技術(shù)得到了飛躍式發(fā)展。常用的電阻率數(shù)值模擬方法有有限差分法、積分方程法、邊界單元法和有限元法,其中有限元法具有較高的求解精度且理論較為成熟,能夠適應(yīng)任意模型形態(tài)。因此,有限元法是目前處理電法正演問題的主流方法。

    ANSYS是一款能實(shí)現(xiàn)多場及多場耦合分析、實(shí)現(xiàn)前后處理、求解及多場分析統(tǒng)一數(shù)據(jù)庫的一體化大型有限元分析軟件,一般而言,ANSYS的基本分析過程可以分為三步,即創(chuàng)建有限元模型、施加載荷進(jìn)行求解和查看分析結(jié)果。對(duì)應(yīng)軟件結(jié)構(gòu)的三個(gè)程序模塊:前處理模塊(PREP7)、分析求解模塊(SOLUTION)和后處理模塊(POST1和POST26)。前處理模塊(PREP7)的主要任務(wù)是建立結(jié)構(gòu)分析的有限單元模型,它是結(jié)構(gòu)分析的開始,一般步驟為分析準(zhǔn)備、設(shè)置單元類型、設(shè)定實(shí)常數(shù)、定義材料屬性、創(chuàng)建模型、劃分網(wǎng)格。分析求解模塊(SOLUTION)的主要任務(wù)是進(jìn)行結(jié)構(gòu)的計(jì)算,一般步驟:定義自由度、施加荷載、設(shè)定分析類型、給定邊界條件及求解等。后處理模塊(POST1和POST26)的主要任務(wù)是進(jìn)行結(jié)構(gòu)的分析,它包括通用后處理器(POST1)和時(shí)程后處理器(PSOT26)兩種,前者主要用于靜態(tài)分析的結(jié)果后處理,而后者主要用于動(dòng)態(tài)分析結(jié)果的后處理。

    2 模擬實(shí)驗(yàn)

    利用ANSYS軟件進(jìn)行地形改正,其一般思路:① 根據(jù)電法工作的實(shí)際情況,依照測點(diǎn)高程建立物理模型,并要根據(jù)計(jì)算精度,對(duì)模型進(jìn)行網(wǎng)格剖分和邊界條件的設(shè)定;② 根據(jù)工作精度的要求,合理選擇單元類型,利用ANSYS程序進(jìn)行純地形條件下各測點(diǎn)的視電阻率值計(jì)算;③ 通過比值法,對(duì)實(shí)測視電阻率數(shù)據(jù)進(jìn)行地形影響的壓制。計(jì)算公式采用ρs改=(ρs實(shí)/ρs地形)×ρ0,通過ANSYS有限元分析得到的純地形視電阻率數(shù)據(jù),對(duì)實(shí)測數(shù)據(jù)進(jìn)行比值壓制,最終得到改正后的ρs改數(shù)據(jù)。總的來說,利用ANSYS軟件進(jìn)行地形改正需要正確記錄實(shí)際地形高程數(shù)據(jù),利用軟件所提供的有限元分析軟件計(jì)算出純地形條件下影響的數(shù)值,再利用實(shí)際測量獲得的電阻率數(shù)值與之進(jìn)行比值法進(jìn)行地形影響消除得到改正后的真實(shí)數(shù)值,基本處理流程見圖1。

    圖1 ANSYS有限元數(shù)據(jù)分析流程圖

    --------------------------------------------------------------------------------------------------------------

    下面利用ANSYS編寫命令流加以說明。

    (1)!進(jìn)入前處理器,建模

    /prep7

    /PNUM,area,1

    et,1,plane230 ! 定義單元類型

    mp,1,resistivity ! 定義材料電阻率參數(shù)

    *ask,KP_N,8,, ! 定義創(chuàng)建地形模型的點(diǎn)數(shù)

    *dim,KP_NX,ARRAY,KP_N,,

    *dim,KP_NY,ARRAY,KP_N,, ! 定義地形節(jié)點(diǎn)的XY數(shù)值

    k,I,KP_NX(I),KP_NY(I) ! 定義關(guān)健節(jié)點(diǎn)

    a,1,2,3,4,5,6,7,8…… ! 由點(diǎn)創(chuàng)建面

    (2)進(jìn)入求解器!將模型網(wǎng)格劃分,定義邊界條件和施加荷載

    asel,s,,,1,2

    aatt,1

    asel,s,,,4

    aatt,2

    asel,s,,,1,4

    esize,,20 ! 定義網(wǎng)絡(luò)劃分尺度

    mshape,1,2d ! 定義剖分類型

    amesh,all ! 開始網(wǎng)絡(luò)劃分

    lsscale,all

    nsel,s,ext ! 選擇邊界上的節(jié)點(diǎn)

    nsel,u,loc,y

    *get,Nnod,node,,count

    *get,Nmin,node,,num,min

    *do,J,1,Nnod,1

    Dist1=SQRT((NX(Nmin)-NAX)**2+(NY(Nmin)-NAY(Nmin))**2)

    Dist2=SQRT((NX(Nmin)-NBX)**2+(NY(Nmin)-NBY(Nmin))**2)

    V0=Currt*Res1*(1/Dist1-1/Dist2)/(2*Pi)

    D,Nmin,VOLT,V0 ! 給邊界上的節(jié)點(diǎn)賦值

    Nmin=ndnext(Nmin)

    *enddo

    fk,2,amps,1

    fk,11,amps,-1 ! 施加電流荷載

    finish

    /sol

    antype,static

    ! /status,solu ! 定義分析類型和求解

    solve

    FINISH

    (3)! 進(jìn)入后處理器,查看結(jié)果

    /post

    *SET,DJ,0.2

    *SET,MN,0.2

    *DIM,V1,ARRAY,20,1,

    *DIM,Ra,TABLE,20,2,

    *DO,I,1,20,1

    V1(I,1)=VOLT(NODE(2+DJ*(I-1),6,0))-VOLT(NODE(2+DJ*I,6,0)) ! 求模擬電位差

    AM=2+DJ*(I-1)

    AN=2+DJ*I

    BM=6-DJ*(I-1)

    BN=6-DJ*I

    KK=2*3.14156/(1/AM-1/AN-1/BM+1/BN)

    Ra(I,1)=2+I*DJ

    Ra(I,2)=V1(I,1)*kk/10

    *enddo

    FINISH

    /GCOL,1,dis

    /GCOL,2,Resistivity

    *VPLOT,Ra(1,1),Ra(1,2),2

    FINISH

    根據(jù)ANSYS軟件進(jìn)行地形改正的思路和分析流程,分別開展了不同坡度角的山谷地形和山脊地形對(duì)直流激發(fā)極化法電阻率影響的數(shù)值模擬(圖2)。圖2a為山脊模型,θ為坡度角,材料電阻率為100 Ω·m,為了研究山脊純地形對(duì)視電阻率的影響程度,分別進(jìn)行了坡度為0°、10°、20°、30°、40°的模型計(jì)算。由圖2a可見,在山脊純地形上得到的視電阻率曲線呈現(xiàn)為低阻特征,并且隨著坡度角的增大,其低阻特征更加明顯。圖2b為山谷模型,θ為坡度角,材料電阻率為100 Ω·m,同樣也進(jìn)行了坡度為0°、10°、20°、30°、40°的模型計(jì)算。由圖2b可見,在山谷純地形上得到的視電阻率曲線呈現(xiàn)為高阻特征,并且隨著坡度角的增大,其高阻特征更加明顯。由圖2可見,純地形異常主要發(fā)生在角域頂點(diǎn)附近,也就是說主要發(fā)生在坡度變化的地方。異常的大小與坡度角的大小有關(guān),即與地形起伏的陡緩程度關(guān)系密切。形態(tài)相同的角域,異常幅度隨坡度的增大而增大。

    圖2 不同地形條件影響下模擬電阻率曲線圖

    為了更好地研究地形對(duì)異常體電阻率曲線的影響,在數(shù)值模擬器中加入異常體進(jìn)行數(shù)值模擬,圖2c為山脊地形下高阻直立薄板電阻率曲線圖,θ為坡度角,背景電阻率為100 Ω·m,直立高阻體電阻率為1000 Ω·m。由圖2c可見,在水平地面時(shí),直立薄板視電阻率異常表現(xiàn)為高阻特征,隨著地形坡度的變大,異常曲線表現(xiàn)為低阻特征,尤其當(dāng)坡度角為30°、40°時(shí),基本為純地形引起的異常曲線特征。圖2d為不同山谷地形下低阻水平薄板電阻率曲線,θ為坡度角,背景電阻率為100 Ω·m,水平薄板的電阻率為10 Ω·m。由圖2d可見,在水平地面時(shí),水平薄板電阻率曲線表現(xiàn)為低阻特征,隨著地形坡度的變大,異常曲線表現(xiàn)為高阻特征,尤其當(dāng)坡度角為30°、40°時(shí),基本為純地形引起的異常曲線特征。由此可以看出,地形起伏對(duì)異常體引起的電阻率曲線影響很大,特別是當(dāng)?shù)匦吻懈顕?yán)重時(shí),對(duì)電阻率資料的影響更大,使異常體電阻率異常發(fā)現(xiàn)較大畸變,甚至出現(xiàn)相反的結(jié)果。圖2e為高阻直立薄板山脊地形改正前后的電阻率曲線。由圖2e可見,藍(lán)色曲線為地改前的電阻率曲線,當(dāng)存在地形影響時(shí),高阻體也表現(xiàn)為低阻特征;紅色曲線為地改后的電阻率曲線,通過比值法改正后,電阻率曲線表現(xiàn)出高阻特征,與實(shí)際地下介質(zhì)情況相符。圖2f為高阻直立薄板山脊地形改正前后的電阻率曲線。由圖2f可見,藍(lán)色曲線為地改前的電阻率曲線,當(dāng)存在地形影響時(shí),由于地形的影響,低阻體在電阻率曲線上也表現(xiàn)為高阻特征;紅色曲線為地改后的電阻率曲線,通過比值法改正后,電阻率曲線表現(xiàn)出低阻特征,與實(shí)際地?cái)嗝娼橘|(zhì)情況相符。通過電流源有限元分析計(jì)算,其結(jié)果充分表明,地形對(duì)電阻率的影響很大,特別是地形起伏較大地區(qū),影響程度更大。通過地形模擬實(shí)驗(yàn)可以總結(jié)得出地形對(duì)電阻率能產(chǎn)生較大影響,具有正地形展現(xiàn)出低阻特性,負(fù)地形展現(xiàn)出高阻特性,這樣的影響隨著地形的復(fù)雜程度不斷提高顯現(xiàn)的越發(fā)明顯,甚至能抵消地下目標(biāo)體所產(chǎn)生的真實(shí)電阻率影響,使得目標(biāo)體的真實(shí)特征不能被探獲。利用ANSYS軟件進(jìn)行快速地形改正有效地解決了地形起伏對(duì)電阻率的影響,使反映的介質(zhì)電阻率回歸為真實(shí)。

    在使用ANSYS軟件進(jìn)行數(shù)值分析的過程中,網(wǎng)格剖分的大小直接影響著計(jì)算結(jié)果的精度,所以網(wǎng)格剖分越細(xì),計(jì)算精度也就有所提高。但是,網(wǎng)格劃分的大小也影響著計(jì)算工作量的大小,直接影響著計(jì)算速度。因此,在應(yīng)用中應(yīng)權(quán)衡這兩種因素,做到綜合考慮。圖3為不同網(wǎng)格剖分密度下的山脊地形電阻率曲線。由圖3可見,隨著剖分尺度的大小變化,曲線的光滑程度也發(fā)生變化,且剖分尺度越大,曲線振蕩程度越小,說明計(jì)算精度越高。在ANSYS穩(wěn)態(tài)電流有限元分析中,程序沒有提供遠(yuǎn)邊界單元類型,在實(shí)際應(yīng)用時(shí),須構(gòu)建一個(gè)遠(yuǎn)邊界模型,對(duì)研究區(qū)和邊界區(qū)進(jìn)行不同的網(wǎng)格劃分。

    圖3 二維純山脊地形電阻率曲線圖(不同剖分尺度)

    3 應(yīng)用實(shí)例

    河南省西部山區(qū)是我國重要的有色金屬成礦區(qū)帶之一,該區(qū)域分布有多個(gè)大大小小的有色貴金屬礦床,為河南省礦產(chǎn)資源需求提供了保障[13-14]。研究區(qū)位于河南省西部洛寧縣,熊耳山北坡,其大地構(gòu)造位置處于中朝準(zhǔn)地臺(tái)南緣、華熊臺(tái)緣凹陷、崤山—魯山拱褶斷束中部(圖4)。地層出露主要為太古宇太華群的角閃斜長片麻巖、熊耳群的安山巖和第四系。區(qū)內(nèi)巖漿活動(dòng)頻繁,巖漿巖種類較多,其中以燕山期酸性巖最為發(fā)育,并與本區(qū)內(nèi)生金屬礦產(chǎn)的形成關(guān)系極為密切。斷裂構(gòu)造十分發(fā)育,并以NE—NNE向斷裂為主,且規(guī)模較大。研究區(qū)圍巖蝕變主要為硅化、絹云母化、綠泥石化、黃鐵礦化、泥化等。硅化表現(xiàn)為石英沿破碎帶充填,形成巨厚層的硅化碎裂巖帶出現(xiàn)。構(gòu)造蝕變帶內(nèi)礦化與蝕變主要有黃鐵礦化、方鉛礦化、碳酸鹽化、高嶺土化、褐鐵礦化、硅化、絹云母化等,金屬礦物表現(xiàn)形式為星點(diǎn)狀、浸染狀、細(xì)脈狀、細(xì)脈侵染狀、團(tuán)塊狀。

    圖4 區(qū)域地質(zhì)礦產(chǎn)圖(a)與研究區(qū)地質(zhì)簡圖(b)

    表1為研究區(qū)電性參數(shù)統(tǒng)計(jì)結(jié)果,由表1可見,研究區(qū)礦石(含礦化體)與圍巖存在著較大的物性差異,目標(biāo)體電性表現(xiàn)特征為低阻高極化特征,為在該研究區(qū)開展物探工作提供物性前提。

    表1 研究區(qū)電性參數(shù)統(tǒng)計(jì)結(jié)果

    本次物探研究工作是在試驗(yàn)區(qū)K9號(hào)礦脈上進(jìn)行,采用時(shí)間域激電偶極-偶極測深裝置,點(diǎn)距20 m,隔離系數(shù)為5,MN極距40 m。儀器為V8,最大供電電流為14 A,6個(gè)電通道方式。圖5為試驗(yàn)剖面測量結(jié)果,由圖5a可見,極化率斷面圖,極化率高值異常與K9號(hào)礦脈基本對(duì)應(yīng),反映出地形對(duì)極化率測量結(jié)果影響較小,而圖5b中所反映電阻率異常與礦脈不對(duì)應(yīng),并且展現(xiàn)出了高電阻率的特征,這一高電阻率特征是地下礦脈真實(shí)電阻率與地形影響產(chǎn)生的電阻率疊加影響產(chǎn)生的綜合反映,其低阻異常向左側(cè)地形相對(duì)低緩的一側(cè)進(jìn)行偏離。這是由于地形切割較大,山谷地形純電阻率異常反映為高阻,致使礦體引起的低電阻異常被掩蓋,發(fā)生偏離現(xiàn)象。圖5c為采用ANSYS軟件計(jì)算出該地形所產(chǎn)生的響應(yīng)值與實(shí)際探測的電阻率值的比值進(jìn)行地形改正計(jì)算,其結(jié)果可見低電阻異常位置向右側(cè)偏移,并與K9號(hào)礦脈基本對(duì)應(yīng),徹底消除了復(fù)雜地形對(duì)探測結(jié)果的影響效應(yīng)。結(jié)合圖5a和圖5c結(jié)果顯示與巖礦石物性結(jié)果一致,表現(xiàn)為低阻高極率特征,證明采用ANSYS軟件進(jìn)行電阻率地形改正使測量結(jié)果回歸真實(shí)反映,效果顯著。

    4 結(jié)論

    (1)激發(fā)極化法是地質(zhì)找礦中應(yīng)用最為廣泛的地球物理方法之一,但是其電阻率數(shù)據(jù)收復(fù)雜地形影響也最為嚴(yán)重,負(fù)地形將形成高電阻率反映,而正地形將產(chǎn)生低電阻率反映,從而實(shí)際物探測量所獲得電阻率帶有地形加載的影響,使得大地介質(zhì)的真實(shí)電阻率得不到有效體現(xiàn)。

    (2)ANSYS軟件是一款能實(shí)現(xiàn)多場及多場耦合分析、實(shí)現(xiàn)前后處理、求解及多場分析統(tǒng)一數(shù)據(jù)庫的一體化大型有限元分析軟件,能快速準(zhǔn)確進(jìn)行多種場景數(shù)值模擬計(jì)算,能夠?qū)?fù)雜地形條件下不同裝置類型的電阻率開展二維、三維地形改正,利用該軟件有限元分析功能模擬計(jì)算出純地形下產(chǎn)生的電阻率響應(yīng),通過比值法將野外所獲取的電阻率值中地形影響消除,獲得大地介質(zhì)真實(shí)電阻率,實(shí)現(xiàn)復(fù)雜地形對(duì)電阻率值的零影響效應(yīng),為物探解釋的可靠性提供保障。

    (3)應(yīng)用ANSYS軟件在研究區(qū)復(fù)雜地形條件下開展電阻率地形改正實(shí)際應(yīng)用,對(duì)比結(jié)果顯示,運(yùn)用該軟件進(jìn)行地形改正方法可行,效果顯著。

    猜你喜歡
    坡度電阻率數(shù)值
    用固定數(shù)值計(jì)算
    數(shù)值大小比較“招招鮮”
    關(guān)于公路超高漸變段合成坡度解析與應(yīng)用
    三維電阻率成像與高聚物注漿在水閘加固中的應(yīng)用
    基于圖像處理的定位器坡度計(jì)算
    電氣化鐵道(2016年4期)2016-04-16 05:59:46
    基于Fluent的GTAW數(shù)值模擬
    焊接(2016年2期)2016-02-27 13:01:02
    隨鉆電阻率測井的固定探測深度合成方法
    坡度在巖石風(fēng)化層解譯中的應(yīng)用
    河北遙感(2015年2期)2015-07-18 11:11:14
    CT和MR對(duì)人上脛腓關(guān)節(jié)面坡度的比較研究
    海洋可控源電磁場視電阻率計(jì)算方法
    又爽又黄a免费视频| 噜噜噜噜噜久久久久久91| 色综合亚洲欧美另类图片| 久久精品国产自在天天线| 亚洲综合色惰| 又粗又硬又长又爽又黄的视频 | 国产一区二区在线观看日韩| 婷婷色综合大香蕉| 高清日韩中文字幕在线| 国产高清视频在线观看网站| 91久久精品国产一区二区三区| 色综合站精品国产| 国产高潮美女av| 秋霞在线观看毛片| 国产色婷婷99| 国产一区二区在线av高清观看| 欧美日韩在线观看h| 日本熟妇午夜| 亚洲最大成人中文| 欧美日韩国产亚洲二区| 日韩国内少妇激情av| 精品久久久久久成人av| 三级毛片av免费| 亚洲av不卡在线观看| 亚洲av免费高清在线观看| 嫩草影院入口| 亚洲av成人av| ponron亚洲| 国产高清有码在线观看视频| 亚洲色图av天堂| 波野结衣二区三区在线| 成人毛片a级毛片在线播放| 亚洲欧美精品专区久久| av天堂在线播放| 久久这里只有精品中国| 在线免费观看的www视频| 欧美日本亚洲视频在线播放| 欧美不卡视频在线免费观看| 精品熟女少妇av免费看| 欧美精品一区二区大全| 亚洲欧美日韩无卡精品| 亚洲国产欧洲综合997久久,| 深爱激情五月婷婷| 欧美色视频一区免费| 亚洲无线观看免费| 成年av动漫网址| 老司机福利观看| 精品无人区乱码1区二区| a级毛色黄片| 成熟少妇高潮喷水视频| 久久99热这里只有精品18| 久久精品夜色国产| 一个人免费在线观看电影| 国产成人a∨麻豆精品| 秋霞在线观看毛片| 99热这里只有是精品50| 日韩,欧美,国产一区二区三区 | 天天一区二区日本电影三级| 色哟哟哟哟哟哟| 看非洲黑人一级黄片| 有码 亚洲区| 久久午夜福利片| 桃色一区二区三区在线观看| 午夜精品国产一区二区电影 | 91狼人影院| 国产探花极品一区二区| 亚洲成a人片在线一区二区| 亚洲精品日韩在线中文字幕 | 国产国拍精品亚洲av在线观看| 2022亚洲国产成人精品| 看十八女毛片水多多多| 在现免费观看毛片| 亚洲av免费在线观看| 激情 狠狠 欧美| 六月丁香七月| 能在线免费观看的黄片| 欧美极品一区二区三区四区| 校园春色视频在线观看| 国产色婷婷99| 精品不卡国产一区二区三区| 男的添女的下面高潮视频| 欧美bdsm另类| 久久久久久大精品| 久久午夜福利片| 国产精品久久久久久精品电影小说 | 菩萨蛮人人尽说江南好唐韦庄 | 男女下面进入的视频免费午夜| 18禁黄网站禁片免费观看直播| 国产成人福利小说| 欧美日韩乱码在线| 一级毛片我不卡| 午夜爱爱视频在线播放| 中文在线观看免费www的网站| 高清在线视频一区二区三区 | 久久久精品大字幕| 人人妻人人看人人澡| 波野结衣二区三区在线| av.在线天堂| 精品不卡国产一区二区三区| 欧美+亚洲+日韩+国产| 国产精品一及| 日本av手机在线免费观看| 国产美女午夜福利| 蜜桃久久精品国产亚洲av| 成人永久免费在线观看视频| 黄色一级大片看看| 波多野结衣高清无吗| 亚洲成人中文字幕在线播放| 久久久久久国产a免费观看| 亚洲aⅴ乱码一区二区在线播放| 久久九九热精品免费| 99热这里只有是精品在线观看| 看片在线看免费视频| 看免费成人av毛片| 久久久久久九九精品二区国产| 亚洲综合色惰| 亚洲七黄色美女视频| 51国产日韩欧美| 观看免费一级毛片| 欧美激情国产日韩精品一区| 国产女主播在线喷水免费视频网站 | 一级av片app| 免费人成在线观看视频色| 国产熟女欧美一区二区| 国产午夜精品久久久久久一区二区三区| 亚洲一级一片aⅴ在线观看| 黄色视频,在线免费观看| 亚洲精品国产av成人精品| 99久国产av精品| 久久久a久久爽久久v久久| 亚洲色图av天堂| 最好的美女福利视频网| 97在线视频观看| 男人的好看免费观看在线视频| 国产一级毛片在线| 国产激情偷乱视频一区二区| 成人永久免费在线观看视频| 欧美+亚洲+日韩+国产| 国产男人的电影天堂91| 色播亚洲综合网| 午夜精品在线福利| 国产 一区 欧美 日韩| 成人高潮视频无遮挡免费网站| 国内少妇人妻偷人精品xxx网站| 国产蜜桃级精品一区二区三区| 欧美激情在线99| 国产亚洲欧美98| 成人欧美大片| 国产精品久久久久久久电影| 18禁黄网站禁片免费观看直播| 亚洲av成人av| 午夜老司机福利剧场| 久久精品国产亚洲av天美| 久久久久免费精品人妻一区二区| 日韩av不卡免费在线播放| 长腿黑丝高跟| 午夜福利在线观看吧| 国产老妇伦熟女老妇高清| 国产亚洲av片在线观看秒播厂 | 少妇的逼好多水| 综合色av麻豆| 日韩亚洲欧美综合| 又黄又爽又刺激的免费视频.| 三级经典国产精品| 国产极品精品免费视频能看的| 国产精品一区二区性色av| 国产麻豆成人av免费视频| 国产真实伦视频高清在线观看| 日本五十路高清| 国产一区二区激情短视频| 好男人视频免费观看在线| 成人亚洲精品av一区二区| 国产午夜精品一二区理论片| 亚洲中文字幕日韩| 高清日韩中文字幕在线| 97热精品久久久久久| 久久人人爽人人片av| 97超碰精品成人国产| 国产亚洲5aaaaa淫片| 丰满乱子伦码专区| 简卡轻食公司| 中文字幕久久专区| 天堂中文最新版在线下载 | 神马国产精品三级电影在线观看| 国产精品日韩av在线免费观看| 亚洲国产色片| 国产色婷婷99| 国产女主播在线喷水免费视频网站 | 久久99热6这里只有精品| 大又大粗又爽又黄少妇毛片口| 色哟哟哟哟哟哟| 天堂√8在线中文| 1024手机看黄色片| 欧美最黄视频在线播放免费| 亚洲国产高清在线一区二区三| 国产精品国产三级国产av玫瑰| 午夜老司机福利剧场| 亚洲国产精品国产精品| 久久鲁丝午夜福利片| 此物有八面人人有两片| 床上黄色一级片| 人妻少妇偷人精品九色| 国产伦精品一区二区三区四那| 熟女电影av网| 日日撸夜夜添| 日韩国内少妇激情av| 91精品国产九色| 伦精品一区二区三区| 国产三级在线视频| 一级毛片久久久久久久久女| 最近的中文字幕免费完整| 日韩大尺度精品在线看网址| 国产精品综合久久久久久久免费| 麻豆国产97在线/欧美| 国产精品久久视频播放| 亚洲国产精品国产精品| 女同久久另类99精品国产91| 国产精品国产三级国产av玫瑰| 午夜老司机福利剧场| 国产高清三级在线| 成人特级黄色片久久久久久久| 一区福利在线观看| 日日摸夜夜添夜夜爱| 精品99又大又爽又粗少妇毛片| 欧美3d第一页| 亚洲无线观看免费| av免费观看日本| 亚洲高清免费不卡视频| 乱码一卡2卡4卡精品| 久久久久久久久久久免费av| 99热网站在线观看| 久久国内精品自在自线图片| 久久精品国产鲁丝片午夜精品| 久久精品国产亚洲av涩爱 | 可以在线观看的亚洲视频| 国产色爽女视频免费观看| 边亲边吃奶的免费视频| 色综合站精品国产| 人妻系列 视频| 国产午夜精品论理片| 国产 一区精品| av视频在线观看入口| 成人亚洲精品av一区二区| 看免费成人av毛片| 久久久成人免费电影| 青青草视频在线视频观看| 国产精品蜜桃在线观看 | 国产一区亚洲一区在线观看| 国产精品久久久久久精品电影小说 | av卡一久久| 男女那种视频在线观看| 天天躁夜夜躁狠狠久久av| 一级二级三级毛片免费看| 日本黄色片子视频| 看免费成人av毛片| 草草在线视频免费看| 国国产精品蜜臀av免费| 欧美xxxx性猛交bbbb| 国产伦理片在线播放av一区 | 夜夜爽天天搞| 日本一本二区三区精品| 欧美日韩一区二区视频在线观看视频在线 | 精品人妻偷拍中文字幕| 91久久精品国产一区二区三区| 婷婷色av中文字幕| 国产黄片美女视频| 丰满的人妻完整版| 最后的刺客免费高清国语| 啦啦啦观看免费观看视频高清| 亚洲激情五月婷婷啪啪| 狠狠狠狠99中文字幕| 在线播放国产精品三级| 悠悠久久av| 综合色丁香网| 亚洲精品乱码久久久v下载方式| av.在线天堂| 狂野欧美激情性xxxx在线观看| 麻豆成人午夜福利视频| 亚洲图色成人| 欧美一区二区国产精品久久精品| 一级毛片久久久久久久久女| 男插女下体视频免费在线播放| 精品免费久久久久久久清纯| 女的被弄到高潮叫床怎么办| 看十八女毛片水多多多| 亚洲国产欧美在线一区| 99热全是精品| 一个人看视频在线观看www免费| 久久这里只有精品中国| 国产真实乱freesex| 麻豆成人av视频| 九九久久精品国产亚洲av麻豆| 精品欧美国产一区二区三| 少妇的逼水好多| 亚洲欧美日韩无卡精品| 午夜老司机福利剧场| 小蜜桃在线观看免费完整版高清| 26uuu在线亚洲综合色| 免费观看的影片在线观看| 99久久精品一区二区三区| 亚洲真实伦在线观看| 久久韩国三级中文字幕| 大又大粗又爽又黄少妇毛片口| 亚洲成人中文字幕在线播放| 久久午夜亚洲精品久久| 亚洲欧美成人精品一区二区| 久久精品国产自在天天线| 好男人在线观看高清免费视频| 欧美+日韩+精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美激情在线99| 成人欧美大片| 一区福利在线观看| 国产精品久久久久久精品电影小说 | 精品人妻熟女av久视频| 久久精品国产亚洲av涩爱 | 国产成人精品一,二区 | 国产精品久久电影中文字幕| 麻豆久久精品国产亚洲av| 精品免费久久久久久久清纯| 最近最新中文字幕大全电影3| 神马国产精品三级电影在线观看| 国产精品久久久久久av不卡| 99精品在免费线老司机午夜| 成人三级黄色视频| 国产精华一区二区三区| 国产极品天堂在线| av在线天堂中文字幕| 久久精品国产鲁丝片午夜精品| 久久精品人妻少妇| 久久韩国三级中文字幕| 老女人水多毛片| 特级一级黄色大片| 免费av毛片视频| 亚洲人与动物交配视频| 国产精品久久久久久精品电影小说 | 国产免费男女视频| 精品人妻一区二区三区麻豆| 大型黄色视频在线免费观看| 性插视频无遮挡在线免费观看| 国产欧美日韩精品一区二区| 少妇被粗大猛烈的视频| 日韩亚洲欧美综合| 日本熟妇午夜| 麻豆国产97在线/欧美| 亚洲av电影不卡..在线观看| 大型黄色视频在线免费观看| a级毛片免费高清观看在线播放| 亚洲图色成人| 欧美又色又爽又黄视频| 嫩草影院新地址| 一级毛片aaaaaa免费看小| 少妇人妻一区二区三区视频| 一本久久中文字幕| 精品少妇黑人巨大在线播放 | 一边亲一边摸免费视频| 午夜福利成人在线免费观看| 国产爱豆传媒在线观看| 一区二区三区免费毛片| 18+在线观看网站| 97超视频在线观看视频| 免费看日本二区| av天堂中文字幕网| 成人漫画全彩无遮挡| 观看免费一级毛片| 久久国产乱子免费精品| 久久久国产成人免费| 亚洲经典国产精华液单| 五月玫瑰六月丁香| 一卡2卡三卡四卡精品乱码亚洲| 一进一出抽搐gif免费好疼| 日本免费a在线| 日本三级黄在线观看| 蜜臀久久99精品久久宅男| 在线观看午夜福利视频| av在线亚洲专区| 97超视频在线观看视频| 亚洲精品乱码久久久v下载方式| 蜜臀久久99精品久久宅男| 午夜福利视频1000在线观看| www.av在线官网国产| 国产精品爽爽va在线观看网站| 观看免费一级毛片| 亚洲自拍偷在线| 亚洲成人久久爱视频| 国产精品1区2区在线观看.| 激情 狠狠 欧美| 啦啦啦观看免费观看视频高清| 国产真实伦视频高清在线观看| 国产精品一区www在线观看| 在线a可以看的网站| 校园人妻丝袜中文字幕| 久久久精品大字幕| 成人国产麻豆网| 日韩一区二区三区影片| 成人亚洲欧美一区二区av| 青春草亚洲视频在线观看| 高清毛片免费看| 97人妻精品一区二区三区麻豆| 国产美女午夜福利| 欧美三级亚洲精品| 亚洲av第一区精品v没综合| 国产午夜精品一二区理论片| 免费不卡的大黄色大毛片视频在线观看 | 亚洲激情五月婷婷啪啪| 能在线免费看毛片的网站| 欧美+日韩+精品| 亚洲内射少妇av| 欧美高清性xxxxhd video| 久久人人精品亚洲av| a级毛片免费高清观看在线播放| 国产精品av视频在线免费观看| 国产老妇女一区| 国产一区二区在线av高清观看| 亚洲国产精品久久男人天堂| 大型黄色视频在线免费观看| 国模一区二区三区四区视频| av黄色大香蕉| 亚洲精品久久久久久婷婷小说 | 国产精品伦人一区二区| 永久网站在线| videossex国产| 看免费成人av毛片| 成人漫画全彩无遮挡| 欧美潮喷喷水| 99热只有精品国产| 亚洲精品456在线播放app| 国产三级中文精品| 老司机影院成人| 国产精品久久久久久久电影| 日日摸夜夜添夜夜爱| 久久亚洲国产成人精品v| 成年女人永久免费观看视频| 晚上一个人看的免费电影| 亚洲欧美日韩无卡精品| 国产中年淑女户外野战色| 国产精品1区2区在线观看.| 99国产精品一区二区蜜桃av| 国产成人一区二区在线| 亚洲精品粉嫩美女一区| 国产成人a区在线观看| 久久久久久伊人网av| 人妻系列 视频| 国产爱豆传媒在线观看| 麻豆国产av国片精品| 欧美日韩一区二区视频在线观看视频在线 | 国产三级在线视频| 亚洲最大成人手机在线| 日本成人三级电影网站| 天美传媒精品一区二区| 免费av观看视频| 又粗又硬又长又爽又黄的视频 | 日韩人妻高清精品专区| 日韩大尺度精品在线看网址| 看黄色毛片网站| 在线免费观看的www视频| 午夜久久久久精精品| 黄片wwwwww| 嘟嘟电影网在线观看| 欧美精品一区二区大全| 国产淫片久久久久久久久| 国产成人freesex在线| 青青草视频在线视频观看| 欧美日本视频| 麻豆国产97在线/欧美| 可以在线观看毛片的网站| 深夜精品福利| 天堂网av新在线| 在线免费观看的www视频| 久久久久久久亚洲中文字幕| 国产国拍精品亚洲av在线观看| 成人欧美大片| 长腿黑丝高跟| 免费一级毛片在线播放高清视频| 一边亲一边摸免费视频| 成人综合一区亚洲| 国产精品一二三区在线看| 日韩国内少妇激情av| 国产爱豆传媒在线观看| 国产人妻一区二区三区在| 91午夜精品亚洲一区二区三区| 熟女人妻精品中文字幕| 永久网站在线| 18禁黄网站禁片免费观看直播| 简卡轻食公司| 国产极品天堂在线| 欧美丝袜亚洲另类| 麻豆久久精品国产亚洲av| 国产精品久久久久久av不卡| 日本熟妇午夜| 亚洲精品国产成人久久av| 免费观看精品视频网站| 国产精品一及| 国产精品永久免费网站| 亚洲一级一片aⅴ在线观看| 国产成人精品一,二区 | 亚洲图色成人| 高清日韩中文字幕在线| 中文在线观看免费www的网站| 色吧在线观看| 亚洲天堂国产精品一区在线| 黑人高潮一二区| 亚洲av成人精品一区久久| 国产亚洲91精品色在线| 精品久久久久久久人妻蜜臀av| 国产视频首页在线观看| 最近最新中文字幕大全电影3| 国产不卡一卡二| 色尼玛亚洲综合影院| 日韩欧美三级三区| 春色校园在线视频观看| 天堂中文最新版在线下载 | 男女下面进入的视频免费午夜| 嘟嘟电影网在线观看| 3wmmmm亚洲av在线观看| 欧洲精品卡2卡3卡4卡5卡区| 久久精品国产清高在天天线| 亚洲在线观看片| 国产精品三级大全| 成人一区二区视频在线观看| 男女做爰动态图高潮gif福利片| 免费观看的影片在线观看| 人妻夜夜爽99麻豆av| 国产精华一区二区三区| 免费av不卡在线播放| 一本久久中文字幕| 12—13女人毛片做爰片一| 国产淫片久久久久久久久| 亚洲精品乱码久久久久久按摩| 干丝袜人妻中文字幕| 婷婷色综合大香蕉| 亚洲久久久久久中文字幕| 丰满人妻一区二区三区视频av| 精品久久久噜噜| 97在线视频观看| 亚洲成av人片在线播放无| 一个人免费在线观看电影| 免费无遮挡裸体视频| 国产精华一区二区三区| 22中文网久久字幕| 久久99热6这里只有精品| 黄片wwwwww| 高清在线视频一区二区三区 | av国产免费在线观看| 性欧美人与动物交配| 黄色配什么色好看| 热99在线观看视频| 国产精品久久电影中文字幕| 日韩精品青青久久久久久| 三级毛片av免费| 91av网一区二区| 尤物成人国产欧美一区二区三区| 亚洲一级一片aⅴ在线观看| 天堂网av新在线| 久久久精品94久久精品| 中文字幕av在线有码专区| 亚洲国产欧美人成| 内射极品少妇av片p| 欧美性感艳星| 亚洲国产精品合色在线| 亚洲欧美日韩卡通动漫| 在线天堂最新版资源| 精品人妻熟女av久视频| 91aial.com中文字幕在线观看| 婷婷色av中文字幕| 永久网站在线| 久久韩国三级中文字幕| 不卡视频在线观看欧美| 久久久久久久久大av| 国产精品久久视频播放| 欧美三级亚洲精品| 亚洲av成人av| 成人午夜精彩视频在线观看| 97在线视频观看| 成人午夜精彩视频在线观看| 欧美日韩乱码在线| 免费av不卡在线播放| av女优亚洲男人天堂| 2022亚洲国产成人精品| 免费看av在线观看网站| www日本黄色视频网| 久久精品国产亚洲av天美| 亚州av有码| 久久久久九九精品影院| 看片在线看免费视频| 日本免费a在线| 欧美日韩综合久久久久久| 亚洲成人久久爱视频| 18禁黄网站禁片免费观看直播| 六月丁香七月| 欧美不卡视频在线免费观看| 久久久久久大精品| 国产成人一区二区在线| 午夜老司机福利剧场| 国产极品天堂在线| 亚洲国产精品国产精品| 国产视频首页在线观看| 高清午夜精品一区二区三区 | 十八禁国产超污无遮挡网站| 人人妻人人看人人澡| 精品国产三级普通话版| 中文精品一卡2卡3卡4更新| 一级黄片播放器| 欧美潮喷喷水| 12—13女人毛片做爰片一| 六月丁香七月| 国产大屁股一区二区在线视频| 啦啦啦韩国在线观看视频| 男人狂女人下面高潮的视频| 国产伦理片在线播放av一区 | 国产精品三级大全| 欧美激情久久久久久爽电影| 好男人在线观看高清免费视频| 亚洲自拍偷在线| 国产精品麻豆人妻色哟哟久久 |