基于Ruger 近似式預(yù)測裂縫型儲層的可行性研究

吳姍姍，高 剛，桂志先，王 鵬，李成毅

(1.長江大學(xué)油氣資源與勘探技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室，湖北 武漢 430100；2.長江大學(xué)地球物理與石油資源學(xué)院，湖北 武漢 430100；3.中國石油集團(tuán)川慶鉆探工程有限公司地球物理勘探公司，四川 成都 610213)

裂縫型儲層的預(yù)測方法根據(jù)資料來源通常分為兩類：一類是通過測井資料與巖心露頭資料定量識別儲層的裂縫，該方法預(yù)測結(jié)果精度較高，但其結(jié)果難以外推至井間[1]；另一類是通過地震資料定性預(yù)測儲層的裂縫，由于地震資料有橫向分辨率高、勘探成本較低、數(shù)據(jù)分布范圍廣等特點(diǎn)，因而運(yùn)用地震資料檢測裂縫的方法廣泛應(yīng)用于裂縫型儲層的預(yù)測[2]。目前運(yùn)用地震資料檢測裂縫的方法主要有：①疊后地震屬性分析法，即利用疊后資料分析其地震屬性，如相干屬性[3]、曲率屬性[4]、螞蟻體屬性[5]等；② 轉(zhuǎn)換波裂縫檢測法，利用橫波分裂產(chǎn)生的快慢橫波時(shí)差來反映裂縫發(fā)育的主方向和裂縫密度，常用的方法有相對時(shí)差梯度法[6]和層剝離法[7]；③疊前縱波方位各向異性檢測法，利用縱波各向異性進(jìn)行裂縫檢測的方法有動校正(Normal Moveout Correction，簡稱NMO)速度方位變化裂縫檢測法、正交地震測線縱波時(shí)差裂縫檢測法、縱波方位AVO(Amplitude Variation with Offset and Azimuth，簡稱AVOZ/AVAZ)和縱波阻抗隨方位角變化(Impedance Variation with Azimuth，簡稱IPVA)裂縫檢測法[8]。

利用疊前縱波資料的方位各向異性特征來識別裂縫型儲層是目前應(yīng)用最為廣泛的方法之一。該技術(shù)的研究始于20 世紀(jì)70 年代，由于受到當(dāng)時(shí)勘探條件的限制，直到20 世紀(jì)90 年代，I.Tsvankin[9]和A.Ruger[10]分別推導(dǎo)了HTI(Horizontal Transverse Isotropy，簡稱HTI)介質(zhì)中純模式波的動校正速度隨方位角變化(Velocity Variation with Azimuth，簡稱VVA)和縱波的反射系數(shù)隨入射角和方位角變化(AVOZ/AVAZ)的近似公式之后，該技術(shù)才引起了廣泛的關(guān)注。隨后S.Malick 等[11]、范國章等[12]、朱成宏等[13]、朱兆林等[14]對地震波方位AVO 特性進(jìn)行了全面的分析，建立了地震波方位AVO 的變化特征與裂縫走向、傾向以及發(fā)育強(qiáng)度的關(guān)系，該方法在實(shí)際裂縫檢測應(yīng)用中也取得了很好的效果[15-19]。該方法在具體應(yīng)用時(shí)僅考慮了各向異性參數(shù)對地震響應(yīng)特征的影響，并未考慮巖性參數(shù)對地震響應(yīng)特征的影響。然而，在實(shí)際資料的應(yīng)用過程中發(fā)現(xiàn)，在儲層與蓋層的彈性參數(shù)存在較為嚴(yán)重疊置現(xiàn)象的情況下，常規(guī)的彈性參數(shù)預(yù)測儲層難度較大，雖然裂縫參數(shù)預(yù)測縱波方位各向異性的相關(guān)技術(shù)方法可以提高儲層的預(yù)測精度，但是存在一定的多解性。因此，筆者以Ruger 近似式為基礎(chǔ)建立模型，將模型中各向異性參數(shù)和巖性參數(shù)設(shè)置為概率密度分布函數(shù)，通過Monte Carlo 隨機(jī)方法進(jìn)行疊前AVAZ正演模擬擬合橢圓，進(jìn)而分析各參數(shù)的隨機(jī)分布對擬合的橢圓信息即裂縫信息的影響。

1 疊前AVAZ 技術(shù)檢測裂縫原理

1.1 常用裂縫模型

在眾多的裂縫巖石物理模型中，目前最為常用的模型是Hudson 的薄硬幣模型和Schoenberg 的線性滑動模型。

Hudson 模型主要是描述彈性固體中包含薄硬幣形狀的橢球裂縫的等效介質(zhì)模型，該模型主要通過裂縫密度e和裂縫的縱橫比χ來表征裂縫系統(tǒng)，裂縫的密度是采用裂縫體積密度的定義方式：

式中a是裂縫的半徑；N/V是單位體積內(nèi)裂縫的數(shù)量；φ是裂縫誘導(dǎo)的孔隙率；χ是裂縫的縱橫比。

Schoenberg 線性滑動模型是基于Backus 細(xì)層層狀介質(zhì)模型提出的。該模型在巖石物理模型的研究中，適用于具有線性連續(xù)邊界的、充滿弱強(qiáng)度(充填物彈性模量小)介質(zhì)的平行層模型。C.J.Hsu 等[20]提出了兩個(gè)無量綱的參數(shù)，用來描述由裂縫引起的巖石剛度變化特征。

式中ΔN和ΔT分別代表由裂縫引起的垂直于裂縫面方向和平行于裂縫面方向的彈性參數(shù)差值，它們的變化范圍為0～1；EN和ET為非負(fù)的、無量綱的與裂縫有關(guān)的柔度參數(shù)；ZN為垂向柔度；ZT為平行于裂縫面的柔度；μ、λ為拉梅參數(shù)。

兩種模型參數(shù)之間具有一定的關(guān)聯(lián)性，ΔN和ΔT與裂縫參數(shù)之間的關(guān)系見式(4)和式(5)。

式中g(shù)=μ/(λ+2μ)；K′和μ′是縫隙中充填物的體積模量和剪切模量。

A.Bakulin 等[21]給出了HTI 介質(zhì)中各向異性參數(shù)(ε(V)、γ、δ(V))與Schoenberg 裂縫模型巖石物理參數(shù)(ΔN和ΔT)之間的關(guān)系，如式(6)—式(8)所示：

將式(4)—式(5)代入式(6)—式(8)可以得到HTI介質(zhì)中各向異性參數(shù)(ε(V)、γ、δ(V))與Hudson 模型所定義的裂縫密度之間的關(guān)系，如式(9)—式(11)所示，揭示了HTI 介質(zhì)中各向異性參數(shù)(ε(V)、γ、δ(V))與裂縫的密度之間的關(guān)系，可以分析L.Thomsen[22]各向異性參數(shù)與裂縫密度及裂縫中流體之間的關(guān)系。

1.2 HTI 介質(zhì)中疊前AVAZ 技術(shù)檢測裂縫原理

A.Ruger 等[10]對HTI 介質(zhì)中地震波傳播理論(圖1)進(jìn)行了深入研究后，在弱各向異性的假設(shè)條件下，

給出了HTI 介質(zhì)中縱波的反射系數(shù)隨入射角和方位角變化的近似公式，即：

式中Z為縱波波阻抗；β為縱波波速；α為橫波波速；(ε(V)、δ(V)為不同于L.Thomsen 裂隙參數(shù)，后續(xù)為了方便，參數(shù)直接寫為ε、δ)，ε、γ分別表示縱波和橫波各向異性；δ表示縱向和橫向各向異性；G=ρβ2表示切向模量；i和φ分別為入射方向和測線方向與裂縫對稱軸的夾角；符號“-”和“Δ”分別表示上下界面參數(shù)的算術(shù)平均值、差值。

HTI 介質(zhì)的各向異性參數(shù)的表達(dá)式：

式中cij為HTI 介質(zhì)彈性常數(shù)。

圖1 HTI 介質(zhì)中縱波傳播示意圖Fig.1 Sketch map of P wave propagation in HIT medium

入射角較小時(shí)，各向異性介質(zhì)中縱波反射振幅與裂縫參數(shù)之間的關(guān)系可以近似表示為：

式中P為縱波垂直入射時(shí)的反射振幅；Giso為反射振幅隨偏移距的變化率，即各向同性梯度；Gani為振幅隨偏移距和方位角的變化率，即各向異性梯度，指示了介質(zhì)受各向異性影響的部分。

在固定入射角的情況下，式(16)可以進(jìn)一步簡化為：

式中A為與炮檢距有關(guān)的偏置因子；B為與炮檢距和裂縫特征有關(guān)的調(diào)制因子；θ是炮檢距方向與裂縫走向的夾角。

式(17)反射系數(shù)R隨著角度θ的變化可用橢圓狀的圖形來表示(圖2)。

參數(shù)A與炮檢距有關(guān)，反映了巖性變化所引起的振幅的變化量；參數(shù)B相當(dāng)于振幅隨炮檢距方位改變的變化量，反映了地層裂縫對反射振幅強(qiáng)度的影響，其大小決定了儲層裂縫的發(fā)育程度，可視為各向異性因子。因此，可以把B/A即橢圓的扁率作為裂縫密度的相對量度，B/A的值越大，說明裂縫發(fā)育越好，反之，裂縫發(fā)育越差，而對于各向同性介質(zhì)，A值就是均勻介質(zhì)下的反射振幅，B值為零。

圖2 地震反射振幅隨方位角變化Fig.2 Variation of seismic reflection amplitude with azimuth

2 Monte Carlo 隨機(jī)正演模擬

疊前AVAZ 技術(shù)檢測裂縫法主要是通過地震資料進(jìn)行橢圓擬合，然后根據(jù)橢圓的信息對裂縫進(jìn)行預(yù)測和評價(jià)。在實(shí)際應(yīng)用中，當(dāng)橫向巖性變化時(shí)，通常是利用平均振幅A的信息，沒有考慮橫向巖性變化對擬合橢圓信息的影響，這在一定程度上造成了多解性。為了了解橫向巖性參數(shù)的變化對橢圓擬合信息的影響，對儲層物性較好，以孔-縫型、裂縫型儲層為主的研究區(qū)域的測井資料進(jìn)行統(tǒng)計(jì)，灰?guī)r儲層的縱波速度為3 000～5 200 m/s，密度為2.40～2.75 g/cm3，蓋層泥巖縱波速度為2 600～4 500 m/s，密度為2.1～2.5 g/cm3，據(jù)此建立含EDA 介質(zhì)的雙層模型(圖3)。模型參數(shù)設(shè)置見表1，入射角為30°，采用Monte Carlo 隨機(jī)方法產(chǎn)生滿足條件的隨機(jī)模型，利用最小二乘法進(jìn)行橢圓擬合，然后分析Ruger公式(式(12))中各個(gè)參數(shù)變化對擬合橢圓信息的影響，進(jìn)而說明該參數(shù)對裂縫評價(jià)的影響。

圖3 含EDA 介質(zhì)的雙層介質(zhì)模型Fig.3 Two-layer model with EDA medium

2.1 各向異性參數(shù)對疊前AVAZ 地震響應(yīng)影響

首先，通過模型分析各向異性參數(shù)單獨(dú)變化對疊前AVAZ 地震響應(yīng)特征的影響。L.Thomsen[23]通過實(shí)驗(yàn)測定了一系列不同巖性的樣品，根據(jù)其測定的不同巖性的各向異性參數(shù)值進(jìn)行統(tǒng)計(jì)，發(fā)現(xiàn)各向異性參數(shù)值基本都小于0.2，所以本文模型各向異性參數(shù)分別設(shè)置為均值0.05、0.10、0.20 的正態(tài)分布，其中標(biāo)準(zhǔn)差設(shè)置為均值的10%，為了保證擬合出的橢圓扁率變化是由各向異性參數(shù)變化所引起，將模型速度和密度設(shè)置為固定的平均值(表1)，利用Ruger 近似式進(jìn)行隨機(jī)正演模擬，分別作出不同各向異性參數(shù)值所對應(yīng)的橢圓扁率B/A與各向異性因子B散點(diǎn)圖(圖4)。

表1 雙層介質(zhì)模型參數(shù)Table 1 Parameters of a two-layer medium model

圖4 不同各向異性參數(shù)的B/A(橢圓扁率)與B(各向異性因子)散點(diǎn)圖Fig.4 The scatter plots of B/A and B with different anisotropic parameters

圖4a—圖4c 分別對應(yīng)各向異性參數(shù)ε(圖4a)、γ(圖4b)、δ(圖4c)為不同均值的正態(tài)分布時(shí)的橢圓的扁率B/A與各向異性因子B散點(diǎn)圖，其中紅色、紫色及藍(lán)色分別為模型的各向異性參數(shù)值設(shè)置為均值0.05、0.10 和0.20 的正態(tài)分布結(jié)果。由圖4 可以看出橢圓的扁率B/A與各向異性因子B的值隨各向異性參數(shù)值增大而增大，說明各向異性越強(qiáng)，其各向異性因子B與橢圓的扁率B/A值越大，證明了利用縱波方位的各向異性特征檢測裂縫是可行的。對比圖4a—圖4c 可以看到，各向異性參數(shù)在同一均值的條件下，參數(shù)γ(圖4b)對應(yīng)的橢圓扁率B/A與各向異性因子B的值最大，其次是參數(shù)δ(圖4c)，參數(shù)ε(圖4a)最小，因此，各向異性參數(shù)γ對橢圓扁率B/A與各向異性因子B影響最大、δ次之、ε最小。

2.2 非各向異性參數(shù)對疊前AVAZ 地震響應(yīng)影響

為了進(jìn)一步分析非各向異性參數(shù)對疊前AVAZ地震響應(yīng)特征的影響，以各向異性參數(shù)γ正態(tài)分布為例(各向異性參數(shù)γ對疊前AVAZ 響應(yīng)特征影響最大，模型參數(shù)γ設(shè)置為均值0.1、標(biāo)準(zhǔn)差0.01 的正態(tài)分布)，然后分別將模型中所研究的巖性參數(shù)根據(jù)表2 設(shè)置為正態(tài)分布，標(biāo)準(zhǔn)差按照模型4 取值，最后通過Monte Carlo 隨機(jī)正演模擬得到對應(yīng)橢圓扁率B/A與各向異性因子B散點(diǎn)圖(圖5)。

圖5a 和圖5d 分別是上層和下層的密度設(shè)置為正態(tài)分布(表2 中模型4)，縱波和橫波速度設(shè)置為均值時(shí)所對應(yīng)的橢圓扁率B/A與各向異性因子B散點(diǎn)圖，可以看出在密度標(biāo)準(zhǔn)差接近均值5%的情況下，其橢圓扁率B/A與各向異性因子B散點(diǎn)圖形態(tài)與圖4b 中γ為均值0.1 模型較為接近，說明密度參數(shù)不是影響橢圓信息的主控因素。圖5b 和圖5e 分別為上層介質(zhì)和下層介質(zhì)的橫波速度設(shè)置為正態(tài)分布(表2 中模型4)，縱波速度和密度設(shè)置為均值時(shí)所對應(yīng)的橢圓扁率B/A與各向異性因子B散點(diǎn)圖，當(dāng)橫波速度的標(biāo)準(zhǔn)差接近均值的10%時(shí)，其圖形的散點(diǎn)離散程度較密度正態(tài)分布時(shí)大，與圖4b 中γ為均值0.1 的模型結(jié)果差異較大，可以看出橫波速度隨機(jī)變化對AVAZ 響應(yīng)特征的影響較大。圖5c 和圖5f 分別為上層介質(zhì)和下層介質(zhì)的縱波速度設(shè)置為正態(tài)分布(表2 中模型4)，橫波速度和密度設(shè)置為均值時(shí)所對應(yīng)的橢圓扁率B/A與各向異性因子B散點(diǎn)圖，綜合來看，當(dāng)縱波速度所取的標(biāo)準(zhǔn)差接近均值的10%時(shí)，其圖形的散點(diǎn)擴(kuò)散程度最大，與圖4b 中γ為均值0.1 的模型結(jié)果差異很大，說明縱波速度的隨機(jī)變化對AVAZ 響應(yīng)特征的影響最大。

表2 不同標(biāo)準(zhǔn)差的巖性參數(shù)表Table 2 Lithology parameters with different standard deviations

圖5 不同非各向異性參數(shù)B/A(橢圓扁率)與B(各向異性因子)散點(diǎn)圖(γ 隨機(jī)變化、ε=0、δ=0)Fig.5 The scatter plots of B/A and B with different non-anisotropic parameters(γ random variation,ε=0,δ=0))

對式(12)進(jìn)行理論分析，可以看出公式中的密度參數(shù)對橢圓的擬合影響較小，速度參數(shù)影響較大。式(12)的理論分析結(jié)果與上述模型分析結(jié)果相符合，因此，模型速度與密度的不確定性會對模型模擬結(jié)果產(chǎn)生影響，其中縱波速度的影響最大，橫波速度的影響次之，密度的影響較小。

由于在實(shí)際地層中速度和密度并不是單獨(dú)變化，為了進(jìn)一步探究非各向異性參數(shù)的隨機(jī)變化對隨機(jī)正演的影響程度，在各向異性參數(shù)ε(圖6a、圖7a、圖8a)、γ(圖6b、圖7b、圖8b)、δ(圖6c、圖7c、圖8c)分別設(shè)置為均值為0.1，標(biāo)準(zhǔn)差為0.01 的正態(tài)分布的基礎(chǔ)上，將地層的速度和密度按照表2(模型1—模型4)設(shè)置為不同標(biāo)準(zhǔn)差的正態(tài)分布，采用Monte Carlo 隨機(jī)方法產(chǎn)生滿足條件的隨機(jī)模型，利用Ruger 近似式正演其疊前AVAZ 地震響應(yīng)，生成橢圓扁率B/A與各向異性因子B散點(diǎn)圖(圖6—圖8)。其中，圖6 為上層介質(zhì)的速度與密度按照模型1—模型4 設(shè)置為不同標(biāo)準(zhǔn)差的正態(tài)分布，下層介質(zhì)的速度與密度設(shè)置為恒定的均值；圖7 為下層介質(zhì)的速度與密度按照模型1—模型4 設(shè)置為不同標(biāo)準(zhǔn)差的正態(tài)分布，上層介質(zhì)的速度與密度設(shè)置為恒定的均值；圖8 為上、下兩層介質(zhì)的速度與密度均按照模型1—模型4 設(shè)置為不同標(biāo)準(zhǔn)差的正態(tài)分布。

圖6 不同標(biāo)準(zhǔn)差的速度與密度模型B/A(橢圓扁率)與B(各向異性因子)散點(diǎn)圖(上層介質(zhì)速度與密度均隨機(jī)變化)Fig.6 B/A and B scatter plots for velocity and density models with different standard deviations (upper layer velocity and density are set to a normal distribution)

圖7 不同標(biāo)準(zhǔn)差的速度與密度模型B/A(橢圓扁率)與B(各向異性因子)散點(diǎn)圖(下層介質(zhì)速度與密度均隨機(jī)變化)Fig.7 B/A and B scatter plots for velocity and density models with different standard deviations(The velocity and density settings of the upper layer are randomly distributed)

圖8 不同標(biāo)準(zhǔn)差的速度與密度模型B/A(橢圓扁率)與B(各向異性因子)散點(diǎn)圖(兩層介質(zhì)速度與密度均隨機(jī)變化)Fig.8 B/A and B scatter plots for velocity and density models with different standard deviations (both the velocity and density of two-layer medium change randomly)

在圖6—圖8 中均可看出：速度與密度設(shè)置為均值時(shí)(模型1)，其橢圓扁率的變化完全是由各向異性參數(shù)的隨機(jī)變化所引起，由于模型2—模型4 的速度與密度設(shè)置為正態(tài)分布，其擬合的橢圓信息的變化是由巖性參數(shù)和各向異性參數(shù)共同變化所引起的結(jié)果。分別將模型2—模型4 與模型1 對比分析：速度與密度在較小的范圍內(nèi)(模型2)隨機(jī)取值時(shí)，模擬出的橢圓信息可以很好地反映裂縫的信息；速度與密度的標(biāo)準(zhǔn)差分別為均值的5%和2.5%時(shí)(模型3)，模擬出的散點(diǎn)圖基本可以反映裂縫的信息；但當(dāng)速度與密度的標(biāo)準(zhǔn)差分別為均值的 10%和 5%時(shí) (模型4)，橢圓信息很大一部分受到速度和密度隨機(jī)波動的影響，其散點(diǎn)分布較為離散。

對比圖6–圖8 可以看出，兩層介質(zhì)(圖8)的速度與密度均為正態(tài)分布，比僅上層介質(zhì)(圖6)或僅下層介質(zhì)(圖7)的速度與密度為正態(tài)分布時(shí)所模擬出的橢圓扁率B/A與各向異性因子B散點(diǎn)圖更離散，由于樣本值滿足正態(tài)分布，其值大部分仍然分布在均值附近，所以其散點(diǎn)密集區(qū)基本一致。圖6—圖8 中間一列的b 圖為各向異性參數(shù)γ設(shè)置為概率密度分布，其B/A與B散點(diǎn)分布范圍最廣；圖6—圖8 左側(cè)一列的a 圖為各向異性參數(shù)ε設(shè)置為概率密度分布，其橢圓扁率B/A與各向異性因子B散點(diǎn)圖分布范圍最??；圖6—圖8 右側(cè)一列c 圖為各向異性參數(shù)δ設(shè)置為概率密度分布時(shí)，其橢圓扁率B/A與各向異性因子B散點(diǎn)圖分布范圍介于a 圖和b 圖之間。分析得出如下結(jié)論：

a.兩層介質(zhì)的速度與密度均設(shè)置為正態(tài)分布，較單層介質(zhì)(上層介質(zhì)或下層介質(zhì))的速度與密度正態(tài)分布對橢圓扁率B/A與各向異性因子B散點(diǎn)圖影響更大。

b.無論是儲層還是蓋層，由橫向巖性變化大而引起速度與密度的隨機(jī)變化時(shí)，橢圓的扁率B/A與各向異性因子B的變化范圍均會變大。當(dāng)橫向巖性變化過大即速度與密度隨機(jī)波動過大(例如模型 4 的速度隨機(jī)變化范圍是其平均值的10%，密度隨機(jī)變化范圍是其平均值的 5%)，會無法準(zhǔn)確識別出由各向異性所引起的變化。在不同巖性的裂縫性儲層中，地層的橫向巖性變化較小時(shí)即速度與密度隨機(jī)變化范圍在較小的范圍內(nèi)(標(biāo)準(zhǔn)差小于均值的10%)才能較好地利用該技術(shù)評價(jià)裂縫。

3 結(jié)論

a.以含EDA 介質(zhì)的雙層模型為基礎(chǔ)，將Ruger近似式中的各向異性參數(shù)和巖性參數(shù)設(shè)置為正態(tài)分布，通過Monte Carlo 隨機(jī)方法進(jìn)行疊前AVAZ 正演模擬擬合橢圓信息，分別分析了各向異性參數(shù)和巖性參數(shù)的隨機(jī)變化對擬合的橢圓的扁率B/A和各向異性因子B產(chǎn)生影響。模擬實(shí)驗(yàn)結(jié)果得知，在利用AVAZ 技術(shù)評價(jià)HTI 介質(zhì)中裂縫發(fā)育程度時(shí)，不僅各向異性參數(shù)ε、γ和δ會對AVAZ 響應(yīng)擬合的橢圓扁率B/A和各向異性因子B產(chǎn)生影響，而且非各向異性參數(shù)如地層的速度和密度(尤其是地層的縱波速度)的隨機(jī)變化也會對擬合的橢圓扁率B/A和各向異性因子B造成一定的影響。

b.各向異性參數(shù)γ對橢圓扁率B/A和各向異性因子B的影響最大，δ次之，ε最小。非各向異性參數(shù)對橢圓扁率B/A和各向異性因子B的影響程度依次為，縱波速度最大、橫波速度次之，密度最小。因此，縱波速度的不確定性是造成AVAZ 的響應(yīng)特征分布的主控因素，橫波速度的影響也較大，密度的影響較小。當(dāng)速度和密度的隨機(jī)波動范圍過大時(shí)，與理想的僅考慮各向異性參數(shù)變化的模型結(jié)果差別較大，所以當(dāng)速度與密度在較小的范圍內(nèi)隨機(jī)取值時(shí)，橢圓扁率B/A和各向異性因子B才可以較好地反映地層各向異性的強(qiáng)弱即裂縫的密度大小。

c.應(yīng)用疊前縱波方位各向異性檢測儲層裂縫時(shí)，應(yīng)當(dāng)充分考慮地層速度和密度橫向的變化對AVAZ 響應(yīng)特征分布的影響。在具體的實(shí)際資料分析中，為了減少地層速度、密度參數(shù)對預(yù)測結(jié)果的影響，在疊前三參數(shù)反演得到空間的縱橫波速度、密度的基礎(chǔ)上，下一步可研究如何利用空間的縱橫波速度、密度來提高橢圓扁率的計(jì)算精度，進(jìn)而提高該參數(shù)預(yù)測裂縫信息的準(zhǔn)確性。