【技术实现步骤摘要】
本专利技术涉及地震数据处理。从一方面来说，本专利技术涉及一种从多个地震波形曲线中产生增强地震显示的方法。绘制地下地质图的地震方法包括地震能源的运用和通过地震检测器(一般称为地震检波器)阵列对它的接收。当在陆地上运用时，该地震能源一般是用电起爆的高爆炸力的炸药，炸药放在孔洞里，孔洞位于地面网格的一个选定格点上，或是一种能够给地面以一系列冲击的能源，比如在Vibroseis振动机中所用的。通过爆破和冲击在地里产生的声波从地层界面传回来，按照穿过的距离和地下的特性以不同的时间达到地面。用地震检波器检测这些返回声波，把这些声波转换成代表它的电信号。把许多地震检波器按照选定的方式排列成阵列，以便最有效地检测这返回的声波，并产生代表它们的电信号，进而从这些信号中可以推导出有关地下地质结构的数据。一般把这些电信号叫做地震波形曲线。图1表示了地震源和检测器的一种典型阵列。把一组地震源(21-25号)和一组地震接收器(31-35号)围绕一个共同的中点来排列。为了便于演示，这里把地震源和地震接收器都等距离地放置，但并不一定要求如此。在地震勘探中，由于障碍物和其它一些因素的影响，通常都不...
【技术保护点】
一种从多个真实地震波形曲线中产生可变振幅／相位地震显示的方法，该方法包括下列步骤：（ａ）对于所述多个真实地震波形曲线的每一个，选择多个数据点，其中所述的多个数据点间距相等；（ｂ）确定在步骤（ａ）中所选择的多个数据点中每个点的真实振幅 ；（ｃ）在步骤（ａ）中所选择的每个数据点上，确定对应于所述多个真实地震波形曲线的多个复地震波形曲线的瞬时相位；（ｄ）根据步骤（ｃ）中在每一数据点所确定的瞬时相位，对步骤（ａ）中所选择的每个数据点赋予彩色一、彩色二、彩色三或彩色四，其 中所述的彩色是根据表Ⅰ指定的，表中的ＢＡ是相角基值。表Ⅰ彩色、范围一、（ＢＡ—９０°...
【技术特征摘要】
所述，所有这些变化都在本发明范围之内。用于可变振幅-相位绘图的源码〔按每秒采样(SPS)〕DIMENSION TARRAY〔4096〕，BUFF〔5000〕，ENV〔4096〕，SLOC1〔2〕，SLOC2〔2〕DIMENSION NUMTRC〔3000〕COMMON/CTRL/NX，NY，NTREC，DELT，NT，TPI，SPI，NREC，SCALE，LINC，ZMAXCOMMON/PARM/XORGIN，YORGIN，IORGIN，PLOTX，PLOTY，XINC，YINC，SIZE，SIZE1COMMON/CTRL2/IFLAG，IXFRST，IXLINE，IP，TAR，IN，NN，IPLTECOMMON/CTRL3/IOP，IOPT1，IOPT2，IOPT3，IOPT4，IOPT5COMMON/CTRL4/ICDPS，ICDPE，ITST，ITED，NSAMP，ICDPNINTEGER＊4 JBUFF〔5000〕，MLINE〔3〕，ICOPON，INKSON，IBLANKINTEGER＊4 IDIR，IDIRR，IDIRL，DDNAME〔2〕，ILOC〔2〕，ILOC2〔2〕INTEGER＊2 IBUFF〔10000〕，ITRHE，IHEAD，IPCODE，IVOL〔3〕，IPROS〔4〕EQUIVALENCE〔BUFF〔1〕，IBUFF〔1〕，JBUFF〔1〕〕EQUIVALENCE〔ILOC〔1〕，SLOC1〔1〕〕，〔ILOC2〔1〕，SLOC2〔1〕〕DATA IDIR/‘ ’/，IDIRR/‘RTOL’/，IDIRL/‘LTOR’/DATA INKSON/‘ ’/，ICDPON/‘ ’/，IBLANK/‘ ’/DATA IHEAD/‘TP’/DATA ITRHE/‘HE’/DATA IPCODE/‘ ’/DATA DDNAME/‘FT29’，‘F001’/WRITE〔6，1075〕检验日时间，并设置间隔定时器CALL TIMINT〔KK〕CALL TIME〔SLOC1〕NX＝0READ〔10，1000〕TPI，SPI，ICDPS，ICDPE，ITST，ITED，IPCODEREAD〔10，1010〕IDIR，ICDPON，INKSON，ZMAXWRITE〔6，1080〕WRITE〔6，1055〕TPI，SPI，ICDPS，ICDPE，ITST，ITED，IDIR，IPCODEIF〔ZMAX.GT.0.0〕IOPT1＝1IF〔INKSON.NE.IBLANK〕IOPT2＝1IF〔ICDPON.NE.IBLANK〕IOPT3＝1IF〔IDIR.EQ.IDIRR〕IOPT4＝1LN＝020 CALL TAPEIO〔‘TAPEIN’，BUFF，LN，‘R’〕IF〔LN.LT.-1〕GO TO 20IF〔LN.EQ.-1〕GO TO 400IF〔IBUFF〔1〕.NE.IHEAD〕GO TO 100IPROS〔1〕＝IBUFF〔5〕IPROS〔2〕＝IBUFF〔6〕IPROS〔3〕＝IBUFF〔7〕IPROS〔4〕＝IBUFF〔8〕MLINE〔1〕＝JBUFF〔5〕MLINE〔2〕＝JBUFF〔6〕MLINE〔3〕＝JBUFF〔7〕IRUN＝IBUFF〔39〕WRITE〔6，1050〕IPROS，MLINE，IRUNWRITE〔6，1060〕NT＝IBUFF〔128〕LT＝1024IF〔NT.GT.1024.AND.NT.LE.2048〕LT＝204812IF〔NT.GT.2048.AND.NT.LE.4096〕LT＝4096IDELT＝IBUFF〔127〕NREC＝IBUFF〔164〕DELT＝IDELTNREC＝ICDPE-ICDPS+1CALL LIMIT100 IF〔IBUFF〔1〕.NE.ITRHE〕GO TO 200DO 50 J＝1，409650 TARRAY〔J〕＝0.0ICDPN＝IBUFF〔63〕IF〔NX.EQ.NREC〕GO TO 400IF〔IBUFF〔6〕.NE.IPCODE〕GO TO 200IF〔ICDPN.LT.ICDPS.OR.ICDPN.GT.ICDPE〕GO TO 200IF〔IBUFF〔6〕.EQ.IPCODE〕NX＝NX+1NTRMP＝IBUFF〔28〕IART＝IBUFF〔4〕WRITE〔6，9998〕DO 250 J＝1，NTRMPN＝IART+JTARRAY〔J〕＝BUFF〔N〕250 CONTINUECALL ENVELP〔TARRAY，LT，ENV〕I＝ITSTK＝ITEDWRITE〔9〕〔TARRAY〔N〕，N＝I，K〕WRITE〔8〕〔ENV〔N〕，N＝I，K〕NUMTRC〔NX〕＝ICDPNWRITE〔6，3004〕ICDPN，NTRMP200 LN＝0GO TO 20绘截面400 WRITE〔6，1065〕CALL IPARM〔IERROR〕IF〔IERROR.EQ.1〕GO TO 5000IF〔IOPT4.EQ.1〕XORGIN＝830CALL IPAGBCALL IPRINT〔1〕CALL ISTART〔IORGIN〕CALL ICHCOL〔1〕CALL APPLECALL TGRID〔NUMTRC〕CALL TITLEWRITE〔6，2040〕CALL ICLB〔0〕CALL TIMINT〔JJ〕CALL TIME〔SLOC2〕INTIME＝KK-JJWRITE〔6，1085〕ILOC，ILOC2，INTIME恢复VDLSERCALL DDCK〔DDNAME，IFG，IVOL，IDEN〕WRITE〔6，1095〕IVOLSTOP5000 WRITE〔6，2000〕STOP1000 FORMAT〔10X，2F10.3，4I5，3 A2〕1010 FORMAT〔16X，A4，6X，A4，6X，A4，F10.2〕1050 FORMAT〔40X，‘FEATURE’，3X，4A2，/40X，‘LINE’，10X，3A4，/，40X，‘RUN＃’，5X，I＊2，///〕1055 FORMAT〔40X，‘TPI’，F11.3，/，40X，‘IPS’，F11.3，/，40X，‘FIRST CDP’，I5，/，40＊X，‘LAST CDP’，I5，/，40X，‘START TIME’，I4，/，40X，‘END TIME’，I4，/，40X＊，‘PLOT DIR’，A4，/，40X，‘PC CODE’，A2〕1060 FORMAT〔///，35X，‘＊＊＊ TRACES PROCESSED ＊＊＊’，///〕1065 FORMAT〔///，35X，‘＊＊＊ PRACE PROCESSING COMPLETED ＊＊＊’，//////〕1075 FORMAT〔1H1，///，100〔‘＊’〕，/////，35X，‘THIS IS THE VARIABLE AMPLITUDE＊PLOTTING PACKAGE FOR SPS’，////，35X，‘FOR PROBLEMS CALL CLIFF SHOCK＊AT X-4559’，/////////〕1080 FORMAT〔35X，‘＊＊＊ INPUT PARAMETERS ＊＊＊＊＊’，///〕1085 FORMAT〔//，35X，‘START TIME’，2A4，//，36X，‘END TIME，’，2A4，＊//，36X，‘CPU IN SECONDS’，I10，//〕1090 FORMAT〔/////////，35X，‘＊＊＊＊＊＊＊＊＊ THIS IS END OF JOB ＊＊＊＊＊＊＊＊＊’〕1095 FORMAT〔35X，‘OUTPUT VOLUME ＊＊’，3A2，‘FILE ＃2 ＊＊’〕2000 FORMAT〔35X，‘＃＃＃＃＃＃ PROGRAM TERMINATED DUE TO CRITICAL ERROR ＃＃＃＃’〕2040 FORMAT〔35X，‘＊＊＊＊ OUTPUT INFORMATION ＊＊＊＊’，/////〕3004 FORMAT〔25X，‘＊＊＊＊ CDP’，I4，‘WITH’，I4，‘SAMPLES IS SUCCESSFULLY PRO＊CESSED ＊＊＊＊＊’〕STOPENDSUBROUTINE APPLECOMMON /CTRL/NX，NY，NTREC，DELT，NT，TPI，SPI，NREC，SCALE，LINC，ZMAXCOMMON /PARM/XORGIN，YORGIN，IORGIN，PLOTX，PLOTY，XINC，YINC，SIZE，SIZE1COMMON /CTRL3/IOP，IOPT1，IOPT2，IOPT3，IOPT4，IOPT5COMMON/CTRL4/ICDPS，ICDPE，ITST，ITED，NSAMP，ICDPNDIMENSION SBUFF〔3000〕，TBUFF〔3000〕，PBUFF〔3000〕，ABUFF〔3000〕DIMENSION U〔4〕，AIN〔100000〕，ZCLASS〔16〕，PIN〔100000〕，IBUFF〔3000〕DIMENSION ISET〔64〕，PCLASS〔5〕DATA PCLASS/-3.1459，-1.57，0.0，1.57，3.1459/DATA ISET/10001，20002，30003，40004，50005，60006，70007，80008，90009，10＊0010，110011，120012，130013，140014，150015，160016，1，2，3，4，5，6，7，8，9，1＊0，11，12，13，14，15，16，100，200，300，400，500，600，700，800，900，1000，1100，＊1200，1300，1400，1500，1600，10100，20200，30300，40400，50500，60600，70700＊，80800，90900，101000，111100，121200，131300，141400，151500，161600/DATA ZCLASS/0.0，400.，800.，1200.，1600.，2000.，2400.，2800.，3200.，3600＊.，4000.，4400.，4800.，5200.，5600.，6000./IF〔IOPT1.EQ.1〕CALL NEWSCL〔ZCLASS，ZMAX〕NSIGN＝1IF〔IOPT4.EQ.1〕NSIGN＝-1NCLZ＝16NCLZ1＝NCLZ-1REWIND 8REWIND 9置绘图坐标XSTART＝YORGINYST＝XORGIN开始内插计算XPIX＝PLOTY/.20+1.5NPIXX＝IFIX〔XPIX〕U〔3〕＝.5NPTS＝NSAMP+1NSAVE＝0NPASS＝0NINT＝100000/NPIXX10 LNT＝NINTU〔2〕＝1./〔〔XPIX-1.5〕/〔NSAMP-1〕〕U〔4〕＝FLOAT〔NPIXX〕NROW＝0NPTS＝NSAMP+1IF〔LNT.GT.〔NX-NSAVE〕〕LNT＝NX-NSAVET方向内插DO 40 I＝1，LNTNSAVE＝NSAVE+1NROW＝NROW+1READ〔8〕〔PBUFF〔K〕，K＝1，NSAMP〕READ〔9〕〔ABUFF〔K〕，K＝1，NSAMP〕PBUFF〔NPTS〕＝PBUFF〔NSAMP〕ABUFF〔NPTS〕＝ABUFF〔NSAMP〕CALL UNWRAP〔PBUFF，NPTS〕CALL APAM〔‘QIN＊’，1，TBUFF，NPIXX，4，0，PBUFF，NPTS，4，0，U，4，4，0〕CALL APAM〔‘QIN＊’，1，SBUFF，NPIXX，4，0，ABUFF，NPTS，4，0，U，4，4，0〕CALL WRAP〔TBUFF，NPIXX〕DO 30 J＝1，NPIXXN＝〔NROW-1〕＊NPIXX+JPIN〔N〕＝TBUFF〔J〕AIN〔N〕＝SBUFF〔J〕30 CONTINUE检验结束的T方向40 CONTINUE开始X方向内插YPIX＝〔NROW-1〕＊XINCYPIX＝YPIX/.2+1.5NPIXY＝IFIX〔YPIX〕U〔2〕＝1./〔〔YPIX-1.5〕/〔NROW-1〕〕U〔4〕＝FLOAT〔NPIXY〕NPTS＝NROW+1开始X内插循环DO 200 J＝1，NPIXXDO 100 I＝1，NROWSBUFF〔I〕＝AIN〔J+〔I-1〕＊NPIXX〕TBUFF〔I〕＝PIN〔J+〔I-1〕＊NPIXX〕100 CONTINUESBUFF〔NPTS〕＝SBUFF〔NROW〕TBUFF〔NPTS〕＝TBUFF〔NROW〕CALL UNWRAP〔TBUFF，NPTS〕CALL APAM〔‘QIN＊’，1，PBUFF，NPIXY，4，0，TBUFF，NPTS，4，0，U，4，4，0〕CALL APAM〔‘QIN＊’，1，ABUFF，NPIXY，4，0，SBUFF，NPTS，4，0，U，4，4，0〕CALL WRAP〔PBUFF，NPIXY〕开始绘X扫描线15A＝PBUFF〔KT〕B＝ABS〔ABUFF〔KT〕〕IF〔A.LT.PCLASS〔1〕〕II＝1IF〔A.GT.PCLASS〔5〕〕II＝4DO 59 I＝1，4IF〔A.GT.PCLASS〔I〕.AND.A.LE.PCLASS〔I+1〕〕II＝I59 CONTINUEIF〔B.GT.ZCLASS〔NCLZ〕〕ICODE＝ISET〔〔II-1〕＊16+16〕IF〔B.LT.ZCLASS〔1〕〕ICODE＝ISET〔〔II-1〕＊16+1〕DO 31 M＝1，NCLZ1IF〔B.GT.ZCLASS〔M〕.AND.B.LE.ZCLASS〔M+1〕〕ICODE＝ISET〔〔II-1〕＊16+M〕31 CONTINUEIF〔B.EQ.0.0〕ICODE＝0IBUFF〔KT〕＝ICODEIF〔ICODE.EQ.0〕KOUNT＝KOUNT+151 CONTINUEDO 71 KT＝1，NPIXYYB＝YST+〔KT-1〕＊.2＊NSIGNYSIDE＝.2CALL IRECT〔YB，XB，YSIDE，.20，IBUFF〔KT〕，0〕71 CONTINUE200 CONTINUE检查、看看是否整个时间片已结束NPASS＝NPASS+1IF〔NSAVE.GE.NX〕GO TO 700YST＝YST+NPIXY＊.2＊NSIGNBACKSPACE 8BACKSPACE 9NSAVE＝NSAVE-1GO TO 10绘图已 完，关闭绘图文件700 CONTINUEXSTART＝XORGIN+PLOTX+25.IF〔IOPT4.EQ.1〕XSTART＝XORGIN-PLOTX-110.IF〔IOPT2.NE.1〕CALL INKS〔XSTART，ISET，PCLASS，ZCLASS〕RETURNENDSUBROUTINE UNWRAP〔SBUFF，K〕DIMENSION SBUFF〔1〕，WRAPN〔8000〕PI＝3.14159TWOPI＝2.＊PINPTS＝K+1WRAPN〔1〕＝0.0DO 100 I＝1，KDIFF＝SBUFF〔I〕-SBUFF〔I+1〕WRAPN〔I+1〕＝WRAPN〔I〕IF〔DIFF.GT.PI〕WRAPN〔I+1〕＝WRAPN〔I〕+1IF〔DIFF.LT.〔-PI〕〕WRAPN〔I+1〕＝WRAPN〔I〕-1100 CONTINUEDO 200 J＝2，NPTSSBUFF〔J〕＝SBUFF〔J〕+WRAPN〔J〕＊TWOPI200 CONTINUERETURN ’/设置波形图例子程序CALL DATE〔IDATE〕JDATE〔1〕＝IDATE〔1〕JDATE〔3〕＝IDATE〔2〕JDATE〔5〕＝IDATE〔3〕C＝ZCLASS〔16〕D＝ZCLASS〔16〕-ZCLASS〔15〕LT＝1024 100 CONTINUEDO 200 J＝1，360TARRAY〔J〕＝TARRAY〔J〕＊〔C+D＊2〕200 CONTINUEDO 300 J＝361，540TARRAY〔J〕＝TARRAY〔J〕＊〔C-D＊3〕300 CONTINUEDO 400 J＝541，720TARRAY〔J〕＝TARRAY〔J〕＊〔C-D＊9〕400 CONTINUEDO 450 J＝1，90ENV〔J〕＝J-90.ENV〔J+90〕＝JENV〔J+180〕＝90.+JENV〔J+270〕＝J-180.ENV〔J+360〕＝J-90.ENV〔J+450〕＝JENV〔J+540〕＝90.+JENV〔J+630〕＝J-180.450 CONTINUEDO 475 J＝1，720ENV〔J〕＝ENV〔J〕/180.＊3.14159475 CONTINUEX1＝XSTART+10.Y＝256.KOUNT＝0DO 502 J＝1，720，4A＝ENV〔J〕B＝ABS〔TARRAY〔J〕〕IF〔A.LT.PCLASS〔1〕〕II＝1IF〔A.GT.PCLASS〔5〕〕II＝4DO 501 I＝1，4IF〔A.GT.PCLASS〔I〕.AND.A.LE.PCLASS〔I+1〕〕II＝I501 CONTINUEIF〔B.GT.ZCLASS〔16〕〕ICODE＝ISET〔〔II-1〕＊16+16〕IF〔B.LT.ZCLASS〔1〕〕ICODE＝ISET〔〔II-1〕＊16+1〕DO 503 K＝1，15IF〔B.GT.ZCLASS〔K〕.AND.B.LT.ZCLASS〔K+1〕〕ICODE＝ISET〔〔II-1〕＊16+K〕503 CONTINUEKOUNT＝KOUNT+1Y1＝Y-〔KOUNT-1〕＊.2CALL IRECT〔X1，Y1，20.，.2，ICODE，0〕502 CONTINUEX1＝XSTART+20.KOUNT＝0DO 600 I＝1，720，4KOUNT＝KOUNT+1YP〔KOUNT〕＝Y-〔KOUNT-1〕＊.2XP〔KOUNT〕＝TARRAY〔I〕/〔C+D＊2〕＊9.+X1600 CONTINUE18CALL ILINE〔X1，256.，.2，7，1〕CALL ILINE〔X1，220.，.2，7，0〕X1＝XSTART+30.CALL ITEXT〔‘+ ’，X1，258.，2.，0，7〕X1＝XSTART+10.CALL ITEXT〔‘- ’，X1，258.，2.，0，7〕X1＝XSTART+20.CALL ITEXT〔‘0 ’，X1，258.，2.，0，7〕X1＝XSTART+10.CALL ITEXT〔JDATE，X1，280.，2.，0，7〕X1＝XSTART-15.CALL ITEXT〔‘PHASE VARIABLE AMPLITUDE ’，X1，270.，3.，0，7〕RETURNENDSUBROUTINE TGRID〔NUMTRC〕COMMON /CTRL/NX，NY，NTREC，DELT，NT，TPI，SPI，NREC，SCALE，LINC，ZMAXCOMMON /PARM/XORGIN，YORGIN，IORGIN，PLOTX，PLOTY，XINC，YINC，SIZE，SIZE1COMMON /CTRL3/ IOP，IOPT1，IOPT2，IOPT3，IOPT4，IOPT5COMMON /C...
【专利技术属性】
技术研发人员：克利佛德兰肖克，科林约翰森，
申请(专利权)人：菲利普石油公司，
类型：发明
国别省市：US[美国]
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。