取自然對數(shù)；系數(shù)Cl、C2、C3、C4、Cs W及Ce分別表示為：
[0035] Cl = d+tan" 0)/2
[0036] = -4化 sin 白）2
[0037] = tan2 日 /2-2 化 sin 目）2
[0038]
[0039]
[0040]
[0041] 其中，0和f分別為入射角和方位角；k為橫波速度與縱波速度的比值。
[0042] 在一個實施例中，第H模塊在建立所述誤差函數(shù)時，
[0043] 先由模擬疊前方位角度道集中的第i個樣點值和實際疊前方位角度道集中的第i 個樣點值建立方位角為取，入射角為0時的模擬疊前方位角度道集和實際疊前方位角度道 集的誤差函數(shù)；再按照先入射角后方位角的順序?qū)煌轿唤呛腿肷浣堑恼`差函數(shù)逐 一累加得到實際疊前方位角度道集和模擬疊前方位角度道集的誤差函數(shù)。
[0044] 在一個實施例中，第H模塊根據(jù)W下表達式建立實際疊前方位角度道集和模擬疊 前方位角度道集的誤差函數(shù)F :
[0045]
[0046] 式中，W為給定的地震子波；達:滅是實際疊前方位角度道集的第i個樣點值；m為界 面數(shù)，n為樣點數(shù)；Pl和P2為方位角的起始值和終止值，Ql和Q2為入射角的起始值和終止 值。J09為方位角為f，入射角為0時的模擬疊前方位角度道集和實際疊前方位角度道集的 誤差函數(shù)，的表達式為：
[0047] 在一個實施例中，第四模塊根據(jù)W下表達式得到所述敏感度矩陣G :
[004引式中，f朽巧;;，郎'說，,，，職娜，&巧布和巧W及0巧巧,分別為所述誤差函數(shù)F 對Zp,, Z,,，P，，5，，e , W及Y ,的偏導數(shù)，具體形式如下：
[005。 在一個實施例中，第四模塊根據(jù)W下表達式求心對Zp，，Z,,，P，，6，，e , W及Y j 的偏導數(shù)：
[0052]
[0053] 式中，,換*表示模擬疊前方位角度道集中第i個樣點值與實際疊前方位角度道集中 第i個樣點值的差值，f，的表達式為；知-;當X依次取Zpi，Zw，P i，6 j，e j '.始巧 扣-:; W及Y i時，每W獄表示腳j分別對Zpj，Zw，P i，6 i，e j W及Y j求偏導；W為給定的地震 子波；n為樣點數(shù)；A和B均為系數(shù)，其中，當X依次為Zp,，Zw，P .，，5，，e , W及Y ,時，A依 次取gpi，gsi，1/P i，1，I W及1，其中g(shù)pi，gsi分別為縱波導納和橫波導納；B是由入射角和方 位角決定的系數(shù)，當X依次為Zpi，Zsi, P i，6 i，e j W及Y j時，B依次取Cl、C2、C3、C4、Cs W 及Ce。
[0054] 與現(xiàn)有技術(shù)相比，上述方案中的一個或多個實施例可W具有如下優(yōu)點或有益效 果：
[0055] 基于實際測井數(shù)據(jù)與井旁地震資料，建立可用于各向異性參數(shù)反演的敏感度矩 陣，得到的敏感度矩陣信息可用于從疊前地震資料中提取各向異性參數(shù)的疊前各向異性參 數(shù)的反演工作，提高算法收斂速度和反演結(jié)果的精度，為尋找裂縫性油藏提供重要參考資 料。
[0056] 本發(fā)明的其它特征和優(yōu)點將在隨后的說明書中闡述，并且部分地從說明書中變得 顯而易見，或者通過實施本發(fā)明而了解。本發(fā)明的目的和其他優(yōu)點可通過在說明書、權(quán)利要 求書W及附圖中所特別指出的結(jié)構(gòu)來實現(xiàn)和獲得。
【附圖說明】
[0057] 附圖用來提供對本發(fā)明的進一步理解，并且構(gòu)成說明書的一部分，與本發(fā)明的實 施例共同用于解釋本發(fā)明，并不構(gòu)成對本發(fā)明的限制。在附圖中：
[0058] 圖1是根據(jù)本發(fā)明第一實施例的用于各向異性參數(shù)反演的敏感度矩陣的計算方 法的流程圖；
[0059] 圖2是某區(qū)域?qū)嶋H的測井數(shù)據(jù)圖；
[0060] 圖3 (a) - (f)是實際疊前方位角度道集；
[0061] 圖4是方位角為180度時的實際疊前角度道集和模擬疊前角度道集W及根據(jù)此時 的角度道集所建立的誤差函數(shù)的各偏導數(shù)曲線圖；
[0062] 圖5是某區(qū)域的敏感度矩陣曲線圖；
[0063] 圖6是根據(jù)本發(fā)明第二實施例的用于各向異性參數(shù)反演的敏感度矩陣的計算系 統(tǒng)的結(jié)構(gòu)示意圖。
【具體實施方式】
[0064] W下將結(jié)合附圖及實施例來詳細說明本發(fā)明的實施方式，借此對本發(fā)明如何應用 技術(shù)手段來解決技術(shù)問題，并達成技術(shù)效果的實現(xiàn)過程能充分理解并據(jù)W實施。需要說明 的是，只要不構(gòu)成沖突，本發(fā)明中的各個實施例W及各實施例中的各個特征可W相互結(jié)合， 所形成的技術(shù)方案均在本發(fā)明的保護范圍之內(nèi)。
[00巧] 第一連施例
[0066] 圖1是根據(jù)本發(fā)明第一實施例的用于各向異性參數(shù)反演的敏感度矩陣的計算方 法的流程圖。下面參考圖1來詳細說明本方法的各個步驟。
[0067] 步驟S110、基于實際測井數(shù)據(jù)建立縱波反射系數(shù)與縱波阻抗、橫波阻抗、密度W及 各向異性參數(shù)的泛函關(guān)系。
[0068] 具體的，基于表達式（1)建立縱波反射系數(shù)與縱波阻抗、橫波阻抗、密度W及各向 異性參數(shù)的泛函關(guān)系：
[0069]
[0070] 式中，咕'(?:?艦堿擇取0，的表示第j個界面的縱波反射系數(shù)，j表示界面的個數(shù)， 取值為j = 0. . . m ;Zp,為第j個界面的縱波阻抗，Zw為第j個界面的橫波阻抗，P，為第j 個界面的密度；5，、e ,和Y ,分別為第j個界面的各向異性參數(shù)；0和f分別為入射角和方 位角。表達式中的In表示取自然對數(shù)。系數(shù)Cl、C2、C3、C4、Cg W及Ce整體示于表達式（2) 中：
[0071] Cl = (l+tan2 日）/2
[0072] 。2 = -4化 sin0)2
[0073] c3 = tan2 日/2-2 化 sin 日）2 似
[0074]
[00 巧]
[0076]
[0077]
[0078] 式中，0和r分別為入射角和方位角，k為橫波速度與縱波速度的比值，通常取為 定值0.5。
[0079] 進一步的，令L = Inx,對表達式（1)進行簡化，得到表達式（3);
[0080]
[008。 式中，ALp、ALs、ALp、AS、AeW及AY分別具有表達式（4)所示的形式：
[0082] A Lp = Lp j"-Lpj = In 狂pj") -In 狂P j)
[0083] A Ls = Ls j"-Lsj = In 狂Si") -In 狂Si)
[0084] ALp = Lp.j"-Lp.j = ln(p .j")-ln(p .j)
[0085] A 5 = 5 j"- S j (4)
[0086] A e = e j"- e j
[0087] A Y = Y jW-Y j
[008引通過步驟SllO,最終得到了如表達式（3)所示的縱波反射系數(shù)rp,的模型，將模型 中的各參數(shù)逐個代入可W求得特定方位角和入射角時的縱波反射系數(shù)r。，的值。
[0089] 步驟S120、基于權(quán)積原理，利用所建立的泛函關(guān)系與給定的地震子波得到模擬疊 前方位角度道集。
[0090] 具體的，基于權(quán)積原理，通過表達式（5)計算得到模擬疊前方位角度道集：
[0091]
[0092] 式中，A鮮為模擬疊前方位角度道集的第i個樣點值，r。，為縱波反射系數(shù)，W為給 定的地震子波，m為界面數(shù)，n為樣點數(shù)。