[發明專利]基于縱橫波速度的脆性礦物指數預測方法

有效

• 申請號: 201710447000.X
• 申請日: 2017-06-14
• 公開（公告）號: CN107203006B
• 公開（公告）日: 2019-01-29
• 發明（設計）人: 鮑祥生;朱立華;周海燕
• 申請（專利權）人: 廣東石油化工學院
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.一種基于縱橫波速度的脆性礦物指數預測方法，其特征在于，包括以下步驟：

寫出三項表示的縱波速度計算方程，將計算孔隙中飽和流體的巖石縱波速度的兩項式Wyllie方程改變為三項表示的方程，所述三項表示方程包括三個部分對巖石縱波速度的貢獻，這三分部為：流體貢獻部分、粘土礦物貢獻部分、脆性礦物貢獻部分；

寫出兩項表示的橫波速度計算方程，將三項表示的縱波速度計算方程的縱波速度換成對應橫波速度，同時將流體貢獻部分給予去除，所述橫波速度計算方程的兩項式Wyllie方程由兩部組成，包括：粘土礦物貢獻部分和脆性礦物貢獻部分；

整理巖石的物理彈性參數，參照巖石物理數據、粘土礦物、脆性礦物的總和定義為1的條件，對三項表示的縱波速度計算方程和兩項表示的橫波速度計算方程進行簡化，得到基于縱橫波速度的脆性礦物指數預測方程，所述巖石的物理彈性參數包括：粘土礦物彈性參數、脆性礦物彈性參數和氣體的彈性參數；

采集應用所需的基礎數據，將采集數據代入基于縱橫波速度的脆性礦物指數預測方程，得到獲得單井或體數據的脆性礦物指數數據，所述基礎數據為：對于單井應用，采集該井的縱橫波速度，對于在有疊前地震數據的地區進行應用，采集基于疊前地震數據進行反演得到縱橫波速度體；

所述寫出三項表示的縱波速度計算方程，將計算孔隙中飽和流體的巖石縱波速度的兩項式Wyllie方程改變為三項表示的方程包括：

兩項式表示的Wyllie方程為：

式中：Φ為巖石的孔隙度，V_fp為流體的縱波速度，V_r為基質的速度，V_p為孔隙中飽和流體巖石的縱波速度；

由于基質是由脆性礦物和粘土礦物組成，因此，兩項式表示的Wyllie方程公式(1)可變形為三項表示Wyllie方程：

式中：Φ為巖石的孔隙度，V_fp為流體的縱波速度，V_m為粘土礦物的體積百分比，V_b為脆性礦物的百分比，即脆性礦物指數，V_mp為地震波在粘土礦物傳播的縱波速度，V_bp為地震波在脆性礦物中傳播的縱波速度；

所述寫出兩項表示的橫波速度計算方程，將三項表示的縱波速度計算方程的縱波速度換成對應橫波速度，同時將流體貢獻部分給予去除為：

由于流體不傳播橫波，所以孔隙中飽和流體巖石的橫波速度可以只用粘土礦物和脆性礦物的橫波貢獻來計算；

式中：V_ms為地震波在粘土礦物中傳播的橫波速度，V_bs為地震波在脆性礦物中傳播的橫波速度，V_s為孔隙中飽含流體巖石的橫波速度。

2.根據權利要求1所述的方法，其特征在于，所述整理巖石的物理彈性參數，參照巖石物理數據、粘土礦物、脆性礦物的總和定義為1的條件，對三項表示的縱波速度計算方程和兩項表示的橫波速度計算方程進行簡化，得到基于縱橫波速度的脆性礦物指數預測方程為：

式中：V_p為孔隙中飽含流體巖石的縱波速度，單位為m/s，V_s為孔隙中飽含流體巖石的橫波速度，單位為m/s，a、b是以單位為m/s的常數，其大小與各個地區獲得巖石物理數據有關。

3.根據權利要求2所述的方法，其特征在于，所述公式(4)中脆性礦物指數V_b的取值范圍為0≤V_b≤1。

4.一種基于縱橫波速度的脆性礦物指數預測方法在基于形成的脆性礦物指數公式和獲得研究區的單井或體數據的脆性礦物指數數據的應用，獲得研究區的單井或體數據的脆性礦物指數數據。