共查询到19条相似文献，搜索用时 218 毫秒
1.
2.
3.
《地球物理学进展》2017,(3)
为了解决三维斜井井间地震数据成像问题,借鉴二维地面地震弯线及宽线处理思路,将三维井间地震成像问题转换成二维成像问题进行研究.利用三维逐段迭代方法对三维斜井井间地质模型进行正演模拟,计算上、下行反射波走时及反射点位置,通过对比正演记录与实际资料的特征,调整地质模型参数获取准确速度场.然后对波场分离后的共炮集数据进行反射波归位,得到共反射点道集数据.最后根据反射点的分布选择一条反射点分布最集中的剖面线,以一定的宽度将剖面线左右的反射点包裹起来;根据叠加宽度划分成大小相等的矩形面元,将面元内的共反射点数据投影到剖面线上进行二维叠加成像.通过模型测试,三维XSP-CDP叠加方法能有效反演地下结构特征,指导三维井间地震采集观测系统的设计并解决不共面斜井叠加成像问题.应用该方法对M地区实际资料进行成像研究,结果表明该方法可有效反演地下三维地质结构特征,为进一步数据处理提供依据. 相似文献
4.
为进一步了解井间地下介质结构及其电性特征,本文用时域有限差分方法(FDTD),实现了井间电磁场响应的数值模拟.激发源是探测深度相对磁偶极子更大一些的电偶极子,首先选择一个位于井轴上的垂直电偶极子作为场源,并假设参与计算的介质相对于发射井井轴是轴向对称的,这样可将研究区域作为二维问题处理.推导了二维井间电磁波传播时域差分公式.该方法适用于任何方向入射的电偶极子源,尤其善于解决频域差分方法所难以描述的宽频脉冲.给出了井间金属圆柱和矿体圆柱数值模拟2个例子,结果表明,时域有限差分方法能有效的模拟井间地下介质中电磁波的传播,揭示电磁波传播规律.该方法速度快、精度高、结果稳定,适合用于井间电磁场反演成像的正演响应计算. 相似文献
5.
井间地震资料测井约束层析成像 总被引：1，自引：0，他引：1
王守东 《CT理论与应用研究》2009,18(4)
利用井间地震资料进行层析成像,由于受观测角的限制,其成像过程具有很强的不稳定性.为了提高成像过程的稳定性,本文将测井资料引入到反演过程中,利用测井资料约束井间地震反演成像过程.采用假设模型进行数值模拟和反演实验,并与无约束反演方法进行对比实验,多种反演结果对比表明,测井约束反演方法改善了层析成像过程的稳定性,提高了成像结果的精度. 相似文献
6.
《地球物理学进展》2017,(1)
随着城市工程勘探及煤矿采矿区勘探要求的不断提高,二维地面及孔间电阻率成像无法确定电性异常沿垂直剖面方向延伸范围,三维地面电阻率成像面临着纵向分辨率小、空间覆盖不均匀等问题.本文提出井地井三维全通道电阻率成像方法,该方法计算时采用除去2个供电电极外其余全部接地电极所采集数据进行反演计算.在成像时首先采用有限差分方法求解3D静电场方程,并基于伴随矩阵方法计算非线性灵敏度矩阵,最后利用牛顿共轭梯度反演方法实现全通道电阻率层析成像.通过理论数据测试表明:全通道3D井地井联合观测方式能有效的对孔间电性结构三维成像,具有较好的纵横向分辨率,可以较好的解决实际工程地质问题及采空区在空间的分布形态难题. 相似文献
7.
8.
《应用地球物理》2019,(3)
为了解决地下复杂微小构造的精细勘探问题,使用基于高斯射线束理论的逆高斯束叠加成像方法对井间地震反射波进行成像研究。井间地震常规射线类叠加成像方法由于覆盖次数受到限制,构造复杂区精细成像质量欠佳,而以波动方程为基础的井间地震偏移方法的成像时效性不高。本文方法借鉴高斯束合成地震记录的思想,将共炮集地震数据逆高斯束分解成共反射点道集数据,然后选取合适的面元进行共反射面元数据叠加,实现了井间地震逆高斯射线束共反射点叠加。与传统的VSP-CDP叠加成像方法相比,成像范围更加广泛,且适应复杂地质构造。该方法不仅能够对二维井间地震勘探资料进行成像,针对三维井间地震资料采用基于宽线处理思路的逆高斯束叠加成像方法仍可以处理复杂构造及斜井成像问题。理论模型及实际资料试算验证了本文研究方法的有效性与稳健性。 相似文献
9.
10.
航空电磁数据的三维解释由于数据量大需要有高效的反演算法作为支撑.本文利用两种目前主流的数值优化技术（非线性共轭梯度和有限内存的BFGS法）实现了三维频率域航空电磁反演,并进一步比较了两种方法的有效性和运算效率.在反演过程中,为了更好地反演异常体的空间位置,模型方差矩阵中的光滑系数在反演起始阶段取值较大;当数据拟合差下降趋于平缓时,再利用较小的光滑因子约束反演过程来实现聚焦和获得精确的反演结果.理论数据反演表明这两种优化策略具有相似的内存需求,但是有限内存的BFGS技术比非线性共轭梯度法在计算时间和模型反演分辨率上具有一定的优越性,因此有限内存BFGS法更适合于求解大规模三维反演问题. 模型试验进一步表明目前主流的迭代法求解技术不适合大规模航空电磁数据反演,未来移动平台多源电磁数据快速正反演可通过引入矩阵分解技术来实现. 相似文献
11.
三维电阻率探测的线性反演和非线性反演中均存在着多解性的固有难题.电阻率线性反演方法的效率较高,但反演结果对初始模型的依赖性较强,易陷入局部极小;而非线性反演方法不依赖初始模型,但搜索效率极低,尚未见到关于三维电阻率非线性反演的文献.针对上述问题,融合线性与非线性反演方法的互补优势,提出了最小二乘法(线性方法)与改进遗传算法(非线性方法)相结合的混合反演方法的概念和思想.首先,提出了将介质电阻率变化范围作为不等式约束引入反演方程的思路,以实现压制多解性、提高可靠性的目标.提出了宽松不等式约束和基于钻孔推断的局部严格不等式约束的获取及定义方法.在此基础上,分别提出了基于不等式约束的最小二乘线性反演方法和遗传算法非线性反演方法.其次,对于遗传算法在变异搜索方向控制、初始群体产生等方面进行了改进,优化了其搜索方向和初始群体多样性.然后,提出了混合反演方法及其实现方案,利用改进遗传算法进行第一阶段反演,发挥其对初始模型的依赖程度低的优势,搜索到最优解附近的空间,输出当前最优个体;利用最小二乘法进行第二阶段反演,将遗传算法得到的当前最优个体作为初始模型,在最优解附近空间执行高效率的局部线性搜索,最终实现地电结构的三维成像.最后,开展了合成数据与实际工程算例验证,与传统最小二乘方法进行了对比,发现混合反演方法在压制多解性、摆脱初始模型依赖和提高反演效果方面有较好效果. 相似文献
12.
A two-dimensional forward and backward algorithm for the controlled-source audio-frequency magnetotelluric (CSAMT) method is developed to invert data in the entire region (near, transition, and far) and deal with the effects of artificial sources. First, a regularization factor is introduced in the 2D magnetic inversion, and the magnetic susceptibility is updated in logarithmic form so that the inversion magnetic susceptibility is always positive. Second, the joint inversion of the CSAMT and magnetic methods is completed with the introduction of the cross gradient. By searching for the weight of the cross-gradient term in the objective function, the mutual influence between two different physical properties at different locations are avoided. Model tests show that the joint inversion based on cross-gradient theory offers better results than the single-method inversion. The 2D forward and inverse algorithm for CSAMT with source can effectively deal with artificial sources and ensures the reliability of the final joint inversion algorithm. 相似文献
13.
3-D simultaneous inversion for velocity and reflector geometry using multiple classes of arrivals in a spherical coordinate frame 总被引：1，自引：0，他引：1
Guo-jiao Huang Chao-ying Bai Xing-wang Li Stewart Greenhalgh 《Journal of Seismology》2014,18(1):123-135
The two key requirements in conducting 3-D simultaneous traveltime tomography on real data at the regional and global scale with multiple classes of arrival time information are (1) it needs an efficient and accurate arrival tracking algorithm for multiply transmitted, reflected (or refracted) and converted waves in a 3-D variable velocity model with embedded velocity discontinuities (or subsurface interfaces), and (2) a subdimensional inversion solver is required which can easily search for different types of model parameters to balance the trade-off between the different types of model parameter updated in the simultaneous inversion process. For these purposes, we first extend a popular grid/cell-based wavefront expanding ray tracing algorithm (the multistage irregular shortest-path ray tracing method), which previously worked only in Cartesian coordinate at the local scale, to spherical coordinates appropriate to the regional or global scale. We then incorporated a fashionable inversion solver (the subspace method) to formulate a simultaneous inversion algorithm, in which the multiple classes of arrivals (including direct and reflected arrivals from different velocity discontinuities) can be used to simultaneously update both the velocity fields and the reflector geometries. Numerical tests indicate that the new inversion method is both applicable and flexible in terms of computational efficiency and solution accuracy, and is not sensitive to a modest level of noise in the traveltime data. It offers several potential benefits over existing schemes for real data seismic imaging. 相似文献
14.
由于地面微地震监测台站布设在地表,会受到地表起伏、低降速带厚度和速度变化的影响,降低了微地震事件的识别准确度和定位精度,限制了地面微地震监测技术在复杂地表地区的应用.因此,将三维地震勘探技术的思路引入到地面微地震监测中,提出了三维地震与地面微地震联合校正方法,将油气勘探和开发技术更加紧密地结合在一起.根据三维地震数据和低降速带测量数据,通过约束层析反演方法建立精确的近地表速度模型,将地面微地震台站从起伏地表校正到高速层中的平滑基准面上,有效消除复杂近地表的影响.其次,根据射孔数据和声波测井速度信息,通过非线性反演方法建立最优速度模型,由于已经消除复杂近地表的影响,在进行速度模型优化时不需要考虑近地表的影响,因而建立的速度模型更加准确.最后,在精确速度模型的基础上,通过互相关方法求取剩余静校正量,进一步消除了复杂近地表和速度模型近似误差的影响.三维地震与地面微地震联合校正方法采用逐步校正的思路,能够有效消除复杂近地表的影响,提高微地震数据的品质和速度模型的精确度,保证了微地震事件的定位精度,具有良好的应用前景. 相似文献
15.
An approximate‐simple method for nonlinear response estimates of reinforced concrete frames subjected to near‐field and far‐field records is presented in this paper. The approximate method is based on equivalent single‐degree‐of‐freedom and linear multi‐degree‐of‐freedom models. In this procedure, the nonlinear maximum roof displacement is estimated using an effective period factor and elastic response spectrum with an equivalent damping. The effective period factor was proposed for far‐field and near‐field ground motion records. For regions of high seismicity, the maximum roof displacement can be estimated by applying an effective period factor of 2.3 and 2.1 for near‐field and far‐field records, respectively, and 9% damped displacement response spectrum. For regions of moderate seismicity, a lower effective period factor of 1.9 and 1.8, for near‐field and far‐field records, respectively, can be applied to estimate the maximum roof displacement. A relationship between linear and nonlinear response of multi‐degree‐of‐freedom systems was also proposed to obtain estimates of the maximum inter‐story drift of nonlinear responding reinforced concrete frames. In addition, the effects of number of ground motion records used in the analyses on the scatters of results were investigated. The required number of ground motions to produce a reliable response was proposed. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
16.
深度速度模型的构建仍然是地震成像中的巨大挑战,获得一个精确的深度速度模型和减少深度成像项目周期都是至关重要的。常规层析反演速度建模每次迭代相当于一次线性反演,且需要重新的拾取工作,导致非常耗时,效率低下。本文提出非线性层析反演速度建模技术来建立速度模型。拾取共成像点道集的RMO量,转换到叠前域,作为运动学不变量,通过层析反演迭代进行模型更新。用一个多次的线性反演来逼近一个非线性的物理过程,避免重复的拾取工作,大大提高项目的运转效率。通过实例证明该方法的有效性。 相似文献
17.
Multistep inversion workflow for 3D long‐offset damped elastic waves in the Fourier domain 
下载免费PDF全文
下载免费PDF全文 We present a new workflow for imaging damped three‐dimensional elastic wavefields in the Fourier domain. The workflow employs a multiscale imaging approach, in which offset lengths are laddered, where frequency content and damping of the data are changed cyclically. Thus, the inversion process is launched using short‐offset and low‐frequency data to recover the long spatial wavelength of the image at a shallow depth. Increasing frequency and offset length leads to the recovery of the fine‐scale features of the model at greater depths. For the fixed offset, we employ (in the imaging process) a few discrete frequencies with a set of Laplace damping parameters. The forward problem is solved with a finite‐difference frequency‐domain method based on a massively parallel iterative solver. The inversion code is based upon the solution of a least squares optimisation problem and is solved using a nonlinear gradient method. It is fully parallelised for distributed memory computational platforms. Our full‐waveform inversion workflow is applied to the 3D Marmousi‐2 and SEG/EAGE Salt models with long‐offset data. The maximum inverted frequencies are 6 Hz for the Marmousi model and 2 Hz for the SEG/EAGE Salt model. The detailed structures are imaged successfully up to the depth approximately equal to one‐third of the maximum offset length at a resolution consistent with the inverted frequencies. 相似文献
18.
球坐标系下多震相走时三参数(速度、震源位置和反射界面)同时反演需要解决两个关键问题:(1)球坐标系下3D速度模型中多次透射、反射(折射)及转换波精确、快速的射线追踪;(2)同时反演时三种不同参数间的强耦合问题.为此,我们将直角坐标系下分区多步不规则最短路径算法推广至球坐标系中,进行区域或者全球尺度的多震相射线追踪.然后将其与适合多参数同时反演的子空间算法相结合,形成一种球坐标系下联合多震相走时三参数同时反演的方法技术.与双参数(速度和反射界面或速度和震源位置)同时反演的数值模拟对比分析显示:三参数与双参数的同时反演结果大体接近,并且它们对到时数据中可容许的随机噪声不太敏感.结果说明本文中的同时反演成像为一种提高成像分辨率,同时反演速度、震源位置和反射界面的有效方法. 相似文献
19.
基于光滑约束的最小二乘法是三维电阻率反演的主要方法,但该方法在某些情况下存在着多解性较强的问题,且普遍耗时较长,严重制约了三维反演方法的推广与发展.为改善上述问题,将表征模型参数变化范围的不等式约束作为先验信息引入最小二乘线性反演方法中,有效地改善了反演结果的精度,降低了反演的多解性问题.为了解决耗时较长的问题,基于预条件共轭梯度(PCG)算法和Cholesky分解法的特点提出了一套优化三维电阻率反演计算效率的计算方案.在该方案中,Cholesky分解法被用来求解敏感度矩阵计算中的多个点源场的正演问题,Cholesky分解法只需对总体系数矩阵进行一次分解,然后对不同的右端向量进行回代即可.将预条件共轭梯度法引入到三维电阻率反演方程的求解中,将雅可比迭代中的对角阵作为预处理矩阵,其具有求逆方便、无需内存空间的特点,有效地加快了收敛速度.对合成数据以及实测数据的反演算例表明,借助不等式约束和反演效率优化方案,最小二乘反演方法可得到较为精确的反演结果,有效地提高了反演计算效率,具有良好的推广前景. 相似文献