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

    三維高斯射線束觀測系統(tǒng)照明及優(yōu)化方法研究

    2015-06-27 05:54:48殷厚成
    石油物探 2015年4期
    關鍵詞:炮點檢波波場

    殷厚成,鄧 飛

    (1.中國石油化工股份有限公司石油物探技術研究院,江蘇南京211103;2.成都理工大學信息科學與技術學院,四川成都610059)

    三維高斯射線束觀測系統(tǒng)照明及優(yōu)化方法研究

    殷厚成1,鄧 飛2

    (1.中國石油化工股份有限公司石油物探技術研究院,江蘇南京211103;2.成都理工大學信息科學與技術學院,四川成都610059)

    基于模型的觀測系統(tǒng)照明分析技術在復雜地區(qū)的觀測系統(tǒng)設計中起著非常重要的作用,但隨著野外勘探復雜程度的增加,三維地質模型也更加復雜,常用的射線方法在精度上已經(jīng)不能很好地滿足實際生產的需求,波動方程類方法則由于計算效率較低在實際使用中存在很大局限性。提出了一種基于高斯射線束的雙向照明方法,給出了方法的實現(xiàn)步驟。該方法將波場分解到具有一定寬度的射線束上實現(xiàn)波場的模擬和延拓,通過檢波點射線坐標反向變換,根據(jù)互易性原理實現(xiàn)觀測系統(tǒng)的雙向照明,不僅大幅提高了計算效率,而且有效改善了照明的精度。根據(jù)高斯射線束照明結果分析和統(tǒng)計炮點對各目的層面元的雙向照明貢獻,優(yōu)選對低能量目的層區(qū)域貢獻大的炮點附近進行補炮,從而實現(xiàn)觀測系統(tǒng)的局部優(yōu)化。采用一個三維地質模型驗證了高斯射線束雙向照明算法的合理性和正確性。

    高斯射線束;射線束照明;雙向照明;觀測系統(tǒng)設計;觀測系統(tǒng)優(yōu)化

    傳統(tǒng)的地震勘探基于水平疊加理論,其前提是假定地下為水平層狀介質;但隨著我國油氣勘探的不斷深入和發(fā)展,探區(qū)的地質目標日益復雜,基于水平層狀假設的常規(guī)觀測系統(tǒng)設計不能滿足勘探需求,尤其是在起伏地表和復雜構造地區(qū)。以地質模型為基礎的照明分析[1]技術是一項對地震波能量分布進行定量分析的技術。通過波場數(shù)值模擬照明分析,可以研究震源及檢波器對目標地質體的照明及接收能量,根據(jù)觀測系統(tǒng)照明能量分布,有效指導觀測系統(tǒng)設計。

    目前常用的觀測系統(tǒng)照明方法主要包括射線法和波動方程法兩類。射線法[2-3]因效率較高很早就被應用于觀測系統(tǒng)設計,其中兩點射線追蹤[4]可用于計算觀測系統(tǒng)對目的層面元的覆蓋次數(shù),以覆蓋次數(shù)近似替代能量來進行照明分析。該方法具有高效、靈活的特點,但存在射線陰影區(qū)、焦散區(qū)等問題,而且以覆蓋次數(shù)近似替代能量也存在一定的精度問題。波動方程類方法包括單程波[5]和全程波[6]兩種,與射線法相比,這類方法通過直接求解波動方程進行波場延拓,解決了多波至以及由速度變化引起的聚焦或焦散效應問題,數(shù)值計算精度高。但是波動方程類方法對硬件設備要求較高,計算效率低,因而在實際應用中受到較多限制。

    射線束方法是一類介于射線和波動方程之間的照明分析方法,高斯射線束[7-9]是射線束方法的代表,它將波場分解到具有一定寬度的射線束上來實現(xiàn)波場的模擬和延拓,不僅具有運動學特征,而且具有動力學特征,可以精確計算中心射線附近波場的振幅和相位,在一定程度上解決了射線的盲區(qū)、焦散問題,并能夠自然地實現(xiàn)多波至。然而,由于高斯射線束自身的特點,直接將高斯射線束方法套用到已有的單程波照明框架上無法將檢波點記錄整體反向延拓,而逐檢波點延拓的效率很低,無法實現(xiàn)高效的雙向照明計算。為了提高高斯射線束的效率,本文在高斯射線束波場延拓的基礎上,通過檢波點射線坐標反向變換,根據(jù)互易性原理實現(xiàn)了快速的高斯射線束雙向照明。

    1 高斯射線束觀測系統(tǒng)照明及優(yōu)化方法

    觀測系統(tǒng)照明實質上就是利用數(shù)值模擬方法研究觀測系統(tǒng)在指定地質模型上的照明能量強度分布,因此現(xiàn)有的地震波場數(shù)值模擬方法均可以應用于照明分析。射線方法是一種最早被人們應用于地震波場數(shù)值模擬的經(jīng)典方法,該方法根據(jù)地質模型的速度變化追蹤地震波在地下介質中的運動軌跡及其在地質界面上的反射和透射變化,算法簡單且效率較高,但其主要反映地震波的運動學特征,不能很好地表現(xiàn)動力學特點。

    傳統(tǒng)的射線方法是將波場分解到單一的射線上來實現(xiàn)波場數(shù)值模擬的,獲取高精度的合成記錄需要大量精確的射線。利用試射法進行的兩點射線追蹤需要耗費大量的時間,而且由于射線的局限性很難有效地解決射線盲區(qū)、多波至等問題。分析認為,傳統(tǒng)射線法的這些缺陷是由于將波場分解到完全理想的高頻射線上造成的,如果將波場分解到一定高頻范圍內的射線能量帶上,就可以從根本上改進射線法的效果。高斯射線束方法就是一種利用“胖射線”思想進行波場數(shù)值模擬的方法,地下介質或檢波點處的波場不再由單一的一條射線決定,而是由多條射線(束)疊加獲得。射線束方法的特性決定了它可以在很大程度上克服傳統(tǒng)射線方法的缺點,因此具有很高的計算效率和較高的精度。

    1.1 三維高斯射線束方法原理

    高斯射線束方法最初由Cerveny等[7-9]在20世紀80年代提出并用于地震記錄正演,后來又被Hill[10-11]用于偏移成像。在高斯射線束數(shù)值模擬中,高斯射線束是波動方程集中于射線附近的高頻漸近解[7],它可以被看作是一條從震源出發(fā)以射線為中心的能量管,射線束的振幅以偏離中心射線的距離呈指數(shù)衰減,因類似于高斯分布而得名,而地下介質中某處或檢波點處的波場則由一定范圍內的多條高斯射線(束)疊加形成。完整的高斯射線束正演算法包括運動學射線追蹤、動力學射線追蹤和波場疊加三個步驟[12]:通過運動學射線追蹤獲得中心射線軌跡和中心射線能量的變化;通過動力學射線追蹤沿中心射線計算高斯射線束的動力學參數(shù),確定高斯射線束的振幅衰減和波前曲率;波場疊加則用于對計算點有貢獻的高斯射線束疊加并形成最終的波場記錄。

    高斯射線束建立在射線坐標系下,三維情況下的射線坐標系如圖1所示,S為中心射線,S附近有一點P。過P點作垂直于射線的平面,與射線交于P′點,P′到起點S0的射線路徑長度為s。n,m為P點垂直于射線的平面內的二維笛卡爾坐標,P點相對于中心射線S的射線坐標記為(s,n,m)。從圖1可以看出,射線坐標系的基矢量由遵循右手坐標系的3個正交的單位矢量et,en,em構成。其中et為射線的單位正切矢量,而en和em垂直于射線,并且保持射線坐標系是正則的。

    圖1 三維射線坐標系

    設在射線起點S0處已知射線出射方向為et0,它與z軸夾角為α,其xoy平面投影與x軸夾角為β,那么et0可表示為:

    (1)

    由于en,em和et是用于確定射線坐標系的單位向量,它們互相垂直,因此可以選擇

    (2)

    當射線在變速介質中前進時,et會發(fā)生變化,en,em也要隨之變化,它們滿足如下微分方程:

    (3)

    在上述射線坐標系下,三維高斯射線束可以表示為:

    (4)

    (5)

    (5)式即為三維高斯射線束動力學方程組,式中v(s)是速度關于局部坐標en,em的2×2二階偏導數(shù)矩陣,矩陣的各個元素表示為:

    (6)

    復值動力學矩陣P(s)/Q(s)決定了高斯射線束的特征,其實部特征值Re[P(s)/Q(s)]決定了高斯射線束在射線上的相前曲率,而其虛部則決定了高斯射線束在垂直于射線截面上的振幅分布,并且三維高斯射線束的半寬度矩陣可以表示為:

    (7)

    L(s)的兩個特征值分別決定了在垂直于射線的平面上呈橢圓分布的振幅長、短軸。

    1.2 基于高斯射線束的雙向照明

    在震源點處將波場分解到一系列的高斯射線束中,利用高斯射線束公式(4),可以將波場延拓到模型的任意位置處,計算地下介質或檢波點處的地震波場并統(tǒng)計該處的照明能量。然而單純通過目的層的入射或者檢波點的接收照明能量并不能直觀地判斷觀測系統(tǒng)的優(yōu)劣,觀測系統(tǒng)照明需要綜合考慮震源和檢波器排列的綜合效應,原因就像文獻[13]中形象比喻的那樣:“在黑暗中用手電筒照亮一個物體,視力好的人可以看到該物體,但即使物體被照得再亮,也無法被盲人看見”。一方面,即使某目的層入射能量再高,若無法被觀測系統(tǒng)接收,那么之后的資料處理也無法使該目的層成像;另一方面,檢波器排列接收的能量再高,若沒有接收到有效目的層的能量,那么成像結果也是無效的。

    評判觀測系統(tǒng)優(yōu)劣的照明能量應能綜合考慮震源和檢波器排列,而雙向照明能量正是這樣一種能量,它是經(jīng)過目的層反射且被檢波點接收的能量。很明顯,這是可以用于目的層成像的有效照明能量。雙向照明最早由Xie等在文獻[5]中提出,朱金平等在文獻[13]中給出了一種效率較高的單程波雙向照明方法DUC-DC。該方法分為兩步:①首先從震源點向下延拓波場至模型底部,然后再將反射波場向上延拓回地表,得到檢波點處的波場,這一步相當于正演;②將檢波點處的波場再次向下延拓至地下介質,獲得雙向照明強度。

    高斯射線束雙向照明的計算,可以采用與單程波方法類似的框架,模仿DUC-DC方法進行。首先利用高斯射線束進行正演模擬計算,獲得各檢波點處的正演記錄;然后逐檢波點使用高斯射線束將波場進行反向延拓,得到地下介質的照明強度。與單程波方法相比不同之處在于,單程波方法的反向延拓可以將檢波點處的記錄作為整體向下延拓,而高斯射線束并不能將多個檢波點處的記錄作為整體一次性向下延拓,必須逐檢波點進行。假設觀測系統(tǒng)中有Ns個震源,每個震源平均影響排列中的Ng個檢波點,那么采用DUC-DC方法的雙向照明框架,需要進行Ns×(1+Ng)次高斯射線束延拓,計算效率很低。

    為了提高高斯射線束雙向照明的效率,應避免從檢波點處進行反向延拓。雙向照明的目的是求出被檢波點接收的面元反射能量,因此,如果能夠充分利用正演時的射線束直接求出雙向照明能量,則將大幅提高計算效率。為此,考查一條從震源出發(fā)的射線束Bi,如圖2所示,該射線束經(jīng)反射面反射回到地表被檢波點G所接收。設反射面上有一點C,求C點關于射線束Bi對檢波點G的雙向照明能量,即求出射線束Bi經(jīng)C點反射被檢波點G所接收的能量。

    圖2 高斯射線束雙向照明計算

    對于給定的射線束Bi,首先根據(jù)高斯射線束公式(4)求出檢波器G點相對于中心射線R的射線坐標(sG,nG,mG),并求出G的振幅A(G)和照明能量。根據(jù)射線束原理,G點接收的能量是由反射面上射線束有效范圍內的反射點共同決定的。由于C點的射線坐標為(sC,nC,mC),而G點的射線坐標是(sG,nG,mG),為了求出反射面上C點對G點的貢獻,可以沿中心射線將G點從sG位置反向延拓至sC位置,記為G′,求出G′點及其射線坐標(sC,nG′,mG′)。進行反向延拓時G′點的振幅等于G點處振幅A(G),如果此時將G′點作為中心射線位置,可求出C點相對于G′點的平面內二維笛卡爾坐標:

    (8)

    將其代入(4)式即可求出G′點對C點的振幅貢獻A(C)。根據(jù)互易性原理,A(C)即為C點對檢波點G的雙向照明振幅。對上述算法思路進行整理,可以得到高斯射線束單炮雙向照明的算法步驟:

    1) 確定射線角度范圍和射線角度間隔;

    2) 遍歷所有的射線角度,重復步驟3)至9);

    3) 按照射線出射角度進行運動學和動力學射線追蹤,確定中心射線路徑R及動力學參數(shù)矩陣P,Q;

    4) 根據(jù)中心射線R,計算反射面元點相對于中心射線的射線坐標,選取高斯射線束半寬度范圍內的面元點,作為計算的有效面元點,記為集合{C};

    5) 根據(jù)中心射線R在地表的出射位置,計算檢波點的射線坐標,選取半寬度范圍內的檢波點,作為有效檢波點,記為集合{G};

    6) 遍歷有效檢波點集合{G},重復步驟7)至9);

    7) 根據(jù)檢波點Gi的射線坐標,按照公式(4)計算檢波點振幅A(Gi);

    8) 遍歷有效面元點集合{C},重復步驟9);

    9) 將檢波點Gi的射線坐標反向延拓到面元點Cj所在的射線坐標系下,得到(SCj,nGi′,mGi′),利用公式(8)和公式(4)求得雙向照明振幅A(i,j)和能量E(i,j),并將求得的照明能量累計到面元Cj上。

    使用上述算法在正演計算后無需逐檢波點進行反向延拓獲取雙向照明強度,只需在射線束正演的過程中增加檢波點射線坐標反向延拓和面元點能量計算的工作。由于正演過程本身就要計算和記錄射線坐標系的變換矩陣,因此求檢波點的反向延拓坐標僅需要乘以變換的逆矩陣即可。使用該方法每炮僅需要做1次射線束延拓運算,Ns個震源僅需進行Ns次射線束延拓運算,計算效率相對于之前的逐檢波點延拓提高了Ns×Ng倍。

    1.3 基于照明結果的觀測系統(tǒng)優(yōu)化

    高斯射線束照明具有很高的計算效率,能夠完成上萬甚至數(shù)十萬炮的雙向照明計算,因此可以利用高斯射線束進行觀測系統(tǒng)論證。首先建立待勘探工區(qū)的三維地質模型,利用常規(guī)觀測系統(tǒng)設計軟件設計幾個備選觀測系統(tǒng),然后利用高斯射線束對不同觀測系統(tǒng)下的地質模型進行雙向照明計算,挑選出目的層照明均勻、能量較高的觀測系統(tǒng)作為優(yōu)選觀測系統(tǒng)。

    由于地表起伏或地下構造復雜,優(yōu)選出的觀測系統(tǒng)對目的層的照明仍然可能存在不均勻和能量較低的暗區(qū),通常需要補充一些炮點以得到優(yōu)化的觀測系統(tǒng),增加低能量區(qū)的照明強度。問題是如何確定合理的加炮區(qū)域,使得增設的炮點對待加強的低能量目的層區(qū)域產生最大化貢獻。為此,可以在照明計算時記錄炮點對各目的層面元的雙向照明貢獻,當選擇出低能量區(qū)域后,查詢并統(tǒng)計出地表炮點對選定區(qū)域的貢獻,優(yōu)先選擇貢獻大的炮點附近作為補炮的候選區(qū)域。在候選區(qū)域補充炮點之后再次計算并分析目的層雙向照明能量,如果達到預期效果則輸出優(yōu)化后的觀測系統(tǒng),否則按照上述方法補充更多的炮點。

    2 模型算例分析

    為了驗證高斯射線束雙向照明算法的合理性和正確性,采用一個三維地質模型進行了試算(圖3)。模型在X,Y,Z三個方向上的尺寸分別為32000,18000,15800m,包含6個層位,4個斷層。各層速度由上至下分別為:5625,6100,6300,4800,5400,6000m/s。

    給定2個備選觀測系統(tǒng)Ⅰ和Ⅱ。觀測系統(tǒng)Ⅰ的基本參數(shù)如下:20線,每線300道,道間距40m,線間距200m,共11500炮。觀測系統(tǒng)Ⅱ的基本參數(shù)如下:28線,每線288道,道間距40m,線間距240m,共10692炮。測試用計算機為聯(lián)想T440p,CPU為i7-4700MQ/2.4GHz,內存4.0GB,64位Win7操作系統(tǒng)。選擇模型第3層為目的層進行高斯射線束雙向照明計算,每炮的單線程平均計算時間為1978ms,使用4線程同時計算,觀測系統(tǒng)Ⅰ照明共耗時87min,觀測系統(tǒng)Ⅱ照明共耗時89min。

    圖4是觀測系統(tǒng)Ⅱ中(13950m,9750m)處炮點的雙向照明效果,可以看到,由于地表起伏和斷層割裂,目的層的照明并沒有呈現(xiàn)出圓形,主要照明能量被分割為明顯的兩塊。圖5對比了觀測系統(tǒng)Ⅰ和觀測系統(tǒng)Ⅱ的雙向照明效果,可以看到觀測系統(tǒng)Ⅱ的照明能量、覆蓋范圍和均勻度均高于觀測系統(tǒng)Ⅰ,因此選擇觀測系統(tǒng)Ⅱ進行進一步優(yōu)化。

    圖3 三維地質模型

    圖4 單炮雙向照明效果

    圖5 觀測系統(tǒng)雙向照明結果對比

    選擇目的層上的低能量區(qū)域,如圖6a所示;統(tǒng)計炮點對選中區(qū)域的照明貢獻并繪制在地表上,如圖6b所示。按照貢獻大小在原炮點附近半炮距位置處增設100個新炮點,增加新炮點后的照明效果如圖6c所示,可以看到低能量區(qū)域的照明效果得到了明顯改善。應用相同方法對觀測系統(tǒng)Ⅱ繼續(xù)進行優(yōu)化,補入650炮點后的最終照明效果如圖6d 所示,與圖5b相比照明效果有顯著提高。

    圖6 觀測系統(tǒng)Ⅱ優(yōu)化

    3 結束語

    照明分析是觀測系統(tǒng)設計與優(yōu)化的重要手段,基于高斯射線束的雙向照明方法通過對原有射線束正演方法進行擴展,僅需要增加少量工作即可完成雙向照明計算,既能保持射線法高效、靈活的特點,又能提高射線法計算的精度,照明能量強弱過渡自然,效果近似于波動方程類方法,能夠很好地應用于復雜三維地質模型的照明分析。

    利用高斯射線束照明結果可以實現(xiàn)觀測系統(tǒng)優(yōu)選。通過對照明結果的統(tǒng)計和分析,可以自動確定目的層局部低能量區(qū),從而確定優(yōu)勢補炮區(qū)域。通過在推薦區(qū)域補炮,可以顯著提高低能量區(qū)的照明效果。

    由于高斯射線束雙向照明方法具有很高的計算效率,因此能夠在普通微機上完成大規(guī)模三維觀測系統(tǒng)的照明計算和分析,滿足野外觀測系統(tǒng)設計的實際需求。

    [1] Berkhout A J,Ongkiehong L,Volker A W F,et al.Comprehensive assessment of seismic acquisition geometries by focal beams,part Ⅰ:theoretical considerations[J].Geophysics,2001,66(3):911-917

    [2] Bear G,Lu C P,Lu R,et al.The construction of subsurface illumination and amplitude maps via ray tracing[J].The Leading Edge,2000,19(7):726-728

    [3] Hoffmann J.Illumination,resolution,and image quality of PP-and PS-waves for survey planning[J].The Leading Edge,2001,20(9):1008-1014

    [4] 徐濤,徐果明,高爾根,等.三維復雜介質的塊狀建模和試射射線追蹤[J].地球物理學報,2004,47(6):1118-1126 Xu T,Xu G M,Gao E G,et al.Block modeling and shooting ray tracing in complex 3D media[J].Chinese Journal of Geophysics,2004,47(6):1118-1126

    [5] Xie X B,Jin S W,Wu R S.Wave-equation-based seismic illumination analysis[J].Geophysics,2006,71(5):169-177

    [6] 陳生昌,馬在田,吳如山.波動方程雙程地下方向照明分析[J].同濟大學學報(自然科學版),2007,35(5):682-684 Chen S C,Ma Z T,Wu R S.Two-way subsurface directional illumination analysis by wave equation[J].Journal of Tongji University(Natural Science),2007,35(5):682-684

    [7] Cerveny V,Psencik I.Gaussian beams and paraxial ray approximation in three-dimensional elastic in homogeneous media[J].Geophysical Journal Royal Astronomical Society,1983,53(1):1-15

    [8] Cerveny V.Gaussian beams synthetic seismograms[J].Journal of Geophysics,1985,58(1):44-72

    [9] Cerveny V.Ray synthetic seismograms for complextwo-dimensional and three-dimensional structures[J].Journal of Geophysics,1985,58(1):2-26

    [10] Hill N R.Gaussian beam migration[J].Geophysics,1990,55(11):1416-1428

    [11] Hill N R.Prestack Gaussian-beam depth migration[J].Geophysics,2001,66(4):1240-1250

    [12] 鄧飛,劉超穎.三維射線快速追蹤及高斯射線束正演[J].石油地球物理勘探,2009,44(2):158-165 Deng F,Liu C Y.3D rapid ray-tracing and Gaussian ray-beam forward simulation[J].Oil Geophysical Prospecting,2009,44(2):158-165

    [13] 朱金平,董良國.地震波雙向照明的概念及計算方法[J].地球物理學報,2011,54(11):2933-2942 Zhu J P,Dong L G.The concept and calculation method of bi-directional seismic illumination[J].Chinese Journal of Geophysics,2011,54(11):2933-2942

    (編輯:戴春秋)

    Research on seismic acquisition geometry illumination and its optimization based on 3D Gaussian beam

    Yin Houcheng1,Deng Fei2

    (1.SinopecGeophysicalResearchInstitute,Nanjing211103,China;2.CollegeofInformationScienceandTechnology,ChengduUniversityofTechnology,Chengdu610059,China)

    The seismic geometry illumination technology based on geological model plays an important role on seismic exploration in complex area.However,since the field exploration becomes more and more complicated,3-D geologic model becomes more complex and the conventional ray method isn’t able to meet the need of actual production.On the other hand,the wave equation method has many limitations in practical application because of low efficiency.Using Gaussian beams to carry out wavefield forward simulation and continuation,Gaussian beam bi-directional illumination method not only keeps the high-efficiency and flexibility but also overcomes some shortcomings of conventional ray method,effectively solves the problem of rapid illumination on complex 3-D geologic model.According to Gaussian beam illumination result we can choose the best seismic geometry,and realize local optimization based on analysis and statistics of the result.

    Gaussian beam,ray beam illumination,bi-directional illumination,seismic geometry design,seismic geometry optimization

    2015-02-27;改回日期:2015-05-15。

    殷厚成(1963—),男,高級工程師,長期從事地震采集技術攻關與研究工作。

    國家科技重大專項“煤層氣地震采集和處理技術”項目(2011ZX05035-002)資助。

    P631

    A

    1000-1441(2015)04-0376-06

    10.3969/j.issn.1000-1441.2015.04.002

    猜你喜歡
    炮點檢波波場
    一種實時頻譜儀中幀檢波器的FPGA 實現(xiàn)
    基于最小炮檢距道快速檢測炮點偏移方法
    彈性波波場分離方法對比及其在逆時偏移成像中的應用
    GSM-R系統(tǒng)場強測試檢波方式對比研究
    對二維地震覆蓋次數(shù)和炮點距的思考
    山西冶金(2018年6期)2018-03-04 12:53:26
    無樁號施工中炮點COG現(xiàn)場快速偏移技術
    科學與財富(2018年1期)2018-03-03 05:12:14
    一二八團開展“夏送清涼”慰問
    兵團工運(2018年8期)2018-01-22 07:50:52
    交錯網(wǎng)格與旋轉交錯網(wǎng)格對VTI介質波場分離的影響分析
    地震學報(2016年1期)2016-11-28 05:38:36
    基于Hilbert變換的全波場分離逆時偏移成像
    旋轉交錯網(wǎng)格VTI介質波場模擬與波場分解
    久久精品国产亚洲av涩爱| 久久99精品国语久久久| av在线老鸭窝| 亚洲不卡免费看| av在线老鸭窝| 亚洲无线观看免费| 91午夜精品亚洲一区二区三区| 九九久久精品国产亚洲av麻豆| 午夜亚洲福利在线播放| 亚洲精品中文字幕在线视频 | 亚洲色图av天堂| 男女边摸边吃奶| 国产国拍精品亚洲av在线观看| 欧美激情在线99| 超碰97精品在线观看| 午夜视频国产福利| freevideosex欧美| 欧美一区二区亚洲| 国产一级毛片在线| 欧美不卡视频在线免费观看| 欧美bdsm另类| 好男人视频免费观看在线| 精品久久久久久成人av| 亚洲在久久综合| 亚洲真实伦在线观看| 精品久久久久久久久av| 国产中年淑女户外野战色| 精品国内亚洲2022精品成人| 在线观看av片永久免费下载| 日韩av免费高清视频| 老司机影院成人| 日本一本二区三区精品| 国内少妇人妻偷人精品xxx网站| 亚州av有码| 一区二区三区乱码不卡18| 内地一区二区视频在线| 精品国产一区二区三区久久久樱花 | 久久久久久国产a免费观看| 免费观看的影片在线观看| 欧美 日韩 精品 国产| 特级一级黄色大片| 最新中文字幕久久久久| 一级毛片我不卡| av线在线观看网站| 欧美三级亚洲精品| 国产成人福利小说| 麻豆精品久久久久久蜜桃| 久久久精品免费免费高清| 国产女主播在线喷水免费视频网站 | 在线免费十八禁| 丝瓜视频免费看黄片| 精品酒店卫生间| 欧美日韩精品成人综合77777| 综合色丁香网| 内地一区二区视频在线| 亚洲欧美一区二区三区黑人 | 亚洲国产日韩欧美精品在线观看| 只有这里有精品99| 国产综合精华液| 麻豆av噜噜一区二区三区| kizo精华| 久久久午夜欧美精品| 亚洲精品国产成人久久av| 97精品久久久久久久久久精品| 成人午夜精彩视频在线观看| 日韩av不卡免费在线播放| 丰满乱子伦码专区| 又爽又黄无遮挡网站| 亚洲,欧美,日韩| 欧美日韩在线观看h| 一级毛片久久久久久久久女| 久久精品久久久久久久性| 日韩欧美国产在线观看| 精品国内亚洲2022精品成人| av线在线观看网站| a级毛色黄片| 国产免费视频播放在线视频 | 日韩国内少妇激情av| 久久久精品免费免费高清| 草草在线视频免费看| 亚洲伊人久久精品综合| 热99在线观看视频| 免费不卡的大黄色大毛片视频在线观看 | 99久国产av精品国产电影| 蜜桃久久精品国产亚洲av| 成人欧美大片| 汤姆久久久久久久影院中文字幕 | 一级爰片在线观看| 纵有疾风起免费观看全集完整版 | 真实男女啪啪啪动态图| 亚洲精品乱码久久久久久按摩| 大话2 男鬼变身卡| 97在线视频观看| 久久精品国产亚洲网站| 日本熟妇午夜| 久久精品国产自在天天线| 少妇猛男粗大的猛烈进出视频 | 一个人看的www免费观看视频| 亚洲伊人久久精品综合| 亚洲欧美一区二区三区国产| 成人高潮视频无遮挡免费网站| 美女主播在线视频| videossex国产| 久久久a久久爽久久v久久| 欧美xxⅹ黑人| 国产综合精华液| 狠狠精品人妻久久久久久综合| 亚洲色图av天堂| 色吧在线观看| 韩国av在线不卡| 精品久久久久久久人妻蜜臀av| av播播在线观看一区| 国产黄色小视频在线观看| 久久久久久久大尺度免费视频| 国产毛片a区久久久久| 菩萨蛮人人尽说江南好唐韦庄| 亚洲va在线va天堂va国产| 夜夜爽夜夜爽视频| 免费观看a级毛片全部| 卡戴珊不雅视频在线播放| 老师上课跳d突然被开到最大视频| 欧美xxⅹ黑人| 日本爱情动作片www.在线观看| 中文字幕免费在线视频6| 波野结衣二区三区在线| 天堂√8在线中文| 久久久久精品性色| 欧美激情国产日韩精品一区| 精品国产一区二区三区久久久樱花 | 日韩电影二区| 婷婷色麻豆天堂久久| 在线观看美女被高潮喷水网站| 成人亚洲精品av一区二区| 成人亚洲精品一区在线观看 | 久久久久久久久久黄片| 2021少妇久久久久久久久久久| 免费观看精品视频网站| 美女大奶头视频| 91久久精品国产一区二区成人| 精品亚洲乱码少妇综合久久| 午夜福利在线观看免费完整高清在| 精品久久国产蜜桃| 女的被弄到高潮叫床怎么办| 国产成人精品久久久久久| 色综合站精品国产| 国产久久久一区二区三区| 亚洲国产高清在线一区二区三| 97在线视频观看| 嘟嘟电影网在线观看| 欧美最新免费一区二区三区| 亚洲精品国产成人久久av| 欧美日韩亚洲高清精品| 国产精品一区二区性色av| 嫩草影院精品99| 亚洲av电影不卡..在线观看| 边亲边吃奶的免费视频| 69av精品久久久久久| videos熟女内射| 男的添女的下面高潮视频| 一级二级三级毛片免费看| 中文字幕免费在线视频6| 国产午夜福利久久久久久| 亚洲精品国产成人久久av| 国产精品蜜桃在线观看| 熟女电影av网| 欧美一区二区亚洲| 国产日韩欧美在线精品| 精品久久久久久久久亚洲| 一区二区三区乱码不卡18| a级一级毛片免费在线观看| 免费观看a级毛片全部| 亚洲综合精品二区| 99re6热这里在线精品视频| 菩萨蛮人人尽说江南好唐韦庄| 国产淫片久久久久久久久| 尾随美女入室| 亚洲精品色激情综合| 久久精品国产亚洲网站| 日本-黄色视频高清免费观看| 精品一区在线观看国产| 老师上课跳d突然被开到最大视频| 男女边摸边吃奶| 卡戴珊不雅视频在线播放| 免费观看精品视频网站| 日本wwww免费看| 麻豆成人午夜福利视频| 99久久九九国产精品国产免费| 国产男人的电影天堂91| 丝瓜视频免费看黄片| 国产永久视频网站| 国产精品久久视频播放| 欧美激情国产日韩精品一区| 97人妻精品一区二区三区麻豆| 久久久久久久午夜电影| 久久人人爽人人爽人人片va| 嫩草影院精品99| 青青草视频在线视频观看| 亚洲国产最新在线播放| 丝袜美腿在线中文| 国产中年淑女户外野战色| 最近最新中文字幕大全电影3| 秋霞伦理黄片| av福利片在线观看| 婷婷色麻豆天堂久久| 日韩大片免费观看网站| 少妇裸体淫交视频免费看高清| av女优亚洲男人天堂| 如何舔出高潮| 夜夜爽夜夜爽视频| 国产av码专区亚洲av| 亚洲精品日韩av片在线观看| 欧美日韩视频高清一区二区三区二| 丝瓜视频免费看黄片| 亚洲精品亚洲一区二区| videossex国产| 欧美高清成人免费视频www| 亚洲国产精品专区欧美| 啦啦啦啦在线视频资源| 午夜福利高清视频| 好男人视频免费观看在线| 亚洲人成网站在线观看播放| 在线观看人妻少妇| 精品久久久久久久久亚洲| 亚洲国产精品成人久久小说| 国产男女超爽视频在线观看| 搡老妇女老女人老熟妇| 久久久欧美国产精品| 99久久精品热视频| 91精品国产九色| av在线蜜桃| 亚洲欧美日韩卡通动漫| 国产伦在线观看视频一区| 免费不卡的大黄色大毛片视频在线观看 | 国产有黄有色有爽视频| 两个人视频免费观看高清| 在线观看人妻少妇| 亚洲av成人精品一二三区| 亚洲精品第二区| 国产精品1区2区在线观看.| 亚洲精品国产成人久久av| 中文在线观看免费www的网站| 人人妻人人澡欧美一区二区| 国产亚洲精品久久久com| 国产一区二区三区综合在线观看 | 一级毛片 在线播放| 少妇熟女欧美另类| 午夜福利高清视频| 欧美成人一区二区免费高清观看| 岛国毛片在线播放| 免费看美女性在线毛片视频| 日韩视频在线欧美| 免费电影在线观看免费观看| 成人二区视频| 2021天堂中文幕一二区在线观| 最后的刺客免费高清国语| 内射极品少妇av片p| 国产黄a三级三级三级人| 丝瓜视频免费看黄片| 只有这里有精品99| 非洲黑人性xxxx精品又粗又长| 夜夜看夜夜爽夜夜摸| 亚洲欧美日韩卡通动漫| 十八禁网站网址无遮挡 | 免费大片黄手机在线观看| 欧美另类一区| 国产亚洲av片在线观看秒播厂 | 人妻系列 视频| 国产免费福利视频在线观看| 午夜免费激情av| 国产成人freesex在线| 一边亲一边摸免费视频| videos熟女内射| 亚洲成人中文字幕在线播放| 国产亚洲一区二区精品| 精品一区二区三区人妻视频| 禁无遮挡网站| 国产有黄有色有爽视频| 丰满少妇做爰视频| 黄色日韩在线| 国产免费视频播放在线视频 | 能在线免费看毛片的网站| 不卡视频在线观看欧美| 中国国产av一级| 免费人成在线观看视频色| 亚洲国产精品专区欧美| 免费黄网站久久成人精品| 亚洲精品,欧美精品| 男人爽女人下面视频在线观看| 波野结衣二区三区在线| 日本爱情动作片www.在线观看| 国产午夜精品一二区理论片| 床上黄色一级片| 中文字幕人妻熟人妻熟丝袜美| 美女大奶头视频| 一级黄片播放器| 午夜福利在线观看免费完整高清在| 男女视频在线观看网站免费| 成人高潮视频无遮挡免费网站| 乱人视频在线观看| 国产亚洲精品av在线| 九色成人免费人妻av| 国产在视频线在精品| 欧美日韩精品成人综合77777| 麻豆精品久久久久久蜜桃| 中文天堂在线官网| 午夜爱爱视频在线播放| ponron亚洲| 亚洲综合精品二区| 全区人妻精品视频| 插阴视频在线观看视频| 精品不卡国产一区二区三区| 嫩草影院入口| 天天一区二区日本电影三级| 啦啦啦中文免费视频观看日本| 日韩在线高清观看一区二区三区| 国产国拍精品亚洲av在线观看| 最近最新中文字幕大全电影3| 亚洲av成人精品一二三区| 久热久热在线精品观看| 国产 一区精品| 亚洲精品第二区| ponron亚洲| 亚洲精品日本国产第一区| 国产熟女欧美一区二区| 精品久久久精品久久久| 久久久久久久久久人人人人人人| 欧美激情在线99| 日韩国内少妇激情av| 18禁在线播放成人免费| 中文资源天堂在线| 欧美性感艳星| 十八禁网站网址无遮挡 | 夫妻性生交免费视频一级片| 亚洲婷婷狠狠爱综合网| 精品不卡国产一区二区三区| 一本久久精品| 日本免费在线观看一区| 午夜精品在线福利| 日韩欧美 国产精品| 禁无遮挡网站| 国产免费一级a男人的天堂| 美女cb高潮喷水在线观看| 欧美日韩综合久久久久久| 色播亚洲综合网| 精品少妇黑人巨大在线播放| 一个人免费在线观看电影| 在线播放无遮挡| 亚洲av成人av| 国产激情偷乱视频一区二区| 久久久久久国产a免费观看| 久久精品熟女亚洲av麻豆精品 | 亚洲av在线观看美女高潮| 日韩精品有码人妻一区| 日本三级黄在线观看| 国产精品av视频在线免费观看| 欧美 日韩 精品 国产| 免费观看a级毛片全部| 女人被狂操c到高潮| 国产精品国产三级专区第一集| 51国产日韩欧美| av网站免费在线观看视频 | av福利片在线观看| 午夜免费激情av| 插逼视频在线观看| 午夜福利视频1000在线观看| 成人毛片60女人毛片免费| 国产片特级美女逼逼视频| 亚洲欧美一区二区三区黑人 | 美女大奶头视频| 久久精品久久久久久久性| 街头女战士在线观看网站| 久久久久性生活片| 久久国产乱子免费精品| 欧美高清性xxxxhd video| 搡女人真爽免费视频火全软件| 中文字幕免费在线视频6| 最新中文字幕久久久久| 伊人久久国产一区二区| 国产一区亚洲一区在线观看| 国产欧美日韩精品一区二区| 综合色av麻豆| 亚洲,欧美,日韩| 一夜夜www| 黄色配什么色好看| 97热精品久久久久久| 久久久久久久久久人人人人人人| 欧美日韩国产mv在线观看视频 | 白带黄色成豆腐渣| 国产亚洲精品av在线| av一本久久久久| 91精品一卡2卡3卡4卡| 欧美性感艳星| 国产毛片a区久久久久| 一个人观看的视频www高清免费观看| 网址你懂的国产日韩在线| 美女内射精品一级片tv| 一级毛片黄色毛片免费观看视频| 日本欧美国产在线视频| av网站免费在线观看视频 | 插逼视频在线观看| 国产老妇伦熟女老妇高清| 纵有疾风起免费观看全集完整版 | 国产男女超爽视频在线观看| 2021少妇久久久久久久久久久| 亚洲av日韩在线播放| av女优亚洲男人天堂| 男女啪啪激烈高潮av片| 欧美丝袜亚洲另类| 亚洲成人精品中文字幕电影| 晚上一个人看的免费电影| 日韩不卡一区二区三区视频在线| 国产亚洲av片在线观看秒播厂 | 日日啪夜夜爽| 亚洲自偷自拍三级| 亚洲人成网站在线观看播放| 国产黄色小视频在线观看| 黄色一级大片看看| 午夜免费激情av| 1000部很黄的大片| 青春草国产在线视频| 亚洲国产精品sss在线观看| 亚洲国产成人一精品久久久| 亚洲精品日韩av片在线观看| 内地一区二区视频在线| 国产成年人精品一区二区| 联通29元200g的流量卡| 综合色av麻豆| 日韩成人伦理影院| 免费观看无遮挡的男女| 国产高清不卡午夜福利| 老司机影院成人| 亚洲精品乱码久久久久久按摩| 成年免费大片在线观看| 久久精品夜夜夜夜夜久久蜜豆| 久久久久久久久久久丰满| 久久久久精品久久久久真实原创| 伊人久久国产一区二区| 九色成人免费人妻av| 亚洲欧美成人精品一区二区| 少妇裸体淫交视频免费看高清| 国产一区二区三区综合在线观看 | 一级二级三级毛片免费看| 国产一区二区三区综合在线观看 | av在线播放精品| 国产男女超爽视频在线观看| 国产精品久久久久久精品电影小说 | 亚洲av中文字字幕乱码综合| 天堂√8在线中文| 97超碰精品成人国产| 一级毛片久久久久久久久女| .国产精品久久| 免费观看的影片在线观看| 国产永久视频网站| 亚洲一区高清亚洲精品| 少妇丰满av| 国产 一区精品| 亚洲精品中文字幕在线视频 | 超碰av人人做人人爽久久| 大香蕉久久网| 亚洲最大成人手机在线| 亚洲真实伦在线观看| 亚洲精品国产av蜜桃| 美女国产视频在线观看| 80岁老熟妇乱子伦牲交| 国产av国产精品国产| 一个人看的www免费观看视频| 毛片女人毛片| av国产久精品久网站免费入址| 高清毛片免费看| 日韩欧美精品v在线| 最近2019中文字幕mv第一页| 精品少妇黑人巨大在线播放| 大香蕉久久网| 国产午夜精品一二区理论片| 人妻夜夜爽99麻豆av| 日韩av免费高清视频| 在线 av 中文字幕| 国产午夜精品一二区理论片| www.色视频.com| 亚洲av电影在线观看一区二区三区 | 99久久精品一区二区三区| 中文欧美无线码| 在现免费观看毛片| 国产伦一二天堂av在线观看| 最近最新中文字幕免费大全7| 少妇熟女欧美另类| 一个人看的www免费观看视频| 亚洲av免费在线观看| 亚洲婷婷狠狠爱综合网| 日日干狠狠操夜夜爽| 九九爱精品视频在线观看| av.在线天堂| 五月玫瑰六月丁香| 色视频www国产| 成人亚洲欧美一区二区av| 亚洲国产欧美在线一区| 成年人午夜在线观看视频 | 国产精品1区2区在线观看.| 美女xxoo啪啪120秒动态图| 亚洲国产精品国产精品| 亚洲真实伦在线观看| 亚洲精品日本国产第一区| 久久久久久久久久久免费av| 久久久久久国产a免费观看| 午夜久久久久精精品| 精华霜和精华液先用哪个| 少妇人妻精品综合一区二区| 插逼视频在线观看| 搡老妇女老女人老熟妇| 免费av观看视频| 少妇的逼好多水| 狠狠精品人妻久久久久久综合| 免费电影在线观看免费观看| 99久久精品一区二区三区| 校园人妻丝袜中文字幕| 高清午夜精品一区二区三区| 免费在线观看成人毛片| 亚洲美女视频黄频| 赤兔流量卡办理| 亚洲欧美日韩卡通动漫| 亚洲av日韩在线播放| 国产精品一二三区在线看| 日本wwww免费看| 成人午夜精彩视频在线观看| 国精品久久久久久国模美| 寂寞人妻少妇视频99o| 国产视频内射| 久久久久久久久大av| av专区在线播放| 国精品久久久久久国模美| 午夜福利网站1000一区二区三区| 男人狂女人下面高潮的视频| 久久久久久久久久黄片| 亚洲性久久影院| 美女cb高潮喷水在线观看| 亚洲精品日韩在线中文字幕| 别揉我奶头 嗯啊视频| 国产亚洲午夜精品一区二区久久 | 国产午夜精品久久久久久一区二区三区| 99久久人妻综合| 国产一区二区亚洲精品在线观看| 成人高潮视频无遮挡免费网站| 老司机影院成人| 久久国产乱子免费精品| 亚洲在线观看片| 校园人妻丝袜中文字幕| 国产免费一级a男人的天堂| 久热久热在线精品观看| 国产男人的电影天堂91| 十八禁网站网址无遮挡 | 亚洲欧美精品专区久久| 日韩欧美国产在线观看| 99九九线精品视频在线观看视频| 精品不卡国产一区二区三区| 色视频www国产| 国产一区亚洲一区在线观看| 搞女人的毛片| 亚洲精品日本国产第一区| 黄片wwwwww| av专区在线播放| 美女内射精品一级片tv| freevideosex欧美| 亚洲天堂国产精品一区在线| 欧美性感艳星| 免费观看av网站的网址| 极品教师在线视频| 亚洲内射少妇av| 人人妻人人澡人人爽人人夜夜 | 国产一级毛片七仙女欲春2| 女人被狂操c到高潮| 一级片'在线观看视频| 国产一区二区三区综合在线观看 | 插阴视频在线观看视频| av播播在线观看一区| av在线老鸭窝| 久久97久久精品| 日韩av不卡免费在线播放| 国产在视频线在精品| 欧美一区二区亚洲| 天堂中文最新版在线下载 | 亚洲在线观看片| 一二三四中文在线观看免费高清| 在线观看av片永久免费下载| 丝瓜视频免费看黄片| 淫秽高清视频在线观看| 亚洲国产高清在线一区二区三| 免费看日本二区| 日本猛色少妇xxxxx猛交久久| 一个人看的www免费观看视频| 三级国产精品片| 亚洲av福利一区| 色综合色国产| 午夜激情福利司机影院| 亚洲av男天堂| 99re6热这里在线精品视频| 亚洲国产色片| 国产成人91sexporn| 国产伦在线观看视频一区| 老司机影院成人| 婷婷六月久久综合丁香| 亚洲精品一二三| 久久久色成人| 最新中文字幕久久久久| 免费av观看视频| 老司机影院成人| 日本wwww免费看| 欧美激情国产日韩精品一区| 久99久视频精品免费| 菩萨蛮人人尽说江南好唐韦庄| 伦精品一区二区三区| 熟女人妻精品中文字幕| 免费av毛片视频| 看免费成人av毛片| 99热全是精品| 亚洲av免费在线观看|