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

        ?

        基于隨機(jī)介質(zhì)正演模擬的塔中順南地區(qū)風(fēng)化殼儲(chǔ)層地震響應(yīng)特征分析

        2013-11-05 06:41:32宋桂橋馬靈偉李宗杰
        石油物探 2013年3期
        關(guān)鍵詞:風(fēng)化殼縫洞孔洞

        宋桂橋,馬靈偉,李宗杰

        (1.中國石油化工股份有限公司油田勘探開發(fā)事業(yè)部,北京 100728;2.中國地質(zhì)大學(xué)(武漢)地球物理與空間信息學(xué)院,湖北武漢 430074;3.中國石油化工股份有限公司西北油田分公司勘探開發(fā)研究院,新疆烏魯木齊 830011)

        塔中順南地區(qū)(順托果勒南區(qū)塊)奧陶系碳酸鹽巖巖溶風(fēng)化殼儲(chǔ)層埋藏深、類型復(fù)雜,非均質(zhì)性強(qiáng),含油氣層段發(fā)育在風(fēng)化殼的頂部。由于風(fēng)化殼界面兩側(cè)存在較大的波阻抗差異,在地震剖面上主要表現(xiàn)為“兩谷夾一峰”的強(qiáng)反射特征,這種強(qiáng)反射掩蓋了風(fēng)化殼頂部儲(chǔ)層發(fā)育的反射特征變化,給風(fēng)化殼儲(chǔ)層地震響應(yīng)特征識(shí)別及儲(chǔ)層預(yù)測帶來很大困難。

        地震波場數(shù)值模擬技術(shù)為研究儲(chǔ)層地震響應(yīng)特征提供了一種方便有效的方法。該方法基于地震、鉆井、測井及地質(zhì)等資料,建立符合地下實(shí)際的地震地質(zhì)模型,通過波動(dòng)方程正演模擬技術(shù)得到復(fù)雜模型條件下的地震波場記錄,在此基礎(chǔ)上,對(duì)地震反射波場進(jìn)行動(dòng)力學(xué)特征及運(yùn)動(dòng)學(xué)特征分析,從而為實(shí)際地震資料波場特征的認(rèn)識(shí)提供一定依據(jù)[1-2]。

        隨機(jī)介質(zhì)建模方法是一種有效刻畫地下介質(zhì)非均質(zhì)性的地質(zhì)建模方法[3-4],能夠?qū)μ妓猁}巖縫洞型儲(chǔ)層的幾何形態(tài)、縫洞的不同尺度及空間變化等統(tǒng)計(jì)特征進(jìn)行很好的描述。我們基于順南地區(qū)測井、巖心測試和地震資料的精細(xì)分析,統(tǒng)計(jì)出奧陶系目標(biāo)儲(chǔ)層的縫洞發(fā)育特征,借助于隨機(jī)介質(zhì)建模方法對(duì)風(fēng)化殼縫洞型儲(chǔ)層進(jìn)行刻畫;在前人研究工作[5-6]的基礎(chǔ)上,重點(diǎn)分析風(fēng)化殼儲(chǔ)層發(fā)育段的地震波場特征變化,優(yōu)選并驗(yàn)證對(duì)風(fēng)化殼儲(chǔ)層特征敏感的地震屬性,為研究區(qū)風(fēng)化殼儲(chǔ)層反射特征識(shí)別及儲(chǔ)層預(yù)測提供參考依據(jù)。

        1 隨機(jī)介質(zhì)建模方法基本原理

        典型的碳酸鹽巖縫洞型儲(chǔ)層由兩部分組成,一是巖石骨架,二是具有一定儲(chǔ)集空間和不同發(fā)育規(guī)模的溶蝕孔洞,如圖1所示。在刻畫此類介質(zhì)時(shí),可將其分為兩部分來描述:一部分是均值不變的地質(zhì)背景,即巖石的骨架;另一部分是具有空間變異特征的隨機(jī)擾動(dòng),即巖石溶蝕孔洞部分,而這個(gè)隨機(jī)擾動(dòng)可以通過隨機(jī)介質(zhì)來描述。

        在二維空間,構(gòu)建某點(diǎn)(x,z)處的巖石物理參數(shù)隨機(jī)模型,如速度v(x,z),并將其分解[7]:

        其中:v0為背景值,設(shè)其為常數(shù)或隨空間坐標(biāo)(x,z)緩慢變化;δv(x,z)是小尺度非均勻介質(zhì)參數(shù),用來刻畫速度空間的變異特征。

        設(shè)

        為空間隨機(jī)擾動(dòng),σ=σ(x,z)是具有零均值、一定自相關(guān)函數(shù)以及方差的空間平穩(wěn)隨機(jī)過程,將(2)式代入(1)式得

        隨機(jī)介質(zhì)在數(shù)學(xué)上可以看作為一個(gè)隨機(jī)序列,一般的隨機(jī)序列不存在理論上的傅里葉變換,這里可由自相關(guān)函數(shù)來描述其傅里葉變換。高斯型、指數(shù)型及Von Karman型自相關(guān)函數(shù)為常用的構(gòu)建隨機(jī)介質(zhì)的自相關(guān)函數(shù)[8]。

        我們基于指數(shù)型橢圓自相關(guān)函數(shù)來構(gòu)建隨機(jī)介質(zhì)自相關(guān)函數(shù),其表達(dá)式為

        其中:a和b分別是自相關(guān)函數(shù)對(duì)應(yīng)的隨機(jī)介質(zhì)在x和z方向上的自相關(guān)長度,描述了隨機(jī)介質(zhì)中擾動(dòng)的平均尺度。

        應(yīng)用截?cái)喾▽?duì)指數(shù)型橢圓自相關(guān)函數(shù)模擬結(jié)果進(jìn)行截取,可以得到滿足一定統(tǒng)計(jì)特征(如孔洞空間尺度、孔洞率φ等)的隨機(jī)孔洞介質(zhì)模型[9-10]。圖2為3個(gè)發(fā)育不同尺度孔洞模型的指數(shù)型橢圓自相關(guān)函數(shù)模擬的結(jié)果,模型長和寬為200m。可以看出,孔洞分布的空間尺度由a,b值的大小決定,孔洞尺度隨a,b值的增大而增大(圖2a,圖2b);孔洞率φ決定儲(chǔ)層內(nèi)孔洞分布的多少,φ值越大,單位面積內(nèi)孔洞的個(gè)數(shù)就越多。

        圖1 典型碳酸鹽巖不同尺度溶蝕孔洞

        圖2 不同空間尺度隨機(jī)介質(zhì)溶蝕孔洞模型

        2 隨機(jī)介質(zhì)模型建立及正演模擬

        2.1 實(shí)際風(fēng)化殼儲(chǔ)層發(fā)育特征

        塔中順南地區(qū)SN1井測井解釋結(jié)果表明,奧陶系風(fēng)化殼表層(一間房組頂部)發(fā)育一段厚度約為8m 的縫洞型儲(chǔ)層,鉆井取心表明裂縫欠發(fā)育,儲(chǔ)集空間以小溶蝕孔洞為主,如表1所示。從時(shí)深關(guān)系上看,儲(chǔ)層發(fā)育位置為3 922ms左右,剛好與研究區(qū)T47反射界面的強(qiáng)波峰反射特征對(duì)應(yīng)(圖3a),而這種強(qiáng)反射界面掩蓋了一間房組頂部儲(chǔ)層的反射特征。

        表1 塔中順南地區(qū)SN1井一間房組(6 522~6 644m)頂部儲(chǔ)層發(fā)育及時(shí)深關(guān)系

        圖3 過SN 1井實(shí)際地震剖面(道號(hào)976)上的 反射界面(a)及其均方根振幅(b)

        2.2 風(fēng)化殼儲(chǔ)層隨機(jī)介質(zhì)模型建立

        縫洞型儲(chǔ)層隨機(jī)介質(zhì)模型的建立大致分為兩步:①依據(jù)地震資料和地質(zhì)沉積相解釋結(jié)果,結(jié)合鉆井、測井等資料統(tǒng)計(jì)研究區(qū)內(nèi)各地層的巖石物理參數(shù),建立符合實(shí)際地下結(jié)構(gòu)的地層框架模型;②在地層框架模型的基礎(chǔ)上,應(yīng)用隨機(jī)建模方法對(duì)目標(biāo)儲(chǔ)層內(nèi)部的溶蝕孔洞進(jìn)行刻畫,即建立隨機(jī)介質(zhì)儲(chǔ)層模型。

        表2給出了研究區(qū)SN1井所揭示的奧陶系各地層的厚度、速度及密度值,恰爾巴克組與一間房組為不整合接觸。隨機(jī)介質(zhì)模型采用圖2c所示的多尺度溶蝕孔洞介質(zhì)模型,圖4為建立的過SN1井縫洞型儲(chǔ)層隨機(jī)介質(zhì)地震地質(zhì)模型。由圖4可見,研究區(qū)風(fēng)化殼表層(一間房組頂部)左半部發(fā)育一套厚度為5~20m 的縫洞型儲(chǔ)層,我們依此來研究風(fēng)化殼表層儲(chǔ)層發(fā)育對(duì)界面反射特征的影響。

        表2 SN1井奧陶系巖石物理參數(shù)

        圖4 過SN 1井地震剖面(道號(hào)976)隨機(jī)介質(zhì)地震地質(zhì)模型

        2.3 隨機(jī)介質(zhì)儲(chǔ)層模型的正演數(shù)值模擬

        常規(guī)的波動(dòng)方程正演數(shù)值模擬方法采用均勻介質(zhì)聲波方程或彈性波方程,難以得到復(fù)雜非均勻孔洞儲(chǔ)層的地震波場響應(yīng)特征。我們采用非均勻介質(zhì)波動(dòng)方程進(jìn)行正演模擬計(jì)算[11-12],從而可以得到接近實(shí)際的地震波場記錄。

        非均勻介質(zhì)波動(dòng)方程正演模擬選用與研究區(qū)實(shí)際地震資料采集接近的觀測系統(tǒng)參數(shù):炮間距50m;道間距50m;排列長度5 950m;滿覆蓋次數(shù)為60次。通過對(duì)實(shí)際地震資料目標(biāo)儲(chǔ)層的頻譜特征進(jìn)行分析,確定數(shù)值模擬所用子波為主頻22Hz的雷克子波。

        3 風(fēng)化殼儲(chǔ)層正演模擬結(jié)果分析

        從厚度為10m 的儲(chǔ)層段與非儲(chǔ)層段正演模擬炮集記錄中各選取一個(gè)地震道進(jìn)行波形對(duì)比可知(圖5a),當(dāng)一間房組頂部發(fā)育儲(chǔ)層時(shí),地震波形的振幅降低,波長變長;對(duì)兩者的頻譜圖(圖5b)進(jìn)行對(duì)比可知,儲(chǔ)層段地震道的頻譜能量減小,主頻降低,高頻部分存在明顯的吸收衰減。

        將隨機(jī)介質(zhì)地震地質(zhì)模型正演模擬炮集記錄進(jìn)行常規(guī)偏移成像處理,結(jié)果如圖6 所示。與圖3a所示實(shí)際地震剖面對(duì)比可見,兩者在T47反射界面均表現(xiàn)為“兩谷夾一峰”的強(qiáng)反射特征,但風(fēng)化殼表層縫洞型儲(chǔ)層的存在使T47界面的反射能量有所變?nèi)酰▓D6中虛線圈所示)。由于儲(chǔ)層發(fā)育的厚度不大,變?nèi)醯姆炔⒉幻黠@,且變?nèi)醯姆扰c儲(chǔ)層發(fā)育情況有關(guān)。

        圖5 儲(chǔ)層段與非儲(chǔ)層段地震道波形(a)及頻譜(b)對(duì)比

        圖6 隨機(jī)介質(zhì)模型正演模擬剖面顯示的T47反射界面

        圖7 基于正演模擬結(jié)果的 界面反射能量隨風(fēng)化殼儲(chǔ)層發(fā)育程度變化定量分析結(jié)果

        圖7為基于正演模擬結(jié)果定量討論的風(fēng)化殼表層儲(chǔ)層發(fā)育(儲(chǔ)層物性、厚度及其空間發(fā)育位置)與T47界面反射能量強(qiáng)弱變化的關(guān)系。圖7a 為10m 厚儲(chǔ)層發(fā)育在距離風(fēng)化殼頂部5m 處,充填不同等效速度時(shí),T47界面反射能量的變化曲線,可以看出,縫洞型儲(chǔ)層等效速度越低(即儲(chǔ)層越發(fā)育),T47界面反射能量減小的幅度越大。圖7b為風(fēng)化殼表層發(fā)育等效速度相同(5 700m/s)而厚度不同(0~80m)的縫洞型儲(chǔ)層時(shí),界面反射能量的變化曲線,可見,隨著風(fēng)化殼頂部儲(chǔ)層厚度的增加,其反射能量逐漸降低,在儲(chǔ)層厚度為50m 左右(λ/4)時(shí)能量減小的幅度最大。圖7c為一間房組頂部以下不同位置發(fā)育等效速度5 700m/s,厚度10m 的儲(chǔ)層時(shí),界面反射能量的變化曲線,可以看出,隨著10m 厚儲(chǔ)層發(fā)育位置與一間房組頂部間的距離增加,界面反射能量先減小后緩慢增加。由此可見,當(dāng)風(fēng)化殼表層發(fā)育儲(chǔ)層時(shí)界面的反射能量降低,且儲(chǔ)層越發(fā)育、儲(chǔ)層的厚度越大,界面反射能量減小的幅度越大。

        4 儲(chǔ)層識(shí)別的敏感地震屬性驗(yàn)證

        基于隨機(jī)介質(zhì)理論建立的風(fēng)化殼縫洞型儲(chǔ)層地震地質(zhì)模型正演模擬結(jié)果表明,當(dāng)風(fēng)化殼表層存在縫洞型儲(chǔ)層時(shí),儲(chǔ)層段相對(duì)于非儲(chǔ)層段使界面的反射能量變?nèi)?,主頻降低,高頻部分存在明顯的吸收衰減。因此,可以借助于對(duì)振幅強(qiáng)度變化敏感的均方根振幅及頻率屬性,對(duì)研究區(qū)界面的反射特征進(jìn)行判識(shí)。

        均方根振幅對(duì)振幅的變化非常敏感,能很好地指示地震剖面上反射波振幅的強(qiáng)弱變化。圖8給出了由實(shí)際地震資料和正演模擬結(jié)果得到的界面均方根振幅剖面,時(shí)窗大小為界面上、下10ms。可以看出,實(shí)際地震資料和正演模擬結(jié)果左側(cè)虛線圈內(nèi)的界面均方根振幅值均明顯降低。由圖4所示的隨機(jī)介質(zhì)地震地質(zhì)模型可知,風(fēng)化殼表層左側(cè)發(fā)育有縫洞型儲(chǔ)層,由此可以證實(shí),圖8a左側(cè)虛線圈內(nèi)反射界面顯示的均方根振幅低異常為縫洞型儲(chǔ)層發(fā)育的標(biāo)記。

        圖9為實(shí)際地震資料與正演模擬結(jié)果單一頻率成分的低頻(18Hz)和高頻(28 Hz)能量體剖面。對(duì)比表明,實(shí)際地震資料及正演模擬結(jié)果的界面低頻(18 Hz)能量體(圖9a,圖9b)強(qiáng)度均高于高頻(28 Hz)能量體(圖9c,圖9d)的強(qiáng)度,整體表現(xiàn)為明顯的低頻特征;就每張剖面上界面的能量體強(qiáng)度變化來說,縫洞型儲(chǔ)層發(fā)育段的能量體強(qiáng)度(圖中虛線圈內(nèi))均低于非儲(chǔ)層發(fā)育段,這種現(xiàn)象在28Hz高頻剖面上更為明顯,可以充分說明儲(chǔ)層段相對(duì)于非儲(chǔ)層段存在更為明顯的高頻吸收衰減。由此可見,除了均方根振幅屬性外,瞬時(shí)頻率屬性亦有助于識(shí)別研究區(qū)奧陶系巖溶風(fēng)化殼縫洞型儲(chǔ)層。

        圖8 實(shí)際地震資料(a)與正演模擬結(jié)果(b)反射界面均方根振幅剖面對(duì)比

        圖9 實(shí)際地震資料與正演模擬結(jié)果的界面分頻能量體剖面對(duì)比

        5 結(jié)束語

        針對(duì)塔中順南地區(qū)奧陶系巖溶風(fēng)化殼儲(chǔ)層地震響應(yīng)特征識(shí)別和儲(chǔ)層預(yù)測的難題,基于隨機(jī)介質(zhì)正演模擬技術(shù),研究了風(fēng)化殼表層儲(chǔ)層發(fā)育時(shí)界面反射特征的變化,定量分析了風(fēng)化殼頂部儲(chǔ)層的物性、厚度及其發(fā)育位置變化對(duì)界面反射能量的影響。結(jié)果表明:當(dāng)風(fēng)化殼頂部發(fā)育儲(chǔ)層時(shí),反射界面相對(duì)于非儲(chǔ)層發(fā)育段表現(xiàn)出不同程度的能量降低、主頻減小,高頻部分存在明顯的吸收衰減;能量降低及頻率吸收衰減的幅度與風(fēng)化殼頂部儲(chǔ)層發(fā)育的程度及厚度有關(guān)。在此基礎(chǔ)上優(yōu)選出均方根振幅及瞬時(shí)頻率為風(fēng)化殼儲(chǔ)層地震響應(yīng)特征識(shí)別相對(duì)敏感的地震屬性,并通過實(shí)際地震資料和正演模擬結(jié)果驗(yàn)證了優(yōu)選屬性的可靠性,為研究區(qū)風(fēng)化殼表層儲(chǔ)層反射特征識(shí)別及儲(chǔ)層預(yù)測提供了一定依據(jù)。

        [1]蔡成國,趙迎月,顧漢明.塔河油田三疊系儲(chǔ)層地震正演模擬及應(yīng)用研究[J].石油天然氣學(xué)報(bào),2008,30(3):92-94 Cai C G,Zhao Y Y,Gu H M.Research and application of seismic forward modeling of Triassic reservoir in Tahe oilfield[J].Journal of Oil and Gas Technology,2008,30(3):92-94

        [2]張雄,石玉,顧漢明.正演模擬技術(shù)在塔河油田于奇地區(qū)的應(yīng)用[J].石油天然氣學(xué)報(bào),2011,33(4):82-87 Zhang X,Shi Y,Gu H M.Application of forward modeling technology in Yuqi area of Tahe oilfield[J].Journal of Oil and Gas Technology,2011,33(4):82-87

        [3]姚姚,奚先.區(qū)域多尺度隨機(jī)介質(zhì)模型及其波場分析[J].石油物探,2004,43(1):1-7 Yao Y,Xi X.Regionalized multi-scale random medium model and its wavefield analysis[J].Geophysical Prospecting for Petroleum,2004,43(1):1-7

        [4]劉春園,朱生旺,魏修成,等.隨機(jī)介質(zhì)地震波正演模擬在碳酸鹽儲(chǔ)層預(yù)測中的應(yīng)用[J].石油物探,2010,49(2):133-139 Liu C Y,Zhu S W,Wei X C,et al.Application of seismic wave forward modeling in random media for carbonate reservoir prediction[J].Geophysical Prospecting for Petroleum,2010,49(2):133-139

        [5]趙迎月,顧漢明,李宗杰,等.塔中地區(qū)奧陶系典型地震體地震識(shí)別模式[J].吉林大學(xué)學(xué)報(bào)(地球科學(xué)版),2010,40(6):1262-1270 Zhao Y Y,Gu H M,Li Z J,et al.The seismic wave field recognition mode of typical geological bodies in the Ordovician in Tazhong area[J].Journal of Jilin University(Earth Seismic Edition),2010,40(6):1262-1270

        [6]張旭光.玉北地區(qū)碳酸鹽巖儲(chǔ)層地震響應(yīng)特征研究[J].石油物探,2012,51(5):493-500 Zhang X G.Study on seismic response characteristics of carbonate reservoir in Yubei area[J].Geophysical prospecting for petroleum,2012,51(5):493-500

        [7]奚先,姚姚.隨機(jī)介質(zhì)模型的模擬與混合型隨機(jī)介質(zhì)[J].地球科學(xué)——中國地質(zhì)大學(xué)學(xué)報(bào),2002,27(1):67-71 Xi X,Yao Y.Simulations of random medium model and intermixed random medium[J].Earth Science(Journal of China University of Geosciences),2002,27(1):67-71

        [8]Ikelle L,Yung S,Daube F.2-D random media with ellipsoidal autocorrelation function[J].Geophysics,1993,58(9):1359-1372

        [9]奚先,姚姚,顧漢明.隨機(jī)溶洞介質(zhì)模型及其波場模擬[J].地球物理學(xué)進(jìn)展,2005,20(2):365-369 Xi X,Yao Y,Gu H M.Random cavity medium model and the wave field simulation[J].Progress in Geophysics,2005,20(2):365-369

        [10]奚先,姚姚,顧漢明.隨機(jī)溶洞介質(zhì)模型的構(gòu)造[J].華中科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2005,33(9):105-108 Xi X,Yao Y,Gu H M.Simulations of random cave medium model[J].Journal of Huazhong University of Science and Technology(Natural Science Edition),2005,33(9):105-108

        [11]姚姚,唐文榜.深層碳酸鹽巖巖溶風(fēng)化殼洞縫型油氣藏可檢測性的理論研究[J].石油地球物理勘探2003,38(6):623-629 Yao Y,Tang W B.Theoretical study of detectable cavern-fractured reservoir in weathered karst of deep carbonatite[J].Oil Geophysical Prospecting,2003,38(6):623-629

        [12]李國平,程利敏,孟令媛,等.非均勻介質(zhì)交錯(cuò)網(wǎng)格高階差分地震波數(shù)值模擬[J].石油天然氣學(xué)報(bào),2010,32(4):220-224 Li G P,Cheng L M,Meng L Y,et al.High-order staggered grid difference seismic wave simulation in heterogeneous media[J].Journal of Oil and Gas Technology,2010,32(4):220-224

        猜你喜歡
        風(fēng)化殼縫洞孔洞
        一種面向孔洞修復(fù)的三角網(wǎng)格復(fù)雜孔洞分割方法
        碳酸鹽巖縫洞儲(chǔ)集體分尺度量化表征
        齊家潛山北段中生界火山巖風(fēng)化殼特征及分布預(yù)測
        錄井工程(2020年2期)2020-07-07 11:40:26
        哈拉哈塘奧陶系縫洞型成巖圈閉及其成因
        孔洞加工工藝的概述及鑒定要點(diǎn)簡析
        收藏界(2019年3期)2019-10-10 03:16:22
        海拉爾盆地W工區(qū)基巖風(fēng)化殼儲(chǔ)層裂縫識(shí)別
        淺析油氣成藏過程中不整合的作用
        玻璃漿料鍵合中的孔洞抑制和微復(fù)合調(diào)控
        云南勐滿紅土風(fēng)化殼特征及其意義
        縫洞型介質(zhì)結(jié)構(gòu)對(duì)非混相氣驅(qū)油采收率的影響
        日本一区二区高清精品| 日本一二三区视频在线| 免费看黑人男阳茎进女阳道视频| 天天爽天天爽天天爽| 国产乱人伦真实精品视频| av资源在线看免费观看| 国产精品爽爽VA吃奶在线观看| 日本在线视频二区一区| 国产亚洲av成人噜噜噜他| 国产三级a三级三级| 无码人妻丰满熟妇区五十路| 少妇丰满大乳被男人揉捏视频| 乱码一二区在线亚洲| 国产艳妇av在线出轨| 一区二区三区婷婷中文字幕| 亚洲av色在线播放一区| 偷拍综合在线视频二区| 国内精品久久久久伊人av| 精品国产午夜理论片不卡| 亚洲中文字幕无码爆乳av| 国产精品丝袜在线不卡| 黄片在线观看大全免费视频| 人妖在线一区二区三区| 亚洲国产婷婷香蕉久久久久久| 丰满少妇呻吟高潮经历| 国产99久久久久久免费看| 在线视频一区色| 亚洲无AV码一区二区三区| 看国产亚洲美女黄色一级片 | 亚洲影院在线观看av| 亚洲av网一区二区三区| 国产欧美一区二区三区在线看| 亚洲欧美综合在线天堂| 无遮挡中文毛片免费观看| 91九色精品日韩内射无| 24小时免费在线观看av| 色欲色欲天天天www亚洲伊| 无码国产精品一区二区免| 中文字幕一区二区三区喷水| 亚洲国产a∨无码中文777| 免费毛片a线观看|