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
- processing system
- data processing
- wave field
- fixed
- 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
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 |
|---|---|---|
| CA2821482C (en) | Systems and methods for two-dimensional domain decomposition during parallel reservoir simulation | |
| CN111650579B (zh) | 一种岩移参数自适应获取的InSAR矿区三维形变估计方法、装置及介质 | |
| Zlatev et al. | Modeling the long-range transport of air pollutants | |
| Plesset et al. | Cavity drag in two and three dimensions | |
| Ammar et al. | Numerical experience with a superfast real Toeplitz solver | |
| CN103699798A (zh) | 一种实现地震波场数值模拟方法 | |
| Brockett | Classification and equivalence in estimation theory | |
| CN110857998A (zh) | 一种基于lowrank有限差分的弹性逆时偏移方法及系统 | |
| Huntington | Frequency distribution of product and quotient | |
| CN1255642A (zh) | 定点lu分解法用于三维地震偏移计算的地震数据处理系统 | |
| Reznik et al. | On the probability of the formation of local groups in random point images | |
| CN109490954B (zh) | 波场正演模拟方法及装置 | |
| CN102411773A (zh) | 面向向量处理器的去均值归一化积相关系数的向量化实现方法 | |
| CN1255643A (zh) | 块追赶法用于三维地震偏移计算的地震数据处理系统 | |
| CN108981673B (zh) | 利用序列星载光学遥感图像计算地面点三维位置的方法 | |
| CN102608658A (zh) | 强能量干扰抑制方法 | |
| CN102831633A (zh) | 一种基于线性纹理滤波的卷积滤波优化方法 | |
| Kytmanov | An algorithm for calculating power sums of roots for a class of systems of nonlinear equations | |
| Dasarathy et al. | New sample complexity bounds for phylogenetic inference from multiple loci | |
| CN107480338A (zh) | 基于时空积分域处理时域边界元奇异积分的方法 | |
| Liu et al. | Uniform $ L_1 $ boundedness of solutions of hyperbolic conservation laws | |
| Lawrie et al. | Fine-grain parallel processing implementations of Kalman filter algorithms | |
| Harris et al. | On approximate integration of Wiener and Poisson process driven stochastic differential equations | |
| Chi | Performance of the SMLR deconvolution algorithm | |
| KR20080068167A (ko) | 윈도우 행렬 분리 연산을 통한 에스에이디 정합 알고리즘기반 스테레오 비젼 시스템 설계 방법 |
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 |