简介：采用声波高阶有限差分算法,克服了有限差分法本身存在的数值频散问题,使用了最佳匹配层吸收边界,较好地消除了边界反射,成功模拟了VSP波场在地下介质中的传播情况。通过三维VSP数值模拟,获得了分辨反射波、绕射波等各种波形的声波VSP剖面,证实了该方法的有效性。

简介：采用高阶交错网格有限差分法算子构建了各向异性介质多分量弹性波动方程的高精度正演模拟和叠前逆时深度偏移的数值离散方程,在正演过程中采用最大绝对振幅能量法实现各向异性介质多分量初至旅行时的计算,并以此作为多分量双程弹性波叠前逆时成像条件,同时给出计算所需的稳定性条件、多分量叠前逆时成像的基本原理、吸收边界条件等,采用凹陷模型合成了共炮点道集,并进行多分量双程各向异性弹性波叠前逆时成像研究,并与对正演炮集进行多分量各向同性叠前逆时成像结果进行对比。计算结果表明,考虑了各向异性弹性波场的矢量特性,能够更为准确地实现叠前多分量弹性波场的成像问题,使地质层位中的断层、断点等复杂目标成像更加清晰准确,并且偏移成像精度较高,因此开展各向异性叠前逆时成像可为当前高精度岩性地震勘探提供方法指导。

简介：在目前常用的周跳探测与修复方法基础上，提出了首先将观测资料按照观测历元不连续分成若干小弧段，然后利用差分法进行周跳探测，根据差分后周跳放大的特性判断周跳和野值，并确定其位置利用宽带组合和电离层组合的方法解算周跳大小。通过实例验证了其有效性。

简介：传统上，有限差分的差分系数一般可以通过泰勒级数展开法或优化方法来极小化频散误差得到。基于泰勒级数展开的差分法在有限的波数范围内精度较高，但在这个范围之外会产生较强的数值频散；基于最小二乘的优化有限差分法能在更大的波数范围内达到较高的精度，并可以在较小的计算需求内获得全局最优解。本文将基于最小二乘的优化有限差分法从二维正演模拟推广到三维，形成了计算效率高、高精度范围宽、适合并行计算的三维声波优化有限差分方法。频散分析及正演模拟表明本文发展的有限差分方法可以很好地压制数值频散。最后，将本文发展的有限差分方法应用到三维逆时偏移的震源波场延拓和检波点波场延拓中，并结合有效边界存储策略与checkpointing技术在GPU集群上实现三维逆时偏移以提高计算效率、减少存储量。三维逆时偏移试算结果表明本文三维优化有限差分方法与传统的有限差分法相比可以获得更高精度的偏移成像结果。

简介：有限差分方法广泛应用于求解许多科技领域所涉及的偏微分方程，高阶显式有限差分方法通常用来提高求解精度，已经提出的高阶隐式有限差分方法和截断高阶显式有限差分方法可用来进一步提高模拟精度而不增加计算量。本文首先计算了针对常规网格上的一阶导数和二阶导数、交错网格上的一阶导数的有限差分系数，发现高阶隐式有限差分系数中存在一些小的系数。频散分析结果表明：忽略这些小的差分系数能够近似维持有限差分的精度，但是显著减小了计算量。然后，引入镜像对称边界条件来提高隐式有限差分方法的精度和稳定性，采用混合吸收边界条件来减小来自模型边界所不需要的反射。最后，给出了针对均匀和非均匀介质模型的弹性波模拟例子，表明了本文方法的优点。

简介：针对有限差分数值模拟的频散问题，本文将交错网格技术和紧致差分格式相结合，推导了横向各向同性介质一阶速度一应力波动方程的紧致交错网格差分格式；对比分析了紧致交错网格差分格式、交错网格差分格式以及紧致差分格式的截断误差主项，并利用Fourier误差分析方法分析了上述三种差分格式的近似精度；在此基础上，分别采用上述三种差分格式进行了波场数值模拟。结果表明，当差分方程阶数相同时，紧致交错网格差分格式截断误差最小，数值频散最弱，差分精度最高，证实了该方法的有效性。

简介：改变航迹线和交会定位公式进行水上固定断面测量，更新了传统水上固定断面复测定位法，新法具有不偏线、点位密、计算简便、精度高等优点，能为生产、设计、研究提供相对精确的地理模型。

简介：在地震波场计算中，计算机只能处理有限的波场区域，导致出现人工边界，而人工边界会影响正常波场的计算和最终成果的解释。简要介绍了声波波动方程模拟中的两大类吸收条件，基于一个简单的二维均匀速度模型，采用有限差分法实现其吸收效果对比。正演模拟结果表明，最佳匹配层吸收边界效果最好。

简介：地震波场正演模拟是地震资料处理、解释中最为重要的技术之一。地震波场正演模拟在大时间步长、长时程的波场延拓中,存在计算不稳定的问题。本文基于声波方程的Hamilton表述,在波动方程求解中用辛差分格式进行时间网格离散,用傅里叶有限差分进行空间网格离散,提出一种新的保结构地震波场正演模拟方法一辛格式傅里叶有限差分法,在保证计算精度的同时提高计算的稳定性。利用声学近似处理空间-波数混合域的积分算子,将该方法推广至各向异性介质。给出各向同性和各向异性条件下的地震正演模拟的计算流程,并将本文方法用于BP盐丘、BPTTI等模型的波场正演模拟。数值算例表明本文开发的方法适用于速度变化剧烈的复杂介质地震波场正演模拟,计算精度高,数值频散小,在各向异性介质正演中能够有效避免qSV波残余,在大时间步长的迭代计算中稳定性好。本文为在辛算法的框架下实现高精度地震正演模拟提供了一种新的选择。

简介：水准网平差及粗差检验计算＠童小华＠张锦￥同济大学测量系￥山西矿业学院地测系基准方程，粗差检验，工程实例水准网平差及粗差检验计算童小华张锦（同济大学测量系，上海，２０００９２）（山西矿业学院地测系，太原，０３０００４）【关键词】基准方程粗差检验工程实例一、水...

简介：Finite-difference(FD)methodsarewidelyusedinseismicforwardmodelingowingtotheircomputationalefficiencybutarenotreadilyapplicabletoirregulartopographies.Thus,severalFDmethodsbasedonthetransformationtocurvilinearcoordinatesusingbody-fittedgridshavebeenproposed,e.g.,standstaggeredgrid(SSG)withinterpolation,nonstaggeredgrid,rotatedstaggeredgrid(RSG),andfullystaggered.TheFDbasedontheRSGissomewhatsuperiortoothersbecauseitsatisfiesthespatialdistributionofthewaveequationwithoutadditionalmemoryandcomputationalrequirements;furthermore,itissimplertoimplement.WeusetheRSGFDmethodtotransformthefirstorderstress–velocityequationinthecurvilinearcoordinatessystemandintroducethehighprecisionadaptive,unilateralmimeticfinite-difference(UMFD)methodtoprocessthefreeboundaryconditionsofanirregularsurface.Thenumericalresultssuggestthattheprecisionofthesolutionishigherthanthatofthevacuumformalism.Whentheminimumwavelengthislow,UMFDavoidsthesurfacewavedispersion.WecompareFDmethodsbasedonRSG,SEM,andnonstaggeredgridandinferthatallsimulationresultsareconsistentbutthecomputationalefficiencyoftheRSGFDmethodishigherthantherest.

简介：由于河口考察，海岸带调查等工作的深入开展，对于各河口和海的分界就需要有明确的划分方法。这对海岸线长度，沿海滩涂面积，内海和内河水域而积的准确量算以及对于科研、教学、经济开发等都有重要意义。

简介：传统卫星测高数据处理中共线平差的精度极大影响最终解算结果,针对经初步处理后卫星测高数据仍可能存在粗差的问题,提出一种融合抗差的共线秩亏平差算法.该算法结合传统共线平差及抗差估计的优点,保留了秩亏自由网平差不需要基准控制点的优良性质,并有效地剔除秩亏自由网中存在的粗差,实现更高精度的解算结果.实验表明,相比较传统共线秩亏平差算法而言,提出的算法解算结果精度更高.

简介：广义S变换可将地震记录分解到时频域,通过时频域上的不同分布,区分有效波和面波。时频域滤波只会消除面波发育时间处的频率成分,而不影响其他时间点的信号,减少了对有效波的影响,同时避免了区域滤波在时窗边缘处的波形畸变作用。理论模型和实际资料测试均取得了良好的处理效果。

简介：利用上海天文台的SHORDE1软件对约3年左右(MJD49001.1-MJD50109.6)的全球SLR资料进行了归算,并将所得的极移序列和相应历元的国际地球自转服务中心(IERS)给出的极移序列作了比较,分析结果表明:上海天文台基于SLR资料的极移序列与IERS极移序列之间存在一定程度的系统偏差。该系统差具有明显的周期特性,周期约为427天,振幅为1mas左右。

简介：用曲线积分代替算术平均法来计算洪水流量演算中的时段平均流量,减少了时段△t对流量值的影响,使计算流量更接近实测值.

简介：在实际地震资料处理解释中,广义S变换的作用举足轻重,其为一种介于连续小波变换和短时傅里叶变换的延伸方法,其突出特点为不同频率处具有不同分辨率.为了对川南丹凤场中联井剖面资料进行储层描述,本文运用了广义S变换分频处理技术,得到地震信号的瞬时参数.通过不同频率剖面的对比与研究,完成储层描述.

简介：波动方程有限差分法是地震数值模拟中的一种重要的方法，对理解和分析地震传播规律、分析地震属性和解释地震资料有着非常重要的意义。但是有限差分法由于其离散化的思想，产生了不稳定性。精细积分法在有限差分法的基础上，在时间域采用解析解的表达形式，在空间域保留任意差分格式，发展成为半解析的数值方法。本文结合并发展了以往学者的成果，推导了任意精细积分法的三维弹性波正演模拟计算公式，并对其稳定性进行了数值分析。在计算实例中，实现了精细积分法二维和三维弹性波模型的地震正演模拟，对计算结果的分析表明，精细积分法反射信号走时准确，稳定性好，弹性波场相较于声波波场，弹性波波场成分更为丰富，包含了更多波型成分（PP-和PS-反射波、透射波和绕射波），这对实际地震资料的解释和储层分析有重要的意义。实践证明，该方法可直接应用到弹性波的地质模型的数值模拟中。

简介：常规多波联合反演采用Zoeppritz方程的近似式构建正演方程,反演过程中需要假定背景纵横波速度比为常数,其反演精度不高,稳定性不好。本文提出了一种基于精确Zoeppritz方程的多波联合反演方法,结合贝叶斯方法进行广义线性反演。本方法基于精确Zoeppritz方程构建正演方程,避免了近似式反演在大角度时引起的误差;利用贝叶斯方法引入模型参数的先验分布信息,作为反演的正则化项,降低了反演的不适定性;反演目标函数中引入低频软约束,稳定了反演低频结果,提高了反演的鲁棒性;在求解反演目标函数时,利用快速算法,降低了反演的运算量。经过模型试算,证明了该方法的优越性和抗噪性;并在实际资料的应用中证明了该方法的实用性和有效性。

简介：频率-波数域单程波算子能高效地模拟地震波在复杂介质中的传播，但是在描述波的大角度传播和速度横向扰动变化较大介质中传播的问题时仍然存在一定误差。这类误差是由于对单平方根算子使用Taylor展开式的近似程度不足所造成。为了进一步提高泰勒展开式的精确性，本文提出一种利用粒子群智能算法优化级数展开系数的高阶广义屏算子对单平方根算子的展开级数进行优化处理。新的偏移算法能在保持单程波偏移算法高效的前提下进一步提高偏移算子在大角度的成像精度和对强横向速度变化介质的适应性。通过脉冲响应实验，验证了基于粒子群算法优化级数的高阶广义屏算子能够提高常规的高阶广义屏算子的成像精度和成像角度。根据对二维SEG/EAGE盐丘模型的成像处理，基于粒子群算法优化级数的高阶广义屏算子对盐丘下面的断层取得了更高质量的成像，说明粒子群优化级数的高阶广义屏算子比常规的高阶广义屏算子具有更好的横向速度适应性。为了检验本文所提算法对实际资料的处理能力，我们利用常规的偏移处理技术和本文所提算法对一条海上二维数据进行了偏移成像处理，对比分析成像剖面发现本文所提算法描述了更加清晰的层位信息和更高质量的偏移剖面。本文所提算法能有效提高高阶广义屏偏移在广角度成像的能力，具有一定实际应用价值。