井控各向異性速度建模技術(shù)在YKL地區(qū)的應(yīng)用

井控各向異性速度建模技術(shù)在YKL地區(qū)的應(yīng)用裴云龍;王立歆;鄔達(dá)理;蔣波;鄧鋒【摘要】Theurgentdemandforhighqualityimagingdata,especiallyhigh-precisiongeologicalstructure,isakeyprobleminthedevelopmentofYKLgasfield.Becauseofmissingnear-offsetdataandlowfoldofseismicdata,well-seismicmisfitforstructuraldepthpredictioninthisregioncan'tmeettheexplorationrequirementifweonlyuseconventionalvelocitymodeling.Therefore,well-controlledanisotropygrid-tomographytechniquewasadoptedtoconductvelocitymodelinglayer-by-layer.Basedonprecisechoosingthereasonablemarkersandbuildingtheinitialvelocitymodelintheselectedarea,weuseknownloggingdatatocalculateanisotropicparameterswithinthisregion.Thenweeliminatetheinfluenceoflargesetforstableformationvelocityongasreservoirdepth,andfromshallowtodeeplayerbylayerinversiononmarkerlayer,updateiterativeanisotropicparameters,velocityanddepthmigrationimaginguntilthepredicteddeptherrorsfromwelllogandseismicdataforeachmarkerlayertobereduced.Inaddition,ithasbeenconfirmedby3drillingwellsthatthemisfitbetweenthepredictedstructuraldepthsrespectivelyfromthefinalseismicmigrationresultsandloggingdepthmeasurementsinthisgasreservoirareallwithin5m.Sowecandrawaconclusionthatourtechniquecansignificantlyimprovetheimagingaccuracyforstructureandfractureofgasreservoirsoastoreducetheriskofproductiondrilling.%YKL氣藏開發(fā)迫切需要高品質(zhì)地震成像資料的支撐,尤其是要求有高精度的構(gòu)造成像資料.由于該區(qū)塊三維地震勘探資料存在缺失近炮檢距、覆蓋次數(shù)低等不利因素,因此,常規(guī)單一的速度建模方法難以解決該區(qū)井震誤差大的問題.為此,針對(duì)性地提出并聯(lián)合應(yīng)用了井控各向異性網(wǎng)格層析和逐層速度建模技術(shù).在選好區(qū)內(nèi)標(biāo)志層、建好初始速度模型的基礎(chǔ)上,利用區(qū)內(nèi)多口已知井信息求取各向異性參數(shù),先消除淺部大套穩(wěn)定地層速度對(duì)氣藏深度的影響,再由淺到深逐層對(duì)標(biāo)志層反演,更新迭代各向異性參數(shù)、速度及深度成像,直至各標(biāo)志層的井震誤差降低到規(guī)定范圍內(nèi).3口驗(yàn)證井的應(yīng)用結(jié)果顯示,實(shí)際資料最終成像成果的氣藏構(gòu)造深度預(yù)測(cè)值與實(shí)鉆井深度誤差均在5m內(nèi),氣藏構(gòu)造成像的精度得到了顯著提高,斷裂成像也得到明顯改善,降低了開發(fā)鉆井的風(fēng)險(xiǎn).期刊名稱】《石油物探》年(卷),期】2017(056)003【總頁數(shù)】10頁(P390-399)【關(guān)鍵詞】疊前深度偏移;井控;各向異性;網(wǎng)格層析;速度建模;井震誤差【作者】裴云龍;王立歆;鄔達(dá)理;蔣波;鄧鋒【作者單位】中國石油化工股份有限公司石油物探技術(shù)研究院,江蘇南京211103;中國石油化工股份有限公司石油物探技術(shù)研究院,江蘇南京211103;中國石油化工股份有限公司石油物探技術(shù)研究院,江蘇南京211103;中國石油化工股份有限公司石油物探技術(shù)研究院,江蘇南京211103;中國石油化工股份有限公司西北油田分公司,新疆烏魯木齊830011【正文語種】中文【中圖分類】P631隨著我國西部YKL氣藏開發(fā)的推進(jìn)，含水量上升成為該氣藏開發(fā)的主要矛盾，氣藏開發(fā)難度越來越大。YKL氣藏歷年實(shí)鉆井?dāng)?shù)據(jù)顯示，氣藏構(gòu)造深度預(yù)測(cè)值與實(shí)鉆數(shù)據(jù)之間誤差較大，尤其在構(gòu)造邊部更是如此，地震勘探成果資料不能精確顯示構(gòu)造和斷裂位置是開發(fā)過程中面臨的主要問題。用時(shí)間域偏移成像數(shù)據(jù)制作構(gòu)造圖，傳統(tǒng)的做法是變速時(shí)深轉(zhuǎn)換成圖［1－4］。利用DIX公式將均方根速度轉(zhuǎn)換成層速度，再將時(shí)間域?qū)游晦D(zhuǎn)換到深度域；或者利用區(qū)域綜合平均速度方法得到時(shí)深關(guān)系，直接將TO圖轉(zhuǎn)換為深度域。受地層傾角和速度橫向變化等因素的影響，時(shí)間域偏移成像數(shù)據(jù)的成像精度不夠準(zhǔn)確。另外，因速度誤差以及變速成圖方法差異的影響，井震誤差往往也較大，準(zhǔn)確預(yù)測(cè)鉆探目標(biāo)深度難度很大［5－6］。當(dāng)?shù)叵麓嬖谒俣犬惓ｓw時(shí)，射線經(jīng)過這些異常體將產(chǎn)生大的偏折，此時(shí)，若按疊前時(shí)間偏移(pre-stacktimemigration,PSTM)的共中心點(diǎn)道集進(jìn)行偏移疊加，顯然得不到真實(shí)的反射。疊前深度偏移的巨大優(yōu)勢(shì)在于，在共深度點(diǎn)道集中，如果速度模型正確，則對(duì)于某個(gè)深度點(diǎn)來說，其道集反射應(yīng)是水平的，由此我們可以判斷模型正確與否［7］。相對(duì)于PSTM，疊前深度偏移(pre-stackdepthmigration,PSDM)在成像精度等方面具有顯著的優(yōu)勢(shì)［8-12］,但因其是一個(gè)對(duì)速度很敏感的處理過程，往往會(huì)遇到建模難題［7］。目前，三維疊前深度偏移速度模型建立的方法有多種，如：測(cè)井約束速度模型優(yōu)化速度建模方法［13］；初始速度模型的建立、模型的修改與驗(yàn)證及模型的修飾性處理和可視化顯示方法［14］；利用地震同相軸的多種信息進(jìn)行層析反演，并在反演中加入含有地質(zhì)含義的傾角信息以增加反演過程的穩(wěn)定性和收斂速度，利用類似于疊加速度譜的斜率速度譜來拾取地震同相軸斜率方法［15］；疊前逐層層速度反演速度建模方法［1617］；以疊前深度偏移為基礎(chǔ)進(jìn)行偏移速度和層速度分析，對(duì)速度模型的層位結(jié)構(gòu)及速度縱、橫向變化，采用三維可視化監(jiān)控和交互方式修改的方法［18］；共聚焦點(diǎn)層析速度建模方法［19］；多域雙重迭代法優(yōu)化更新層速度模型的方法［20］；聯(lián)合應(yīng)用基于層位的層析反演和基于網(wǎng)格的層析反演兩種速度建模方法［21］；VTI介質(zhì)多參數(shù)聯(lián)合走時(shí)層析成像的方法［22］等目前應(yīng)用最廣、最為成熟的速度建模方法是網(wǎng)格層析速度建模方法［23］。各種速度建模方法都有其適應(yīng)性，由于YKL氣藏區(qū)三維地震勘探資料存在缺失近炮檢距資料等不利因素，單獨(dú)采用某一種速度建模方法難以解決該區(qū)塊三維資料的深度偏移準(zhǔn)確成像問題。針對(duì)上述難點(diǎn)問題，綜合采用了井控［24］各向異性網(wǎng)格層析逐層反演多種速度建模技術(shù)，建立了合理的深度偏移速度模型，采用Kirchhoff方法進(jìn)行疊前深度偏移成像，最終在深度域?qū)⒛康膶臃謱訑?shù)據(jù)直接生成構(gòu)造圖。對(duì)3口井進(jìn)行驗(yàn)證，實(shí)鉆井深與地震預(yù)測(cè)深度誤差都在5m內(nèi)，較好地解決了YKL氣藏井震誤差大的問題。同時(shí)，通過重新處理，成像的品質(zhì)也得到明顯提高，為氣藏開發(fā)提供了可靠的基礎(chǔ)資料。1．1井震誤差分析Kirchhoff疊前深度偏移層速度往往要高于測(cè)井的層速度，若將深度偏移成像深度直接進(jìn)行井深預(yù)測(cè)，絕對(duì)誤差會(huì)很大（見表1中的各向同性速度建模誤差）。在探區(qū)速度橫向變化很小的情況下，利用實(shí)際測(cè)井的層速度計(jì)算深度，將時(shí)間域成像數(shù)據(jù)通過井校變速成圖能夠得到較準(zhǔn)確的預(yù)測(cè)深度。YKL探區(qū)的地層產(chǎn)狀較為平坦，屬弱各向異性地區(qū)，但由于橫向速度存在一定的變化，時(shí)間域成像數(shù)據(jù)通過井校變速成圖存在較大誤差（見表1中的井校變速成圖誤差，19m以內(nèi)），尤其是井少的情況下，井間數(shù)據(jù)通過插值連接，難以得到精確的構(gòu)造形態(tài)。而采用井控各向異性網(wǎng)格層析逐層速度建模技術(shù)則不同，井間數(shù)據(jù)連接是靠數(shù)據(jù)本身驅(qū)動(dòng)的，在深度偏移成像條件約束下形成的構(gòu)造形態(tài)更為可靠，井震誤差?。ㄒ姳?中的井控VTI速度建模誤差，3米以內(nèi)）。1．2誤差原因分析分析認(rèn)為，影響氣藏井震誤差的主要原因有以下幾個(gè)方面。1） YKL氣藏探區(qū)三維地震資料采集參數(shù)為：8線8炮，接收線距200m,接收點(diǎn)距50m，炮線距300m，炮點(diǎn)距100m，滿覆蓋次數(shù)24次，面元25mx50m,最小偏移距600m，單邊接收。原始資料覆蓋次數(shù)太低，且缺失近偏移足巨數(shù)據(jù)，導(dǎo)致偏移后的道集在0~600ms基本沒有數(shù)據(jù)（圖1）,不利于建立近地表速度模型；淺層0~2500m的共成像點(diǎn)（CIP）道集中有效反射同相軸信息很少（圖），淺層剩余延遲（RMO）無法準(zhǔn)確求取，嚴(yán)重影響速度建模精度。2） 淺層的測(cè)井曲線數(shù)據(jù)不可靠（圖3）,標(biāo)定困難,無法準(zhǔn)確獲取淺層地質(zhì)分層,影響淺層速度估算的精度。3） 存在一定的速度各向異性問題,常規(guī)的各向同性偏移不能滿足該區(qū)的深度偏移需要,需采用各向異性參數(shù)建立速度模型。2．1井控各向異性網(wǎng)格層析速度建模技術(shù)思路及流程首先,利用疊前時(shí)間偏移得到的成像剖面及均方根速度,結(jié)合區(qū)內(nèi)地質(zhì)認(rèn)識(shí),建立初始速度模型。在此基礎(chǔ)上進(jìn)行各向同性網(wǎng)格層析速度建模及疊前深度偏移,對(duì)偏移后CIP道集進(jìn)行剩余延遲分析，逐步拉平CIP道集，優(yōu)化速度模型，將每個(gè)深度偏移CIP道集的同相軸拉平。其次,生成各向異性初始速度模型。選好區(qū)內(nèi)標(biāo)志層,利用區(qū)內(nèi)多口已知井信息求取各向異性參數(shù),先消除大套穩(wěn)定地層速度對(duì)氣藏深度的影響,將第1個(gè)淺部穩(wěn)定標(biāo)志層的井震誤差消除。然后,采用逐層建模技術(shù),從淺到深在已知井的控制下,經(jīng)過多次更新各向異性參數(shù)和深度域?qū)铀俣饶Ｐ?逐層做好中深部標(biāo)志層的偏移成像,直至各標(biāo)志層的井震誤差降低到規(guī)定范圍內(nèi)，確保最終YKL氣藏構(gòu)造成像準(zhǔn)確。井控各向異性網(wǎng)格層析逐層速度建模技術(shù)流程如圖4所示。2．2初始速度模型建立實(shí)際處理YKL地區(qū)資料，建立各向異性參數(shù)的初始模型包括以下3方面工作。1）預(yù)處理，特別是做好靜校正處理。本區(qū)采用低速帶微測(cè)井成果約束下的折射層析靜校正及地表一致性剩余靜校正聯(lián)合應(yīng)用，避免靜校正不當(dāng)引起不利影響，生成可靠的疊前道集，再通過多次Kirchhoff疊前時(shí)間偏移均方根速度優(yōu)化迭代得到偏移成像數(shù)據(jù)體，然后對(duì)偏移體進(jìn)行井震標(biāo)定落實(shí)大套穩(wěn)定地層，并對(duì)其進(jìn)行構(gòu)造解釋。2）利用解釋的大套層位按照空間插值方式約束均方根速度。通過DIX公式將均方根速度轉(zhuǎn)換為深度域?qū)铀俣?，再?duì)層速度進(jìn)行編輯、沿層外推和沿層平滑，使層速度與大套地層趨勢(shì)更吻合，得到初始深度域?qū)铀俣饶Ｐ?，進(jìn)行目標(biāo)線Kirchhoff疊前深度偏移。3）從淺至深優(yōu)化初始速度模型。經(jīng)過多輪各向同性網(wǎng)格層析速度模型更新優(yōu)化及目標(biāo)線Kirchhoff疊前深度偏移迭代，確保主要地層CIP道集拉平以及偏移成像效果最佳，將最終迭代更新的各向同性網(wǎng)格層析速度模型V作為井控各向異性（VTI）網(wǎng)格層析的初始速度模型V0。另外，通過各向異性參數(shù)掃描，獲得￡與6的比例關(guān)系，當(dāng)設(shè)定￡=6時(shí)，道集同相軸更平，模型速度與井速度擬合更好，層速度網(wǎng)格層析反演后，得到的構(gòu)造成像效果更好，且井震誤差小。各向異性參數(shù)的初始值設(shè)為：60=0，￡0=0。2．3井控各向異性網(wǎng)格層析速度建模概念各向異性網(wǎng)格層析速度反演原理是從地下成像點(diǎn)向地表進(jìn)行射線追蹤，將之前在成像道集上拾取的剩余延遲沿著該射線路經(jīng)進(jìn)行反投影從而獲得各向異性參數(shù)和速度的更新量。需要注意的是，速度及各向異性參數(shù)是在規(guī)則的稀疏網(wǎng)格上更新的。各向異性網(wǎng)格層析方程為：式中：Atray表示深度偏移道集的延遲剩余時(shí)（RMOs）;flin表示各向異性參數(shù)的線性函數(shù)；z,dk表示解層析成像方程得到的各向異性中間參數(shù)剩余量［25］。當(dāng)?shù)卣鸩ㄍㄟ^射線追蹤從地下構(gòu)造反射到地表時(shí)，時(shí)差是沿著追蹤射線的各向異性參數(shù)的一個(gè)線性函數(shù)，那么可以首先進(jìn)行射線追蹤，然后通過反演計(jì)算來更新各向異性參數(shù)，其中這些參數(shù)描述了以相速度波傳播的彈性參數(shù)，也就是湯姆森(Thomsen)各向異性介質(zhì)參數(shù)。在各向異性介質(zhì)中，Atray/AWmisstie(旅行時(shí)／井震誤差)也可以認(rèn)為是各向異性參數(shù)誤差的線性函數(shù)［26］，因此，公式(1)也可表示為：公式(2)定義為井控各向異性(VTI+WellTie)網(wǎng)格層析速度迭代方程。井控各向異性參數(shù)計(jì)算公式為：式中：AzI為由疊前深度偏移剖面獲得的各向同性層厚度；AzA為由測(cè)井時(shí)深關(guān)系獲得的各向異性層厚度；VaO為各向異性層速度；V0為各向異性初始層速度；6為P波在垂直方向上的差異［27－30］。2．4井控各向異性參數(shù)模型修正VTI各向異性介質(zhì)中的P波橫向速度比垂向速度快，測(cè)井速度是垂向分量，地震波速度包括垂向和橫向分量。在各向異性介質(zhì)中測(cè)井速度通常比地震各向同性的層速度小，如圖5所示。Thomsen參數(shù)8代表P波水平方向和垂直方向的關(guān)系，其大小體現(xiàn)了由于各向異性而導(dǎo)致的同相軸彎曲程度。￡的求取通常依靠拾取剩余延遲譜，經(jīng)過井控層析校正后，剩余延遲譜基本歸于零值附近，中遠(yuǎn)偏移距處的同相軸得到拉平，從而進(jìn)一步提高偏移成像的效果。鑒于該區(qū)資料缺失近偏移距信息且覆蓋次數(shù)低，過少的CIP無法形成可靠的剩余延遲譜，為此，采取井控各向異性網(wǎng)格層析速度建模方法，借助已知井信息，在偏移成像數(shù)據(jù)體上逐層統(tǒng)計(jì)井震誤差，并根據(jù)公式(3)求取6,且令8=6,獲得初始8,通過逐層多次迭代，直到誤差完全消除。井控各向異性參數(shù)模型修正的步驟如下。1）在時(shí)間域偏移成果數(shù)據(jù)上，通過井震標(biāo)定精確解釋主要標(biāo)志層（如：T42,T23，T43，其中T43是YKL氣藏目的層）。）將時(shí)間域T42轉(zhuǎn)到深度域，得到該層位的深度值，依據(jù)27口井信息求得的井震誤差，并計(jì)算A80,然后將初始速度模型V0和△60,Az0共3個(gè)參數(shù)代入各向異性（VTI）網(wǎng)格層析迭代方程（2）,求解該方程獲得速度更新量AV1,從而得到新的速度模型V1和新的偏移成像數(shù)據(jù)體。）將時(shí)間域按照新速度模型V1轉(zhuǎn)到深度域，得到新的地震深度值，統(tǒng)計(jì)27口井井震誤差可得到新的A81,并將新的速度模型V1,A81和As13個(gè)參數(shù)通過各向異性（VTI）網(wǎng)格層析迭代方程，獲得速度更新量AV2,得到新的速度模型V2和新的偏移成像數(shù)據(jù)體。）按照步驟3）進(jìn)行網(wǎng)格層析迭代，經(jīng)過4~6次迭代，最后A8n更新量趨于0,可認(rèn)為井震誤差基本消除，即消除了以上淺層速度對(duì)下覆地層的影響。同時(shí)和也隨著誤差的減小而減??；同理，按照成像過程地層的井震誤差經(jīng)過3~5次迭代可消除，最后對(duì)目的層進(jìn)行2~4次迭代消除井震誤差。2．5各向異性參數(shù)模型經(jīng)過逐層井控,采用各向異性網(wǎng)格層析速度建模對(duì)速度進(jìn)行更新,最終使各向異性的層速度與實(shí)鉆井分層速度一致，求取的3層的各向異性參數(shù)S分布在-0.12~0．07（圖6）。3．1速度模型與氣藏分層的對(duì)比YKL氣藏區(qū)不同速度建模方法得到的速度剖面結(jié)果如圖7所示。從圖7可以看出,不同速度建模方法得到的速度剖面存在明顯差異,井控各向異性網(wǎng)格層析速度建模方法得到的中間層速度變小了（圖7b）。YKL氣藏區(qū)不同速度建模方法主要標(biāo)志層與層速度的匹配對(duì)比結(jié)果如圖8所示。從圖8可以看出，井控各向異性網(wǎng)格層析速度場(chǎng)分層界面與氣藏分層界面得到較好的匹配，說明速度得到了修正。3．2井震誤差情況YKL氣藏區(qū)不同速度建模方法過部分井的疊前深度偏移剖面對(duì)比結(jié)果如圖9所示從圖9可以看出，各向同性網(wǎng)格層析速度建模法存在明顯的深度誤差問題；應(yīng)用井控各向異性網(wǎng)格層析逐層速度建模后，成像精度得到明顯改善，氣藏構(gòu)造地震分層深度與測(cè)井實(shí)鉆分層深度基本閉合。3．3實(shí)際井的驗(yàn)證YKL氣藏區(qū)不同速度建模方法3口驗(yàn)證井的井震誤差情況對(duì)比結(jié)果如圖10所示。由圖10a可以看出，各向同性網(wǎng)格層析成像速度建模方法的井震誤差大；由圖10b可以看出，井控各向異性網(wǎng)格層析成像速度建模的井震誤差小，目的層氣藏深度預(yù)測(cè)值與實(shí)鉆井分層深度誤差在5m內(nèi)。3．4剖面品質(zhì)及構(gòu)造圖YKL氣藏區(qū)新、老偏移剖面對(duì)比結(jié)果如圖11所示。從圖11可以看出，經(jīng)過井控各向異性網(wǎng)格層析速度建模和疊前深度偏移等技術(shù)的應(yīng)用，獲得的處理成像剖面品質(zhì)明顯提高，構(gòu)造成像精度更高，波組特征更清楚，斷裂清晰。YKL氣藏區(qū)不同方法得到的構(gòu)造圖對(duì)比結(jié)果如圖12所示。從圖12可以看出，兩種方法得到的構(gòu)造形態(tài)有一定差異，經(jīng)已有鉆井資料證實(shí)白堊系構(gòu)造成像精度得到明顯提高，為YKL氣藏主體區(qū)井位加密調(diào)整以及外圍滾動(dòng)外擴(kuò)提供可靠數(shù)據(jù)支撐。井控各向異性網(wǎng)格層析速度建模和逐層成像速度建模方法技術(shù)的聯(lián)合應(yīng)用，是目前深度域提高速度建模精度和成像精度較好的方法，也是目前工業(yè)界較為普遍的做法類似YKL氣藏區(qū)地震勘探覆蓋次數(shù)低，缺失近偏移距資料，近地表速度無法準(zhǔn)確獲取的難題，采用速度掃描、井控各向異性網(wǎng)格層析迭代及標(biāo)志反射層控制的手段解決淺層地層速度變化及各向異性等問題，是獲得高精度成像結(jié)果的關(guān)鍵，較好地解決了YKL氣藏井震誤差大的問題，得到的成果資料更好地滿足了氣藏開發(fā)對(duì)地震資料的高精度要求。在速度建模中重點(diǎn)考慮以下幾點(diǎn)：1） 處理解釋人員緊密結(jié)合，落實(shí)地震波組與測(cè)井分層的對(duì)應(yīng)關(guān)系，層位解釋、標(biāo)定必須準(zhǔn)確；2） 消除井震誤差首先需消除淺層大套穩(wěn)定地層速度對(duì)氣藏深度的影響，為后續(xù)深度成像迭代打好基礎(chǔ)；3） 對(duì)井震誤差要進(jìn)行合理性分析，做好井控約束外推的邊界處理，逐步降低井震誤差以提高成像精度；4） 盡可能充分利用多口已知井資料，提高約束的可靠性。致謝：本研究得到了西北油田分公司勘探開發(fā)研究院及項(xiàng)目組等有關(guān)專家的大力支持，在此致以衷心的感謝?！鞠嚓P(guān)文獻(xiàn)】孫開峰，葉勇?泛克里格法在三維地震資料變速成圖中的應(yīng)用[J]?石油物探，1998,37（2）：112－117SUNKF，YEY．ApplyingtheuniversalKirgingtothevariable－velocitymappingof3－Dseismicdata[J]．GeophysicalProspectingforPetroleum，1998，37（2）：112－117楊江峰，王咸彬，洪太遠(yuǎn)，等．準(zhǔn)噶爾盆地南緣山前帶構(gòu)造變速成圖[J]．石油物探，2008，47（2）：179－182YANGJF，WANGXB，HONGTY，etal．ApplicationofvaryingvelocitymappingtechniqueformountainfrontinsourthernJunggarbasin[J]．GeophysicalProspectingforPetroleum，2008，47（2）：179－182葉勇，孫武亮，孫開峰．塔中圍斜區(qū)連片精細(xì)速度建模與變速成圖[J]．石油物探，2010，49（4）：401－406YEY，SUNWL，SUNKF．Detailedvelocitymodelingandvelocity－variablestructuremappingofeight－piecingprocessinginTazhongarea[J]．GeophysicalProspectingforPetroleum，2010，49（4）：401－406[4]楊勤林，王彥春，張靜，等?濱里海盆地M區(qū)塊鹽下構(gòu)造變速成圖[J].石油物探，2012,51（4）：377－382YANGQL，WANGYC，ZHANGJ，etal．Variablevelocitymappingforsub－saltstructureinBlockM，Pre－CaspianBasin[J]．GeophysicalProspectingforPetroleum，2012，51(4)：377－382[5]JONESIF．3－Dprestackdepthmigrationandvelocitymodelbuilding[J]．TheLeadingEdge，1998，19(7)：897－906[6]RESHELM．Theuseof3－Dprestackdepthimagingtoestimatelayervelocitiesandreflectorposition[J]．Geophysics，1997，62(1)：206－210[7]李家康，樊佳芳，邊西燕，等?大規(guī)模并行處理機(jī)三維疊前深度偏移[J].石油物探，1998，37(2)：77－87LIJK，F(xiàn)ANJF，BIANXY，etal．3－Dprestackdepthmigrationonmassiveparallelprocessingsystem[J]．GeophysicalProspectingforPetroleum，1998，37(2)：77－87齊中山，王汝珍，敬朋貴，等．川東北地區(qū)膏鹽巖影響下的儲(chǔ)層成像技術(shù)與效果分析[J].石油物探，2011,50(6):607-611QIZS，WANGRZ，JINGPG，etal．TheimagingtechniquesofreservoirsinfluencedbysalinedepositintheNortheasternSichuanBasin[J].GeophysicalProspectingforPetroleum，2011，50(6):607－611劉西寧，劉司紅，馬秀國，等.疊前深度偏移技術(shù)在焉耆盆地寶南地區(qū)的應(yīng)用[J].石油物探，2006，45(2):197－201LIUXN，LIUSH，MAXG，etal.Theapplicationandinvestigationof3－DprestackdepthmigrationinBaonanareaofYanqibasin[J].GeophysicalProspectingforPetroleum，2006，45(2):197－201郝守玲，趙群.橫向速度變化對(duì)構(gòu)造成像影響的物理模擬研究[J].石油物探，2008，47(1):49－54HAOSL，ZHAOQ.Physicalmodelingresearchfortheeffectoflateralvelocityonstruc－turalimaging[J].GeophysicalProspectingforPetroleum，2008，47(1):49－54李鵬，劉志成，楊楠，等.玉北地區(qū)碳酸鹽巖縫洞型儲(chǔ)層精細(xì)成像方法應(yīng)用研究[J].石油物探，2015，54(4):443－451LIP，LIUZC，YANGN，etal.Applicationoffinetargetreservoirimaginginfractured－vuggycarbonatereservoirYubeiarea[J].GeophysicalProspectingforPetroleum，2015，54(4):443－451王麗，屠世杰，趙傳雪，等.疊前深度偏移技術(shù)在江蘇火成巖地區(qū)的應(yīng)用[J].石油物探，2006，45(1):79－82WANGL，TUSJ，ZHAOCX，etal.TheapplicationofprestackdepthmigrationinigneousrocksareaoftheJiangsuOilfield[J].GeophysicalProspectingforPetroleum，2006，45(1):79－82胡英，姚逢昌?用于疊前深度偏移速度建模的一種新方法[J].石油勘探與開發(fā)，2000,27(2)：62－64HUY，YAOFC．Anewapproachtoprestackdepthmigrationmodelbuilding[J]．PetroleumExplorationandDevelopment，2000，27(2)：62－64[14]方伍寶?三維疊前深度偏移的建模技術(shù)[J].石油物探，2002,41(2):132-135FANGWB．Themodelbuildingtechnologyin3－Dprestackdepthmigration[J]．GeophysicalProspectingforPetroleum，2002，41(2):132－135[15]井西利，楊長(zhǎng)春，李幼銘?建立速度模型的層析成像方法研究[J].石油物探，2002,41(1):72－76JINGXL，YANGCC，LIYM．Studyoftomographicinversionintheconstructionofvelocitymodel[J]．GeophysicalProspectingforPetroleum，2002，41(1):72－76[16]陳遵德，徐廣民，賀振華，等．三維逐層層速度反演方法及誤差分析[J]．石油物探，2000，39(4):12－18CHENZD，XUGM，HEZH，etal．The3Dinversionmethodoflayer－by－layerintervalvelocityanderroranalysis[J]．GeophysicalProspectingforPetroleum，2000，39(4):12－18[17]王者順，樊佳芳，王衛(wèi)江?疊前逐層成像速度建模方法及其應(yīng)用[J].石油物探，2004,43(4):384－387WANGZS，F(xiàn)ANJF，WANGWJ．Layer－by－layerprestackimagingvelocitymodelbuildingandapplication[J]．GeophysicalProspectingforPetroleum，2004，43(4):384－387[18]羅省賢，李錄明．三維疊前深度偏移速度模型建立方法[J]．石油物探，1999，38(4):1－6LUOXX，LiLM．Velocitymodelconstructionmethodsfor3－Dprestackdepthmigration[J]．GeophysicalProspectingforPetroleum，1999，38(4):1－6[19]辛可鋒，王華忠，馬在田，等．共聚焦點(diǎn)層析速度建模方法[J]．石油物探，2005，44:329－333XINKF，WANGHZ，MAZT，etal．TomographicvelocitymodelingbasedonCFPtechnology[J]．GeophysicalProspectingforPetroleum，2005，44(4):329－333[20]胡英，張研，陳立康，等．速度建模的影響因素與技術(shù)對(duì)策[J]．石油物探，2006，45:503－507HUY，ZHANGY，CHENLK，etal．Influencingfactorsandtechnicalstrategiesforvelocitymodeling[J]．GeophysicalProspectingforPetroleum，2006，45(5):503－507[21]馬彥彥，李國發(fā)，張星宇，等．疊前深度偏移速度建模方法分析[J]．石油地球物理勘探，2014，49(4):687－693MAYY，LIGF，ZHANGXY，etal．Strategyofvelocitymodelbuildinginpre－stackdepthmigration[J]．OilGeophysicalProspecting，2014，49(4):687－693[22]劉玉柱，王光銀，董良國，等.VTI介質(zhì)多參數(shù)聯(lián)合走時(shí)層析成像方法[J].地球物理學(xué)報(bào)，2014，57（10）：3402－3410LIUYZ，WANGGY，DONGLG，etal．JointinversionofVTIparametersusingnonlineartraveltimetomography[J]．ChineseJournalofGeophysics，2014，57（10）：3402－3410[23]李慧，成德安，金婧?網(wǎng)格層析成像速度建模方法與應(yīng)用[J].石油地球物理勘探，2013，48（增刊1）：12－16LIH，CHENGDA，JINJ．Velocitymodelbuildingbasedongridtomography[J]．OilGeophysicalProspecting，2013，48（S1）：12－16[24]康有元，馬豐臣，滿紅霞，等．面向油田開發(fā)的井控地震資料處理技術(shù)及應(yīng)用[