【導讀】雙狐變速成圖系統(tǒng)是在多年的塔里木變速成圖攻關(guān)和生產(chǎn)過程中逐漸發(fā)展和完善起來的，和模型層析法在內(nèi)的四種成圖方法。進的設(shè)備、技術(shù)、數(shù)學理論都應用于此。雙狐軟件的主要開發(fā)者是在多年的油田科研生產(chǎn)中脫。內(nèi)外油田專家的認可。

　　

【正文】 d，輸出 (x,y)文件。在填寫參數(shù)時，第 3 行是時間的范圍，通常根據(jù)作圖的深度初步給個范圍，大于最大時間 500 毫秒足夠了（如最大時間 4500 毫秒，我們給 5000） ，最小通常給 0，可以減少數(shù)據(jù)量，第 1 行是 x 坐標，第 2 行是 y 坐標，根據(jù)坐標范圍給定 deltp 的值，使初始網(wǎng)格數(shù)目 number 合適，通常 10number30 即可，然后加密次數(shù)給 3（2 的 3 次方） ，這樣總共的網(wǎng)格點數(shù)目就是252222522227222＝8640000，可信度通常設(shè)置在 10 到 20 左右，根據(jù)資料情況適當調(diào)整，不可給太大，否則平滑太厲害，把速度變化的真實規(guī)律丟掉了。生成的文件大小為 ，因此網(wǎng)格數(shù)如果給的太大，數(shù)據(jù)量很大，計算忙，網(wǎng)格數(shù)要給的合適。雙狐軟件 71 圖 421 三維網(wǎng)格化參數(shù)第 112 步：調(diào)用 dfCreatef(x,y).exe 對 數(shù)據(jù)進行網(wǎng)格華，看偏移 T0 圖，在此只需修改網(wǎng)格數(shù)、加密次數(shù)、平滑參數(shù)，其它不用修改，具體使用前面已有詳細說明，生成網(wǎng)格文件 (x,y)，等值線文件 ，文件路徑和文件名都自動傳遞，不用修改。第 13 步：調(diào)用 DFVf(xyz) 模塊，輸入平均速度體 (x,y)，由 T0 時間 計算深度，輸出初始深度文件 ，文件名也是自動傳遞，不需要做任何修改。圖 422 速度體分析第 115 步：對計算出的初始深度進行網(wǎng)格華，注意在網(wǎng)格華時平滑大一點，獲得光滑的深度曲面。輸出曲面網(wǎng)格文件 (x,y)，等值線文件 ，第 16 步：調(diào)用 dfDrawf(x,y).exe 模塊，在 file\open 中打開 (x,y)，然后 file\save differential，選擇 4D，輸出 x、y 方向的偏導數(shù)，文件名自動傳遞，生成 (x,y)、 (x,y) 兩個文件。第 118 步：調(diào)用 模塊對偏導數(shù)做 3d 體網(wǎng)格華，在此對偏導數(shù)的平滑要大一點，如 50。只需要選擇 3d、修改網(wǎng)格數(shù)目、加密次數(shù)、給平滑系數(shù)，其它不用修，文件名自動傳遞。第 120 步：調(diào)用 模塊，對初始深度進行反偏移，輸出偏移后的深度散點，每個點的偏移量圖 ，也就是時間的反偏移量。雙狐軟件 72 圖 423 空間曲面偏移 第 21 步：根據(jù)反偏移量對時間 進行反偏移，輸出反偏移時間 。第 22 步：調(diào)用 DFVf(xyz)，用反偏移時間到疊加的平均速度場中提取平均速度，輸出 。第 23 步：由平均速度 和時間偏移量計算偏移后的平均速度 。第 225 步：調(diào)用 dfCreatef(x,y).exe 對偏移平均速度進行適當平滑，生成曲面網(wǎng)格文件(x,y)，等值線文件 。只需修改網(wǎng)格參數(shù)、加密次數(shù)、平滑系數(shù)，文件名自動傳遞。第 26 步：調(diào)用 z=f(x,y).exe 模塊，用時間 乘以速度 (x,y)，即（h=1/2*t*v）得到深度 。第 27－37 步：基準面校正，如果基準面為統(tǒng)一基準面（如 1300） ，就用深度 直接減去 1300。如果是浮動基準面，就有垂直校正和傾斜校正兩種方式，垂直校正是將基準面網(wǎng)格化以后，用深度直接減去基準面，傾斜校正需要求浮動基準面的偏導數(shù)，傾斜校正時深度要校正，平面方向坐標位置也要變，同事斷層也要校正，使用的是 ：輸出校正后的深度文件 。雙狐軟件 73 圖 424 浮動基準面傾斜校正第 339 步：對校正后的深度做網(wǎng)格化，只需修改網(wǎng)格化參數(shù)、平滑系數(shù)、文件名自動傳遞，輸出網(wǎng)格文件 (x,y)，輸出等值線文件 。第 40 步：井誤差校正，調(diào)用 模塊，調(diào)用時自動讀入井分層文件，如果數(shù)據(jù)格式不對，可以編輯，然后點“輸入井分層”重新讀入文件，然后根據(jù)分層文件的格式選擇 x坐標所在的列，y 坐標所在的列，井口高程所在的列 h0，目的層所在的列 dh，計算方式和文件名不用修改，然后點擊“計算誤差” ，生成的 是井點誤差對比值， 是井點誤差和邊界誤差值，邊界誤差是根據(jù)井點誤差通過反距離加權(quán)計算出來的，然后生成校正量曲面，如果只有一口井，直接將誤差值傳遞給校正量，如果多口井，調(diào)用網(wǎng)格化模塊對誤差進行網(wǎng)格化，在網(wǎng)格化時根據(jù)誤差點的分布有時需要人為控制，比如加控制點等，保證校正時不出現(xiàn)畸變，生成的校正量曲面 (x,y)，然后用校正量對深度進行校正，如果計算方式不改，那么校正方式也不需要修改，文件名自動傳遞，輸出校正后的文件 。雙狐軟件 74 圖 425 井誤差校正第 41 步：將深度值乘以－1，變?yōu)楹０?0 以下，輸出文件 。第 42 步：網(wǎng)格化結(jié)果，得到構(gòu)造圖，然后在 dfdraw 圖形平臺中進行修飾輸出。層位控制法成圖：層位控制法也是基于 Dix 公式計算，先計算出地層傾角，然后將通過疊加速度傾角校正將其轉(zhuǎn)換成均方根速度，然后通過公式 Dix 公式計算出層速度，最后轉(zhuǎn)換成平均速度。下面來介紹量板法數(shù)據(jù)準備和變速流程：數(shù)據(jù)準備：（如圖 426 所示）輸入速度場文件（ *.dfv）和層位斷層文件名創(chuàng)建列表（可參考前面的量板法） ，輸入 Vsp 數(shù)據(jù)可以選擇（可是看幫助，一般不建議使用） ，網(wǎng)格參數(shù)以及其他參考量板法和平均速度法。雙狐軟件 75 圖 426 層位控制法流程生成流程生成：（如圖 427 所示）層位控制法流程。圖 427 層位控制法流程步驟詳解：第一步：T0 時間控制層生成（同量板法） ，如果速度沒有進行傾角校正，請在“計算傾角”上打勾，還可以每層 T0 劃分小層。第二步：疊加速度數(shù)據(jù)準備， （如圖 428 所示）由第一步計算出的傾角進行 X 和 Y 方向的偏導三維網(wǎng)格化（注意光滑程度） ，然后疊加速度傾角校正（如圖 429 所示） ，轉(zhuǎn)換成均方根速度，最后速度體網(wǎng)格。雙狐軟件 76 圖 428 網(wǎng)格化 X 方向偏導圖 429 疊加速度傾角校正第三步：層速度計算（如圖 430 所示）：運用 Dix 公式計算出各層的層速度。雙狐軟件 77 圖 430 計算層速度第四步：層速度平滑（如圖 431 所示）：一般選擇分步執(zhí)行，以便觀測每一層的層速度是否有異常點，單獨給出每一層的平滑參數(shù)。圖 431 層速度平滑流程第五步：速度場的建立（如圖 432 所示）：通過各層 T0 和斷層以及各層層速度建立平均速度場。雙狐軟件 78 圖 432 平均速度場建立第 6 步：平均速度體網(wǎng)格化（同量板法和常規(guī)變速） 。第 7 步：應用：對任意給定的 t0 時間層位，要做變速，就可以用時間層位到上面形成的數(shù)據(jù)體中去插值求對應深度。基準面校正和井誤差校正可以參看常規(guī)變速成圖。模型層析法成圖：前面介紹的常規(guī)變速空校是基于迪克斯（Dix）公式進行速度計算的，它是假設(shè)在地層傾角不大的情況下，疊加速度近似于均方根速度。但在實際生產(chǎn)中，可能地層傾角較大，如山前地區(qū)，由于構(gòu)造運動，地層傾角很大，有的甚至在 600到 900 之間，這時迪克斯（Dix）公式已經(jīng)不實用，必須尋求新的速度計算方法。模型層析法是解決高陡構(gòu)造速度計算的很好的解決方法。模型層析法變速成圖思路： 當?shù)谝粚拥膶铀俣纫阎?，?t0 時間圖和第一層層速度可以計算射線的初始入射角和第一層的初始深度，由入射角等于反射角的原理，修改第二層層速度，利用曲射線追蹤迭代和非雙曲線擬合，當擬合的疊加速度與已知速度相等時停止跌代，即求出第二層的層速度，同時求出射線路徑，有射線路徑就可以計算反射界面和反射點偏離入射點的位置（即偏移量） 。以此類推，逐層計算 。下面以實例來介紹模型層析法的流程和步驟：數(shù)據(jù)準備：在做之前首先準備數(shù)據(jù)，數(shù)據(jù)格式與平均速度法使用數(shù)據(jù)格式相同，不同的地方是 T0 數(shù)據(jù)用疊加的，如果提供的數(shù)據(jù)是偏移的，那么在之前要做反偏移，得到疊加的時間，而且要有多層 T0 數(shù)據(jù)，上面要有控制層，最后是將有傾角變化的層解釋出來，一層數(shù)據(jù)是不能使用該方法的，數(shù)據(jù)單位統(tǒng)一成米單位。生成執(zhí)行步驟：調(diào)用 模塊，如下：圖 433雙狐軟件 79 選擇準備好的疊加速度場，設(shè)置工區(qū)的參數(shù)，可以用 dfCreatef(x,y).exe 網(wǎng)格化一個范圍比較大的 T0 文件，然后在文件\另存為下給文件名，擴展名用“ini”,保持網(wǎng)格化參數(shù)，再通過讀參數(shù)文件來獲取參數(shù)。然后給輸出文件名（如 ） ，保持路徑自動生成，然后生成處理流程，生成處理流程時自動建立 t0、refraction、TH 子目錄，分別存放 T0 時間數(shù)據(jù)、射線追蹤數(shù)據(jù)和結(jié)果數(shù)據(jù)。生成的流程如下圖：圖 434第一步：獲得時間傾角在射線追蹤時，首先要確定射線的初始入射角，這個入射角可以通過時間傾角和層速度來獲得，時間傾角的獲得方式也有多種：第一種：由解釋層位來獲得：由解釋層位獲得時間傾角的思路是將已知的解釋層位加密，按等間隔劃分出多個小層，然后對各個小層時間作平面網(wǎng)格化，在網(wǎng)格化時適當光滑，獲得光滑曲面，然后再每層曲面上沿X 方向和 Y 方向分別求導數(shù)，所有層位的 X 方向?qū)?shù)放一個文件，格式為（x,y,t,dx） ，Y 方向?qū)?shù)放一個文件，格式為（x,y,t,dy ） ，這個導數(shù)文件就是時間傾角。將準備好的 T0 和斷層文件放到一個目錄中，如放到 T0 目錄。創(chuàng)建解釋層位文件名列表：這步調(diào)用 模塊，然后按從上到下的順序分別選擇層位的時間散點和斷層文件，如果沒斷層就不選，選好一層，然后點“添加文件”添加到下面的列表中，如果選錯可以刪除重選。選好所有的后直接確定， “保存文件名列表”的文件名程序自動傳遞，不需要修改。計算時間傾角：雙狐軟件 80 這步調(diào)用 模塊，這步中只需修改“每層劃分的小層數(shù)” ，如果用原始 T0 時間層位計算傾角，就改成“1” ，如果層位少，可以加密，如給“2” ，就是把每層加密一層，就是在原層位之間多插值一層，如給“3” ，表示每層之間插入 2 層。其它文件名自動傳遞，不用修改。生成一個執(zhí)行步驟 ，圖 435圖 436第二種：由疊加速度體獲得時間傾角在無法采取第一種方式獲得時間傾角的時候，我們可以通過從疊加速度體中來獲取時間傾雙狐軟件 81 角，這種方式獲得時間傾角的精度完全取決于疊加速度的質(zhì)量和譜點的密度（建議少用） 。模型層析法步驟詳解：在應用模型層析法射線追蹤時要準備數(shù)據(jù)，首先是疊加速度，然后是控制射線入射角的時間傾角及控制射線追蹤的控制層位。下面詳細介紹模型層析法的各執(zhí)行步驟：第一步：疊加速度數(shù)據(jù)準備由于在射線追蹤過程中，擬合的疊加速度要與已知疊加速度比較，因此在數(shù)據(jù)準備時要輸入疊加速度，為了保證疊加速度的密度和規(guī)律性，我們將雙狐格式的疊加速度轉(zhuǎn)換成（x,y,t,va）格式。對疊加速度體做三維空間立體網(wǎng)格化，保證速度在橫向和縱向的規(guī)律性。第三步：時間傾角網(wǎng)格化時間傾角確定了射線的入射角，由于地層的傾角變化是比較穩(wěn)定的，同時也為了保證射線追蹤的穩(wěn)定性，那么就必須獲得傾角的光滑曲面，在這里通過對時間傾角體即（x,y,t,dx）和（x,y,t,dy）做三維立體空間網(wǎng)格化來獲得，在網(wǎng)格化的時候進行大平滑，但不能太大，否則反映不出真正的傾角變化情況。第四步：獲得射線追蹤的控制層位在射線追蹤時是從上到下逐層追蹤，那么必須有進行追蹤的控制層位，這些控制層可以是解釋的已知層位，為了更加接近真實情況，可以對已知層位劃分小層，如把求傾角的層位當控制層。還有就是在工區(qū)范圍內(nèi)，劃分等時間間隔的水平層位來作為追蹤的控制層位，比如時間間隔為：100ms，要做到 5000ms，需要生成 50 個間隔為 100ms 的控制層文件。第五步：射線追蹤由于射線追蹤的步驟比較多，因此采用批處理來執(zhí)行，在射線追蹤之前首先生成射線追蹤的處理流程，在本程序中每一層的追蹤過程用一個批處理來實現(xiàn)，在每一層的追蹤過程中需要輸入該層的控制 t0 時間、疊加速度、時間傾角來計算出這層的層速度和射線路徑，輸入輸出文件都是通過參數(shù)文件方式來自動傳遞。輸入信息的參數(shù)文件 ，生成的追蹤步驟放在 refractio\ ，執(zhí)行追蹤時用 打開 進行追蹤計算。第六步：生成結(jié)果數(shù)據(jù)體通過前面的各層射線追蹤，在空間上就可以形成時間－深度、時間－x 方向偏移量、時間－y 方向偏移量的數(shù)據(jù)體，通過參數(shù)文件提取這些數(shù)據(jù)并分別對這三個數(shù)據(jù)體做三維空間網(wǎng)格化，生成數(shù)據(jù)體時可以生成疊加的和偏移的，如果生成疊加的，那么偏移量單獨生成 2 個文件，時深轉(zhuǎn)換時先用疊加時間時深轉(zhuǎn)換然后由偏移量進行偏移，如果生成偏移數(shù)據(jù)體，那么就將偏移量考慮在深度數(shù)據(jù)體中，成圖時直接進行時深轉(zhuǎn)換。第七步：應用對任意給定的 t0 時間層位，要做變速，就可以用時間層位到上面形成的三個數(shù)據(jù)體中去插值求對應深度和 x、y 反向的偏移量，并進行偏移歸位，如果是用疊加的 T0 數(shù)據(jù)成圖，選擇疊加剖面成圖，如果用偏移剖面成圖選擇偏移剖面。注意：在射線追蹤過程中輸入的控制層位要求是沒有做過偏移歸位的水平疊加 t0 數(shù)據(jù)。如果在計算時間傾角的時候，如果自動網(wǎng)格化時，出現(xiàn)不收斂的情況，主要是網(wǎng)格化參數(shù)給的不合適，這時停止計算，修改 ini 參數(shù)文件，然后重新執(zhí)行。雙狐軟件 82 （三） 、其它偏移反偏移：（圖 437 所示）偏移與反偏移一般采用直射線，可參考常規(guī)變速和模型層析生成的批處理流程。圖 437 偏移反偏移浮動基準面校正：（圖 438 所示）基準面校正有垂直和傾斜校正兩種，傾斜校正能更好的校正由于基準面傾斜引起的傾斜效應，保證有傾角地層的高點歸位。