本公開涉及地球物理,尤其涉及一種地層模型構(gòu)建方法、裝置和電子設(shè)備及存儲介質(zhì)。
背景技術(shù):
1、克里金插值方法是一種利用空間變量的空間相關(guān)信息對空間已知點賦予一定的權(quán)值,然后通過加權(quán)平均法對空間未知點進行估計的方法。其首先由測井數(shù)據(jù)計算變差函數(shù)得到地下介質(zhì)的空間變程,以表征地下介質(zhì)的空間相關(guān)信息。隨后將空間變程代入空間相關(guān)函數(shù)中得到空間協(xié)方差函數(shù),并由此建立克里金矩陣和克里金向量。最后由克里金矩陣和克里金向量計算得到空間已知點的克里金權(quán)重系數(shù),對空間已知點按其各自的克里金權(quán)重系數(shù)對其進行加權(quán)平均即可得到克里金插值結(jié)果。
2、在克里金插值過程中,為了簡化插值過程,通常假設(shè)空間中所有點在不同方向上的空間相關(guān)信息是相同的,利用空間變程來表征地下介質(zhì)的各向同性空間結(jié)構(gòu)。但是地下介質(zhì)由于沉積環(huán)境、構(gòu)造運動以及上覆壓力等地質(zhì)因素的影響,其空間結(jié)構(gòu)通常表現(xiàn)出各向異性的特征,不同方向上的空間相關(guān)長度是不同的。因此基于空間變程的各向同性空間結(jié)構(gòu)克里金插值方法得到的克里金插值結(jié)果與地下介質(zhì)的真實情況相差較大,難以符合地下介質(zhì)的實際空間結(jié)構(gòu)。
技術(shù)實現(xiàn)思路
1、根據(jù)本公開實施例的一方面,提供了一種地層模型構(gòu)建方法,包括:獲取地層空間的地震數(shù)據(jù)和該地層空間中未知空間點和已知空間點的位置坐標、以及已知空間點的觀測值;根據(jù)該地層空間的地震數(shù)據(jù)估計各個未知空間點和各個已知空間點的隨機介質(zhì)統(tǒng)計特征參數(shù);根據(jù)各個未知空間點和各個已知空間點的位置坐標和隨機介質(zhì)統(tǒng)計特征參數(shù),確定每個已知空間點與其他已知空間點及各個未知空間點連線方向上的各向異性空間相關(guān)長度;將確定得到的各向異性空間相關(guān)長度代入空間相關(guān)函數(shù),以獲得每個已知空間點與其他已知空間點及各個未知空間點的各向異性空間協(xié)方差函數(shù);針對每個未知空間點:根據(jù)確定得到的各向異性空間協(xié)方差函數(shù),確定每個已知空間點的各向異性克里金權(quán)重系數(shù);根據(jù)確定得到的各向異性克里金權(quán)重系數(shù)和各個已知空間點的觀測值,確定該未知空間點的預測值。
2、在一些可能的實施方式中,隨機介質(zhì)統(tǒng)計特征參數(shù)包括:橫向自相關(guān)長度、縱向自相關(guān)長度以及自相關(guān)角度。
3、在一些可能的實施方式中,按照以下公式空間點i與空間點j連線方向上的各向異性空間相關(guān)長度lφij:
4、
5、式中,ai為空間點i的橫向自相關(guān)長度、bi為空間點i的縱向自相關(guān)長度,θi為空間點i的自相關(guān)角度,δxij和δyij分別為空間點i和空間點j之間的橫縱坐標差。
6、在一些可能的實施方式中,根據(jù)地層空間的地震數(shù)據(jù)估計各個未知空間點和各個已知空間點的隨機介質(zhì)統(tǒng)計特征參數(shù),包括:針對每個空間點:以該空間點為滑動窗中心點,計算滑動窗范圍內(nèi)地震數(shù)據(jù)的功率譜和地震子波的功率譜;基于該功率譜法得到隨機介質(zhì)功率譜;基于維納-辛欽定理對隨機介質(zhì)功率譜做傅里葉反變換得到隨機介質(zhì)自相關(guān)函數(shù);將隨機介質(zhì)自相關(guān)函數(shù)二值化后擬合為橢圓;估計所述橢圓的長軸,獲得橫向自相關(guān)長度;估計所述橢圓的短軸,獲得縱向自相關(guān)長度;確定所述長軸與水平軸形成的夾角,獲得自相關(guān)角度。
7、根據(jù)本公開實施例的另一方面,提供了一種地層模型構(gòu)建裝置,包括:獲取模塊,用于獲取地層空間的地震數(shù)據(jù)和該地層空間中未知空間點和已知空間點的位置坐標、以及已知空間點的觀測值;處理模塊,用于根據(jù)該地層空間的地震數(shù)據(jù)估計各個未知空間點和各個已知空間點的隨機介質(zhì)統(tǒng)計特征參數(shù);根據(jù)各個未知空間點和各個已知空間點的位置坐標和隨機介質(zhì)統(tǒng)計特征參數(shù),確定每個已知空間點與其他已知空間點及各個未知空間點連線方向上的各向異性空間相關(guān)長度;將確定得到的各向異性空間相關(guān)長度代入空間相關(guān)函數(shù),以獲得每個已知空間點與其他已知空間點及各個未知空間點的各向異性空間協(xié)方差函數(shù);針對每個未知空間點:根據(jù)確定得到的各向異性空間協(xié)方差函數(shù),確定每個已知空間點的各向異性克里金權(quán)重系數(shù);根據(jù)確定得到的各向異性克里金權(quán)重系數(shù)和各個已知空間點的觀測值,確定該未知空間點的預測值。
8、在一些可能的實施方式中,隨機介質(zhì)統(tǒng)計特征參數(shù)包括:橫向自相關(guān)長度、縱向自相關(guān)長度以及自相關(guān)角度。
9、在一些可能的實施方式中,處理模塊,用于按照以下公式空間點i與空間點j連線方向上的各向異性空間相關(guān)長度lφij:
10、
11、式中,ai為空間點i的橫向自相關(guān)長度、bi為空間點i的縱向自相關(guān)長度,θi為空間點i的自相關(guān)角度,δxij和δyij分別為空間點i和空間點j之間的橫縱坐標差。
12、在一些可能的實施方式中,處理模塊,用于針對每個空間點:以該空間點為滑動窗中心點,計算滑動窗范圍內(nèi)地震數(shù)據(jù)的功率譜和地震子波的功率譜;基于功率譜法得到隨機介質(zhì)功率譜;基于維納-辛欽定理對隨機介質(zhì)功率譜做傅里葉反變換得到隨機介質(zhì)自相關(guān)函數(shù);將隨機介質(zhì)自相關(guān)函數(shù)二值化后擬合為橢圓;估計所述橢圓的長軸,獲得橫向自相關(guān)長度;估計所述橢圓的短軸,獲得縱向自相關(guān)長度;確定所述長軸與水平軸形成的夾角,獲得自相關(guān)角度。
13、根據(jù)本公開實施例的另一方面,提供了一種電子設(shè)備,包括:處理器;以及存儲程序的存儲器,其中,所述程序包括指令,所述指令在由所述處理器執(zhí)行時使所述處理器執(zhí)行前述提及的方法。
14、根據(jù)本公開實施例的另一方面,提供了一種存儲有計算機指令的非瞬時計算機可讀存儲介質(zhì),其中,所述計算機指令用于使計算機執(zhí)行前述提及的方法。
15、本公開實施例中提供的一個或多個技術(shù)方案,根據(jù)地震數(shù)據(jù)估計隨機介質(zhì)統(tǒng)計特征參數(shù),通過隨機介質(zhì)統(tǒng)計特征參數(shù)來表征地下介質(zhì)的各向異性空間結(jié)構(gòu)特征,并由此計算地下介質(zhì)中空間點在不同方向上的空間相關(guān)長度,克服了傳統(tǒng)克里金插值方法假設(shè)地下介質(zhì)空間結(jié)構(gòu)是各向同性的假設(shè),能夠建立更為誤差更小、更符合地下介質(zhì)實際空間結(jié)構(gòu)特征以及插值過程不確定性更小的克里金插值結(jié)果。并且,在已知空間點的數(shù)量適中、較少的情形,仍然具有更好的應用效果。尤其是在已知空間點的數(shù)量較少的情形,該方法對地層模型中薄層的刻畫是出色的,無論是薄層的邊界還是其連續(xù)性都具有很好的效果,薄層細節(jié)信息豐富。
1.一種地層模型構(gòu)建方法,其特征在于,包括:
2.如權(quán)利要求1所述的方法,其特征在于,所述隨機介質(zhì)統(tǒng)計特征參數(shù)包括:橫向自相關(guān)長度、縱向自相關(guān)長度以及自相關(guān)角度。
3.如權(quán)利要求2所述的方法,其特征在于,按照以下公式空間點i與空間點j連線方向上的各向異性空間相關(guān)長度lφij:
4.如權(quán)利要求2所述的方法,其特征在于,所述根據(jù)所述地層空間的地震數(shù)據(jù)估計各個未知空間點和各個已知空間點的隨機介質(zhì)統(tǒng)計特征參數(shù),包括:
5.一種地層模型構(gòu)建裝置,其特征在于,包括:
6.如權(quán)利要求5所述的裝置,其特征在于,所述隨機介質(zhì)統(tǒng)計特征參數(shù)包括:橫向自相關(guān)長度、縱向自相關(guān)長度以及自相關(guān)角度。
7.如權(quán)利要求6所述的裝置,其特征在于,所述處理模塊,用于按照以下公式空間點i與空間點j連線方向上的各向異性空間相關(guān)長度lφij:
8.如權(quán)利要求6所述的裝置,其特征在于,所述處理模塊,用于針對每個空間點:
9.一種電子設(shè)備,包括:
10.一種存儲有計算機指令的非瞬時計算機可讀存儲介質(zhì),其中,所述計算機指令用于使計算機執(zhí)行根據(jù)權(quán)利要求1-4中任一項所述的方法。