共查询到20条相似文献，搜索用时 171 毫秒
1.
为对复杂声速剖面进行反演,笔者采用速度分析建模与伴随状态法层析成像反演建模相结合的方法获得声速剖面沿垂向速度。采用Munk公式描述水层速度垂向变化,并引入随机声速剖面复杂速度建模。在此基础上,使用海底反射走时伴随状态法层析成像反演复杂声速剖面。将声速剖面较准确的速度分析结果作为初始迭代,通过求解伴随偏微分方程得到非线性目标函数的梯度,避免了对Fréchet导数的大量计算。较准确叠加速度初始模型的引入提高了迭代效率与计算精度。最终得到了与真实模型速度吻合的结果。 相似文献
2.
Chirp浅剖图像可以显示东沙群岛西南海域存在泥火山.通过基于Chirp浅剖记录数据的反演方法研究,定量地分析泥火山区海底物性参数.首先由Chirp浅剖记录的反射振幅计算了海底反射系数,结果表明泥火山区的反射变化很大(0.1~0.8).对于海底弱反射(＜0.45) 区域,用Schock-Stoll模型以及物理参数间的经验公式反演海底浅表层物理性质参数(速度、密度、孔隙度等参数);对于海底强反射(＞0.45) 区域,用Gardner经验公式反演海底浅表层物理性质参数.反演结果表明研究区海底物性变化大,声速最高达5 237 m/s,密度最大为2.673 g/cm3.反演的物性参数与实验室测量值吻合度较高.Chirp浅剖数据反演是一种有效可行的遥测海底声学参数方法. 相似文献
3.
介绍了在倾斜均匀各向同性层状介质情况下,根据地震反射走时反演二维地层结构和地层速度的地震射线层析成像方法——走时反演。正演模型射线追踪是根据Fermat原理实现的,即求解满足该原理的非线性方程组得到射线与界面的交点,进而计算相应的走时。反演是先假设一初始模型,用最优化方法使射线追踪走时与观察走时的残差极小。最后计算了分辨矩阵和信息密度矩阵,以评价反演结果。对有噪情形也进行了反演。 相似文献
4.
二维地层结构和速度的走时反演 总被引：1，自引：0，他引：1
介绍了在倾斜均匀各向同性层状介质情况下，根据地震反射走时反演二维地层结构和地层一地震射线层析成像方法－走时反演。正演模型射线追踪是根据Ｆｅｒｍａｔ原理实现的，即求解满足该原理的非线性方程组得到射线与界面的交点，进而计算相应的走时。反演是先假设一初始模型，用最优化方法使射线追踪走时的残差极小。最后计算了分辨矩阵和信息密度矩阵，以评价反演结果。对有噪情形也进行了反演。 相似文献
5.
6.
浅海海底起伏和速度变化对OBC资料成像质量产生较大影响。建立浅海海底表层速度模型不仅能够解决OBC资料的静校正问题,也可用于海底反射系数计算、双检资料合并、多次波压制等,但是,目前针对浅海海底表层速度建模的研究还不多。文中提出了针对浅海地区OBC资料的海底表层速度建模的三维初至走时反演技术,主要包括:(1)震源位置校正技术。根据地震波在海水中传播特征,把在海水中激发的震源位置校正至海底,使震源和接收点都位于海底,利于初至走时反演计算;(2)快速三维初至走时反演方法。利用回折波走时和射线方程,形成了高效率初至走时反演方法。将该技术应用于胜利油田莱州湾浅海区海底OBC资料的处理中,建立了三维海底表层速度模型,用此速度模型进行静校正,取得了良好的应用效果。 相似文献
7.
作为一种模拟地震波走时的传统途径，射线追踪不仅对地震波理论研究具有重要意义。而且也可直接应用于地震波反演及偏移成像等过程。在一般条件下，射线追踪只能用数值方法来实现，但当介质上具有垂向不均匀性时，射线路径及其走时可以由解析式给出。根据适当的解析公式，可建立反射点和射线路径以及走时之间的关系，并用VC^ 6．0编制了基于图形界面的计算程序。这些工作对于逐步开发地震波的偏移和反偏移程序具有基础性的意义。 相似文献
8.
在利用有限差分等基于网格的数值分析方法解地震波走时所满足的程函方程时,由于速度模型的网格化离散等原因,会使走时在各网格节点之间不具有计算射线路径所要求的光滑性,即走时在邻近网格节点之间不具有连续的一阶导数。因此,直接利用网格节点走时计算射线路径会使最终的射线路径不光滑。为解决这个问题,已有研究者提出了基于B样条插值的逆向梯度方案（法）。然而,在速度发生突变时,B样条逆向梯度法所计算出的射线路径会具有较大的误差。针对这个问题,首先采用适合于解最小零偏差逼近及最佳平方逼近问题的Chebyshev多项式取代B样条对来自于分区多级计算方案的网格节点走时进行最佳逼近,得到在最小平方意义下的最优走时公式;然后采用与B样条逆向梯度法类似的计算过程得到光滑的射线路径。数值实验表明,利用Chebyshev多项式逼近走时可以得到具有很高精度的多次反射射线路径,在多次波偏移成像研究中具有潜在的价值。 相似文献
9.
跨孔雷达走时层析成像主要利用雷达波的走时进行反演,走时提取的正确与否将直接影响到层析成像的效果。数字图像分割法基于凸集投影（POCS）方法,使用能量比彩色图像分割技术准确提取走时。数字图像分割法提取走时首先应用在折射地震波的数据处理中。笔者首次将数字图像分割法提取走时的方法应用到跨孔雷达走时层析成像中,使用迭代线性反演算法重建了雷达波速度场。反演过程中,使用最小二乘QR分解法（LSQR）求解线性方程组,利用弯曲射线追踪技术构建雅可比矩阵,走时的计算值则由多模板快速推进算法（MSFM）得到。为了验证数字图像分割法在走时层析成像中的效果,使用一组合成数据及一组实测数据分别对比了基于数字图像分割法和基于传统能量比法获得的层析成像反演结果。对比结果表明,使用数字图像分割法得到的层析成像结果更为精确,误差更小,能为判断地下雷达波速度场提供更为有力的帮助。 相似文献
10.
由于地球物理反演存在固有的多解性问题,如何减小多解性,提高地球物理方法反演结果的可靠性和准确性是当今地球物理领域研究的热点。为了减小大地电磁与地震反射波走时资料反演的多解性,在Gallardo和Meju研究的交叉梯度理论的基础上,实现大地电磁与地震反射波走时数据的二维联合反演算法;设计理论地电模型和速度模型进行合成数据的单独反演和联合反演试算,并验证大地电磁和地震反射波走时数据二维联合反演算法的正确性和有效性。反演结果表明,利用交叉梯度方法耦合大地电磁电阻率参数和地震波速度参数能使联合反演过程中的这两个模型在更新的时候相互约束,联合反演结果比单一反演结果更接近真实模型,从而使反演结果更可靠,减少了反演解的多解性。 相似文献
11.
分离P-P波和P-SV波最常用的是τ-p变换法, 但是在τ-p域中, P-P波和P-SV波常有重叠部分, 不能很好分离。地震勘探中, 当炮检距小于界面深度时, P-P波和P-SV波的时距曲线方程都可近似为抛物线方程。针对这种情况, 提出了一种基于抛物线方程的坐标拉伸τ-p变换的P-P、P-SV 波分离方法, 该方法首先对地震记录沿空间方向上进行坐标拉伸, 再进行τ-p正变换, 在τ-p域中分离P-P波和P-SV波, 其次对P-P波和P-SV波分别进行τ-p反变换, 实现P-P波和P-SV波的分离。与常规τ-p变换法对比分析表明, 利用该方法能很好地将合成地震记录中的P-P波和P-SV波分离。 相似文献
12.
This paper presents a semianalytical approach for solving first-order perturbation (FOP) equations, which are used to describe dissolution-timescale reactive infiltration instability (RII) problems in fluid-saturated rocks. The proposed approach contains two parts because the chemical dissolution reaction divides the whole problem domain into two subdomains. In the first part, the interface-condition substitution strategy is used to derive the analytical expressions of purely mathematical solutions for the FOP equations in the upstream subdomain, where the dissolution chemical reaction is ceased and the FOP equations are weakly coupled. In the second part, the finite element method (FEM) is used to derive the analytical expressions of numerical solutions for the FOP equations in the downstream subdomain, where the dissolution chemical reaction needs to be considered and the FOP equations are strongly coupled so that it is impossible to derive purely mathematical solutions for them. Particular attention is paid to the development of the element-by-element forward marching strategy, which is associated with the use of the FEM for solving this new kind of scientific problem. The related analytical results demonstrated that (1) both the dynamic characteristic of a reactive infiltration system and the dimensionless wavenumber can have pronounced influences on the distribution of the FOP dimensionless acid concentration within the entire domain of the dissolution-timescale RII problems in fluid-saturated rocks and (2) the FOP dimensionless acid concentration distribution exhibits two significantly different patterns in the upstream and downstream subdomains of the dissolution-timescale RII system. 相似文献
13.
地震AVO （amplitude versus offset）技术是一项利用振幅信息研究岩性、检测油气的地震勘探技术。常规方法基于Zoeppritz方程计算模型界面处的反射系数,而实际地下非均匀介质中地震反射特征不仅与入射角度、物性差异有关,还与入射波频率、地层厚度、薄互层结构等因素有关。为此,应用传播矩阵理论充分考虑与这些因素有关的调谐干涉等传播效应,针对实际数据计算非均匀地下介质的高精度合成地震记录,对比Zoeppritz方程、Shuey二项近似方程、Shuey三项近似方程以及传播矩阵算法的模拟效果。研究发现：在小角度入射时Zoeppritz方程、Shuey二项近似方程和Shuey三项近似方程的反射振幅和波形基本一致,大角度入射时Zoeppritz方程与Shuey三项近似方程接近;Zoeppritz算法的模拟结果在小角度入射和浅层情况下与传播矩阵算法差别较小,而在大角度入射和深层情况下与传播矩阵差别较大,说明层间多次波的调谐干涉等传播效应不可忽略。 相似文献
14.
从Aki和Richards在弱反差条件下给出的P-S波反射系数近似公式出发,提出了一种在小角度入射范围内的P-S波反射系数近似公式.为了得到该近似公式,对原公式中的角度项进行了二项式展开,直接利用入射角代替入射角的正弦.不同模型的计算表明,该近似在小角度入射范围内相当于Aki和Richards的原近似公式,在大角度时有一定改善.新的近似公式形式简单、涉及的弹性参数较少,容易应用到反演当中.另外,利用一阶小角度近似结果可以对横波阻抗进行反演. 相似文献
15.
Two types of garnet porphyroblast occur in the Schneeberg Complex of the Italian Alps. Type 1 porphyroblasts form ellipsoidal pods with a centre consisting of unstrained quartz, decussate mica and small garnet grains, and a margin containing large garnet grains. Orientation contrast imaging using the scanning electron microscope shows that the larger marginal garnet grains comprise a number of orientation subdomains. Individual garnet grains without subdomains are small (< 50 µm), faceted and idioblastic, and have simple zoning profiles with Ca‐rich cores and Ca‐poor rims. Subdomains of larger garnet grains are similar in size to the individual, small garnet grains. Type 2 porphyroblasts comprise only ellipsoidal garnet, with small subdomains in the centre and larger subdomains at the margin. Each subdomain has its own Ca high, Ca dropping towards subdomain boundaries. Garnet grains, with or without subdomains, all have the same Ca‐poor composition at rims in contact with other minerals. The compositional zonation patterns are best explained by simultaneous, multiple nucleation, followed by growth and amalgamation of individual garnet grains. The range of individual garnet and garnet subdomain sizes can be explained by a faster growth rate at the porphyroblast margin than in the centre. The difference between Type 1 and Type 2 porphyroblasts is probably related to the growth rate differential across the porphyroblast. Electron backscatter diffraction shows that small, individual garnet grains are randomly oriented. Large marginal garnet grains and subdomain‐bearing garnet grains have a strong preferred orientation, clustering around a single garnet orientation. Misorientations across subdomain boundaries are small and misorientation axes are randomly oriented with respect to crystallographic orientations. The only explanation that fits the observational data is that individual garnet grains rotated towards coincident orientations once they came into contact with each other. This process was driven by the reduction of subdomain boundary energy associated with misorientation loss. Rotation of garnet grains was accommodated by diffusion in the subdomain boundary and diffusional creep and rigid body rotation of other minerals (quartz and mica) around the garnet. An analytical model, in which the kinetics of garnet rotation are controlled by the rheology of surrounding quartz, suggests that, at the conditions of metamorphism, the rotation required to give a strong preferred orientation can occur on a similar time‐scale to that of porphyroblast growth. 相似文献
16.
Acoustic imaging and sensor modeling are processes that require repeated solution of the acoustic wave equation. Solution of the wave equation can be computationally expensive and memory intensive for large simulation domains. One scheme for speeding up solution of the wave equation is the operator-based upscaling method. The algorithm proceeds in two steps. First, the wave equation is solved for fine grid unknowns internal to coarse blocks assuming the coarse blocks do not need to communicate with neighboring blocks in parallel. Second, these fine grid solutions are used to form a new problem which is solved on the coarse grid. Accurate and efficient wave propagation schemes also must avoid artificial reflections off of the computational domain edges. One popular method for preventing artificial reflections is the nearly perfectly matched layer (NPML) method. In this paper, we discuss applying NPML to operator upscaling for the wave equation. We show that although we only apply NPML to the first step of this two step algorithm (directly affecting the fine grid unknowns only), we still see a significant reduction of reflections back into the domain. We describe three numerical experiments (one homogeneous medium experiment and two heterogeneous media examples) in which we validate that the solution of the wave equation exponentially decays in the NPML regions. Numerical experiments of acoustic wave propagation in two dimensions with a reasonable absorbing layer thickness resulted in a maximum pressure reflection of 3–8%. While the coarse grid acceleration is not explicitly damped in our algorithm, the tight coupling between the two steps of the algorithm results in only 0.1–1% of acceleration reflecting back into the computational domain. 相似文献
17.
18.
Pirnia Pouyan Duhaime François Ethier Yannic Dubé Jean-Sébastien 《Acta Geotechnica》2020,15(10):2877-2889
This paper presents a coupled finite and discrete-element model (FEM and DEM) to simulate internal erosion. The model is based on ICY, an interface between COMSOL, an FEM engine, and YADE, a DEM code. With this model, smaller DEM subdomains are generated to simulate particle displacements at the grain scale. Particles in these small subdomains are subjected to buoyancy, gravity, drag and contact forces for short time steps (0.1 s). The DEM subdomains provide the macroscale (continuum) model with a particle flux distribution. Through a mass conservation equation, the flux distribution allows changes in porosity, hydraulic conductivity and hydraulic gradient to be evaluated for the same time steps at a larger, continuum scale. The updated hydraulic gradients from the continuum model provide the DEM subdomains with updated hydrodynamic forces based on a coarse-grid method. The number of particles in the DEM subdomains is also updated based on the new porosity distribution. The hierarchical multiscale model (HMM) was validated with the simulation of suffusion. Results for the proposed HMM algorithm are consistent with results based on a DEM model incorporating the full sample and simulation duration. The proposed HMM algorithm could enable the modelling of internal erosion for soil volumes that are too large to be modelled with a single DEM subdomain.
相似文献19.
近些年来地下水中多组分反应-运移模型在地球科学及环境领域开始得到应用,但其求解较为复杂,为了提高计算效率,可以采用去耦合化方法处理,从而使模型求解得到简化。针对自然界中广泛存在的非均质地质体,提出该类条件下的去耦合化方法,即根据水-岩间、水溶组分间反应的不同,将整个研究区划分为若干子区域,获得对应的去耦合化矩阵。对化学场中各子区域间相邻边界进行设定,达到简化模型求解的目的。最后,以一维非均质介质中基于热力学平衡的反应-运移问题作为算例,基于以上方法进行求解,并与该算例经PHAST软件所示的结果较为一致。结果表明,基于去耦合化方法获得的各离子浓度随时间演变和沿空间分布特征与PHAST所示的结果较为一致,显示该方法在非均质区域模拟溶质运移等方面具有较好的适用性。 相似文献
20.
介绍了联合运用积分方程法与有限元法(简称混合法)来计算三维非均匀介质中电磁场分布的理论方法,并进行了数值模拟实验。混合法的原理是引入一个包围非均匀目标体的虚构边界,在边界内部的场用有限元法模拟,在边界外部(包括边界)的场用积分方程表达,二者在边界上通过场的连续性耦合起来。数值实验结果表明,混合法既能显著地减小网格规模,又能灵活地模拟复杂的介质情况,且计算精度较高。 相似文献