亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        采煤工作面底板水雙頻激電法數(shù)值仿真與探測(cè)試驗(yàn)

        2022-09-20 02:33:00胡雄武張平松吳榮新付茂如
        煤炭學(xué)報(bào) 2022年8期
        關(guān)鍵詞:激電富水電阻率

        胡雄武,徐 標(biāo),張平松,吳榮新,付茂如,程 林

        (1.安徽理工大學(xué) 地球與環(huán)境學(xué)院,安徽 淮南 232001; 2.礦山地質(zhì)災(zāi)害防治與環(huán)境保護(hù)安徽普通高校重點(diǎn)實(shí)驗(yàn)室,安徽 淮南 232001;3.深部煤礦采動(dòng)響應(yīng)與災(zāi)害防控國家重點(diǎn)實(shí)驗(yàn)室,安徽 淮南 232001)

        水害是我國華北地區(qū)深部煤層開采主要災(zāi)害源之一,實(shí)施采場(chǎng)底板巖層富水性探測(cè)是我國水害防治工作中的一項(xiàng)重點(diǎn)內(nèi)容,礦井電法是采場(chǎng)底板巖層富水性探查的主要方法。近40 a來,圍繞采場(chǎng)空間開展的礦井電法先后形成了礦井直流電法、礦井音頻電透視和礦井瞬變電磁法等多個(gè)分支,各分支方法在探測(cè)理論、儀器設(shè)備和數(shù)據(jù)解譯等方面發(fā)展迅速,應(yīng)用效果得到較大程度改善。但實(shí)踐表明,現(xiàn)有礦井電法勘探結(jié)果對(duì)巖層富水性判定尚存在2個(gè)問題。其一,現(xiàn)階段圍繞采場(chǎng)空間開展的礦井電法是通過觀測(cè)電阻率在地層縱橫向上的變化,并以低電阻率區(qū)域?qū)r層富水性進(jìn)行定性或半定量判定。而實(shí)測(cè)電阻率既受巖石礦物成分及結(jié)構(gòu)、孔隙度、含水率、壓力、溫度以及層理等諸多內(nèi)在因素影響,又與采場(chǎng)巷道起伏以及巷道底板導(dǎo)電不均勻性等外在因素有關(guān)。因此,實(shí)測(cè)電阻率影響因素多,對(duì)巖層富水區(qū)判定存在顯著的多解性。其二,煤層采動(dòng)裂隙的發(fā)育范圍一般大于煤層的開采范圍,對(duì)采動(dòng)裂隙發(fā)育范圍內(nèi)的巖層實(shí)施富水性探測(cè)是礦井電法需要滿足的勘探要求。礦井雙巷多角度瞬變電磁測(cè)深法、礦井單巷電測(cè)深法(包括直流電測(cè)深、高密度直流電阻率法等)和雙巷電透視法(包括直流電透視、音頻電透視、雙巷網(wǎng)絡(luò)并行電法等)是采場(chǎng)底板巖層富水性的主要觀測(cè)方法。其中,礦井瞬變電磁法一般將實(shí)測(cè)電阻率數(shù)據(jù)人為放置在線圈法向,不同角度的數(shù)據(jù)通過空間插值進(jìn)行成像,其結(jié)果對(duì)中遠(yuǎn)距離目標(biāo)體有不同程度的放大;礦井單巷電法通常將實(shí)測(cè)電阻率放置在巷道正下方進(jìn)行反演,而事實(shí)上實(shí)測(cè)電阻率是測(cè)量裝置體積范圍內(nèi)的綜合響應(yīng),其電阻率成像所判定的巖層富水區(qū)不能明確是位于巷道沿采煤工作面內(nèi)側(cè)、外側(cè)或是巷道底板,且當(dāng)工作面傾斜寬度較大時(shí),該觀測(cè)方式對(duì)工作面內(nèi)底板富水區(qū)的探測(cè)能力較差;雙巷電透視法是將電阻率反演空間人為限制在工作面內(nèi)部底板進(jìn)行成像,其結(jié)果一般對(duì)采煤工作面內(nèi)部底板巖層的富水區(qū)判定效果較好,但對(duì)巷道底板及其外側(cè)富水區(qū)的判定能力不足。為提高勘探效果,現(xiàn)場(chǎng)常采用單巷測(cè)深和雙巷透視進(jìn)行聯(lián)合觀測(cè),但數(shù)據(jù)卻是分別進(jìn)行處理和解釋,尚未實(shí)現(xiàn)不同觀測(cè)數(shù)據(jù)的有效融合,對(duì)富水區(qū)空間定位精度的提高作用較小。以上反映礦井電法觀測(cè)方式多樣且不同觀測(cè)方式下的數(shù)據(jù)體未有效融合,對(duì)探測(cè)范圍內(nèi)巖層富水區(qū)的空間定位精度低。由此可見,采用以電阻率差異為物理基礎(chǔ)的礦井電法對(duì)采場(chǎng)底板巖層富水性進(jìn)行探測(cè)尚有較大不足,難以實(shí)現(xiàn)精準(zhǔn)探測(cè)目標(biāo),必須深入挖掘巖層富水條件下的電學(xué)特性并給出精準(zhǔn)定位方法。

        激電效應(yīng)是伴隨電阻率測(cè)量過程中附加的一種電學(xué)特性,是由巖礦石及其所含水溶液在電流作用下所發(fā)生的復(fù)雜電化學(xué)活動(dòng)性所致。近年來該方法在礦井地質(zhì)探測(cè)領(lǐng)域發(fā)揮作用。圍繞煤體結(jié)構(gòu),柳蘇、劉耀寧等研究了不同煤體的激電特性。圍繞礦井水害,文獻(xiàn)[21-24]研究了巷道掘進(jìn)前方含水地電異常的雙頻激電場(chǎng)響應(yīng)特征。劉盛東團(tuán)隊(duì)開展了應(yīng)力加載過程中煤系地層不同巖石的時(shí)頻電性特征等基礎(chǔ)研究工作;劉樹才、劉志新等團(tuán)隊(duì)研究了礦井地質(zhì)條件下典型含導(dǎo)水構(gòu)造異常體的時(shí)間域激電響應(yīng)特征,取得了良好的實(shí)際應(yīng)用效果。因激勵(lì)場(chǎng)源性質(zhì)的差異,激電方法發(fā)展過程中先后形成時(shí)間域和頻率域激電法。前期針對(duì)礦井采場(chǎng)空間的大量直流(時(shí)間域)電法探測(cè)實(shí)踐表明,采場(chǎng)巷道地電場(chǎng)噪聲干擾大、電極接地條件差、場(chǎng)源激勵(lì)弱等問題突出,而時(shí)域激電的充放電過程往往需要在較長(zhǎng)時(shí)間內(nèi)向大地供入很強(qiáng)的電流才能明顯觀測(cè)到激電效應(yīng),且其觀測(cè)需要不極化電極,同時(shí)還面臨采場(chǎng)空間工業(yè)游散電流的強(qiáng)干擾問題。因此,受礦井物探裝備本安或防爆要求限制,以及采場(chǎng)大寬面的實(shí)際條件,時(shí)域激電法的適用性受到制約;而頻率域激電法可在供電電流較小情況下,通過選頻和濾波功能克服電極極化不穩(wěn)定、不良接地條件以及工業(yè)游散電流和天然大地電流場(chǎng)等因素的干擾,能夠在較長(zhǎng)傳播距離范圍內(nèi)有效捕捉電性信息且可提供多種激電參數(shù)。以上反映頻率域激電方法在礦井采場(chǎng)空間實(shí)施探測(cè)應(yīng)用相對(duì)有優(yōu)勢(shì)。

        為此,筆者在視電阻率全方位探測(cè)方法基礎(chǔ)上,擬提出煤層底板富水區(qū)雙頻激電全方位探測(cè)方法。結(jié)合淮南礦區(qū)深部煤層頂?shù)装宓刭|(zhì)與地電條件,以采煤面底板灰?guī)r富水區(qū)為研究對(duì)象,通過構(gòu)建不同異常地電模型,采用數(shù)值仿真技術(shù)揭示全方位觀測(cè)方法的激電場(chǎng)響應(yīng)特征,并通過理論與現(xiàn)場(chǎng)實(shí)踐驗(yàn)證該方法的優(yōu)越性。

        1 雙頻激電全方位探測(cè)方法

        1.1 測(cè)試原理

        頻率域激電的物理化學(xué)依據(jù)是巖石電阻率的頻率特性差異,即巖石電阻率隨頻率改變而改變。如果向巖石內(nèi)部分別供幅值相等的高、低頻交變電流,則低頻電流周期長(zhǎng),激電效應(yīng)大;高頻電流周期短,激電效應(yīng)小,2者的差異反映了巖石激電效應(yīng)的強(qiáng)度(圖1)。實(shí)際測(cè)量中,可獲得視復(fù)電阻率和視頻散率2個(gè)激電視參數(shù)。一般認(rèn)為,巖石含水率越大,其視復(fù)電阻率越小,視頻散率越大;反之,視復(fù)電阻率越大,視頻散率越小。這種激電視參數(shù)隨巖石含水率的定性變化關(guān)系是利用雙頻激電法解釋地層富水性的主要依據(jù)。

        圖1 復(fù)電阻率振幅隨頻率的變化關(guān)系

        1.2 觀測(cè)方法

        采煤工作面激電觀測(cè)方法與視電阻率全方位觀測(cè)方法一致,即在工作面上、下巷道底面分別布置沿巷道走向的測(cè)線SL1和SL2,測(cè)線內(nèi)等間距布置64或96個(gè)電極,并分別由交流電法儀控制,儀器之間具有時(shí)間同步功能。形成雙巷電法測(cè)試系統(tǒng)后,依次向SL1和SL2測(cè)線內(nèi)每個(gè)電極供入高低頻組合的交變電流,同時(shí)測(cè)量SL1和SL2測(cè)線內(nèi)非供電電極的交流電位差。

        圖2 單-偶極測(cè)量裝置

        圖3 單-偶極透視測(cè)量裝置

        1.3 激電視參數(shù)反演

        激電反演過程與電阻率反演類似,可先構(gòu)建圓滑模型反演的目標(biāo)函數(shù)()為

        ()=[-()][-()]+α

        (1)

        式中,為實(shí)測(cè)視電阻率(可用低頻視復(fù)電阻率代替)和視頻散率;()為對(duì)應(yīng)的理論值;為權(quán)值矩陣;為模型向量;為圓滑因子;為粗糙度矩陣。

        不同的是,()需對(duì)模型電阻率和頻散率進(jìn)行共同求導(dǎo)。限于篇幅,此處僅闡述模型的頻散率求導(dǎo)過程。一般將地質(zhì)體視為體極化介質(zhì),對(duì)空間進(jìn)行網(wǎng)格劃分后,第點(diǎn)的視頻散率為

        (2)

        S對(duì)第個(gè)模型塊的頻散率求導(dǎo),得

        (3)

        根據(jù)等效電阻率=(1-),可知

        (4)

        式中,為真電阻率;為等效電阻率;為頻散率。

        將式(3)代入式(4),整理可得偏導(dǎo)數(shù)為

        (5)

        由此,可據(jù)式(6)得到模型的修改量,即可得到新的預(yù)測(cè)模型:

        (+α=[-()]

        (6)

        式中,=?()?為Jacobian矩陣。

        經(jīng)過循環(huán)迭代求解,可計(jì)算出新模型理論值,不斷與實(shí)測(cè)數(shù)據(jù)擬合對(duì)比,直到精度滿足要求。

        2 數(shù)值仿真

        2.1 全空間激電視參數(shù)求解

        全空間條件下,穩(wěn)定電流場(chǎng)的邊值問題為

        (7)

        式中,為待求電位;,,為點(diǎn)電流源空間位置;為供電電流;為無窮遠(yuǎn)邊界;為供電電極點(diǎn);為電源點(diǎn)至邊界的距離;為介質(zhì)分界面法線方向。

        對(duì)式(7)進(jìn)行有限元變分求解,并對(duì)模擬空間進(jìn)行網(wǎng)格劃分,可得大型線性方程組

        [][]=[]

        (8)

        式中,為系數(shù)矩陣;為待求電位矩陣;為與場(chǎng)源位置有關(guān)的向量。

        求解式(8),可得各網(wǎng)格節(jié)點(diǎn)的電位。若用Cole-Cole模型來描述巖石的激電效應(yīng),則交流電場(chǎng)作用下巖石的電阻率隨頻率的改變過程可寫為

        (9)

        式中,(i)為介質(zhì)復(fù)電阻率;為零頻電阻率;為充電率;為頻率相關(guān)系數(shù);為時(shí)間常數(shù);為角頻率,=2π。

        將式(9)代入式(7),則式(8)計(jì)算的電位Δ(i)為復(fù)量。全空間PDP裝置形式下,可據(jù)式(10)和(11)分別求得視復(fù)電阻率()和視頻散率(,),其中和分別為高、低頻率,,分別為供電電極至接收電極和之間的距離。

        經(jīng)產(chǎn)蛋白酶活性篩選,共有7株菌有產(chǎn)蛋白酶能力,其中有2株乳酸菌,2株真菌,3株芽孢桿菌,分別編號(hào)為R1,R2,Z1,Z2,Y1,Y2,Y3;經(jīng)產(chǎn)淀粉酶活性篩選,共有3株菌有產(chǎn)淀粉酶能力,且均為具有產(chǎn)蛋白酶能力的芽孢桿菌,即為Y1,Y2,Y3,詳細(xì)結(jié)果見表1。

        (10)

        (11)

        2.2 模型構(gòu)建與參數(shù)設(shè)置

        圖4 研究區(qū)地質(zhì)與地電特征

        由各地層電性對(duì)比可見,煤層頂板為低阻層,底板則呈現(xiàn)出“低阻-高阻”交替互層??紤]到砂巖和灰?guī)r等地層在不含或弱含水條件下,地層對(duì)應(yīng)的充電率、時(shí)間常數(shù)等激電效應(yīng)參數(shù)取值較小,泥巖則相對(duì)略大。結(jié)合鉆孔揭露的地層厚度,在充分考慮不同地層的激電效應(yīng)強(qiáng)度基礎(chǔ)上,合理設(shè)置Cole-Cole模型中描述激電效應(yīng)的相關(guān)參數(shù)(表1),據(jù)此建立背景模型(圖5,其中,模型層序從上至下與表1一致)。

        表1 地層電性參數(shù)

        圖5 研究區(qū)地電模型及觀測(cè)布置

        為模擬煤層底板灰?guī)r富水區(qū)的激電視參數(shù)響應(yīng),令工作面左下角為坐標(biāo)原點(diǎn),沿巷道走向?yàn)檎S,沿工作面傾向和頂板分別為和正軸。設(shè)工作面走向長(zhǎng)度480 m,寬度200 m。分別在巷道1和巷道2底面布設(shè)測(cè)線SL1和SL2,電極分別編為1~96號(hào)和97~192號(hào),相鄰電極距5 m。按前文所述,激電數(shù)據(jù)將由單巷測(cè)深和雙巷透視數(shù)據(jù)合成,其中前者包含不同極距的正向和反向PDP數(shù)據(jù),后者則包含SL1線內(nèi)電極發(fā)射SL2線內(nèi)相鄰電極組接收的數(shù)據(jù)和SL2線內(nèi)電極發(fā)射SL1線內(nèi)相鄰電極組接收的數(shù)據(jù),如以1~96號(hào)電極按順序供高低頻交流電,分別獲得97~192號(hào)電極中相鄰2個(gè)電極間對(duì)應(yīng)頻率的電位差,即Δ(),Δ(),Δ(),…,Δ()和Δ();反之,以97~192號(hào)電極按順序供高低頻交流電,分別獲得1~96號(hào)電極中相鄰2個(gè)電極間對(duì)應(yīng)頻率的電位差,即Δ(),Δ(),Δ(),…,Δ()和Δ()。

        表2 不同圓柱體的幾何參數(shù)

        2.3 激電視參數(shù)響應(yīng)特征

        (1)單巷激電視參數(shù)曲線特征。

        Model-1~4對(duì)應(yīng)不同極距的視復(fù)電阻率及其背景值曲線如圖6所示。

        圖6 不同極距的單巷視復(fù)電阻率剖面曲線

        圖7 不同極距的單巷視頻散率剖面曲線

        (2)雙巷激電視參數(shù)曲線特征。

        透視數(shù)據(jù)中部分發(fā)射點(diǎn)對(duì)應(yīng)的視復(fù)電阻率曲線如圖8所示。

        圖8 不同發(fā)射點(diǎn)的雙巷透視視復(fù)電阻率曲線

        圖9 不同發(fā)射點(diǎn)的雙巷透視視頻散率曲線

        以上表明,各發(fā)射點(diǎn)激電視參數(shù)曲線在圓柱體附近聚焦,總體表現(xiàn)為低阻高頻散特征,發(fā)射點(diǎn)與圓柱體的位置關(guān)系決定視參數(shù)曲線形態(tài)及性質(zhì),接收點(diǎn)與圓柱體的位置關(guān)系決定視參數(shù)異常幅度。雙巷透視數(shù)據(jù)對(duì)面內(nèi)及巷道底板異常的響應(yīng)強(qiáng),對(duì)面外異常的響應(yīng)弱,捕捉能力差。

        2.4 全方位反演分析

        圖10 不同位置模型的全方位數(shù)據(jù)反演結(jié)果

        3 工程實(shí)踐

        某礦1234采煤工作面走向長(zhǎng)800 m,傾斜長(zhǎng)195 m。工作面內(nèi)煤厚平均7.8 m,煤層底界面距下伏厚層灰?guī)r頂界面約28 m,據(jù)已知資料,巖層內(nèi)部裂隙發(fā)育,為主要含水層,嚴(yán)重威脅工作面安全生產(chǎn)?,F(xiàn)場(chǎng)采用雙頻激電全方位探測(cè)方法,測(cè)試時(shí)選用了礦用并行交直流電法儀,沿采煤工作面退遲方向布置了3個(gè)測(cè)站,相鄰測(cè)站之間重疊75 m,共控制工作面走向長(zhǎng)度800 m。每個(gè)測(cè)站由2條測(cè)線構(gòu)成,分別布置在上、下巷道底面,2測(cè)線走向起點(diǎn)和終點(diǎn)相同。其中每測(cè)線內(nèi)布置電極64個(gè),電極間距5 m,控制測(cè)線315 m。在每站數(shù)據(jù)采集時(shí),上、下巷道測(cè)線分別由1臺(tái)電法儀通過時(shí)間約定的方式同步實(shí)施16 Hz和128 Hz的交流供電與電位數(shù)據(jù)采集,即上、下巷道內(nèi)每個(gè)電極供電時(shí),同時(shí)完成上、下巷道所有電極的數(shù)據(jù)接收(供電電極除外)。

        對(duì)上、下巷道3個(gè)測(cè)站數(shù)據(jù)進(jìn)行解編,可提取單巷測(cè)深和雙巷透視激電數(shù)據(jù),將2者合成全方位數(shù)據(jù)并通過三維反演,進(jìn)一步獲得工作面底板巖層的反演電阻率和頻散率立體結(jié)果。考慮到該工作面主要灰?guī)r含水層位于煤層底板30 m以下,此處從全方位反演結(jié)果中提取深度35 m的切片數(shù)據(jù)進(jìn)行成像。

        圖12(a)為煤層底板35 m深度反演電阻率切片。可見,反演電阻率存在3處較顯著的低阻異常區(qū),分別為巷道2一側(cè)的0~450 m段、550~780 m段以及巷道1外側(cè)的200~550 m段。圖12(b)為煤底板35 m深度反演頻散率切片。圖12(b)顯示存在1處高頻散異常區(qū),主要分布在巷道2一側(cè)的80~260 m段,沿工作面內(nèi)外有半徑約為60 m的影響范圍;其他區(qū)域頻散率值則相對(duì)偏低。對(duì)照?qǐng)D12,并結(jié)合已知資料可知,低值異常區(qū)內(nèi)均分布有落差不等的斷層,分析斷層帶內(nèi)裂隙發(fā)育并伴隨巖性改變,是引起低阻異常的主要原因。相比而言,頻散率對(duì)巖性改變不敏感,僅在其中的一處低阻異常區(qū)內(nèi)有明顯的高值異常,且相對(duì)聚焦。按照巖層富水區(qū)同時(shí)滿足低阻高頻散的判定原則,解釋高頻散異常區(qū)為該工作面底板灰?guī)r層位的核心富水區(qū)。為驗(yàn)證探測(cè)結(jié)果的準(zhǔn)確性,礦方針對(duì)工作面不同區(qū)段實(shí)施了多個(gè)鉆孔,其探查結(jié)果反映在解釋的低阻高頻散異常區(qū)內(nèi)鉆孔出水量約18 m/h,其他區(qū)域鉆孔出水量明顯偏小,是局部巖層裂隙發(fā)育并少量含水引起。綜上可知,通過對(duì)采煤工作面雙頻激電場(chǎng)的全方位觀測(cè)與反演,基于激電視參數(shù)反演結(jié)果的綜合解釋,可以有效避免電阻率的多解性問題,并實(shí)現(xiàn)對(duì)工作面底板富水異常區(qū)的高精度定位。

        圖11 不同位置模型的雙巷透視數(shù)據(jù)反演結(jié)果

        圖12 工作面底板35 m深度激電視參數(shù)擬斷面

        4 結(jié) 論

        (1)鑒于現(xiàn)有以電阻率為物理基礎(chǔ)的礦井電法對(duì)煤層工作面底板水害探測(cè)存在多解性強(qiáng)及定位不準(zhǔn)的問題,提出了采煤工作面雙頻激電全方位探測(cè)方法。

        (2)基于不同異常地電模型的數(shù)值模擬與理論分析,揭示了激電視參數(shù)響應(yīng)特征:① 單巷測(cè)深激電視參數(shù)表現(xiàn)低阻高頻散異常,利用視參數(shù)曲線“交點(diǎn)”可準(zhǔn)確定位異常體的走向位置。② 雙巷透視激電視參數(shù)曲線在異常體附近聚焦,表現(xiàn)低阻高頻散異常;發(fā)射點(diǎn)與異常體的位置關(guān)系決定視參數(shù)曲線形態(tài)及性質(zhì);接收點(diǎn)與異常體的位置關(guān)系決定異常幅度;雙巷透視數(shù)據(jù)對(duì)面內(nèi)及巷道底板異常的響應(yīng)強(qiáng),對(duì)面外異常的響應(yīng)弱,捕捉能力差。

        (3)數(shù)值模型激電數(shù)據(jù)的反演結(jié)果表明,與雙巷透視數(shù)據(jù)相比,全方位數(shù)據(jù)能實(shí)現(xiàn)單巷測(cè)深和雙巷透視數(shù)據(jù)的約束反演,使異常體得以準(zhǔn)確歸位,體現(xiàn)了全方位激電探測(cè)方法的優(yōu)越性。工程實(shí)踐進(jìn)一步表明,該方法能有效降低電性異常的多解性并提高空間定位精度,改善了煤層底板巖層富水區(qū)的判定效果。

        猜你喜歡
        激電富水電阻率
        大功率激電測(cè)深方法在豫西董家埝銀礦床勘查中的應(yīng)用
        高頻大地電磁測(cè)深與激電中梯在金礦勘查中的應(yīng)用研究
        大功率激電測(cè)量在冀北溫家營(yíng)—馬家溝銀多金屬礦勘查中的應(yīng)用
        激電聯(lián)合剖面在判斷矽卡巖型礦床礦體產(chǎn)狀中的應(yīng)用
        富水砂卵石地層錨索施工工藝工法
        超級(jí)高密度電法在新疆某礦區(qū)富水性評(píng)價(jià)應(yīng)用研究
        三維電阻率成像與高聚物注漿在水閘加固中的應(yīng)用
        掘進(jìn)巷道構(gòu)造富水性電法探測(cè)綜合應(yīng)用研究
        隨鉆電阻率測(cè)井的固定探測(cè)深度合成方法
        海洋可控源電磁場(chǎng)視電阻率計(jì)算方法
        国产边摸边吃奶叫床视频| 亚洲深深色噜噜狠狠网站| 亚洲高清中文字幕视频| 制服丝袜一区二区三区| 无遮挡18禁啪啪羞羞漫画| 欧美寡妇xxxx黑人猛交| 成av人片一区二区三区久久| 婷婷四房播播| 日韩国产成人精品视频| 国产成年无码v片在线| 香蕉久久人人97超碰caoproen| 在线视频你懂的国产福利| 精品人伦一区二区三区蜜桃麻豆| 白白视频在线免费观看| 久久久人妻丰满熟妇av蜜臀| 精品私密av一区二区三区| 成年人观看视频在线播放| 免费av一区二区三区无码| 亚洲色欲色欲大片www无码| 欧美激情区| 久久久亚洲精品午夜福利| 亚洲精品一区二区三区国产| 偷拍视频这里只有精品| 白白发在线视频免费观看2| 青青草原亚洲| 欧美精品中文字幕亚洲专区| 亚洲tv精品一区二区三区| 亚洲国产精品av麻豆一区| 蜜桃av噜噜一区二区三区9| 无码一区二区三区亚洲人妻| 99久久综合狠狠综合久久 | 国产情侣一区在线| 丰满少妇一区二区三区专区| 痴汉电车中文字幕在线| 无码一区二区三区在线| 自拍偷自拍亚洲精品播放| 日本a在线播放| 精品国产亚洲av高清日韩专区| 亚洲中文字幕久久精品品| 亚洲乱色伦图片区小说| 亚洲av无码精品色午夜蛋壳|