多波资料速度分析（spswva）

《多波资料速度分析（spswva）》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、多波资料速度分析(SPSWVA)1功能描述1.1功能概述速度参数是转换波资料处理中最重要而又最难掌握的参数之一。因此，如何才能从转换波中求取精确的合理的横波速度VS，是整个转换波处理中的重点。由于转换波传播时间中含有P波速度VP和横波速度VS两个因素，两者互相制约。使转换波速度分析比P波速度分析复杂得多。本模块的计算方法完全从转换波的时距方程出发，在事先求取到高精度的P波速度的基础上，直接从转换波中求取横波速度VS。该方法比目前广泛采用的先求VPS，再由公式求取VS的方法精度高得多。此外，本模块还有以下几方面的功

2、能：·对CSP道集中的P波求VP，对SV波求VS。·对CDP道集中的P波求VP，对SV波求VS。·从绘图仪输出等值线形式的速度谱图。·从打印机输出能量极大值的时间及速度。·为交互解释软件输出道集及能量矩阵文件。·输出的转换波速度的坐标可选择（tOPS，VS），（tOP,VS）,(tOS,VS)中任一种。1.2功能细述（1）在CSP道对P波及SV波进行速度分析。（2）在CDP道对P波及SV波进行速度分析。（3）对P波可求（top,Vprms）。（4）对SV波可求（tops,Vsrms）,(top,Vsrms),(t

3、os,Vsrms)及（tops,Vps）。（5）打印机输出极大值时间及速度。（6）绘图仪输出等值线形式的速度谱。1.3数学与物理方法速度参数是转换波资料处理中最重要而又最难掌握的参数之一。因此，如何才能从转换波中求取精确、合理的横波速度VS，是整个转换波处理的重点。由于转换波传播时间中含有P波速度VP和横波速度VS两个因素，且两者互相制约，使转换波速度分析比P波速度分析复杂得多。本节讨论的计算方法完全从转换波的精确时距方程出发，在事先求取到高精度P波速度的基础上，直接从转换波中求取横波速度VS。该方法的精度比采用

4、近似公式先求VPS，再由公式求取VS的方法高得多。此外，本方法不仅可在CCL道集应用，也可推广到CSP道集速度分析。由转换波运动学原理的讨论已知，转换波时距方程为(1)式中，纵波垂直往返时间；横波垂直往返时间；--纵波均方根速度；9--横波均方根速度；H--界面深度；XP--转换点到炮点的距离；XS=X-XP，X为炮检距，XS为转换点到接收点的距离。由(2)有(3)式中，为转换波垂直往返时间，即纵波垂直入射到界面，转换为横波再垂直反射到地面的时间，为纵横波速度比。从而，可用下面的动校正量公式求取速度：(4)此外，

5、还可采取以下技术措施提高速度的精度：(1)采用无畸变动校正方法，使速度精度和分辨率不受动校拉伸畸变的影响。(2)速度谱能量矩阵采用相似性加权迭加：(5)式中，(6)为振幅值，N为总道数，)为叠加道，为扫描速度。(3)用二次曲线求极值的方法寻找能量团的极值点，提高速度精度。92输入/输出数据说明2.1输入数据a.GRISYS格式炮集或CDP道集记录。b.参考速度文件。2.2输出数据a.绘图仪输出等值线形式的速度谱图。b.输出内部格式速度谱能量矩阵及CDP道集文件（交互用）。c.输出自动检测的极大值速度文件。3外部接

6、口描述3.1地震子系统子程序AMINIT，GOSPAR，GETHED，GOSCAR，GOSFIL，ERRMSI，ERRMSS，GOSABT，WTRACE1，RTRACE1，RESFLG，GOSDFL，CHAFLG3.2地震数据库子程序3.3科学库子程序SIN，COS，EXP，LN3.4绘图子程序INITPL，PLOT，SYSMBLL，NUMBER，LINEWT4．模块结构描述译码及参数准备输入一炮或几个CDP道集炮道集？调CDP速度分析子程序炮集速度分析调绘图子程序绘速度谱速度点完否？输出速度谱及结果NYNY94

7、．1主程序与子程序调用关系SPSWVA→SPSWV01→SPSWV02→SPSWV034．2程序逻辑结构描述（1）SPSWVA子程序完成译码，准备参数并调用SPSWV01子程序。（2）SPSWV01子程序完成炮集多波速度分析及调用SPSWV02、SPSWV03子程序。（3）SPSWV02子程序完成CDP道集多波速度分析。（4）SPSWV03子程序完成绘制等值线形式的速度谱。4.3编码参数设计（1）控制卡列号区位内容1标识：3-8名字SPSWVA10-80第一选件空或SV,PY空——P波、SH波速度分析SV——P-

8、SV波速度分析PY——P-SH波速度分析第二选件空或CP空——CSP道集速度分析CP——CDP道集速度分析（2）参数（10-80列）DVn:扫描速度间隔缺省：n=20DTOn:扫描t0间隔（样点个数）缺省：n=2VWn:扫描速度范围（m/s）缺省：n=500LRn:Vs/Vp低限（×100）缺省：n=50GRn:Vs/Vp高限（×100）缺省：n=90WLn:子波长度（m