獲取風化層的地質參數的方法技術

本發明專利技術提供一種獲取風化層的地質參數的方法，所述方法包括：(a)獲取目標區域的標準的初至數據；(b)確定最小偏移距；(c)確定每個點的風化層速度值；(d)計算選擇的m個點中的任意兩點之間的風化層速度值的變差函數值；(e)計算選擇的n個點中的任意兩點之間的微測井資料的速度值的變差函數值；(f)計算m個點中的任意一點與n個點中的任意一點之間的關于速度的交叉編譯函數；(g)獲得與m個點中的每個點的風化層速度值對應的m個加權系數和與n個點中的每個點的微測井數據的速度值對應的n個加權系數；(h)計算風化層速度的估計值。采用本發明專利技術獲取風化層的地質參數的方法可得到更符合地質分布規律的淺層速度結構。

【技術實現步驟摘要】
獲取風化層的地質參數的方法
本專利技術總體來說涉及石油地震資料分析處理領域，更具體地講，涉及一種基于協同克里金技術融合微測井資料的獲取風化層的地質參數的方法。
技術介紹
在地震勘探中，近地表的低降速帶(即，風化層)的速度和厚度不均勻變化會引起地震波傳播旅行時的不均勻延遲，致使實際觀測得到的反射波到達時間發生畸變，極大地影響著地震剖面的成像和構造形態。靜校正就是消除地表低降速帶對勘探資料的影響，對于復雜地表條件下地震資料的處理，求取正確的靜校正量是非常關鍵和必要的。如果不能很好的解決嚴重的靜校正問題，就會在地震剖面上產生假象或不成像，就會給地質解釋帶來嚴重困難，甚至錯誤的地質構造認識。目前，在解決靜校正問題過程中，對于存在穩定折射層的區域，一般采用折射靜校正方法(例如，基于迭代法的折射靜校正方法和廣義互換靜校正方法)進行計算。但這兩類方法在使用時，對于風化層速度，一般參照工作經驗，給定常速度值進行計算。然而，很多地區的風化層速度橫向變化劇烈，實際風化層速度與前述參照工作經驗給定的常速度值存在較大差距。此外，有的方法還會對實際風化層速度進行擬合計算或者結合微測井資料進行速度反演，其效果優于給定常速度值計算的結果。但采用速度擬合的方法仍存在地震剖面的成像橫向不連續現象，需要進一步處理；而進行速度反演，則計算時間長，微測井資料的好壞將直接影響速度反演的結果，不確定程度高，從而影響處理后疊加剖面成像效果，使得構造形態不準確，信噪比大大降低，不能達到資料精細處理的目的。
技術實現思路
本專利技術的目的是提出一種基于協同克里金技術融合微測井資料的獲取風化層的地質參數的方法，通過變異函數這一數學工具，分別刻畫出微測井資料的速度值、厚度值和由初至數據擬合出的風化層速度、厚度的變化趨勢，再運用交叉編譯函數得到兩類數據的整體變化趨勢，并加以它們的相關系數約束，最終得到符合地質分布規律的淺層速度結構。本專利技術的一方面提供一種獲取風化層的地質參數的方法，所述方法包括：(a)獲取目標區域的標準的初至數據；(b)確定最小偏移距；(c)在所述目標區域內選擇m個點，從獲取的初至數據中選擇m個點中的每個點在最小偏移距范圍內的初至數據，并對選擇的m個點中的每個點在最小偏移距范圍內對應的初至數據進行擬合，確定每個點的風化層速度值，m為大于等于1的自然數；(d)根據確定的每個點的風化層速度值，計算選擇的m個點中的任意兩點之間的風化層速度值的變差函數值；(e)在所述目標區域內選擇n個點，從所述目標區域的微測井資料中提取n個點中的每個點的微測井資料的速度值，并根據提取的每個點的微測井資料的速度值，計算選擇的n個點中的任意兩點之間的微測井資料的速度值的變差函數值，n為大于等于1的自然數；(f)根據確定的m個點中的每個點的風化層速度值和n個點中的每個點的微測井資料的速度值，計算m個點中的任意一點與n個點中的任意一點之間的關于速度的交叉編譯函數；(g)根據步驟(d)～步驟(f)的結果獲得與m個點中的每個點的風化層速度值對應的m個加權系數和與n個點中的每個點的微測井數據的速度值對應的n個加權系數；(h)根據獲得的與m個點中的每個點的風化層速度值對應的m個加權系數、與n個點中的每個點的微測井數據的速度值對應的n個加權系數、確定的m個點中的每個點的風化層速度值和n個點中的每個點的微測井資料的速度值，計算風化層速度的估計值。可選地，步驟(g)可包括：根據步驟(d)～步驟(f)的結果利用協同克里金方程組獲得與m個點中的每個點的風化層速度值對應的m個加權系數和與n個點中的每個點的微測井數據的速度值對應的n個加權系數。可選地，步驟(g)可包括：利用下面的協同克里金方程組計算與m個點中的每個點的風化層速度值對應的m個加權系數和與n個點中的每個點的微測井數據的速度值對應的n個加權系數，其中，Cx(xi,xp)為n個點中的第i點的坐標位置xi與n個點中的第p點的坐標位置xp對應的微測井資料的速度值的變差函數值，Cx(0,xi)為待求位置0與n個點中的第i點的坐標位置xi對應的微測井資料的速度值的變差函數值，Cy(yq,yj)為m個點中的第q點的坐標位置yq與m個點中的第j點的坐標位置yj對應的風化層速度值的變差函數值，Cy(0,yj)為待求位置0與m個點中的第j點的坐標位置yj對應的風化層速度值的變差函數值，Cc(xi,yj)和Cc(yj,xi)為關于n個點中的第i點的坐標位置xi對應的微測井資料的速度值與m個點中的第j點的坐標位置yj對應的風化層速度值的交叉編譯函數，為協同克里金方程組中與n個點中的第i點的微測井數據的速度值對應的加權系數，為協同克里金方程組中與m個點中的第j點的風化層速度值對應的加權系數，μx、μy為計算方程組的解與原始值的誤差的參數。可選地，步驟(h)可包括：利用下面的公式計算風化層速度的估計值Z*(0)，其中，Z*(0)為風化層速度在待求位置0處的估計值，Zx(xi)為n個點中的第i點的坐標位置xi對應的微測井資料的速度值，為與n個點中的第i點的微測井數據的速度值對應的加權系數，Zy(yj)為m個點中的第j點的坐標位置yj對應的風化層速度值，為與m個點中的第j點的風化層速度值對應的加權系數。可選地，步驟(c)可還包括：對選擇的m個點中的每個點在最小偏移距范圍內對應的初至數據進行擬合，確定每個點的炮點延遲時的初始值，其中，所述方法可還包括：(o)確定最大偏移距；(p)從獲取的初至數據中選擇m個點中的每個點在最小偏移距至最大偏移距范圍內的初至數據，對選擇的m個點中的每個點在最小偏移距至最大偏移距范圍內對應的初至數據進行擬合，確定每個點的高速層速度的初始值；(q)根據確定的m個點中的每個點的高速層速度的初始值和m個點中的每個炮點延遲時的初始值，利用折射初至時間公式依次迭代計算m個點中的每個點的檢波點延遲時、炮點延遲時和高速層速度。可選地，步驟(q)可包括：(q1)根據m個點中的每個點的高速層速度、炮點延遲時、初至數據中每個點對應的炮檢對從激發到接收產生的時間、每個點對應的炮點與檢波點之間的距離，利用折射初至時間公式計算m個點中的每個點的檢波點延遲時，其中，將m個點中的每個點的高速層速度的初始值和m個點中的每個點的炮點延遲時的初始值作為初值代入；(q2)根據m個點中的每個點的高速層速度、步驟(q1)得到的m個點中的每個點的檢波點延遲時、初至數據中每個點對應的炮檢對從激發到接收產生的時間、每個點對應的炮點與檢波點之間的距離，利用折射初至時間公式計算m個點中的每個點的炮點延遲時；(q3)根據步驟(q1)得到的m個點中的每個點的檢波點延遲時、步驟(q2)得到的m個點中的每個點的炮點延遲時、初至數據中每個點對應的炮檢對從激發到接收產生的時間、每個點對應的炮點與檢波點之間的距離，利用折射初至時間公式計算m個點中的每個點的高速層速度；(q4)根據步驟(q1)得到的m個點中的每個點的檢波點延遲時、步驟(q2)得到的m個點中的每個點的炮點延遲時、步驟(q3)得到的m個點中的每個點的高速層速度、每個點對應的炮點與檢波點之間的距離，利用折射初至時間公式計算對應的初至數據中每個點對應的炮檢對從激發到接收產生的時間；(q5)當所述時間小于等于設定的時間范圍時，將步驟(q1)～步驟本文檔來自技高網...

【技術保護點】
一種獲取風化層的地質參數的方法，所述方法包括：(a)獲取目標區域的標準的初至數據；(b)確定最小偏移距；(c)在所述目標區域內選擇m個點，從獲取的初至數據中選擇m個點中的每個點在最小偏移距范圍內的初至數據，并對選擇的m個點中的每個點在最小偏移距范圍內對應的初至數據進行擬合，確定每個點的風化層速度值，m為大于等于1的自然數；(d)根據確定的每個點的風化層速度值，計算選擇的m個點中的任意兩點之間的風化層速度值的變差函數值；(e)在所述目標區域內選擇n個點，從所述目標區域的微測井資料中提取n個點中的每個點的微測井資料的速度值，并根據提取的每個點的微測井資料的速度值，計算選擇的n個點中的任意兩點之間的微測井資料的速度值的變差函數值，n為大于等于1的自然數；(f)根據確定的m個點中的每個點的風化層速度值和n個點中的每個點的微測井資料的速度值，計算m個點中的任意一點與n個點中的任意一點之間的關于速度的交叉編譯函數；(g)根據步驟(d)～步驟(f)的結果獲得與m個點中的每個點的風化層速度值對應的m個加權系數和與n個點中的每個點的微測井數據的速度值對應的n個加權系數；(h)根據獲得的與m個點中的每個點的風化層速度值對應的m個加權系數、與n個點中的每個點的微測井數據的速度值對應的n個加權系數、確定的m個點中的每個點的風化層速度值和n個點中的每個點的微測井資料的速度值，計算風化層速度的估計值。...

【技術特征摘要】
1.一種獲取風化層的地質參數的方法，所述方法包括：(a)獲取目標區域的標準的初至數據；(b)確定最小偏移距；(c)在所述目標區域內選擇m個點，從獲取的初至數據中選擇m個點中的每個點在最小偏移距范圍內的初至數據，并對選擇的m個點中的每個點在最小偏移距范圍內對應的初至數據進行擬合，確定每個點的風化層速度值，m為大于等于1的自然數；(d)根據確定的每個點的風化層速度值，計算選擇的m個點中的任意兩點之間的風化層速度值的變差函數值；(e)在所述目標區域內選擇n個點，從所述目標區域的微測井資料中提取n個點中的每個點的微測井資料的速度值，并根據提取的每個點的微測井資料的速度值，計算選擇的n個點中的任意兩點之間的微測井資料的速度值的變差函數值，n為大于等于1的自然數；(f)根據確定的m個點中的每個點的風化層速度值和n個點中的每個點的微測井資料的速度值，計算m個點中的任意一點與n個點中的任意一點之間的關于速度的交叉編譯函數；(g)根據步驟(d)～步驟(f)的結果獲得與m個點中的每個點的風化層速度值對應的m個加權系數和與n個點中的每個點的微測井數據的速度值對應的n個加權系數；(h)根據獲得的與m個點中的每個點的風化層速度值對應的m個加權系數、與n個點中的每個點的微測井數據的速度值對應的n個加權系數、確定的m個點中的每個點的風化層速度值和n個點中的每個點的微測井資料的速度值，計算風化層速度的估計值，其中，步驟(c)還包括：對選擇的m個點中的每個點在最小偏移距范圍內對應的初至數據進行擬合，確定每個點的炮點延遲時的初始值，其中，所述方法還包括：(o)確定最大偏移距；(p)從獲取的初至數據中選擇m個點中的每個點在最小偏移距至最大偏移距范圍內的初至數據，對選擇的m個點中的每個點在最小偏移距至最大偏移距范圍內對應的初至數據進行擬合，確定每個點的高速層速度的初始值；(q)根據確定的m個點中的每個點的高速層速度的初始值和m個點中的每個點的炮點延遲時的初始值，利用折射初至時間公式計算m個點中的每個點的檢波點延遲時、炮點延遲時和高速層速度，以根據m個點中的每個點的風化層速度值、m個點中的每個點的高速層速度、m個點中的每個點的炮點延遲時或者檢波點延遲時，計算m個點中的每個點的風化層厚度值，其中，步驟(d)還包括：根據確定的每個點的風化層厚度值，計算選擇的m個點中的任意兩點之間的風化層厚度值的變差函數值；步驟(e)還包括：從所述目標區域的微測井資料中提取n個點中的每個點的微測井資料的厚度值，并根據提取的每個點的微測井資料的厚度值，計算選擇的n個點中的任意兩點之間的微測井資料的厚度值的變差函數值；步驟(f)還包括：根據確定的m個點中的每個點的風化層厚度值和n個點中的每個點的微測井資料的厚度值，計算m個點中的任意一點與n個點中的任意一點之間的關于厚度的交叉編譯函數；步驟(g)還包括：根據步驟(d)～步驟(f)的結果獲得與m個點中的每個點的風化層厚度值對應的m個加權系數和與n個點中的每個點的微測井數據的厚度值對應的n個加權系數；步驟(h)還包括：根據獲得的與m個點中的每個點的風化層厚度值對應的m個加權系數、與n個點中的每個點的微測井數據的厚度值對應的n個加權系數、確定的m個點中的每個點的風化層厚度值和n個點中的每個點的微測井資料的厚度值，計算風化層厚度的估計值，其中，步驟(q)包括：(q1)根據m個點中的每個點的高速層速度、炮點延遲時、初至數據中每個點對應的炮檢對從激發到接收產生的時間、每個點對應的炮點與檢波點之間的距離，利用折射初至時間公式計算m個點中的每個點的檢波點延遲時，其中，將m個點中的每個點的高速層速度的初始值和m個點中的每個點的炮點延遲時的初始值作為初值代入；(q2)根據m個點中的每個點的高速層速度、步驟(q1)得到的m個點中的每個點的檢波點延遲時、初至數據中每個點對應的炮檢對從激發到接收產生的時間、每個點對應的炮點與檢波點之間的距離，利用折射初至時間公式計算m個點中的每個點的炮點延遲時；(q3)根據步驟(q1)得到的m個點中的每個點的檢波點延遲時、步驟(q2)得到的m個點中的每個點的炮點延遲時、初至數據中每個點對應的炮檢對從激發到接收產生的時間、每個點對應的炮點與檢波點之間的距離，利用折射初至時間公式計算m個點中的每個點的高速層速度；(q4)根據步驟(q1)得到的m個點中的每個點的檢波點延遲時、步驟(q2)得到的m個點中的每個點的炮點延遲時、步驟(q3)得到的m個點中的每個點的高速層速度、每個點對應的炮點與檢波點之間的距離，利用折射初至時間公式計算對應的初至數據中每個點對應的炮檢對從激發到接收產生的...

【專利技術屬性】
技術研發人員：張亨，陳愛萍，何光明，劉奇琳，曹中林，周強，漆韜，李若禹，
申請(專利權)人：中國石油集團川慶鉆探工程有限公司地球物理勘探公司，
類型：發明
國別省市：四川;51