本發(fā)明涉及油田開發(fā)技術(shù)領(lǐng)域,特別是涉及到一種自動修正構(gòu)造層面的方法。
背景技術(shù):
隨著油藏描述的發(fā)展,三維地質(zhì)建模技術(shù)日益普及,在地質(zhì)研究和數(shù)值模擬中的應(yīng)用也越來越廣泛。地質(zhì)建模是綜合應(yīng)用各種資料、多學(xué)科技術(shù)的過程,構(gòu)造層面是地質(zhì)建模必備的內(nèi)容,也是一切模型的基礎(chǔ)。然而受不同專業(yè)、不同軟件和研究者水平的差異,從工作站或者地震數(shù)據(jù)體中導(dǎo)出的構(gòu)造層面往往存在奇異點(diǎn),有些是工作的疏忽、有些是軟件本身的限制、有些可能是目前手段無法處理的。對這些奇異點(diǎn)的處理非常困難,而且對于兩個層面之間的修正更加復(fù)雜,一直困擾科研人員。為此我們發(fā)明了一種新的自動修正構(gòu)造層面的方法,解決了以上技術(shù)問題。
技術(shù)實(shí)現(xiàn)要素:
本發(fā)明的目的是提供一種簡單、易于操作并且非常實(shí)用的自動修正構(gòu)造層面的方法。
本發(fā)明的目的可通過如下技術(shù)措施來實(shí)現(xiàn):自動修正構(gòu)造層面的方法,該自動修正構(gòu)造層面的方法包括:步驟1,通過精細(xì)構(gòu)造解釋確定構(gòu)造層面,導(dǎo)入構(gòu)造層面;步驟2,對奇異點(diǎn)進(jìn)行識別;步驟3,根據(jù)閥值對奇異點(diǎn)進(jìn)行修正,得到新的構(gòu)造層面;步驟4,進(jìn)行相鄰層面對比,將相鄰層面相減,得到厚度面;步驟5,進(jìn)行厚度差異修正;步驟6,輸出修正好的構(gòu)造層面。
本發(fā)明的目的還可通過如下技術(shù)措施來實(shí)現(xiàn):
在步驟1中,通過地震資料,利用合成記錄及反演技術(shù),確定區(qū)域構(gòu)造特征及構(gòu)造層面,或根據(jù)井點(diǎn)插值得到的構(gòu)造層面。
在步驟2中,通過巖心、錄井、測井資料,確定閥值Φ;在構(gòu)造層面上任取一點(diǎn)A,給定一半徑R的圓,計算點(diǎn)A與圓內(nèi)任意一點(diǎn)Bn之間的高度差ΔZn,若ΔZn≤Φ,則A不是奇異點(diǎn);若ΔZn>Φ,則A是奇異點(diǎn)。
在步驟3中,令ΔZn=Φ,對奇異點(diǎn)進(jìn)行校正,則點(diǎn)A的高度ZA等于點(diǎn)Bn的 高度與ΔZn值之和,即ZA=ZBn+Φ。
步驟3還包括,當(dāng)?shù)玫叫碌臉?gòu)造層面有斷層發(fā)育時,給定斷層邊界之后,依據(jù)斷層的屬性和斷層斷距這些參數(shù)整體調(diào)整層面,將上盤和下盤按照斷距斷開。
在步驟4中,為了確保每個構(gòu)造層面都沒有異常點(diǎn),對任意上下相鄰的層面Cn和Cn+1,分別取構(gòu)造層面任意相同橫坐標(biāo)xi和縱坐標(biāo)yi的點(diǎn),將這些點(diǎn)構(gòu)造深度Zi相減,得到厚度值ΔHi;若ΔHi小于0,說明構(gòu)造層面Cn和Cn+1存在異常點(diǎn)和則重復(fù)步驟2和步驟3,反之,流程進(jìn)入到步驟5。
在步驟5中,通過井間地層對比得到各個相鄰構(gòu)造層面間的理論下限值,將該下限值確定為當(dāng)則保留原值;當(dāng)則將該厚度賦值為那么構(gòu)造層面Cn+1的點(diǎn)的Z值等于構(gòu)造層面Cn的點(diǎn)的Z值與之和,即
本發(fā)明中的自動修正構(gòu)造層面的方法,包括:構(gòu)造層面導(dǎo)入、奇異點(diǎn)識別、奇異點(diǎn)修正、邊界或斷層條件約束、相鄰層面對比、厚度差異修正、輸出層面7個步驟。通過精細(xì)構(gòu)造解釋確定構(gòu)造層面,導(dǎo)入;運(yùn)用設(shè)計好的算法對奇異點(diǎn)進(jìn)行識別;根據(jù)閥值對奇異點(diǎn)進(jìn)行修正;相鄰層面相減,得到厚度面;對負(fù)值區(qū)域進(jìn)行校正;導(dǎo)出修正好的構(gòu)造層面。該方法簡單、易于操作并且非常實(shí)用,能夠快速識別構(gòu)造層面存在的奇異點(diǎn)并通過電腦編程實(shí)現(xiàn)自動校正,可以有效地提高工作效率。
附圖說明
圖1為本發(fā)明的自動修正構(gòu)造層面的方法的一具體實(shí)施例的流程圖;
圖2為本發(fā)明的一具體實(shí)施例中任意點(diǎn)A同相鄰數(shù)據(jù)位置示意圖;
圖3為本發(fā)明的一具體實(shí)施例中差值ΔZ求取示意圖。
具體實(shí)施方式
為使本發(fā)明的上述和其他目的、特征和優(yōu)點(diǎn)能更明顯易懂,下文特舉出較佳實(shí)施例,并配合附圖所示,作詳細(xì)說明如下。
如圖1所示,圖1為本發(fā)明的自動修正構(gòu)造層面的方法的流程圖。
步驟101,構(gòu)造層面導(dǎo)入。通過地震資料,利用合成記錄及反演技術(shù), 確定區(qū)域構(gòu)造特征及構(gòu)造層面,也可以是根據(jù)井點(diǎn)插值得到的構(gòu)造層面;將構(gòu)造層面的數(shù)據(jù)按(x,y,z)格式導(dǎo)入。流程進(jìn)入到步驟102。
步驟102,奇異點(diǎn)識別。通過巖心、錄井、測井資料,確定閥值Φ;在構(gòu)造層面上任取一點(diǎn)A,給定一半徑R的圓,計算點(diǎn)A與圓內(nèi)任意一點(diǎn)Bn之間的高度差ΔZn,見圖2、圖3,圖2中點(diǎn)A為構(gòu)造層面任意一點(diǎn),點(diǎn)B1、B2...Bn為以點(diǎn)A為中心,半徑為R的圓范圍內(nèi)任意數(shù)據(jù)點(diǎn);圖3中ΔZ=ZA-ZBn。若ΔZn≤Φ,則A不是奇異點(diǎn);若ΔZn>Φ,則A是奇異點(diǎn)。流程進(jìn)入到步驟103。
步驟103,奇異點(diǎn)修正。令ΔZn=Φ,對奇異點(diǎn)進(jìn)行校正,則點(diǎn)A的高度等于點(diǎn)Bn的高度與ΔZn值之和,即ZA=ZBn+Φ。流程進(jìn)入到步驟104。
步驟104,邊界或斷層條件約束。通過對構(gòu)造層面的奇異點(diǎn)進(jìn)行校正后,得到新的構(gòu)造層面,分別命名為C1、C2...Cn,如果構(gòu)造層面沒有斷層發(fā)育,則直接進(jìn)入步驟105;如果構(gòu)造層面有斷層發(fā)育,給定斷層邊界之后,可依據(jù)斷層的屬性(正斷層或逆斷層)、斷層斷距等參數(shù)整體調(diào)整層面,將上盤和下盤按照斷距斷開。流程進(jìn)入到步驟105。
步驟105,相鄰層面對比。為了確保每個構(gòu)造層面都沒有異常點(diǎn),對任意上下相鄰的層面Cn和Cn+1,分別取構(gòu)造層面任意相同橫坐標(biāo)xi和縱坐標(biāo)yi的點(diǎn),將這些點(diǎn)構(gòu)造深度Zi相減,得到厚度值ΔHi;若ΔHi小于0,說明構(gòu)造層面Cn和Cn+1存在異常點(diǎn)和則重復(fù)步驟102、103、104。反之,流程進(jìn)入到步驟106。
步驟106,厚度差異修正。通過井間地層對比得到各個相鄰構(gòu)造層面間的理論下限值,將該下限值確定為當(dāng)則保留原值;當(dāng)則將該厚度賦值為那么構(gòu)造層面Cn+1的點(diǎn)的Z值等于構(gòu)造層面Cn的點(diǎn)的Z值與之和,即流程進(jìn)入到步驟107。
步驟107,輸出層面。根據(jù)以上步驟,得到了每個層面數(shù)據(jù)校正后的結(jié)果(x,y,z校正),將最終修正好的各個構(gòu)造層面按照相應(yīng)的格式導(dǎo)出。流程結(jié)束。