• <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麻豆| 男女国产视频网站| 亚洲国产毛片av蜜桃av| av播播在线观看一区| 国产高清三级在线| 美女内射精品一级片tv| 中文字幕人妻丝袜制服| 国产精品久久久久久精品古装| 熟女av电影| 一二三四中文在线观看免费高清| 精品国产露脸久久av麻豆| 国产日韩欧美在线精品| 中文字幕精品免费在线观看视频 | 成年人免费黄色播放视频 | 亚洲精品视频女| 嫩草影院入口| 国产精品秋霞免费鲁丝片| 国产伦精品一区二区三区四那| 午夜91福利影院| 欧美性感艳星| 亚洲综合色惰| 中文欧美无线码| 最近中文字幕2019免费版| 久久亚洲国产成人精品v| 乱码一卡2卡4卡精品| 国产亚洲一区二区精品| av福利片在线观看| 六月丁香七月| 色婷婷av一区二区三区视频| 欧美精品一区二区大全| 亚洲av成人精品一区久久| 毛片一级片免费看久久久久| xxx大片免费视频| 精品久久久久久久久亚洲| 精品国产一区二区三区久久久樱花| 日韩精品免费视频一区二区三区 | 国产av国产精品国产| 亚洲国产精品专区欧美| 久久久亚洲精品成人影院| 成年人免费黄色播放视频 | 免费黄频网站在线观看国产| 亚洲精品,欧美精品| 亚洲国产欧美日韩在线播放 | 偷拍熟女少妇极品色| 精品亚洲成国产av| 亚洲成人手机| 国产熟女午夜一区二区三区 | 国产91av在线免费观看| 亚洲av免费高清在线观看| 99久久综合免费| 视频中文字幕在线观看| 久久久久人妻精品一区果冻| 久久午夜综合久久蜜桃| 丰满乱子伦码专区| 少妇的逼好多水| 我的女老师完整版在线观看| 视频区图区小说| 丝袜喷水一区| 日韩免费高清中文字幕av| 久久久午夜欧美精品| 国产色爽女视频免费观看| 日韩成人av中文字幕在线观看| 人人妻人人澡人人看| 亚洲电影在线观看av| 亚洲精华国产精华液的使用体验| 纯流量卡能插随身wifi吗| 欧美 亚洲 国产 日韩一| 国产精品三级大全| 美女国产视频在线观看| 男女啪啪激烈高潮av片| 人人妻人人澡人人爽人人夜夜| 国产高清三级在线| 在线免费观看不下载黄p国产| 国产男女内射视频| 婷婷色麻豆天堂久久| 久久鲁丝午夜福利片| 在线播放无遮挡| 免费看av在线观看网站| 青青草视频在线视频观看| 国产高清三级在线| 97精品久久久久久久久久精品| 青春草国产在线视频| 国产精品蜜桃在线观看| 曰老女人黄片| 国产美女午夜福利| 毛片一级片免费看久久久久| 搡女人真爽免费视频火全软件| 国产黄片美女视频| 免费观看av网站的网址| 国产免费一级a男人的天堂| 亚洲精品乱码久久久久久按摩| 日日摸夜夜添夜夜爱| 精品一区在线观看国产| 成人特级av手机在线观看| 91精品国产国语对白视频| 亚洲婷婷狠狠爱综合网| 乱系列少妇在线播放| 免费观看的影片在线观看| 久久国内精品自在自线图片| 亚洲四区av| 插阴视频在线观看视频| 最近中文字幕高清免费大全6| 桃花免费在线播放| 亚洲欧美精品专区久久| 老司机影院毛片| 91精品一卡2卡3卡4卡| 久久久a久久爽久久v久久| 三级国产精品片| 五月天丁香电影| 在线观看免费高清a一片| 亚洲欧洲日产国产| 99热网站在线观看| 国产精品成人在线| 我的老师免费观看完整版| 国产成人freesex在线| 国产又色又爽无遮挡免| 精品久久久噜噜| 极品少妇高潮喷水抽搐| 亚洲人成网站在线播| 精品亚洲成国产av| 亚洲图色成人| 国产日韩欧美亚洲二区| 又黄又爽又刺激的免费视频.| 97在线人人人人妻| 国产成人精品一,二区| 国产精品一二三区在线看| 2021少妇久久久久久久久久久| 久久久精品94久久精品| 一级毛片 在线播放| 亚洲欧美成人综合另类久久久| 日韩av不卡免费在线播放| 永久网站在线| 久久av网站| 亚洲人成网站在线观看播放| 国产精品99久久99久久久不卡 | 最近手机中文字幕大全| 一级片'在线观看视频| 亚洲美女黄色视频免费看| 欧美bdsm另类| 久久99一区二区三区| 99国产精品免费福利视频| 国产成人午夜福利电影在线观看| 日韩在线高清观看一区二区三区| 亚洲av二区三区四区| kizo精华| 久久精品国产亚洲av涩爱| 久久久久国产网址| 国产高清不卡午夜福利| 男女免费视频国产| 日韩欧美 国产精品| 2022亚洲国产成人精品| 国产 一区精品| 97超视频在线观看视频| 多毛熟女@视频| 亚洲欧洲精品一区二区精品久久久 | 精品少妇黑人巨大在线播放| 80岁老熟妇乱子伦牲交| 国产91av在线免费观看| 亚洲精品成人av观看孕妇| 欧美 亚洲 国产 日韩一| 美女内射精品一级片tv| 免费大片18禁| 91久久精品国产一区二区三区| 日韩中文字幕视频在线看片| 九九久久精品国产亚洲av麻豆| 少妇人妻一区二区三区视频| 午夜日本视频在线| 两个人的视频大全免费| 中文字幕人妻丝袜制服| 国产男女超爽视频在线观看| 99re6热这里在线精品视频| 人人妻人人看人人澡| 亚洲丝袜综合中文字幕| 一个人免费看片子| 午夜视频国产福利| 两个人免费观看高清视频 | 极品少妇高潮喷水抽搐| 日韩亚洲欧美综合| 色视频在线一区二区三区| 精品国产国语对白av| 久久久久人妻精品一区果冻| √禁漫天堂资源中文www| 国产成人精品婷婷| 欧美精品国产亚洲| 国产一区亚洲一区在线观看| 中国三级夫妇交换| 亚洲国产av新网站| 日韩中文字幕视频在线看片| 国产精品人妻久久久影院| 久久人人爽人人爽人人片va| 激情五月婷婷亚洲| 嫩草影院新地址| 天堂8中文在线网| 欧美日韩视频高清一区二区三区二| 日韩亚洲欧美综合| 国产美女午夜福利| 麻豆成人午夜福利视频| 免费看av在线观看网站| 免费少妇av软件| 亚洲精品一二三| 久久99一区二区三区| 久久久久久久久久久久大奶| 亚洲精品亚洲一区二区| 九九久久精品国产亚洲av麻豆| 久久99精品国语久久久| 七月丁香在线播放| 久久ye,这里只有精品| 在线观看免费高清a一片| 国内精品宾馆在线| 三级国产精品片| 丁香六月天网| 女人久久www免费人成看片| 亚洲精品视频女| 国产永久视频网站| 最新的欧美精品一区二区| 日韩av不卡免费在线播放| 欧美三级亚洲精品| 99国产精品免费福利视频| 亚洲成色77777| 久久久亚洲精品成人影院| 欧美精品人与动牲交sv欧美| 91午夜精品亚洲一区二区三区| 人体艺术视频欧美日本| 久久女婷五月综合色啪小说| 永久免费av网站大全| 日产精品乱码卡一卡2卡三| 女人精品久久久久毛片| 午夜免费鲁丝| 国产精品熟女久久久久浪| 国产无遮挡羞羞视频在线观看| 极品少妇高潮喷水抽搐| 亚洲欧美成人精品一区二区| 高清午夜精品一区二区三区| 亚洲欧美成人综合另类久久久| 国产亚洲最大av| 成人无遮挡网站| 美女视频免费永久观看网站| av专区在线播放| 成人免费观看视频高清| 久久久亚洲精品成人影院| 人体艺术视频欧美日本| 伊人久久精品亚洲午夜| 日韩成人伦理影院| 岛国毛片在线播放| a级一级毛片免费在线观看| 女人久久www免费人成看片| 国产伦理片在线播放av一区| 欧美最新免费一区二区三区| 亚洲国产精品999| 午夜福利,免费看| 久久久久久久久久成人| 国产欧美亚洲国产| 边亲边吃奶的免费视频| 中文欧美无线码| 黑人高潮一二区| 少妇裸体淫交视频免费看高清| 亚洲综合精品二区| 国产无遮挡羞羞视频在线观看| 在线观看一区二区三区激情| 波野结衣二区三区在线| 亚洲欧美精品专区久久| av一本久久久久| 高清不卡的av网站| 国产成人aa在线观看| 亚洲国产成人一精品久久久| 久久久欧美国产精品| 99re6热这里在线精品视频| 国产男女超爽视频在线观看| 成人特级av手机在线观看| 人妻人人澡人人爽人人| 伊人久久国产一区二区| 亚洲成人av在线免费| 最新中文字幕久久久久| 啦啦啦中文免费视频观看日本| 美女视频免费永久观看网站| 国产伦理片在线播放av一区| 久久99精品国语久久久| 亚洲精品国产成人久久av| 在线免费观看不下载黄p国产| 国产精品国产av在线观看| 婷婷色麻豆天堂久久| 青青草视频在线视频观看| 久热久热在线精品观看| 国产在线视频一区二区| av福利片在线| 亚洲婷婷狠狠爱综合网| 十八禁高潮呻吟视频 | 国产永久视频网站| 亚洲国产最新在线播放| 亚洲性久久影院| 亚洲自偷自拍三级| 亚洲图色成人| 亚洲精品,欧美精品| 亚洲av免费高清在线观看| 99re6热这里在线精品视频| 爱豆传媒免费全集在线观看| av线在线观看网站| av福利片在线| 99热全是精品| 肉色欧美久久久久久久蜜桃| 亚洲成人手机| 亚洲在久久综合| 色哟哟·www| 精品亚洲乱码少妇综合久久| 最黄视频免费看| 日韩 亚洲 欧美在线| 九九在线视频观看精品| 91久久精品国产一区二区成人| 久热这里只有精品99| 七月丁香在线播放| av线在线观看网站| 欧美日韩视频精品一区| 亚洲精品成人av观看孕妇| 天堂中文最新版在线下载| av在线播放精品| 99久久精品一区二区三区| 亚洲伊人久久精品综合| 国产欧美日韩一区二区三区在线 | 亚洲精品色激情综合| 在线观看www视频免费| 天堂中文最新版在线下载| 男男h啪啪无遮挡| 人妻 亚洲 视频| 各种免费的搞黄视频| 蜜桃久久精品国产亚洲av| 久久久亚洲精品成人影院| 最近手机中文字幕大全| 亚洲成人av在线免费| 亚洲av福利一区| 建设人人有责人人尽责人人享有的| 中文字幕人妻丝袜制服| 插阴视频在线观看视频| 国产精品.久久久| 青春草亚洲视频在线观看| 精品久久久噜噜| 少妇人妻久久综合中文| 国产在线一区二区三区精| 欧美日韩视频精品一区| 超碰97精品在线观看| 日韩av在线免费看完整版不卡| 亚洲av免费高清在线观看| 国产精品一区二区在线观看99| 日本av手机在线免费观看| 人妻 亚洲 视频| 大码成人一级视频| 日本欧美视频一区| 亚洲自偷自拍三级| 国产日韩欧美在线精品| 五月玫瑰六月丁香| 日韩视频在线欧美| 一边亲一边摸免费视频| 免费黄频网站在线观看国产| 色吧在线观看| 国产在线一区二区三区精| 少妇被粗大猛烈的视频| 久热久热在线精品观看| 街头女战士在线观看网站| 国产中年淑女户外野战色| 夜夜看夜夜爽夜夜摸| 国产免费一区二区三区四区乱码| 免费av不卡在线播放| 黄色毛片三级朝国网站 | 99九九在线精品视频 | 大香蕉97超碰在线| 最近中文字幕2019免费版| 亚洲色图综合在线观看| 九草在线视频观看| 久久国产亚洲av麻豆专区| 麻豆成人av视频| 欧美国产精品一级二级三级 | 久久综合国产亚洲精品| 亚洲久久久国产精品| 丰满迷人的少妇在线观看| 午夜免费男女啪啪视频观看| 亚洲精品国产成人久久av| 人妻 亚洲 视频| 精品久久久久久久久av| 久久这里有精品视频免费| 啦啦啦在线观看免费高清www| 亚洲人成网站在线播| 美女主播在线视频| 国产爽快片一区二区三区| 3wmmmm亚洲av在线观看| 国产精品蜜桃在线观看| 国产极品粉嫩免费观看在线 | av在线app专区| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品国产av蜜桃| 午夜av观看不卡| 三级国产精品片| 这个男人来自地球电影免费观看 | 建设人人有责人人尽责人人享有的| 欧美精品一区二区大全| 一本大道久久a久久精品| 九九在线视频观看精品| 大香蕉97超碰在线| 成年美女黄网站色视频大全免费 | 亚洲精品一区蜜桃| 99久久中文字幕三级久久日本| 97在线人人人人妻| 久久久久视频综合| 欧美97在线视频| 国产亚洲欧美精品永久| 两个人免费观看高清视频 | 国产一区二区三区综合在线观看 | 精品久久久久久久久av| 成人午夜精彩视频在线观看| 午夜久久久在线观看| 国产淫语在线视频| 黑人猛操日本美女一级片| 国产在线男女| 中国国产av一级| 国产毛片在线视频| 黄色怎么调成土黄色| 亚洲在久久综合| 天堂8中文在线网| 国产探花极品一区二区| 最近中文字幕高清免费大全6| 色5月婷婷丁香| 亚洲av福利一区| 精品亚洲乱码少妇综合久久| 草草在线视频免费看| 日日啪夜夜爽| 一本色道久久久久久精品综合| 97在线人人人人妻| 国产极品粉嫩免费观看在线 | 自线自在国产av| av播播在线观看一区| 国产美女午夜福利| 99久国产av精品国产电影| 你懂的网址亚洲精品在线观看| 春色校园在线视频观看| 你懂的网址亚洲精品在线观看| 99久久精品国产国产毛片| 三级国产精品欧美在线观看| 日韩 亚洲 欧美在线| 国产精品国产三级专区第一集| 久久久久网色| freevideosex欧美| 超碰97精品在线观看| 五月伊人婷婷丁香| 97在线人人人人妻| 一级片'在线观看视频| a级一级毛片免费在线观看| 国产亚洲最大av| 国产成人一区二区在线| 在线观看免费高清a一片| 久久国产精品大桥未久av | 亚洲av中文av极速乱| 在线免费观看不下载黄p国产| 在线观看免费日韩欧美大片 | 国产毛片在线视频| 少妇猛男粗大的猛烈进出视频| 国产白丝娇喘喷水9色精品| 99九九在线精品视频 | 国产精品女同一区二区软件| 国产免费视频播放在线视频| 大陆偷拍与自拍| av不卡在线播放| 22中文网久久字幕| 国产精品国产三级国产专区5o| 国产亚洲5aaaaa淫片| 91精品国产九色| 精品国产一区二区久久| 精品久久久噜噜| 成人特级av手机在线观看| 亚洲,欧美,日韩| 午夜福利视频精品| 寂寞人妻少妇视频99o| 狂野欧美白嫩少妇大欣赏| 中文欧美无线码| 熟妇人妻不卡中文字幕| 国产精品人妻久久久影院| 日韩一区二区视频免费看| 国产高清国产精品国产三级| 欧美人与善性xxx| 久久韩国三级中文字幕| 久久人妻熟女aⅴ| 最近2019中文字幕mv第一页| 国产 一区精品| 久久久久久久久久久丰满| 亚洲av福利一区| kizo精华| 国产精品久久久久久av不卡| 永久网站在线| 一级黄片播放器| av福利片在线| 国产成人精品一,二区| 欧美国产精品一级二级三级 | 丰满少妇做爰视频| 久久久久久久久久久久大奶| 免费看av在线观看网站| 国产成人91sexporn| 久久精品熟女亚洲av麻豆精品| 亚洲av电影在线观看一区二区三区| 精品视频人人做人人爽| 久久人人爽av亚洲精品天堂| 国产一区二区三区综合在线观看 | 国产有黄有色有爽视频| 国产精品嫩草影院av在线观看| 午夜福利网站1000一区二区三区| 秋霞在线观看毛片| 亚洲欧美成人综合另类久久久| 在线观看一区二区三区激情| 丰满乱子伦码专区| 免费大片黄手机在线观看| av卡一久久| 午夜福利在线观看免费完整高清在| 亚州av有码| 男人添女人高潮全过程视频| 欧美日韩视频高清一区二区三区二| 波野结衣二区三区在线| 精品人妻熟女av久视频| 五月伊人婷婷丁香| 免费少妇av软件| 99九九线精品视频在线观看视频| 成人亚洲精品一区在线观看| 男人爽女人下面视频在线观看| 国产精品福利在线免费观看| 又爽又黄a免费视频| 一边亲一边摸免费视频| 国产在线视频一区二区| 插阴视频在线观看视频| 日日啪夜夜撸| 亚洲国产精品成人久久小说| 少妇人妻一区二区三区视频| 99热6这里只有精品| 亚洲精品久久久久久婷婷小说| 在线亚洲精品国产二区图片欧美 | 午夜91福利影院| 伊人亚洲综合成人网| 国产精品久久久久久精品古装| av在线观看视频网站免费| 中文字幕制服av| 色吧在线观看| 亚洲无线观看免费| 国产高清有码在线观看视频| 黑丝袜美女国产一区| 久久久久国产精品人妻一区二区| 久久精品国产自在天天线| 午夜影院在线不卡| 观看美女的网站| 肉色欧美久久久久久久蜜桃| 色视频www国产| 我的女老师完整版在线观看| 18+在线观看网站| 热re99久久国产66热| 午夜福利网站1000一区二区三区| 嘟嘟电影网在线观看| 国产精品人妻久久久影院| 欧美人与善性xxx| 美女cb高潮喷水在线观看| 亚洲精品久久久久久婷婷小说| 乱人伦中国视频| 男女国产视频网站| 国产精品一区二区在线观看99| 最近2019中文字幕mv第一页| 一区二区av电影网| 另类亚洲欧美激情| videossex国产| 亚洲精品久久久久久婷婷小说| 久久精品国产鲁丝片午夜精品| 精品一区在线观看国产| 看十八女毛片水多多多| 久久精品夜色国产| 精品一区二区免费观看| 少妇被粗大的猛进出69影院 | 哪个播放器可以免费观看大片| 日本欧美视频一区| 啦啦啦中文免费视频观看日本| 精品国产露脸久久av麻豆| 欧美xxxx性猛交bbbb| 九九久久精品国产亚洲av麻豆| 国产午夜精品一二区理论片| av卡一久久| 婷婷色麻豆天堂久久| 中文字幕人妻熟人妻熟丝袜美| 伦理电影免费视频| 九九在线视频观看精品| 久久精品久久久久久久性| 爱豆传媒免费全集在线观看| 99re6热这里在线精品视频| 色婷婷久久久亚洲欧美| 日韩伦理黄色片| 99久久精品一区二区三区| 亚洲成人手机| 国产精品国产三级专区第一集| 91成人精品电影| 午夜免费鲁丝| 交换朋友夫妻互换小说| 丝瓜视频免费看黄片| 国产精品一区二区在线不卡| 欧美日韩综合久久久久久| 天堂俺去俺来也www色官网| 国产欧美另类精品又又久久亚洲欧美| 免费看光身美女| 国产 精品1| 色婷婷av一区二区三区视频| 日韩成人伦理影院| 国产精品熟女久久久久浪| 最黄视频免费看| av视频免费观看在线观看| videossex国产| av播播在线观看一区| 亚洲人与动物交配视频| 国产精品秋霞免费鲁丝片| 中国美白少妇内射xxxbb| 中文字幕精品免费在线观看视频 | 永久免费av网站大全| 高清毛片免费看| 9色porny在线观看| 午夜福利在线观看免费完整高清在| 中文字幕制服av| 日韩强制内射视频|