一階彈性波方程的變網(wǎng)格高階有限差分數(shù)值模擬

1、一階彈性波方程的變網(wǎng)格高階有限差分數(shù)值模擬2019年12月第43卷第6期正演技術(shù)一階彈性波方程的變網(wǎng)格高階有限差分數(shù)值模擬李振春3張慧張華（中國石油大學(xué)地球資源與信息學(xué)院,山東東營257061）李振春,張慧,張華.一階彈性波方程的變網(wǎng)格高階有限差分數(shù)值模擬.石油地球物理勘探,2019,43（6）:711716摘要使用可變網(wǎng)格的有限差分法進行地震模擬有許多獨特的優(yōu)點,主要表現(xiàn)為對地質(zhì)模型的離散化更為合理,在低速帶和復(fù)雜構(gòu)造區(qū)域,可將局部網(wǎng)格劃分得相對精細些,不僅提高了模擬精度,消除了因采樣不足導(dǎo)致的頻散現(xiàn)象,而且可以減少計算機內(nèi)存需求,保持模型計算的靈活性。本文提出一種新的基于高階交錯網(wǎng)格技術(shù)的

2、彈性波數(shù)值模擬方法,通過改變網(wǎng)格的空間步長實現(xiàn)了局部網(wǎng)格加密技術(shù),彌補了常規(guī)網(wǎng)格的缺陷和不足。試算結(jié)果表明,本文提出的算法穩(wěn)定性較好,且能夠提高模擬精度,減少計算時間,提高計算效率。關(guān)鍵詞變網(wǎng)格交錯網(wǎng)格高階有限差分彈性波1引言法技術(shù)（地震采集。為適應(yīng)復(fù)雜勘探條件的需要,、高效率的地震數(shù)值模擬技術(shù)變得越來越重要。自Virieux1率先采用交錯網(wǎng)格有限差分法模擬P2SV波在非均勻介質(zhì)中的傳播以來,交錯網(wǎng)格成為了有限差分法正演模擬中普遍使用的網(wǎng)格剖分技術(shù)。同二階位移方程相比,使用交錯網(wǎng)格的一階速度應(yīng)力方程,無論是模擬精度還是模擬的穩(wěn)定性都有了較大提高。但是如果局部地區(qū)構(gòu)造較為復(fù)雜,或者存在低速層時,

3、為了提高復(fù)雜構(gòu)造的正演精度,通常對整個模型使用較小的采樣網(wǎng)格,以保證對頻散現(xiàn)象的抑制和對穩(wěn)定性的保障。這對保持模型計算的靈活性非常不利,還會對模型的其他部分造成過采樣,增加驚人的計算量。所以,設(shè)計一種更合理的離散化方式進行正演模擬是非常必要的。從空間采樣的角度考慮,既要提高模擬精度、又要降低計算機內(nèi)存需求的最有效方法就是在模型的不同區(qū)域采用不同步長的網(wǎng)格，即變網(wǎng)格（也稱為不規(guī)則網(wǎng)格）。為此Jastram等提出了基于二維聲波方程針對某一深度變網(wǎng)格步長的算法2;Jastram等=vx3;4Pitarka實現(xiàn)了低階精度的不5。本文在前人研究的基礎(chǔ)上提出一種高階變網(wǎng)格數(shù)值模擬方法,實現(xiàn)了高精度的局部網(wǎng)

4、格加密技術(shù),另外還改進了前人對變化網(wǎng)格過渡帶采用的波場插值算法,避免了在過渡帶區(qū)域進行大量的插值計算,也不需要平滑函數(shù)進行平滑,同時還保持了算法的穩(wěn)定性。2方法原理文中所提出的算法是基于Virieux1和Graves6提出的一階速度應(yīng)力方程,對于各向同性介質(zhì),有txx+5t5tvx=b(5xzxz)+fxtxz+5t5tvz=b(5xzzz)+fz入+2u)5xvx+入5t5zvztxx=(入+2u)5zvz+入5t5xvxtzz=(5zvx+5xvz)5Ttxz二口(1)tt式中:vx,vz分別表示速度；txx,zz,xz分別表示應(yīng)力向量；5x、5z、5t分別表示一階偏微分算子5/5x、5/

5、5z、5/5t;b表示密度的倒數(shù)；fx,fz為體力向量；入及卩為拉梅常數(shù)。3山東省東營市中國石油大學(xué)地球資源與信息學(xué)院,257061本文于2019年1月8日收到,修改稿于同年7月30日收到。本課題由國家973專題(2019CB209605)、國家863專題(2019AA06Z206)和CNPC物探重點實驗室中國石油大學(xué)(華東)研究室資助。712石油地球物理勘探2019年通過對時間導(dǎo)數(shù)采用二階近似,得到彈性波一階速度應(yīng)力方程的離散形式為vxvn+2i,jn-2xxi+2,jni,jtxx+Dt+Atbx(Dxzxz)|ni,jnA1=dxi/2,A2=dxi/2n-1n-1k=1n+2zi+2,

6、j+2=vn-2zi+2,j+2t+Atbz(Dxxz+Dtzzz)|i+,j+22n+2i+2,jn+2i+2,jAn-1=k=1dxi+k+dxi/2,An=Edxi-k+dxi/2Ettn+1xxi+2,jn+1zzi+2,j=tnzzi+2,jDxg(x,z)入+2口)Dxvx+入+At(Dzvz圖1計算對稱點為i+的差分系數(shù)所需的變網(wǎng)格節(jié)點2n-1=T)Dzvz+入+At(入+2uDxvx同理，以i為對稱點的差分算子為An-1=n-1k=1n+1nTxz(Dzvx+Dxvz)xzi,j+2=txzi,j+2+At(n+2i,j+2dxEi+(k-1)+dxi+n-1/2,An=k=1

7、dxEi-k+dxi-n/2。(2)式中:Dx,Dz分別表示對x,z的一階微分算子；bx=(bi,j+bi+1,j)/2,bz=(bi,j+bi,j+1)/2;uxz=(1/ui,j+1/ui+1,j+1/ui,j+1+1/ui+1,j+1)/4。用變量g來表示速度VX、vz或應(yīng)力向量TXX、ttzz、xz,則在常步長的交錯網(wǎng)格技術(shù)中,用下式計算式(2)中的g關(guān)于x的2n階空間偏導(dǎo)數(shù)7Cg(x+AxAxEnii=1下面給出變網(wǎng)格算法的任意偶數(shù)階精度差分近似式及差分系數(shù)的求取方法,即nDxg(x,z)=i-1cE2i-1g(x+A2i-1,z)+(4)+c2i-A2i,z)n,)-i；i是空間差

8、分算子，它dg關(guān)于x。令g(x,z)=gzexp(ikx),則式(4)可寫成n(3)-g(x-Axi-lzx2n式中，系數(shù)Cin由x(x2iT,z)和g(x-Ax2i-1,z)。與常網(wǎng)格不同,變網(wǎng)格的差分算子是空間變化ik=i-lcE2i-lA2iT)+c2iexp(-ikA2i)exp(+ik(5)的。給定地質(zhì)模型,在確定網(wǎng)格剖分后就可求出各網(wǎng)格節(jié)點的差分算子供遞推計算時調(diào)用。在交錯網(wǎng)格技術(shù)中差分算子有兩個對稱點，即i+或i,不2同的對稱點對應(yīng)的Ai也不同。下面以對稱點i+2為例求取Ai,如圖1所示。對式(5)中的指數(shù)項進行泰勒展開，它的2n階泰勒展開式為2Ai)1+ikAi+(ik)2Ae

9、xp(iki+22n-12n(kAi)2n-1+o(A(6)+i)(2n-1)!把式(6)代入式(5)并整理,可以得到2222AA)()(c1AA2ik=(c1+c2+-+c2n)+ik(c1A-c+c+i1222n2n1+c2A2+c2n2n)+2+(i)336332n-1(c1AA31+C2A2+c2n2n)+(i)n-12n-12n-1(clA-c2A+c2nA2n)12(2n-l)!(7)2n-1寫成矩陣的形式1iA11-iA231iA2n-11-iA2nc1c2二2n-22n-10i(i)2n-2A12n-22n-1(i)2n-1A12n-2(i)2n-2A22n-12n-lA-(i

10、)2(i)2n-2A2n2n-1A-(i)2n2n-2(i)2n-2A2n-12n-1-(i)2n-1A2n-1c2n-1c2n00(8)第43卷第6期李振春等:一階彈性波方程的變網(wǎng)格高階有限差分數(shù)值模擬713通過解方程組(8),可以得到有限差分算子Dx的系數(shù)ci。同樣地，差分算子Dz的差分系數(shù)也是通過這種方法求得。對于交錯的網(wǎng)格，由于變量所定義的網(wǎng)格點的位置不同，對于每一個差分算子Dx或Dz都有兩組不同的差分系數(shù)。在使用變網(wǎng)格時，網(wǎng)格步長的變化可能會導(dǎo)致數(shù)值反射的出現(xiàn)。究其原因是因為波場離散后，相速度是網(wǎng)格步長的函數(shù)。當(dāng)相速度梯度較大時，即使速度和密度都沒有變化，入射波的能量也會部分反射回來

11、，導(dǎo)致了數(shù)值反射現(xiàn)象。因為頻散程度與頻率有關(guān)，所以當(dāng)入射波的頻率較低時，產(chǎn)生的數(shù)值反射較高頻時弱。在網(wǎng)格變化的區(qū)域?qū)Σ▓鲞M行插值計算雖可以在某種程度上壓制數(shù)值反射，但是插值算法計算量大，且效果并不令人滿意，所以本文沒有采用插值計算。本文提出的算法是根據(jù)各點所對應(yīng)的步長的變化計算其差分系數(shù)，在一個循環(huán)內(nèi)完成所有點的波場計算，不會產(chǎn)生數(shù)值反射，并且減小了計算量。式中:h=max(dxi,dzi)表示最大的空間網(wǎng)格步長；vmin是地震波速的最小值。本文所提算法需要滿足的穩(wěn)定性條件為Wmin(dxi，dzi)ni=1(-1)i-1Ci式中:Vmax為最大速度值；P為密度；Ci為式(3)中做2n階空間展

12、開時的系數(shù)。4數(shù)值模擬實例為了驗證提出的變網(wǎng)格差分算法的穩(wěn)定性及有效性,設(shè)計兩個模型正演實例來說明文中所述方法的實際效果。4.1。模X3000m,縱波速度為3500m/s,橫0125來計算，密度為常數(shù)。P波震源激發(fā)，震源主頻為30Hz,位于(1500m,1500m)處。變網(wǎng)格步長為10m,5m。與其進行對比的常規(guī)交錯網(wǎng)格步長為5m。數(shù)值模擬中，差分解法的精度為o(At2,Ax10)。彈性波波場快照和單炮記錄分別記錄于圖2和圖3中。從上面的彈性波波場快照和單炮記錄中可以看出,變網(wǎng)格的過渡區(qū)沒有產(chǎn)生數(shù)值邊界反射,說明本文提出的算法成功地消除了空間網(wǎng)格變化造成的數(shù)值反射,進而說明了該算法的合理性和良

13、好的穩(wěn)定性。3法較好地解決了邊界反射問題。使用透射邊界能夠吸收絕大部分邊界反射的能量,剩余能量通過衰減邊界再次減小,這樣計算量增加不大,吸收也較徹底,可以較好地消除人為邊界反射。穩(wěn)定性條件是有限差分數(shù)值模擬中十分重要的問題9。經(jīng)過驗證,利用本文所提出的算法能夠精確模擬的最大的地震波頻率為fmax圖2a變網(wǎng)格波場(左)和常規(guī)網(wǎng)格波場(右)水平分量快照(t=360ms)714石油地球物理勘探2019年第43卷第6期李振春等:一階彈性波方程的變網(wǎng)格高階有限差分數(shù)值模擬715412背斜模型模型大小為3000mX2000m,模型中橫波速度取泊松比為0.25來計算；P波震源激發(fā)，震源子波主頻為20Hz。采

14、用精度為。(At2,Ax10)的常規(guī)交錯網(wǎng)格有限差分法和變網(wǎng)格有限差分法進行模擬對比。模型參數(shù)及單炮記錄如圖4和圖5所示。從上面的模擬效果對比分析得知:在用常規(guī)交錯網(wǎng)格有限差分法模擬的炮記錄中,由于背斜的影響,自背斜頂界面反射以下出現(xiàn)了嚴重的數(shù)值頻散;而用變網(wǎng)格算法進行模擬時,由于在背斜低速區(qū)域采用小網(wǎng)格步長計算,故由背斜影響導(dǎo)致的頻散得到了較好地壓制,充分顯示了此算法較高的模擬精度。圖4背斜模型716石油地球物理勘探2019年5結(jié)論本文在網(wǎng)格過渡區(qū)域,沒有使用插值計算,而是根據(jù)各點所對應(yīng)的步長的變化計算其差分系數(shù),在一個循環(huán)內(nèi)完成對所有點的波場計算,成功地消除了數(shù)值反射,且同波場插值計算相比

15、減少了計算量。對于常規(guī)網(wǎng)格算法難以解決的地下復(fù)雜構(gòu)造的問題,采用變網(wǎng)格算法可較好地解決,且能較靈活地適應(yīng)地層模型的變化,對低速體的適應(yīng)性較常規(guī)網(wǎng)格算法有了較大改善另外,相對常規(guī)網(wǎng)格算法而言,變網(wǎng)格算法節(jié)省了計算機內(nèi)存,減少了計算時間,且模擬的區(qū)域越大,節(jié)省的時間就會越多,避免了對整個模型的精細采樣造成的資源浪費。本文提出的算法精度可達任意偶數(shù)階精度,高階變網(wǎng)格算法的提出,尤其是與交錯網(wǎng)格有限差分法相結(jié)合,可以有效地提高模擬精度和計算效率。參考文獻1VirieuxJ.P2media:velocity2stressfinite2differencemethod.Geo2physics,1986,5

16、1,889012JastramCandBehleA.Acousticmodelingonaverti2callyveryinggrid.GeophysProsp,1992,40,1571693JastramCandTessmerE.Elasticmodelingonagridofverticallyvaryingspacing.GeophysProsp,1994,42,357370FalkJ,TessmerEandGajewskiD.Tubewavemodel2ingbythefinite2differencemethodwithvaryinggridspacing.Prosp,1996,14

17、8,7793PitarkaA.3Delasticfinite2differencemodelingofseis2micmotionusingstaggeredgridswithnonuniformspacing.BulletinofSeismologicalofAmerica,1999,89,5468GravesRW.Simulatingseismicwavepropagationin3Delasticusingstaggered2gridfinite2differences.BullSeismSocAm,1996,86,109111077董良國,馬在田,曹景忠.一階彈性波方程交錯網(wǎng)格高階差分

18、解法.地球物理學(xué)報,2000,43(3):4114198廖振鵬,周正華,張艷紅.,2002,45(4):534545,地球物理學(xué)報,2000,43(6):(本文編輯:馮杏芝)(上接第688頁)虛線框在橫向上的分布與圖8中虛線框在橫向上的分布比較可以得知,阻抗低值處并非都富集天然氣。通過時頻分析技術(shù)和阻抗反演技術(shù)相結(jié)合,最終確定儲層的分布范圍為圖8中白色實線框內(nèi)。從對實際資料的應(yīng)用和分析中發(fā)現(xiàn),將時頻分析技術(shù)和波阻抗反演技術(shù)相結(jié)合可以更加精細地刻畫礁灘相儲集層的空間分布、油氣的富集范圍,計算油氣儲層的厚度,估算油氣的儲量,從而有效地對礁灘相儲層進行預(yù)測,降低鉆探風(fēng)險。只對垂直剖面進行了探討,建議

19、在實際的生產(chǎn)研究中也應(yīng)該考慮對沿層切片、水平時間切片及波阻抗剖面和單頻剖面疊合圖進行分析,從而更加準確地對礁灘相油氣藏進行預(yù)測。參考文獻1ThomasEnglehart;石天譯.墨西哥灣中部氣云四分量地震數(shù)據(jù)的解釋.石油勘探開發(fā)情報,2001(6):33432李等華,唐躍,殷積峰等.川東黃龍場構(gòu)造上二疊統(tǒng)長興組生物礁特征與潛伏礁預(yù)測.中國地質(zhì),2019,33(2):4274343CastagnaJP,SunShengjie,SiegfriedRW.Instantane2ousspectralanalysis:Detectionoflow2frequencyshad2owsassociatedw

20、ithhydrocarbons.TheLeadingEdge,2019,1201274董臣強,王軍,張金偉.時頻分析技術(shù)在三角洲層序分析中的應(yīng)用.斷塊油氣,2002,9(2):18205張宏彬,何樵登等.寬帶約束反演.石油物探,1995,34(1):1106張霖斌，姚振興等快速模擬退火算法及應(yīng)用.石油地球物理勘探,1997,32(5):6546597王赟，王妙月，彭蘇萍.地球物理隨機聯(lián)合反演.地球物理學(xué)報,1999,42(Supp.1):1421518寇曉麗，劉三陽.基于模擬退火的粒子群算法求解約束優(yōu)化問題.吉林大學(xué)學(xué)報(工學(xué)版),2019,37(1):1361405結(jié)束語通過對礁灘相油氣儲集

21、層的特征分析和研究,提出了將時頻分析技術(shù)和波阻抗反演技術(shù)相結(jié)合對礁灘相油氣藏進行預(yù)測的方法。文中分析了單獨應(yīng)用兩種技術(shù)對礁灘相油氣藏預(yù)測的優(yōu)缺點。綜合兩種技術(shù)的優(yōu)點,可以更加有效地對礁灘相油氣藏進行預(yù)測。通過對實際資料的應(yīng)用發(fā)現(xiàn),采用文中提出的礁灘相預(yù)測研究方法和思路,可以有效識別儲集層的位置,精細刻畫儲集層的范圍,計算儲集層的厚度,減少鉆井的風(fēng)險,更好地指導(dǎo)油氣開發(fā)。文中(本文編輯:馮小球)Vol.43No.6Abstracts町IVariable2gridhigh2orderfinite2differencenumericsimulationoffirst2orderelasticwave

22、equation.LiZhen2chun1,ZhangHui1andZhangHua1.OGP,2019,43(6):711716Usingvariable2gridfinite2differencealgorithmtocarryouttheseismicsimulationhasmanyu2niqueadvantagesthatmainlyaremorereasonablefordiscretizationofgeologicmodel,andrelativelyfinedivisionoflocalgridcanbeimplementedinlow2velocityzoneandcomp

23、lexstructuralarea,whichnotonlyimprovesthesimulatedprecisionandeliminatesthedispersionresultedfrominsuffi2cientsampling,butalsocanreducetherequire2mentofcomputermemory,keepingtheflexibilityofmodelcomputation.Thepaperpresentedanewelasticwavenumericsimulationmethodbasedonhigh2orderstaggered2girdtechniq

24、ue,whichreal2izedthetechniquedensifyinglocalgridbychangingthespatialstepofgrid,remedyingtshortcomingofordinaryshowedthatthethepaperischaracterizedby,improvedpreci2sionofsimulation,reducingcomputationaltimeandimprovingcomputationalefficiency.Keywords:variablegrid,staggeredgrid,high2or2derfinitediffer

25、ence,elasticwave1.CollegeofEarthResourceandInformation,Chi2naUniversityofPetroleum,DongyingCity,Shan2dongProvince,257061,ChinaElasticwavesimulationofP2andS2wavesseparationby22Dstaggeredgridandapplication.ZhangJian2lei1,TianZhen2ping1andWangCheng2xiang1.OGP,2019,43(6):717722Throughsystematicallystudy

26、ingonthemeth2odseparatingP2wavefromS2wavebysecond2orderelasticwaveequationandfirst2orderstaggeredgridalgorithm,thepaperpresentedtheformularepre2sentationsimulatingP2andS2wavesseparationinfirst2orderstaggeredgrid.TheP2andS2wavesarenaturallyseparatedbutdonotbecoupledtogetherinasimulatingprocessoftheme

27、thod,producingfullP2waveandfullS2wavesimulatedrecords.Theresultsofnumericsimulationshowedthecorrect2nessofthemethod.ThereisimportantpracticalmeaningforknowingpropagationruleofelasticwaveandguidingnumericsimulationofreservoirbystudyoffullP2waveandfullS2wavesimulation.Keywords:elasticwavesimulation,P2

28、andS2waveseparation,staggeredgrid1,TechnicalResearchCenterofGeophysicalPros2pecting,BGP,ZhuozhouCity,HebeiProvince,072751,ChinaStaggered2gridpseudo2spectrumsimulationoffrac2tureandcavereservoir.LiuJiong,MaJian2wei1,YangHui2zhu1andBaJing1.OGP,2019,43(6):723727Thepapercarriedouttheanalysisandcorrela2t

29、ionoftraditionalpseudo2spectrumandimprovedstaggered2gridpseudo2spectrumfromtheoreticalmethodcases,thesuperiorityof2method.Themoreofdifferentiallocalizationofpseudo2spectrumisshowedbyfirst2orderdifferentialcomputationalcaseofpulsese2quence;thecorrectnessofstaggered2gridpseudo2spectrummethodisprovedby

30、22Dacousticwavecomputationalcase;theprecisionofthestaggered2gridpseudo2spectrummethodbeinghigherthanthatoftraditionalpseudo2spectrumisprovedbythecomputationalcaseofelasticwaveinnon2hom2ogeneousmedium.Finally,thepropagationissueofelasticwaveinnon2homogeneousfractureandcavereservoirissimulatedbystagge

31、red2gridpseu2do2spectrummethod.Keywords:numericsimulationofwaveequation,traditionalpseudo2spectrum,staggered2gridpseu2do2spectrum,fractureandcavereservoir1.InstituteofSeismicWaveExplorationandDe2velopment,TsinghuaUniversity,BeijingCity,100084,ChinaMethodforidentifyingfluidpropertyinvolcanitereservoi

32、r.ZhangLi2hua1,PanBao2zhi1,ShanGang2yilandHeLi2.OGP,2019,43（6）:728730Thevulcanitereservoirischaracterizedbycomplexlithologyandvarietyofporestructure,whichbroughtaboutdifficultyinidentifyingande2valuatingthelithologyofvolcanitereservoir.Tak2ingdeepvolcanitereservoirinDaqingOilfieldas2019年12月第43卷第6期作者

33、介紹呂公河教授級高工,1963年生;1984年畢業(yè)于長春地質(zhì)學(xué)院物探專業(yè),1996年獲得青島海洋大學(xué)海洋地質(zhì)專業(yè)碩士學(xué)位,2019年畢業(yè)于西安交通大學(xué)工程力學(xué)專業(yè),獲得博士學(xué)位;現(xiàn)任勝利油田地球物理勘探高級專家,長期從事灘海以及高分辨率地震勘探技術(shù)研究工作。陳浩林高級工程師,1969年生;1991年畢業(yè)于中國地質(zhì)大學(xué)（武漢）地球物理專業(yè),獲學(xué)士學(xué)位;2001年畢業(yè)于中國石油大學(xué)（北京）地球物理專業(yè),獲工程碩士學(xué)位;2019年畢業(yè)于中國礦業(yè)大學(xué)（北京）地球物理專業(yè),獲博士學(xué)位。長期從事地震資料采集、處理方法及理論研究。王振華碩士,1958年生;1982年畢業(yè)于華東石油學(xué)院勘探地球物理專業(yè),198

34、5年畢業(yè)于華東石油學(xué)院北京研究生部應(yīng)用地球物理專業(yè)。主要從事地球物理勘探技術(shù)的開發(fā)與應(yīng)用?，F(xiàn)在北京華陽捷泰技術(shù)有限公司工作。單聯(lián)瑜高級工程師,1959年生;1982年畢業(yè)于華東石油學(xué)院地球物理勘探專業(yè),究和地震資料處理研究所在讀博士王玉貴高級工程師;年畢業(yè)于石油大學(xué)（華東）,2019年獲中國石油大學(xué)應(yīng)用地球物理專業(yè)碩士學(xué)位,現(xiàn)為中國地質(zhì)大學(xué)地球物理與信息技術(shù)學(xué)院在讀博士。一直從事VSP資料處理、解釋及方法研究工作。裴正林博士,副研究員,1962年生;1984年畢業(yè)于西安地質(zhì)學(xué)院物探專業(yè),2000年畢業(yè)于中國地質(zhì)大學(xué)（北京），獲博士學(xué)位,20002019年在中國石油大學(xué)（北京）從事博士后研究。主要從事地震波傳播理論及數(shù)值模擬,井間地震層析成像,小波域地震屬性提取和油氣識別,以及小波理論應(yīng)用研究。獲省部級科技進步一等獎1項,二等獎1項,三等獎2項。合著1部,發(fā)表論文50余篇。段云卿副教授,1958年生;1982年畢業(yè)于華東石油學(xué)院物探專業(yè),現(xiàn)在中國地質(zhì)大學(xué)（北京）從事教學(xué)和地震資料處理及解釋方法研究。已發(fā)表論文多篇。何兵壽副教授,1973年生;1996年畢業(yè)于中國礦業(yè)大學(xué)（徐州）物探專業(yè),2002年獲石油大學(xué)（北京）應(yīng)用地球物理專業(yè)