本發(fā)明涉及油田開發(fā)，特別是涉及到一種基于鉆井約束的速度模型校正方法。

背景技術：

1、隨著我國油氣勘探開發(fā)不斷推進，三維地震技術已經(jīng)成為勘探、開發(fā)、生產(chǎn)整個流程中每個環(huán)節(jié)必不可缺的關鍵內(nèi)容，通過三維地震資料不僅要研究地下的形態(tài)構造，而且越來越關注三維空間的儲層預測問題。而隨著各類地震屬性分析技術以及反演技術的快速發(fā)展，通過三維地震資料可獲得的信息也來越來多，比如豐富的地下三維空間的有關巖性、物性、流體的分布信息。然而，由于地震勘探得到的地震資料往往都是時間域的，即使進行三維地震資料深度偏移處理，也很難得到準確的深度信息，特別是與鉆井深度的對應尤為困難，無法滿足研究人員的精度要求。

2、目前，前人對于二維構造層面的時深轉(zhuǎn)換做了大量的研究，已經(jīng)比較成熟。然后，這些方法并未對三維地震資料中針對儲層內(nèi)部的小層進行精細的時深轉(zhuǎn)換以及與鉆井深度、巖性對應等問題，轉(zhuǎn)換后得到的地震資料往往容易出現(xiàn)變形和穿層的現(xiàn)象。

3、速度建模的準確性和合理性將決定地震資料時深轉(zhuǎn)換后的精度，隨著研究的不斷深入，目前的速度建模方法主要分為以地震為主的速度建模的方法和以測井為主的速度建模方法，然而這兩種方法都有著局限性，以地震為主的速度建模方法往往在縱向上準確度較差，以測井為主的速度建模方法則是在橫向上準確度較低。

4、在申請?zhí)枺篶n201510080055.2的中國專利申請中，涉及到一種井約束二維地震變速度場非線性誤差校正方法，該誤差校正方法包括步驟：s1、通過插值法推演換算，以改善、補正二維工區(qū)資料品質(zhì)差區(qū)段及無地震數(shù)據(jù)區(qū)的疊加速度，對網(wǎng)格點上所有疊加速度進行多項式擬合，建立初始疊加速度場；s2、生成井口平均速度；s3、構造誤差函數(shù)，即分別制作井旁疊加速度與井口平均速度、以及與非井旁地震道疊加速度的誤差曲線；s4、將步驟s3得到的非線性變化函數(shù)值作為校正量，對網(wǎng)格點內(nèi)地震疊加速度進行誤差校正，建立更合理精度更高的變速度場；s5、利用校正后的變速度平均速度模型，對解釋的時間層位進行時深轉(zhuǎn)換，即實現(xiàn)變速度成圖。該發(fā)明方法有助于降低變速度場誤差，從而提高變速度構造成圖精度。

5、在申請?zhí)枺篶n201210462223.0的中國專利申請中，涉及到一種地層的層速度模型的建模方法。所述方法包括：采集工區(qū)的地震數(shù)據(jù)；利用地震數(shù)據(jù)以及工區(qū)的鉆井的數(shù)據(jù)，產(chǎn)生層速度場；對層速度場執(zhí)行反射層析成像法速度迭代處理，從而生成優(yōu)化的層速度場；從優(yōu)化的層速度場提取各個地層的層速度平面圖；利用鉆井的數(shù)據(jù)對提取的各個地層的層速度平面圖進行校正；基于校正后的層速度平面圖生成層速度場；當用于進行校正的校正量超出預定范圍時，針對生成的層速度場進行到執(zhí)行反射層析成像法速度迭代處理的步驟。

6、在申請?zhí)枺篶n201610367256.5的中國專利申請中，涉及到一種構建地層速度模型的方法，包括：(a)獲取工區(qū)現(xiàn)有的地震資料；(b)獲得工區(qū)的疊前時間偏移地震數(shù)據(jù)體和疊前時間偏移成像速度；(c)獲得工區(qū)中各個地層的時間域的地震反射層位數(shù)據(jù)和斷層數(shù)據(jù)；(d)獲取工區(qū)的疊前深度偏移速度體；(e)根據(jù)聲波測井曲線和疊前時間偏移地震數(shù)據(jù)體確定工區(qū)中鉆井的時深曲線；(f)獲取工區(qū)中鉆井的井震偽速度；(g)建立工區(qū)的格架模型；(h)根據(jù)疊前深度偏移速度體、鉆井的時深曲線、鉆井的井震偽速度和所述工區(qū)的格架模型建立地層速度模型；(i)利用隨鉆資料對地層速度模型進行校正，以獲得隨鉆地層速度模型。根據(jù)所述方法，能夠有效提高地層速度模型的精度。

7、在申請?zhí)枺篶n201911325343.4的中國專利申請中，涉及到一種井控速度場建模方法，步驟如下：s1、建立微測井淺地表低速體模型，并在低速體模型上標定高速頂界面；s2、利用時距曲線對低速體模型進行時差校正，獲得淺層區(qū)域速度模型；s3、建立深層區(qū)域的初始速度模型；s4、采集測井信息建立井速度模型，并利用井速度模型校正初始速度模型；s5、對校正后的初始速度模型進行偏移處理，獲取連井剖面；s6、基于連井剖面查找速度異常點，并校正速度異常點；s7、重復步驟s5、s6，直到?jīng)]有速度異常點，獲得深層區(qū)域速度模型，進而獲得最終的井控速度模型。該發(fā)明方法充分利用了微測井和測井資料，在實際使用過程中可以有效提高速度建模的精度，減少建模誤差，適應于多種地形。

8、以上現(xiàn)有技術均與本發(fā)明有較大區(qū)別，未能解決我們想要解決的技術問題，為此我們發(fā)明了一種新的基于鉆井約束的速度模型校正方法。

技術實現(xiàn)思路

1、本發(fā)明的目的是提供一種通過反距離加權創(chuàng)建速度模型功能計算得到層速度的基于鉆井約束的速度模型校正方法。

2、本發(fā)明的目的可通過如下技術措施來實現(xiàn)：基于鉆井約束的速度模型校正方法，該基于鉆井約束的速度模型校正方法包括：

3、步驟1，構建初始速度模型；

4、步驟2，通過速度譜數(shù)據(jù)建立層速度模型；

5、步驟3，井校時間域模型；

6、步驟4，進行模型約束的速度反演；

7、步驟5，進行速度轉(zhuǎn)換，并計算井校平均速度；

8、步驟6，通過平均速度進行時深轉(zhuǎn)換，得到深度域的反演層速度、地震屬性體、層位數(shù)據(jù)。

9、本發(fā)明的目的還可通過如下技術措施來實現(xiàn)：

10、在步驟1，在層序?qū)涌丶s束下，結(jié)合速度譜數(shù)據(jù)或者精細的井震標定時深關系數(shù)據(jù)，建立初始速度模型。

11、步驟1包括：

12、步驟11、將速度譜插值為地震三維網(wǎng)格下的疊加速度體；

13、步驟12、井震標定得到時深關系，該時深關系是一維的時間和深度對應的數(shù)據(jù)；

14、步驟13、按照時間順序調(diào)整層位的層序順序，并建立層位構造框架結(jié)構；

15、步驟14、將解釋的抽稀層位插值為滿覆蓋層位；

16、步驟15、使用井軌跡和解釋層位數(shù)據(jù)參與求交點計算，得到層位處的時間和深度；將其保留作為約束點數(shù)據(jù)。

17、在步驟2，進行基于層控約束的地震層速度建模，首先依據(jù)地震解釋的層位，建立本地區(qū)的構造格架，在構造格架的基礎上考慮速度橫向變化和縱向變化的速度場建模。

18、步驟2包括：

19、步驟21、依據(jù)dix公式，將第一步插值得到的疊加速度體，轉(zhuǎn)換為層速度體；

20、步驟22、依次沿著層位框架，從淺部層位到深部層位，提取沿層的層速度，作為層位的屬性；

21、步驟23、提取層位上的層速度趨勢面，獲得沿層的層速度變化趨勢；

22、步驟24、對層位上的層速度進行平滑處理，使得沿層的層速度是一個較為平滑的速度場，去除初始層速度中的畸變數(shù)據(jù)；

23、步驟25、依據(jù)層速度趨勢面，對光滑后的層速度進行校正；

24、步驟26、層間按照沉積模式包括同沉積、頂超、削截、底超，對層間速度進行小層插值。

25、在步驟3，利用井旁測井數(shù)據(jù)或者隨鉆獲得的聲波測井數(shù)據(jù)，進行測井數(shù)據(jù)歸一化，并對目的層建立的精細的層速度模型進行校正，實現(xiàn)基于鉆井約束的實時速度模型建模。

26、步驟3包括：

27、步驟31、對聲波測井進行去野值、缺失值插值預處理；

28、步驟32、提取聲波測井的壓實趨勢，依據(jù)該壓實趨勢對聲波測井的頂部缺失段進行延拓；

29、步驟33、分析測井數(shù)據(jù)目的層段的值域分布規(guī)律，對聲波測井進行標準化校正；

30、步驟34、沿著井軌跡，提取初始層速度模型的井旁道；

31、步驟35、對聲波測井進行低頻提取，獲得聲波測井的低頻曲線，依據(jù)該低頻曲線，對井旁道初始層速度進行校正；

32、步驟36、獲得所有井的校正數(shù)據(jù)，對初始層速度模型進行體校正；

33、步驟37、井校正后就可以得到比較準確的時間域?qū)铀俣饶Ｐ汀?/p>

34、在步驟4，將井校之后的速度體作為初始模型，用模型進行約束，加入地震信息，提取模型與地震數(shù)據(jù)的權重系數(shù)，進行速度反演，得到精細層速度體、密度體及阻抗體。

35、步驟4包括：

36、步驟41、依據(jù)疊后地震數(shù)據(jù)，提取該地震數(shù)據(jù)的子波；

37、步驟42、依據(jù)該子波，對已知井的井旁速度模型道進行反演參數(shù)提??；

38、步驟43、確定模型修改量和地震數(shù)據(jù)影響權重系數(shù)，確保反演得到速度的分辨率達到要求；

39、步驟44、在確定反演參數(shù)后，對校正后的層速度模型進行廣義線性速度反演。

40、在步驟5，利用層速度體、均方根速度體進行轉(zhuǎn)換獲得平均速度體；通過井時深數(shù)據(jù)，選擇地震層位的沉積模式，對平均速度體進行速度場空間校正，獲得更好、更為準確的平均速度體。

41、步驟5包括：

42、步驟51、基于高精度的層速度，進行平均速度轉(zhuǎn)換；

43、步驟52、提取已知井的時深關系，并進行平均速度提??；

44、步驟53、沿著井軌跡，將井上的平均速度和初始轉(zhuǎn)換的平均速度，進行比較，得到平均速度的修正量，沿著層位對該修正量進行插值，獲得層位上的平均速度修正量；

45、步驟54、依據(jù)層位框架模型，依次對層位上的平均速度修正量進行三維空間體插值，獲得三維平均速度修正體；

46、步驟55、將三維平均速度修正體應用到初始的平均速度體，得到高準確度的平均速度體。

47、在步驟6，在建立的高精度層速度場的基礎上可以進行平均速度場轉(zhuǎn)換，利用平均速度體，將層速度體、地震數(shù)據(jù)體、屬性體、層位、斷層這些數(shù)據(jù)在時間域和深度域之間快速、準確的實現(xiàn)轉(zhuǎn)換。

48、步驟6包括：

49、步驟61、在高精度、高準確度的平均速度體控制下，依次將時間域的速度體，地震屬性體轉(zhuǎn)換到深度域；

50、步驟62、在高精度、高準確度的平均速度體控制下，依次將時間域的層位轉(zhuǎn)換到深度域；

51、步驟63、在高精度、高準確度的平均速度體控制下，依次將時間域的斷層轉(zhuǎn)換到深度域。

52、本發(fā)明中的基于鉆井約束的速度模型校正方法，基于速度模型校正技術，在構造、層序模型的約束下，將速度譜等資料與鉆井的聲波、vsp等數(shù)據(jù)結(jié)合，對速度模型中的平均速度場、層速度場進行校正。本發(fā)明則將地震資料與測井資料相結(jié)合，在建立井區(qū)初始三維時間域?qū)铀俣饶Ｐ突A上，結(jié)合隨鉆測井、錄井信息，獲得更為準確的三維時間域?qū)铀俣饶Ｐ?，并通過精細的井震標定技術，獲得準確的聲波測井時深關系，將三維時間域速度模型轉(zhuǎn)換為三維深度域速度模型，通過將地震速度場與測井相結(jié)合，從而實現(xiàn)橫向、縱向準確性較高的時深轉(zhuǎn)換。