CN1255642A - 定点lu分解法用于三维地震偏移计算的地震数据处理系统 - Google Patents
定点lu分解法用于三维地震偏移计算的地震数据处理系统 Download PDFInfo
- Publication number
- CN1255642A CN1255642A CN 99122562 CN99122562A CN1255642A CN 1255642 A CN1255642 A CN 1255642A CN 99122562 CN99122562 CN 99122562 CN 99122562 A CN99122562 A CN 99122562A CN 1255642 A CN1255642 A CN 1255642A
- Authority
- CN
- China
- Prior art keywords
- earthquake
- calculation
- data processing
- processing system
- wave field
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000004364 calculation method Methods 0.000 title claims description 25
- 238000000034 method Methods 0.000 title claims description 24
- 238000012545 processing Methods 0.000 title claims description 9
- 238000000354 decomposition reaction Methods 0.000 claims description 16
- 238000013213 extrapolation Methods 0.000 claims description 16
- 230000005012 migration Effects 0.000 claims description 11
- 238000013508 migration Methods 0.000 claims description 11
- 238000005516 engineering process Methods 0.000 claims description 5
- 239000011159 matrix material Substances 0.000 description 15
- 230000008569 process Effects 0.000 description 4
- 230000009467 reduction Effects 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000000739 chaotic effect Effects 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种定点LU分解法用于“三维地震偏移”计算的地震数据处理系统,包括该系统中的“三维地震偏移”操作程序,实施该操作采用“波场外推”技术,其特征在于“波场外推”的计算方法为定点LU分解逼近法。本发明实现后,可使其“三维地震偏移”操作的计算过程更加方便、快捷、并提高运算精度,从而进一步提高地震勘探工作的总体质量水平。
Description
本发明涉及地震勘探数据处理技术，特别是“三维地震偏移”的定点LU分解逼近技术。
地震勘探工作主要由两大工作系统组成，一是地震数据采集系统，二是地震数据处理系统。其中，地震数据处理系统包括大型运算设备、数据处理程序、适用于相应程序的运算方法等。在数据处理程序中，一个重要的组成部分是“三维地震偏移”，即通过一定的技术手段来最大限度地消除在地面上测试得到的用来描述工区地质图像的三维地震数据与工区实际地质结构之间的偏差，以便得到能够更为准确地反映工区实际地质结构的勘探结果图像。到目前为止，用来实现“三维地震偏移”操作的主要途径是“波场外推”技术，即通过一定的运算方法，使得在地面上取得的地震波场数据，转化成地下震点处波场数据，再利用该数据制得反映上述工区实际地质结构情况的勘探结果图像。现有的利用“波场外推”技术来实现“三维地震偏移”操作的计算方法较多，主要是“ADI交替方向隐式法”和尚未推广应用的“降维LU分解逼近法”。前者的缺欠是计算误差明显地与地下反射层倾角、方位及速度变化有关，这不仅增加了运算难度，也使波场外推算子的圆对称性不能保证，因而计算速度与计算精度都不够理想；后者对“波场外推”方程的算子矩阵采用非线性的方法作LU分解、而后作降维处理，即把远离对角线的元素取为零，再经两次递归多项式除法后完成求解计算。这种方法的局限性首先在于引入非线性方法后，使计算过程出现混沌状态，造成计算困难。其次，到目前为止，还没有足够的理论依据来证明该方法能够保证外推算子的圆对称性，缺少定量分析，计算的可靠性尚未得到证实。
本发明的目的是为人们提供一种采用定点LU分解逼近法实施“三维地震偏移”操作计算任务的地震数据处理系统，以便使“三维地震偏移”操作的计算过程更加方便、快捷，并提高运算精度，从而进一步提高地震勘探工作的总体质量水平。
本发明的目的是这样实现的，它包括地震数据处理系统中的“三维地震偏移”操作程序，实施该操作采用“波场外推”技术，其特征在于“波场外推”的计算方法为定点LU分解逼近法。
本发明实施后与现有技术相比具有以下优点：1、相对于现有的降维LU分解逼近法而言，由于其运算中不引入非线性方法，因而运算过程筒单、明了、易于操作。其次，能较好地保持波场外推算子的圆对称性，提高波场外推算子的逼近精度。再次，能够取得很高的运算效率。2、与现有的ADI方法相比，由于其运算过程与反射层倾角、方位、及其速度变化等复杂的附加因素无关，因而，一方面可大幅度提高运算速度，另一方面也有利于保证波场外推算子的圆对称性，进而提高运算精度。
下面结合实施例对本发明做进一步详细地说明。
利用本发明提供的地震数据处理系统，进行“地震偏移”操作时，也要通过“波场外推”途径，这与已有技术是相同的，而其“波场外推”的计算方法，则采用定点LU分解逼近法，这与现有的降维LU分解逼近方法有着明显的区别，这种区别在于它不引入非线性。定点LU分解逼近的计算步骤及其算式如下：
波场外推方程的系数矩阵为
其中，
k＝1，2...，Nx；这里，Nx为每个子块矩阵的行数，即拉直前数据矩阵的行数，Ny为每个子块矩阵的列数即拉直前数据矩阵的列数；矩阵1的行数和列数均为N＝Nx×Ny；矩阵1的阶数为(Nx×Ny)×((Nx×Ny))。我们规定
A＝R+LLT (4)式4中R是块不完全cholesky分解(BIC)的剩余矩阵，矩阵L是个块两对角型的下三角阵
L的块对角线上的每个子块矩阵为两对角型
L的块超对角线上的每个子块矩阵为对角阵
由A-8表达的就是波场外推方程的系数矩阵的定点的LU分解算法。
该分解得到的因子矩阵是特别稀疏的，而且分解式的计算并不涉及非线性或者引入边界条件之类的人为的假设，因此，计算效率及精度均很高。
Claims (1)
1、定点LU分解法用于三维地震偏移计算的地震数据处理系统，包括该系统中的“三维地震偏移”操作程序，实施该操作采用“波场外推”技术，其特征在于“波场外推”的计算方法为定点LU分解逼近法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 99122562 CN1255642A (zh) | 1999-11-23 | 1999-11-23 | 定点lu分解法用于三维地震偏移计算的地震数据处理系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 99122562 CN1255642A (zh) | 1999-11-23 | 1999-11-23 | 定点lu分解法用于三维地震偏移计算的地震数据处理系统 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN1255642A true CN1255642A (zh) | 2000-06-07 |
Family
ID=5282559
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 99122562 Pending CN1255642A (zh) | 1999-11-23 | 1999-11-23 | 定点lu分解法用于三维地震偏移计算的地震数据处理系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN1255642A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100344993C (zh) * | 2003-10-23 | 2007-10-24 | Pgs美洲公司 | 叠加前地震数据的向下外推方法 |
-
1999
- 1999-11-23 CN CN 99122562 patent/CN1255642A/zh active Pending
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN100344993C (zh) * | 2003-10-23 | 2007-10-24 | Pgs美洲公司 | 叠加前地震数据的向下外推方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114117906A (zh) | 基于观测数据自编码的多尺度无监督地震波速反演方法 | |
| Michael et al. | Coseismic stress changes induced by the 1989 Loma Prieta, California earthquake | |
| CN108181653B (zh) | 针对vti介质逆时偏移方法、设备及介质 | |
| CN110857998B (zh) | 一种基于lowrank有限差分的弹性逆时偏移方法及系统 | |
| CN112632837B (zh) | 一种地下结构纵向抗震数值确定方法 | |
| CN101950034B (zh) | 一种用于地震反演的自然电位曲线基线漂移校正方法 | |
| CN105139339A (zh) | 基于多级滤波和样本匹配的偏振图像超分辨率重建方法 | |
| Tang et al. | Generalized pulse-spectrum technique for 2-D and 2-phase history matching | |
| CN106959467A (zh) | 地震波速度反演方法和装置 | |
| Nasiri et al. | Seismic response analysis of embankment dams under decomposed earthquakes | |
| CN111123374A (zh) | 一种基于匹配滤波的探地雷达全波形反演方法 | |
| CN118311653A (zh) | 基于矩阵信息降维的线性卡尔曼地震随机反演方法及系统 | |
| CN1255642A (zh) | 定点lu分解法用于三维地震偏移计算的地震数据处理系统 | |
| CN111175816A (zh) | 油藏改造实时构建微地震裂缝网络的方法及装置 | |
| CN117706623A (zh) | 一种深度学习和波动方程联合驱动的微地震速度反演方法 | |
| CN113568033B (zh) | 一种三维非规则采样地震采集观测系统设计方法和装置 | |
| CN100349005C (zh) | 地震数据多炮组合叠前深度偏移方法 | |
| CN118941246A (zh) | 一种施工现场预制梁钢筋生产管理系统 | |
| Farmer et al. | A computer-aided method for calculating the distributions of strain-rate and strain from an experimental flow field | |
| Liu et al. | A generalized pulse-spectrum technique (GPST) for determining time-dependent coefficients of one-dimensional diffusion equations | |
| Augarde et al. | Numerical modelling of tunnel installation procedures | |
| CN114428307A (zh) | 一种地震多数据联合的三维速度建模方法及系统 | |
| CN117272731B (zh) | 基于卷积核分解-Winograd快速卷积算法的地震波场模拟方法 | |
| Sandhu | Finite element analysis of coupled deformation and fluid flow in porous media | |
| CN120742414B (zh) | 一种基于点扩散函数校正的隧道超前探测方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C12 | Rejection of a patent application after its publication | ||
| RJ01 | Rejection of invention patent application after publication |