基于地質(zhì)模型的全方位觀測(cè)系統(tǒng)參數(shù)論證方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明涉及地震勘探觀測(cè)系統(tǒng)參數(shù)論證，特別是涉及到一種基于地質(zhì)模型的全方 位觀測(cè)系統(tǒng)參數(shù)論證方法。
【背景技術(shù)】
[0002] 以往觀測(cè)系統(tǒng)參數(shù)論證只是在全工區(qū)選擇2-3個(gè)參數(shù)論證點(diǎn)進(jìn)行論證，所使用的 方法也比較簡(jiǎn)單，都是基于各向同性介質(zhì)的假設(shè)下計(jì)算的，而在西部非常復(fù)雜的山前帶區(qū) 域，由于地質(zhì)構(gòu)造已經(jīng)不滿足各向同性介質(zhì)假設(shè)的條件，因此這種方法已經(jīng)不能滿足技術(shù) 設(shè)計(jì)的需要，首先2-3個(gè)論證點(diǎn)不能滿足需要，所計(jì)算出來(lái)的結(jié)果具有偏頗性，另外計(jì)算方 法不合適，計(jì)算結(jié)果不滿足實(shí)際生產(chǎn)需要。為此我們發(fā)明了一種新的基于地質(zhì)模型的全方 位觀測(cè)系統(tǒng)參數(shù)論證方法，解決了以上技術(shù)問(wèn)題。

【發(fā)明內(nèi)容】

[0003] 本發(fā)明的目的是提供一種以地質(zhì)模型為基礎(chǔ)，在模型的基礎(chǔ)上計(jì)算論證每個(gè)物理 點(diǎn)采集參數(shù)需求的基于地質(zhì)模型的全方位觀測(cè)系統(tǒng)參數(shù)論證方法。
[0004] 本發(fā)明的目的可通過(guò)如下技術(shù)措施來(lái)實(shí)現(xiàn)：基于地質(zhì)模型的全方位觀測(cè)系統(tǒng)參數(shù) 論證方法，該基于地質(zhì)模型的全方位觀測(cè)系統(tǒng)參數(shù)論證方法包括：步驟1，采集多種地質(zhì)數(shù) 據(jù)，建立工區(qū)地質(zhì)模型；步驟2,選擇地質(zhì)模型中的一個(gè)層位，采用公式計(jì)算層位上每個(gè)物 理點(diǎn)的面元的需求；步驟3,將每個(gè)物理點(diǎn)的面元需求顯示在模型中對(duì)應(yīng)的層面上；步驟4， 選擇地質(zhì)模型中的一個(gè)層位，采用公式計(jì)算層位上每個(gè)物理點(diǎn)的最大炮檢距的需求；步驟 5,將每個(gè)物理點(diǎn)的最大炮檢距需求顯示在模型中對(duì)應(yīng)的層面上；以及步驟6,確定滿足生 產(chǎn)需要的觀測(cè)系統(tǒng)參數(shù)。
[0005] 本發(fā)明的目的還可通過(guò)如下技術(shù)措施來(lái)實(shí)現(xiàn)：
[0006] 在步驟1中，采集的地質(zhì)數(shù)據(jù)包括層位數(shù)據(jù)、速度資料、密度資料、深度資料、斷層 資料、頻率資料這些資料。
[0007] 在步驟1中，層位數(shù)據(jù)為txt格式，包含各個(gè)層位的每個(gè)物理點(diǎn)的X橫坐標(biāo)、Y縱坐 標(biāo)、Z深度信息，各個(gè)層位中每個(gè)離散物理點(diǎn)的密度小于該區(qū)地震勘探所需成像面元大小， 速度資料、頻率資料、密度資料、斷層資料與層位數(shù)據(jù)是對(duì)應(yīng)的。
[0008] 在步驟1中，采用層位數(shù)據(jù)建立工區(qū)地質(zhì)模型，并將各層位對(duì)應(yīng)的速度資料、密度 資料、斷層資料導(dǎo)入，每個(gè)物理點(diǎn)賦予其地質(zhì)屬性。
[0009]，在步驟2中，在采用公式計(jì)算層位上每個(gè)物理點(diǎn)的面元的需求時(shí)，采用兩種公式 分別進(jìn)行計(jì)算，分別是滿足最高無(wú)混疊頻率法則和滿足橫向分辨率法則公式；
[0010] 滿足最高無(wú)混疊頻率法則的公式為：b=Vint/(4XFmaxXSin0 )
[0011] 式中：b-面兀大小，Vint-上一層層速度，F(xiàn)max-最商無(wú)混疊頻率，Q地層 傾角；
[0012] 滿足橫向分辨率要求，依據(jù)保證良好橫向分辨率面元大小的經(jīng)驗(yàn)公式為：
[0013] b=Vint/(2XFdom)
[0014] 式中：b-面元大小，Vint-上一層層速度，F(xiàn)dom-反射層視主頻，采用這兩種公 式，分別進(jìn)行計(jì)算。
[0015] 在步驟3中，原有地質(zhì)模型不變,將每個(gè)物理點(diǎn)的面元需求疊合顯示在層面上,采 用顏色代表面元需求的大小，根據(jù)步驟2所采用的兩種公式計(jì)算出的兩種面元大小，分別 顯示這兩種面元大小。
[0016] 在步驟4中，在采用公式計(jì)算層位上每個(gè)物理點(diǎn)的最大炮檢距的需求時(shí)，采用兩 種公式分別進(jìn)行計(jì)算，分別是滿足動(dòng)校正拉伸需要法則和滿足速度分析精度法則公式；
[0017] 滿足動(dòng)校正拉伸需要法則的公式為：
[0019]式中：
[0020] Vrms--均方根速度；
[0021]Xnax-最大炮檢距；
[0022] t〇--雙程旅行時(shí)；
[0023]D--拉伸系數(shù)；
[0024] 滿足速度分析精度法則的公式為：
[0026] 式中：Av-待測(cè)速度增量，
[0027]Xmax-最大炮檢距
[0028]At--速度分析中可以檢測(cè)出的最小正常時(shí)差值
[0029]t〇--雙程旅行時(shí)；
[0030] fp反射波主頻；
[0031]vRMS--均方根速度。
[0032] 在步驟5中，原有地質(zhì)模型不變，將每個(gè)物理點(diǎn)的最大炮檢距需求疊合顯示在層 面上，采用顏色代表最大炮檢距需求的大小，根據(jù)步驟4所采用的兩種公式計(jì)算出的兩種 最大炮檢距大小，分別顯示這兩種最大炮檢距大小。
[0033] 在步驟6中，通過(guò)對(duì)比分析各個(gè)層位的參數(shù)需求大小，并結(jié)合施工成本需求，合理 選擇滿足全區(qū)生產(chǎn)需要的觀測(cè)系統(tǒng)參數(shù)。
[0034] 本發(fā)明中的基于地質(zhì)模型的全方位觀測(cè)系統(tǒng)參數(shù)論證方法，通過(guò)收集工區(qū)的以往 地質(zhì)資料，建立精細(xì)的地質(zhì)模型，在該模型的基礎(chǔ)上采用公式計(jì)算每個(gè)物理點(diǎn)觀測(cè)系統(tǒng)參 數(shù)需求，然后將參數(shù)需求疊合顯示在地質(zhì)模型上，這種顯示方式比較直觀的顯示出全區(qū)的 參數(shù)需求，因此更加容易確定出了適合于該層位地震采集需求的參數(shù)，可以較為直觀、全面 的分析全區(qū)的參數(shù)需求，更加科學(xué)合理。該發(fā)明采用地質(zhì)模型來(lái)計(jì)算每個(gè)物理點(diǎn)的參數(shù)需 求，克服了以往參數(shù)論證方法中，只選擇幾個(gè)（一般為3-4個(gè)）物理點(diǎn)進(jìn)行論證，無(wú)法滿足 全區(qū)所有物理點(diǎn)參數(shù)需求的矛盾，論證結(jié)果更加全面，可以彌補(bǔ)以往方法中論證點(diǎn)少的不 足，和彌補(bǔ)以往參數(shù)論證方法無(wú)法了解全區(qū)參數(shù)需求的缺點(diǎn)，該方法適合于所有二維、三維 地震勘探的觀測(cè)系統(tǒng)參數(shù)論證。
【附圖說(shuō)明】
[0035] 圖1為本發(fā)明的基于地質(zhì)模型的全方位觀測(cè)系統(tǒng)參數(shù)論證方法的一具體實(shí)施例 的流程圖；
[0036] 圖2是本發(fā)明的一具體實(shí)施例中某一層位地質(zhì)模型；
[0037] 圖3是本發(fā)明的一具體實(shí)施例中某一層位中考慮動(dòng)校拉伸的最大炮檢距分析展 示圖；
[0038] 圖4是本發(fā)明的一具體實(shí)施例中某一層位中考慮速度分析精度的最大炮檢距分 析展示圖。
【具體實(shí)施方式】
[0039] 為使本發(fā)明的上述和其他目的、特征和優(yōu)點(diǎn)能更明顯易懂，下文特舉出較佳實(shí)施 例，并配合所附圖式，作詳細(xì)說(shuō)明如下。
[0040] 如圖1所示，圖1為本發(fā)明的基于地質(zhì)模型的全方位觀測(cè)系統(tǒng)參數(shù)論證方法的流 程圖。
[0041] 在步驟101，采集多種地質(zhì)數(shù)據(jù)，建立精確的工區(qū)地質(zhì)模型。采集的地質(zhì)數(shù)據(jù)包括 層位數(shù)據(jù)、速度、密度、深度、斷層、頻率等資料，采用層位數(shù)據(jù)建立工區(qū)地質(zhì)模型，并將各 層位對(duì)應(yīng)的速度、密度、斷層信息導(dǎo)入，每個(gè)物理點(diǎn)賦予其地質(zhì)屬性。如圖2所示，圖2是本 發(fā)明的一具體實(shí)施例中某一層位地質(zhì)模型。模型是計(jì)算的基礎(chǔ)，模型要求精確才能保證參 數(shù)計(jì)算的正確性。所述的收集工區(qū)資料中，層位數(shù)據(jù)的特征是txt格式，包含各個(gè)層位的每 個(gè)物理點(diǎn)的X(橫坐標(biāo)）、Y(縱坐標(biāo)）、Z(深度）信息，各個(gè)層位中每個(gè)離散物理點(diǎn)的密度應(yīng) 該小于該區(qū)地震勘探所需成像面元大小，速度資料、頻率資料、密度資料、斷層資料應(yīng)該與 層位數(shù)據(jù)是對(duì)應(yīng)的。所述的建立模型中，各個(gè)層位中的每個(gè)物理點(diǎn)需要賦予對(duì)應(yīng)的速度、密 度，在后續(xù)的公式計(jì)算中需要調(diào)用這些參數(shù)進(jìn)行計(jì)算，為使模型建立的精細(xì)，需要導(dǎo)入斷層 信息。流程進(jìn)入到步驟102。
[0042] 在步驟102,選擇其中的一個(gè)層位，采用公式計(jì)算層位上每個(gè)物理點(diǎn)的面元的需 求，采用兩種公式分別進(jìn)行計(jì)算，分別是滿足最高無(wú)混疊頻率法則和滿足橫向分辨率法則 公式。所述的計(jì)算每個(gè)物理點(diǎn)的參數(shù)需求中，面元的選擇主要是依據(jù)兩種法則，
[0043] 第一是滿足最高無(wú)混疊頻率需求，滿足偏移成像時(shí)不產(chǎn)生偏移噪音，即滿足最高 無(wú)混疊頻率法則，依據(jù)公式：b=Vint/(4XFmaxXSin0 )
[0044] 式中：b-面兀大小，Vint-上一層層速度，F(xiàn)max-最商無(wú)混疊頻率，9 地層 傾角。
[0045] 根據(jù)以上公式，計(jì)算各層位上每個(gè)物理點(diǎn)的面元需求。
[0046] 第二是滿足橫向分辨率要求，依據(jù)保證良好橫向分辨率面元大小的經(jīng)驗(yàn)公式：
[0047]b=Vint/