反射系數(shù)與縱波阻抗、橫波阻抗、密度以及各向異性參數(shù) 的泛函關(guān)系； 基于褶積原理，利用所建立的泛函關(guān)系與給定的地震子波得到模擬疊前方位角度道 集； 根據(jù)實(shí)際疊前方位角度道集和模擬疊前方位角度道集建立誤差函數(shù)； 對(duì)誤差函數(shù)求偏導(dǎo)數(shù)得到用于各向異性參數(shù)反演的敏感度矩陣。2. 根據(jù)權(quán)利要求1所述的方法，其特征在于，根據(jù)以下表達(dá)式建立所述泛函關(guān)系： rp j = C1 (lnZpj+「lnZpj) +c2 (lnZsj+「lnZsj) +c3 (In P j+「ln P j) +C4 ( δ j+「δ j) +。5 ( ε j+「ε j )+c6 ( γ .j+「γ .j) 式中，rp]表示第j個(gè)界面的縱波反射系數(shù)，j表示界面的個(gè)數(shù)；Zw為第j個(gè)界面的縱 波阻抗，Zs]為第j個(gè)界面的橫波阻抗，P ,為第j個(gè)界面的密度；δ ρ ε ,和Y ,分別為第j 個(gè)界面的各向異性參數(shù)；In表示取自然對(duì)數(shù)；系數(shù)Cl、c2、c3、c 4、C5以及C6分別為： C1 = (1+tan2 θ ) /2 c2 = -4(k sin θ)2 c3 = tan2 θ /2-2 (k sin θ )2其中，Θ和f分別為入射角和方位角；k為橫波速度與縱波速度的比值。3. 根據(jù)權(quán)利要求2所述的方法，其特征在于，在建立所述誤差函數(shù)的步驟中， 先由模擬疊前方位角度道集中的第i個(gè)樣點(diǎn)值和實(shí)際疊前方位角度道集中的第i個(gè)樣 點(diǎn)值建立方位角為，入射角為Θ時(shí)的模擬疊前方位角度道集和實(shí)際疊前方位角度道集的 誤差函數(shù)； 再按照先入射角后方位角的順序?qū)?duì)應(yīng)不同方位角和入射角的誤差函數(shù)逐一累加得 到實(shí)際疊前方位角度道集和模擬疊前方位角度道集的誤差函數(shù)。4. 根據(jù)權(quán)利要求3所述的方法，其特征在于，根據(jù)以下表達(dá)式建立實(shí)際疊前方位角度 道集和模擬疊前方位角度道集的誤差函數(shù)F :式中，w為給定的地震子波；是實(shí)際疊前方位角度道集的第i個(gè)樣點(diǎn)值；m為界面 數(shù)，η為樣點(diǎn)數(shù)；Pl和P2為方位角的起始值和終止值，Ql和Q2為入射角的起始值和終止 值；>V?為方位角為-入射角為Θ時(shí)的模擬疊前方位角度道集和實(shí)際疊前方位角度道集的 誤差函數(shù)，iw)i的表達(dá)式為：5. 根據(jù)權(quán)利要求4所述的方法，其特征在于，根據(jù)以下表達(dá)式得到所述敏感度矩陣G :式中以及分別為所述誤差函數(shù)F對(duì) Zw，Z、i，P p s P ε li以及Y j的偏導(dǎo)數(shù)，具體形式如下：6. 根據(jù)權(quán)利要求5所述的方法，其特征在于，根據(jù)以下表達(dá)式求》對(duì)ZP]，ZS]，P ,，δ ,， ε _j以及Y _j的偏導(dǎo)數(shù)：式中，爲(wèi)a表示模擬疊前方位角度道集中第i個(gè)樣點(diǎn)值與實(shí)際疊前方位角度道集中第i 個(gè)樣點(diǎn)值的差值，^的表達(dá)式為：；當(dāng)X依次取zpj，Zsj，P j，δ j，ε j以 及Y j時(shí)，dVAiW表示9分別對(duì)Zpj，Zsj，P j，δ y ε ;以及γ ;求偏導(dǎo)數(shù)；w為給定的地震 子波；η為樣點(diǎn)數(shù)；A和B均為系數(shù),其中，當(dāng)X依次為Zw, Zs.j, P .j, δ ^ ε 以及γ 時(shí)，A依 次取gw，gsj，Ι/P y 1，1以及1，其中g(shù)w，gsj分別為縱波導(dǎo)納和橫波導(dǎo)納；B是由入射角和方 位角決定的系數(shù)，當(dāng)X依次為Z pj，Zsj，P j，δ j，ε j以及Y j時(shí)，B依次取Cl、c2、c3、c4、C 5以 及 Cf;。7. -種用于各向異性參數(shù)反演的敏感度矩陣的計(jì)算系統(tǒng)，包括以下模塊： 第一模塊，其基于實(shí)際測(cè)井?dāng)?shù)據(jù)建立縱波反射系數(shù)與縱波阻抗、橫波阻抗、密度以及各 向異性參數(shù)的泛函關(guān)系； 第二模塊，其基于褶積原理，利用所建立的泛函關(guān)系與給定的地震子波得到模擬疊前 方位角度道集； 第三模塊，其根據(jù)實(shí)際疊前方位角度道集和模擬疊前方位角度道集建立誤差函數(shù)； 第四模塊，其對(duì)誤差函數(shù)求偏導(dǎo)數(shù)得到用于各向異性參數(shù)反演的敏感度矩陣。8. 根據(jù)權(quán)利要求7所述的系統(tǒng)，其特征在于，所述第一模塊根據(jù)以下表達(dá)式建立所述 泛函關(guān)系： rp j = C1 (lnZpj+「lnZpj) +c2 (lnZsj+「lnZsj) +c3 (In P j+「ln P j) +C4 ( δ j+「δ j) +。5 ( ε j+「ε j )+c6 ( γ .j+「γ .j) 式中，rp]表示第j個(gè)界面的縱波反射系數(shù)，j表示界面的個(gè)數(shù)；Zw為第j個(gè)界面的縱 波阻抗，Zs]為第j個(gè)界面的橫波阻抗，P ,為第j個(gè)界面的密度；δ ρ ε ,和Y ,分別為第j 個(gè)界面的各向異性參數(shù)；In表示取自然對(duì)數(shù)；系數(shù)CpCpCyCzp C5以及C6分別表示為： C1 = (1+tan2 Θ ) /2 c2 = -4(k sin θ)2其中，Θ和P分別為入射角和方位角；k為橫波速度與縱波速度的比值。9. 根據(jù)權(quán)利要求8所述的系統(tǒng)，其特征在于，所述第三模塊在建立所述誤差函數(shù)時(shí)， 先由模擬疊前方位角度道集中的第i個(gè)樣點(diǎn)值和實(shí)際疊前方位角度道集中的第i個(gè)樣 點(diǎn)值建立方位角為心入射角為Θ時(shí)的模擬疊前方位角度道集和實(shí)際疊前方位角度道集的 誤差函數(shù)； 再按照先入射角后方位角的順序?qū)?duì)應(yīng)不同方位角和入射角的誤差函數(shù)逐一累加得 到實(shí)際疊前方位角度道集和模擬疊前方位角度道集的誤差函數(shù)。10. 根據(jù)權(quán)利要求9所述的系統(tǒng)，其特征在于，所述第三模塊根據(jù)以下表達(dá)式建立實(shí)際 疊前方位角度道集和模擬疊前方位角度道集的誤差函數(shù)F :式中，w為給定的地震子波；是實(shí)際疊前方位角度道集的第i個(gè)樣點(diǎn)值；m為界面 數(shù)，η為樣點(diǎn)數(shù)；Pl和P2為方位角的起始值和終止值，Ql和Q2為入射角的起始值和終止 值；.V#為方位角為，入射角為Θ時(shí)的模擬疊前方位角度道集和實(shí)際疊前方位角度道集的 誤差函數(shù)，.Si 5的表達(dá)式為：11. 根據(jù)權(quán)利要求10所述的系統(tǒng)，其特征在于，所述第四模塊根據(jù)以下表達(dá)式得到所 述敏感度矩陣G :Zw，Z、i，μ，s p ε i以及的偏導(dǎo)數(shù)，具體形式如下：12. 根據(jù)權(quán)利要求11所述的系統(tǒng)，其特征在于，所述第四模塊根據(jù)以下表達(dá)式求對(duì) zw，Zs.j，p .j，δ .j，ε .j 以及 Y .j 的偏導(dǎo)數(shù)：式中，表示模擬疊前方位角度道集中第i個(gè)樣點(diǎn)值與實(shí)際疊前方位角度道集中第i 個(gè)樣點(diǎn)值的差值，的表達(dá)式為：：當(dāng)X依次取Zpj，Zsj，P j，δ j，ε j以 及Y i時(shí)，表示鄭分別對(duì)Zpj，Zsj，P j，δ y ε j以及Y j求偏導(dǎo)；w為給定的地震子 波；η為樣點(diǎn)數(shù)；A和B均為系數(shù),其中，當(dāng)X依次為Zw, Zs.j, P .j, δ ^ ε 以及γ 時(shí)，A依次 取gw，gsj，1，1以及1，其中g(shù)w，g sj分別為縱波導(dǎo)納和橫波導(dǎo)納；B是由入射角和方位 角決定的系數(shù)，當(dāng)X依次為Zpj，Z sj，P p δ p ε j以及Y j時(shí)，B依次取Cl、c2、c3、c4、C 5以及 C6O
【專利摘要】本發(fā)明公開了一種用于各向異性參數(shù)反演的敏感度矩陣的計(jì)算方法及系統(tǒng)，該方法包括以下步驟：基于實(shí)際測(cè)井?dāng)?shù)據(jù)建立縱波反射系數(shù)與縱波阻抗、橫波阻抗、密度以及各向異性參數(shù)的泛函關(guān)系；基于褶積原理，利用所建立的泛函關(guān)系與給定的地震子波得到模擬疊前方位角度道集；根據(jù)實(shí)際疊前方位角度道集和模擬疊前方位角度道集建立誤差函數(shù)；對(duì)誤差函數(shù)求偏導(dǎo)數(shù)得到用于各向異性參數(shù)反演的敏感度矩陣。得到的敏感度矩陣信息可用于從疊前地震資料中提取各向異性參數(shù)的疊前各向異性參數(shù)反演工作，提高算法收斂速度和反演結(jié)果的精度，為尋找裂縫性油藏提供重要參考資料。
【IPC分類】G01V1/28
【公開號(hào)】CN105527648
【申請(qǐng)?zhí)枴緾N201410520158
【發(fā)明人】白俊雨
【申請(qǐng)人】中國石油化工股份有限公司, 中國石油化工股份有限公司石油物探技術(shù)研究院
【公開日】2016年4月27日
【申請(qǐng)日】2014年9月30日