本發(fā)明涉及裂縫參數(shù)預(yù)測(cè)的領(lǐng)域，尤其涉及一種基于橫波速度變化率的裂縫檢測(cè)方法、存儲(chǔ)介質(zhì)及設(shè)備。
背景技術(shù)：
1、非常規(guī)(包括頁(yè)巖氣、致密砂巖氣等)天然氣勘探開發(fā)中還有一系列難題亟需解決，其中的裂縫地震預(yù)測(cè)技術(shù)就是其中一項(xiàng)極其關(guān)鍵的技術(shù)，比如在塔河的奧陶系碳酸鹽巖油藏儲(chǔ)集體以巖溶縫洞為主，對(duì)該地區(qū)進(jìn)行裂縫檢測(cè)或預(yù)測(cè)的需求也比較大。
2、裂縫地震預(yù)測(cè)主要是基于研究寬方位地震數(shù)據(jù)的方位各向異性，即彈性參數(shù)在不同方向存在差異，通過(guò)橢圓擬合技術(shù)找出彈性參數(shù)的快方向和慢方向(例如快橫波、慢橫波方向和大小)，從而確定彈性參數(shù)各向異性主方向(快方向即橢圓長(zhǎng)軸方向)和各向異性強(qiáng)度(橢圓長(zhǎng)軸和短軸之比)。基于方位各向異性反演需要利用寬方位地震資料，但目前寬方位地震資料并不普遍，在沒(méi)有寬方位地震資料的地區(qū)，如何提高裂縫檢測(cè)的精度和可靠性就是一個(gè)亟需解決的難題。以往基于地震屬性如相干、構(gòu)造曲率和螞蟻體等方法主要檢測(cè)較大尺度的小斷層，并不是裂縫或裂隙，裂縫檢測(cè)精度低。
3、專利號(hào)為cn104678434b的專利公開了一種預(yù)測(cè)儲(chǔ)層裂縫發(fā)育參數(shù)的方法，利用不同方位角的地震數(shù)據(jù)提取目的層雙程走時(shí)值，進(jìn)行歸一化處理，得到經(jīng)歸一化處理后的各個(gè)目的層的雙程走時(shí)值；然后對(duì)設(shè)定網(wǎng)格上各點(diǎn)分別生成地震玫瑰圖和綜合裂縫發(fā)育強(qiáng)度平面圖，將兩圖疊合后進(jìn)行顯示，獲得裂縫參數(shù)平面圖。
4、專利號(hào)為cn108562936a的專利公開了一種裂縫預(yù)測(cè)方法、系統(tǒng)、存儲(chǔ)介質(zhì)及終端，根據(jù)ovt道集和地震速度數(shù)據(jù)得到n個(gè)方位入射角道集；再?gòu)念A(yù)設(shè)的地震屬性集中選取任一地震屬性，再根據(jù)地震屬性分別對(duì)n個(gè)方位入射角道集進(jìn)行計(jì)算,得到地震屬性在n個(gè)方位入射角道集對(duì)應(yīng)的方位的強(qiáng)度；對(duì)n個(gè)方位的屬性強(qiáng)度進(jìn)行各向異性計(jì)算,得到計(jì)算結(jié)果；根據(jù)計(jì)算結(jié)果對(duì)待預(yù)測(cè)區(qū)域的裂縫進(jìn)行處理，得到裂縫的發(fā)育密度和優(yōu)勢(shì)方位方向。
5、上述公開的兩個(gè)方法在應(yīng)用時(shí)均需要寬方位地震數(shù)據(jù)，在沒(méi)有寬方位地震資料的地區(qū)進(jìn)行檢測(cè)時(shí)，裂縫檢測(cè)的精度和可靠性差。
技術(shù)實(shí)現(xiàn)思路
1、本發(fā)明解決在沒(méi)有寬方位地震資料的地區(qū)，裂縫地震預(yù)測(cè)的精度和可靠性差的問(wèn)題，提供一種基于橫波速度變化率的裂縫檢測(cè)方法、存儲(chǔ)介質(zhì)及設(shè)備，以解決該技術(shù)問(wèn)題，利用avo三參數(shù)反演得到橫波速度反射率，并通過(guò)尋找δvs/vs的強(qiáng)振幅區(qū)或其道積分?jǐn)?shù)據(jù)低值區(qū)異常，準(zhǔn)確高效地檢測(cè)裂縫。
2、為解決上述技術(shù)問(wèn)題，本發(fā)明的技術(shù)方案如下：
3、一種基于橫波速度變化率的裂縫檢測(cè)方法，包括以下步驟：
4、s1.根據(jù)疊前crp道集數(shù)據(jù)和地震速度，以zoeppritz?aki&
5、rechards線性方程開展三參數(shù)avo反演，得到待檢測(cè)區(qū)域橫波速度變化率δvs/vs數(shù)據(jù)；
6、s2.對(duì)待檢測(cè)區(qū)域奧陶系地層的橫波速度變化率δvs/vs數(shù)據(jù)進(jìn)行處理，得到：
7、橫波速度變化率δvs/vs的最大振幅圖、rms振幅圖，和/或
8、對(duì)所述橫波速度變化率δvs/vs進(jìn)行道積分處理，得到橫波速度vs的相對(duì)值rvs，并處理得到相對(duì)值rvs最小值的平面圖；
9、s3.在所述橫波速度變化率δvs/vs的最大振幅圖、rms振幅圖上圈定強(qiáng)振幅異常區(qū)，和/或
10、在相對(duì)值rvs最小值的平面圖上圈定低值異常區(qū)；
11、s4.所述異常區(qū)為裂縫發(fā)育區(qū)或縫洞發(fā)育區(qū)。
12、優(yōu)選的，在步驟s1中，計(jì)算橫波速度變化率的公式為：
13、
14、式中，r為反射系數(shù)；θ為入射角，單位：度；vp為縱波速度，單位：m/s；vs為橫波速度,單位：m/s；δvs/vs為橫波速度變化率，單位：無(wú)；ρ為密度，單位：g/cc。
15、優(yōu)選的，在步驟s3中，預(yù)先根據(jù)裂縫發(fā)育井段標(biāo)定橫波速度變化率δvs/vs數(shù)據(jù)的強(qiáng)振幅范圍以及橫波速度相對(duì)值rvs低值的范圍，再以標(biāo)定的范圍在所述橫波速度變化率δvs/vs的最大振幅圖、rms振幅圖以及相對(duì)值rvs最小值的平面圖上標(biāo)定異常區(qū)。
16、優(yōu)選的，在步驟s3中，同時(shí)圈定所述強(qiáng)振幅異常區(qū)與低值異常區(qū)的條件下，將所述強(qiáng)振幅異常區(qū)與所述低值異常區(qū)相互驗(yàn)證，以提高對(duì)裂縫發(fā)育區(qū)/縫洞發(fā)育區(qū)的判斷精度，可通過(guò)實(shí)鉆井進(jìn)行驗(yàn)證。
17、優(yōu)選的，所述裂縫檢測(cè)方法的適用條件包括：地震采集inline方向與待檢測(cè)區(qū)域的主斷裂走向垂直。
18、一種計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)，所述計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)中存儲(chǔ)有計(jì)算機(jī)執(zhí)行指令，所述計(jì)算機(jī)執(zhí)行指令被處理器執(zhí)行時(shí)用于實(shí)現(xiàn)上述的基于橫波速度變化率的裂縫檢測(cè)方法。
19、一種計(jì)算機(jī)設(shè)備，所述計(jì)算機(jī)設(shè)備包括處理器和存儲(chǔ)器，所述存儲(chǔ)器存儲(chǔ)有計(jì)算機(jī)程序，所述計(jì)算機(jī)程序由所述處理器加載并執(zhí)行以實(shí)現(xiàn)上述的基于橫波速度變化率的裂縫檢測(cè)方法。
20、本發(fā)明技術(shù)方案的有益技術(shù)效果：
21、本方案利用avo三參數(shù)反演得到橫波速度反射率，再通過(guò)尋找δvs/vs的強(qiáng)振幅區(qū)或其道積分?jǐn)?shù)據(jù)低值區(qū)異常來(lái)實(shí)現(xiàn)對(duì)裂縫進(jìn)行直接有效檢測(cè)，由于該方法是建立在垂直裂縫面方向傳播的橫波速度衰減明顯大于縱波速度衰減理論基礎(chǔ)之上，因此裂縫檢測(cè)精度較疊后地震預(yù)測(cè)方法高，尤其適用于在缺少寬方位地震資料的地區(qū)檢測(cè)裂縫。
22、同時(shí)本方案也解決了目前在沒(méi)有寬方位地震資采集的條件下進(jìn)行裂縫高精度有效檢測(cè)的難題。使用本發(fā)明的技術(shù)方法將顯著提高常規(guī)地震屬性預(yù)測(cè)裂縫的可靠性和精度，較常規(guī)地震屬性裂縫檢測(cè)方法更具優(yōu)勢(shì)。
1.一種基于橫波速度變化率的裂縫檢測(cè)方法，其特征在于，包括以下步驟：
2.如權(quán)利要求1所述的一種基于橫波速度變化率的裂縫檢測(cè)方法，其特征在于，在步驟s1中，計(jì)算橫波速度變化率的公式為：
3.如權(quán)利要求1所述的一種基于橫波速度變化率的裂縫檢測(cè)方法，其特征在于，在步驟s3中，預(yù)先根據(jù)裂縫發(fā)育井段標(biāo)定橫波速度變化率δvs/vs數(shù)據(jù)的強(qiáng)振幅范圍以及橫波速度相對(duì)值rvs低值的范圍，再以標(biāo)定的范圍在所述橫波速度變化率δvs/vs的最大振幅圖、rms振幅圖以及相對(duì)值rvs最小值的平面圖上標(biāo)定異常區(qū)。
4.如權(quán)利要求1所述的一種基于橫波速度變化率的裂縫檢測(cè)方法，其特征在于，在步驟s3中，同時(shí)圈定所述強(qiáng)振幅異常區(qū)與低值異常區(qū)的條件下，將所述強(qiáng)振幅異常區(qū)與所述低值異常區(qū)相互驗(yàn)證，以提高對(duì)裂縫發(fā)育區(qū)/縫洞發(fā)育區(qū)的判斷精度，可通過(guò)實(shí)鉆井進(jìn)行驗(yàn)證。
5.如權(quán)利要求1所述的一種基于橫波速度變化率的裂縫檢測(cè)方法，其特征在于，所述裂縫檢測(cè)方法的適用條件包括：地震采集inline方向與待檢測(cè)區(qū)域的主斷裂走向垂直。
6.一種計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)，其特征在于，所述計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)中存儲(chǔ)有計(jì)算機(jī)執(zhí)行指令，所述計(jì)算機(jī)執(zhí)行指令被處理器執(zhí)行時(shí)用于實(shí)現(xiàn)如權(quán)利要求1至5任一項(xiàng)所述的基于橫波速度變化率的裂縫檢測(cè)方法。
7.一種計(jì)算機(jī)設(shè)備，其特征在于，所述計(jì)算機(jī)設(shè)備包括處理器和存儲(chǔ)器，所述存儲(chǔ)器存儲(chǔ)有計(jì)算機(jī)程序，所述計(jì)算機(jī)程序由所述處理器加載并執(zhí)行以實(shí)現(xiàn)如權(quán)利要求1至5任一所述的基于橫波速度變化率的裂縫檢測(cè)方法。