一種基于巖石微結構及介質(zhì)彈性參數(shù)的地震激發(fā)模擬方法
【專利摘要】本發(fā)明提供了一種基于巖石微結構及介質(zhì)彈性參數(shù)的地震激發(fā)模擬方法，屬于地球物理勘探領域。本發(fā)明方法包括以下步驟：(1)巖石微結構調(diào)查和求取介質(zhì)彈性參數(shù)；(2)建立研究區(qū)近地表速度模型、深層速度模型和構造模型；(3)數(shù)值模擬算法實現(xiàn)及建立包含彈性介質(zhì)的近震源波場函數(shù)；(4)進行基于巖石微結構及介質(zhì)彈性參數(shù)的地震激發(fā)模擬。利用本發(fā)明方法可以進行基于巖石微結構和介質(zhì)彈性參數(shù)的地震激發(fā)模擬，并建立圍巖激發(fā)特征數(shù)據(jù)庫，最大程度在室內(nèi)完成不同地區(qū)激發(fā)參數(shù)的選擇。
【專利說明】一種基于巖石微結構及介質(zhì)彈性參數(shù)的地震激發(fā)模擬方法
【技術領域】
[0001] 本發(fā)明屬于地球物理勘探領域，具體涉及一種基于巖石微結構及介質(zhì)彈性參數(shù)的 地震激發(fā)模擬方法。
【背景技術】
[0002] 對于特殊構造地區(qū)（如山前帶灰?guī)r出露區(qū)）來說，提高資料信噪比是必須要解決 的首要問題，而激發(fā)是影響單炮品質(zhì)的主要因素之一。彈性介質(zhì)近震源波場研究表明：地震 激發(fā)與介質(zhì)的彈性參數(shù)和結構相關，激發(fā)介質(zhì)的巖石力學性質(zhì)和物理結構決定激發(fā)的振幅 特性和頻率特性，所以震源激發(fā)地震子波的特性如何，很大程度影響地震資料的品質(zhì)。目前 波場模擬技術都是運用雷克子波作為震源子波進行模擬，再對模擬結果進行分析，這樣只 是簡單地分析地質(zhì)構造對波場特征的影響，沒有從本質(zhì)上討論不同巖性地層對激發(fā)子波的 影響，進而對地震資料的影響，導致野外實際采集施工中，炸藥激發(fā)參數(shù)（藥量、井深、組合 和爆速等）的選擇必須進行大量的試驗來確定，勢必造成人力、財力和物力的損失，還不安 全。因此有必要開展基于巖石微結構和介質(zhì)彈性參數(shù)的地震激發(fā)模擬技術，并建立模塊化 的圍巖激發(fā)特征數(shù)據(jù)庫，最大程度在室內(nèi)完成不同地區(qū)激發(fā)參數(shù)的選擇，最終有效指導野 外地震生產(chǎn)。
【發(fā)明內(nèi)容】
[0003] 本發(fā)明的目的在于解決上述現(xiàn)有技術中存在的難題，提供一種基于巖石微結構及 介質(zhì)彈性參數(shù)的地震激發(fā)模擬方法，通過在室內(nèi)開展不同特殊構造地區(qū)基于巖石微結構和 介質(zhì)彈性參數(shù)的地震激發(fā)模擬，進行巖石微結構調(diào)查和應用巖石物理實驗來測定介質(zhì)的彈 性參數(shù)，進而進行基于巖石微結構及介質(zhì)彈性參數(shù)的地震激發(fā)模擬，由模擬結果未定量 化地確定野外實際地層中的各地區(qū)地震采集炸藥震源激發(fā)參數(shù)（藥量、井深、組合和爆速 等），使激發(fā)的地震子波振幅和頻率特性最優(yōu)化，最終使采集到的地震數(shù)據(jù)的品質(zhì)得到極大 改善。
[0004] 本發(fā)明是通過以下技術方案實現(xiàn)的：
[0005] -種基于巖石微結構及介質(zhì)彈性參數(shù)的地震激發(fā)模擬方法，包括：
[0006] (1)巖石微結構調(diào)查和求取介質(zhì)彈性參數(shù)；
[0007] (2)建立研究區(qū)近地表速度模型、深層速度模型和構造模型；
[0008] (3)數(shù)值模擬算法實現(xiàn)及建立包含彈性介質(zhì)的近震源波場函數(shù)；
[0009] (4)進行基于巖石微結構及介質(zhì)彈性參數(shù)的地震激發(fā)模擬。
[0010] 所述步驟⑴包括：
[0011] (11)對從研究區(qū)收集來的實際巖樣進行電鏡掃描完成巖樣的微結構調(diào)查；
[0012] (12)不同圍壓條件下介質(zhì)彈性參數(shù)關系的實驗和分析以及不同飽和度條件下介 質(zhì)彈性參數(shù)關系的實驗和分析。
[0013] 所述步驟（12)是這樣實現(xiàn)的：
[0014] 通過實驗室測量不同圍壓條件下所取巖芯的縱橫波速度和泊松比，以及介質(zhì)的彈 性參數(shù)，獲得縱橫波速度和泊松比與介質(zhì)彈性參數(shù)之間的關系；
[0015] 通過實驗室測量不同飽和度條件下所取巖芯的縱橫波速度和泊松比，以及介質(zhì)的 彈性參數(shù)，獲得縱橫波速度和泊松比與介質(zhì)彈性參數(shù)之間的關系；
[0016] 進而建立適合該地區(qū)近地表條件下的介質(zhì)速度和彈性參數(shù)的數(shù)學模型。
[0017] 所述步驟（2)是這樣實現(xiàn)的：通過研究區(qū)的地震疊后資料建立構造模型，通過研 究區(qū)的鉆井、VSP測井資料建立深層速度模型，通過步驟（1)得到的介質(zhì)速度和彈性參數(shù)的 數(shù)學模型以及步驟（1)得到的巖樣的微結構調(diào)查建立近地表速度模型；最后將近地表層速 度和深層速度模型合并得到最終的速度模型。
[0018] 所述步驟（3)中數(shù)值模擬算法實現(xiàn)包括：
[0019] (31)建立波動方程的一階速度-應力方程；
[0020] (32)實現(xiàn)交錯網(wǎng)格有限差分法的數(shù)值模擬算法：在步驟（31)中的一階速度-應 力方程基礎上得到相應的交錯網(wǎng)格差分格式，進而形成數(shù)值模擬算法。
[0021] 所述步驟（3)建立的包含彈性介質(zhì)的近震源波場函數(shù)為：
【權利要求】
1. 一種基于巖石微結構及介質(zhì)彈性參數(shù)的地震激發(fā)模擬方法，其特征在于：所述方法 包括： (1) 巖石微結構調(diào)查和求取介質(zhì)彈性參數(shù)； (2) 建立研究區(qū)近地表速度模型、深層速度模型和構造模型； (3) 數(shù)值模擬算法實現(xiàn)及建立包含彈性介質(zhì)的近震源波場函數(shù)； (4) 進行基于巖石微結構及介質(zhì)彈性參數(shù)的地震激發(fā)模擬。
2. 根據(jù)權利要求1所述的基于巖石微結構及介質(zhì)彈性參數(shù)的地震激發(fā)模擬方法，其特 征在于：所述步驟（1)包括： (11) 對從研究區(qū)收集來的實際巖樣進行電鏡掃描完成巖樣的微結構調(diào)查； (12) 不同圍壓條件下介質(zhì)彈性參數(shù)關系的實驗和分析以及不同飽和度條件下介質(zhì)彈 性參數(shù)關系的實驗和分析。
3. 根據(jù)權利要求2所述的基于巖石微結構及介質(zhì)彈性參數(shù)的地震激發(fā)模擬方法，其特 征在于：所述步驟（12)是這樣實現(xiàn)的： 通過實驗室測量不同圍壓條件下所取巖芯的縱橫波速度和泊松比，以及介質(zhì)的彈性參 數(shù)，獲得縱橫波速度和泊松比與介質(zhì)彈性參數(shù)之間的關系； 通過實驗室測量不同飽和度條件下所取巖芯的縱橫波速度和泊松比，以及介質(zhì)的彈性 參數(shù)，獲得縱橫波速度和泊松比與介質(zhì)彈性參數(shù)之間的關系； 進而建立適合該地區(qū)近地表條件下的介質(zhì)速度和彈性參數(shù)的數(shù)學模型。
4. 根據(jù)權利要求3所述的基于巖石微結構及介質(zhì)彈性參數(shù)的地震激發(fā)模擬方法，其特 征在于：所述步驟（2)是這樣實現(xiàn)的：通過研究區(qū)的地震疊后資料建立構造模型，通過研究 區(qū)的鉆井、VSP測井資料建立深層速度模型，通過步驟（1)得到的介質(zhì)速度和彈性參數(shù)的數(shù) 學模型以及步驟（1)得到的巖樣的微結構調(diào)查建立近地表速度模型；最后將近地表層速度 和深層速度模型合并得到最終的速度模型。
5. 根據(jù)權利要求4所述的基于巖石微結構及介質(zhì)彈性參數(shù)的地震激發(fā)模擬方法，其特 征在于：所述步驟（3)中數(shù)值模擬算法實現(xiàn)包括： (31) 建立波動方程的一階速度-應力方程； (32) 實現(xiàn)交錯網(wǎng)格有限差分法的數(shù)值模擬算法：在步驟（31)中的一階速度-應力方程 基礎上得到相應的交錯網(wǎng)格差分格式，進而形成數(shù)值模擬算法。
7. 根據(jù)權利要求6所述的基于巖石微結構及介質(zhì)彈性參數(shù)的地震激發(fā)模擬方法，其特 征在于：所述步驟（4)包括： (41) 完成包含彈性介質(zhì)的近震源波場函數(shù)的加載； (42) 進行基于巖石微結構及介質(zhì)彈性參數(shù)的地震激發(fā)模擬：運用步驟（3)得到的數(shù)值 模擬算法得到基于介質(zhì)彈性參數(shù)的激發(fā)模擬記錄。
8. 根據(jù)權利要求7所述的基于巖石微結構及介質(zhì)彈性參數(shù)的地震激發(fā)模擬方法，其特 征在于：所述步驟（41)包括： (411) :對步驟（2)得到的最終的速度模型進行網(wǎng)格剖分，根據(jù)水平和垂直網(wǎng)格點坐標 確定震源位置（NXs，NZs); (412) :在以震源位置（NXs，NZs)為中心的周圍四個網(wǎng)格點上加載用所述彈性介質(zhì)的近 震源波場函數(shù)表示的震源子波。
【文檔編號】G01V1/104GK104122578SQ201310142832
【公開日】2014年10月29日 申請日期:2013年4月23日 優(yōu)先權日:2013年4月23日
【發(fā)明者】肖云飛, 殷厚成, 王躍 申請人:中國石油化工股份有限公司, 中國石油化工股份有限公司石油物探技術研究院