本發(fā)明屬于火山巖儲層含氣飽和度和滲透率預測，尤其涉及一種火山巖儲層含氣飽和度與滲透率地震反演方法及系統。
背景技術：
1、隨著非常規(guī)油氣資源開發(fā)的不斷推進，火山巖儲層因其資源潛力大而日益受到關注。然而，火山巖儲層成藏機制特殊，巖性成分多樣、構造復雜、成巖演化劇烈，形成了低孔隙度、低滲透率、流體非均勻分布等儲層特征；同時，火山巖儲層在空間形態(tài)上表現出橫向和縱向快速變化的特征，且其空間展布具有較強的離散性和不連續(xù)性；在微觀尺度上，孔隙、溶蝕孔洞、微裂縫及節(jié)理等多種孔滲結構相互交織，形成了復雜的儲集空間與滲流通道?；鹕綆r儲層中，含氣飽和度是儲層含氣富集程度的重要指標，滲透率控制致密火山巖地層中的流體運移，二者是影響儲層開發(fā)潛力的重要評價參數，因此，應用地震方法預測含氣飽和度和滲透率的空間分布度對儲層識別與優(yōu)化開發(fā)具有重要意義。
2、根據巖石物理理論，巖石孔隙中的流體能夠顯著影響地震波的傳播特性。傳統gassmann理論可用于描述單一流體飽和巖石的彈性響應，但對多種不相溶流體的復雜情況適用性有限，為克服該局限，現有技術中將gassmann理論與等效介質方法結合，提升流體響應的模擬精度，在低頻地震條件下，當流體壓力平衡，流體體積模量與wood理論一致；在無壓力平衡下，voigt理論變得適用；介于二者之間的狀態(tài)通常采用brie指數模型近似；papageorgiou?et?al.進一步引入靜態(tài)毛細管壓力參數來描述孔隙空間內不相溶流體之間的壓力差，計算混合流體的等效體積模量。盡管上述研究在流體響應模擬方法上取得了重要進展，但目前仍缺乏對與流體混合狀態(tài)有關的地層滲透性彈性及地震響應的有效描述方法。因此，亟需開展火山巖儲層巖石物理建模研究，為流體和滲透率巖石物理響應模擬提供有效方法。
3、地震巖石物理反演作為連接地震數據與巖石物理參數的重要手段，已成為通過地震數據提取儲層參數的重要工具，其核心在于通過巖石物理模型，應用測井與地震數據實現儲層參數的定量預測。目前，常用的反演策略主要包括兩類：一是先由疊前地震數據反演彈性參數，在此基礎上應用巖石物理關系開展儲層巖石物理參數的定量解釋；二是通過巖石物理約束構建關于儲層物性參數的反射系數公式，并開發(fā)相應的物性參數的直接反演方法。直接反演方法可減少間接解釋過程在誤差積累的影響，提升預測結果的可靠性，但是目前針對非均質性強、孔滲特征復雜的火山巖儲層，仍缺乏專門針對含氣飽和度與滲透率參數的地震直接反演方法。因此，亟需建立一種以巖石物理為驅動的疊前地震反演方法，用于火山巖儲層含氣飽和度與滲透率的直接反演與定量預測。
技術實現思路
1、本發(fā)明實施例的目的在于提供一種火山巖儲層含氣飽和度與滲透率地震反演方法，旨在解決上述背景技術中提出的問題。
2、本發(fā)明實施例是這樣實現的，一種火山巖儲層含氣飽和度與滲透率地震反演方法，包括以下步驟：
3、采集數據，建立火山巖儲層的巖石物理模型；
4、建立基于模型驅動的巖石物理參數井中反演方法，計算滲透率指示參數和微裂縫孔隙度參數；
5、基于參數和，結合測井數據構建含氣性地震指示因子和滲透率地震指示因子；通過井中實測數據建立與含氣飽和度、與滲透率之間的定量關系；
6、結合巖石物理與地震反射理論，構建關于和的pp波反射系數公式；
7、基于構建的pp波反射系數公式開發(fā)適用于火山巖儲層的疊前地震反演方法，計算地震指示因子和，量化為含氣飽和度和滲透率。
8、本發(fā)明實施例提供的一種火山巖儲層含氣飽和度與滲透率地震反演方法，在現有技術中的基礎上，建立了適用于火山巖儲層的巖石物理模型，引入滲透率指示參數，有效實現流體混合狀態(tài)與孔隙連通性等因素影響的滲透率特征的定量描述，同時，提出了含氣性地震指示因子和滲透率地震指示因子，并建立以和參數化表達的pp波反射系數公式；在此基礎上，開發(fā)適用于火山巖復雜地層結構的疊前地震反演方法，由疊前地震數據計算指示因子和；最終，綜合井中巖石物理分析和地震反演結果，實現火山巖儲層含氣飽和度與滲透率的空間分布的定量預測，為火山巖儲層含氣飽和度與滲透率的定量預測提供了高精度反演方法與新的技術思路。
1.一種火山巖儲層含氣飽和度與滲透率地震反演方法，其特征在于，包括以下步驟：
2.根據權利要求1所述的火山巖儲層含氣飽和度與滲透率地震反演方法，其特征在于，所述采集數據，建立火山巖儲層的巖石物理模型的步驟中，火山巖儲層的巖石物理模型的公式為：
3.根據權利要求2所述的火山巖儲層含氣飽和度與滲透率地震反演方法，其特征在于，所述建立基于模型驅動的巖石物理參數井中反演方法，計算滲透率指示參數和微裂縫孔隙度參數的步驟中，反演的目標函數為：
4.根據權利要求3所述的火山巖儲層含氣飽和度與滲透率地震反演方法，其特征在于，所述基于參數和，結合測井數據構建含氣性地震指示因子和滲透率地震指示因子的步驟中，所述含氣性地震指示因子的計算公式為：
5.根據權利要求4所述的火山巖儲層含氣飽和度與滲透率地震反演方法，其特征在于，所述結合巖石物理與地震反射理論，構建關于和的pp波反射系數公式的步驟中，所述pp波反射系數公式為：
6.根據權利要求1所述的火山巖儲層含氣飽和度與滲透率地震反演方法，其特征在于，所述基于構建的pp波反射系數公式開發(fā)適用于火山巖儲層的疊前地震反演方法，計算地震指示因子和的步驟，具體包括：
7.一種火山巖儲層含氣飽和度與滲透率地震反演系統，用于實現如權利要求1-6任一所述的火山巖儲層含氣飽和度與滲透率地震反演方法，其特征在于，包括：