[發明專利]一種基于混沌屬性的砂礫巖體空間展布預測方法及裝置

有效

• 申請號: 201710113208.8
• 申請日: 2017-02-28
• 公開（公告）號: CN106896405B
• 公開（公告）日: 2019-06-21
• 發明（設計）人: 張永華;陳祥;杜偉;羅家群;朱景修;嚴移勝;楊云飛;鄭凱文;劉俊青;朱軍;邱雯;張帆
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司河南油田分公司勘探開發研究院
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 鄭州睿信知識產權代理有限公司 41119
• 代理人: 陳浩

1.一種基于混沌屬性的砂礫巖體空間展布預測方法，其特征在于，該預測方法包括以下步驟：

1)將三維地震數據體中的地震反射振幅作為三維空間變量，計算任意一點P(x、y、t)地震反射振幅A(x、y、t)的梯度；任意一點P(x、y、t)地震反射振幅A(x、y、t)的梯度為：

▽A＝(▽x,▽y,▽t)^T

其中A表示三維地震數據、▽表示求梯度，上標T表示矩陣轉置；

2)在給定范圍內，計算出每個點的梯度向量，并建立相應的協方差矩陣；在定給范圍內各點的梯度向量為：

其中N為時窗長度的樣點數；

3)計算協方差矩陣的特征值，根據協方差矩陣的特征值確定混沌屬性值J，根據混沌屬性值J確定砂礫巖體的大小，以實現對砂礫巖體空間展布的預測；根據混沌屬性值J確定砂礫巖體的大小的步驟為：

A)在相應的地質層位上、統計研究區內鉆遇砂礫巖與未鉆遇砂礫巖的井段，利用合成記錄對目的層進行標定，確定砂礫巖體與非砂礫巖體在混沌屬性體上對應的屬性值；B)通過多口井的統計，將統計井的混沌屬性值利用直方圖進行分析，對于未能完全區分的混沌屬性值范圍，即重疊區，根據重疊區間求取二者的均值作為門檻值；C)根據求取的門檻值大小及混沌屬性含義，大于門檻值的混沌屬性部分為預測的砂礫巖體。

2.根據權利要求1所述的基于混沌屬性的砂礫巖體空間展布預測方法，其特征在于，該方法還包括利用合成記錄對目的層進行標定，對屬性體資料進行解釋，通過巖層屬性提取對砂礫巖體進行描述。

3.根據權利要求1或2所述的基于混沌屬性的砂礫巖體空間展布預測方法，其特征在于，混沌屬性值J的計算公式為：

其中λ_max為協方差矩陣的最大特征值，λ_min為協方差矩陣的最小特征值，λ_mid為協方差矩陣的中間特征值。

4.根據權利要求1所述的基于混沌屬性的砂礫巖體空間展布預測方法，其特征在于，橫向上J值較大的區域反映砂體的平面展布，縱向上J值較小的區間反映砂礫巖體形成的間歇期。

5.根據權利要求1所述的基于混沌屬性的砂礫巖體空間展布預測方法，其特征在于，在利用三維地震數據計算地震反射振幅梯度時需對三維地震數據進行疊后去噪處理。

6.一種基于混沌屬性的砂礫巖體空間展布預測裝置，其特征在于，該裝置包括梯度計算模塊、協方差矩陣構建模塊和混沌屬性值確定模塊，

所述的梯度計算模塊用于將三維地震數據體中的地震反射振幅作為三維空間變量，計算任意一點P(x、y、t)地震反射振幅A(x、y、t)的梯度；任意一點P(x、y、t)地震反射振幅A(x、y、t)的梯度為：

▽A＝(▽x,▽y,▽t)^T

其中A表示三維地震數據、▽表示求梯度，上標T表示矩陣轉置；

所述的協方差矩陣構建模塊用于在給定范圍內，計算出每個點的梯度向量，并建立相應的協方差矩陣；在定給范圍內各點的梯度向量為：

其中N為時窗長度的樣點數；

所述的混沌屬性值確定模塊用于計算協方差矩陣的特征值，根據協方差矩陣的特征值確定混沌屬性值J，并根據混沌屬性值J確定砂礫巖體的大小；根據混沌屬性值J確定砂礫巖體的大小的步驟為：

A)在相應的地質層位上、統計研究區內鉆遇砂礫巖與未鉆遇砂礫巖的井段，利用合成記錄對目的層進行標定，確定砂礫巖體與非砂礫巖體在混沌屬性體上對應的屬性值；B)通過多口井的統計，將統計井的混沌屬性值利用直方圖進行分析，對于未能完全區分的混沌屬性值范圍，即重疊區，根據重疊區間求取二者的均值作為門檻值；C)根據求取的門檻值大小及混沌屬性含義，大于門檻值的混沌屬性部分為預測的砂礫巖體。

7.根據權利要求6所述的基于混沌屬性的砂礫巖體空間展布預測裝置，其特征在于，該裝置還包括用于利用合成記錄對目的層進行標定，對屬性體資料進行解釋，通過巖層屬性提取對砂礫巖體進行描述的模塊。

8.根據權利要求6或7所述的基于混沌屬性的砂礫巖體空間展布預測裝置，其特征在于，所述混沌屬性值確定模塊采用的混沌屬性值計算公式為：

其中λ_max為協方差矩陣的最大特征值，λ_min為協方差矩陣的最小特征值，λ_mid為協方差矩陣的中間特征值。