地震波速

摘要： 超深层裂缝性非常规天然气藏的开发在能源领域所占比重越来越大,但目前关于该类型储层在原位应力下各向异性地震波速及渗透率的基础研究仍鲜有报道,一定程度上制约了其高效开发策略的制定进程。以塔里木油田克深2号裂缝性致密砂岩天然气储层为工程背景,利用先进的地球物理成像真三轴测试系统,开展了储层衰竭期间原位应力下完整砂岩与不同倾角裂缝岩样的各向异性地震波速以及各向异性渗透率演化规律实验研究,探讨了压缩波、剪切波、纵横波速比、3个主应力方向渗透率各向异性特性与内在机制,研究了地震波速与渗透率的内在关联特征。研究结果表明:首先,裂缝引起的结构非均质性以及三向不等应力作用下岩石各主应力方向微观结构的变形差异导致了地震波速的各向异性;储层衰竭期间地震波速度均随有效应力的增大近似线性增加,压缩波速度明显大于剪切波速度;纵横波速比随有效应力增大而单调减小,表明剪切波对外部荷载敏感性更强;通过压缩波速和纵横波波速比数据对比,发现剪切波S2速度对裂缝倾角更为敏感。其次,原位应力下各主方向渗透率均表现出显著的应力敏感性和结构(裂缝)敏感性,早期应力大小对储层衰竭中后期渗透率变化具有重要影响,渗透率各向异性随裂缝倾角减小而更加显著;与裂缝平面正(斜)交方向的渗透率相对比较平滑,而沿裂缝平面方向渗透率呈现出一定的波动特性,这主要是由于在压力流体冲刷和岩石骨架进一步受压变形双重作用下,破碎的微颗粒会在孔隙/裂隙中不断运移,造成岩石内部孔隙/裂隙周期性堵塞和疏通引起的。此外,在实验应力条件下(走滑应力状态),裂缝倾角越小,水平主应力方向上的地震波速度-渗透率相关性越不明显;最大与最小水平主应力方向上的地震波速度-渗透率关系应区分对待,最大水平主应力方向地震波速与渗透率的相关性明显优于最小水平主应力方向。

摘要： 结合新疆某水利工程实例介绍综合物探方法在中浅覆盖层,特别是覆盖层厚度在100m左右的地层中的运用.通过综合运用电测深法和浅层地震反射波法,能解决电性差异小,浅层折射测试深度浅等问题,能有效地提高物探解释的准确性和可靠性,提高物探解释精度.是比较经济、普遍的方法.

摘要： 地震波泊松比蕴涵着丰富的地质信息,根据泊松比可以识别地壳岩石的基性程度,进而分析地球动力学与成矿背景.本文搜集整理了国家地震科学数据共享中心等数据平台自20世纪80年代至2017年底的天然地震震相资料,经过数据筛选,得到华南研究区内M≥2.0级地震事件11410个.通过最小二乘法拟合数据的走时-震中距关系直线方程,得到Pg波平均波速6.11 km/s,Sg波平均波速3.60 km/s,首波Pn波平均波速8.06 km/s,Sn波平均波速4.59 km/s.在此基础上,求得华南地壳平均波速比Vp/Vs为1.697,由公式计算得到华南地壳的平均泊松比为0.234,壳幔界面附近的平均波速比Vp/Vs为1.756,计算得到壳幔界面附近的平均泊松比为0.260.走时反演给出的泊松比结果显示扬子地块、华夏地块以及东南沿海地壳的物性差异.华夏地块与扬子地块之间存在大面积的低泊松比异常区,与江南造山带存在良好的对应关系;而东南沿海、钦杭成矿带泊松比明显较高,解释为幔源物质注入地壳的结果.

摘要： The high-temperature hydrothermal activity region in the Western Sichuan Plateau is located in the northeast margin of the eastern syntaxis of the Tibetan Plateau.This area is the easternmost part of the Mediterranean-Himalayan geothermal belt,and also a notable region of geothermal research.The hot springs in this region are distributed extensively along the Jinshajiang fault,Dege-Xiangcheng fault,Garzê-Litang fault and Xianshuihe fault,forming the NW-SE strip-like high-temperature hydrothermal activity area.The geothermal analysis of this high temperature hydrothermal activity area is an important foundation for studying the deep geodynamics process of the eastern syntaxis of the Tibetan Plateau,which also contributes to exploiting and utilizing the geothermal resources in the western Sichuan Plateau.In this paper,we compare the crust-mantle thermal structure acquired from the helium isotopes measurement with geophysical results of the aeromagnetic,gravity data and shear wave data.The Moho and Curie point depth is obtained by inversion.Besides,the deep geothermal structure and geophysical background are discussed in the typical anomaly area of Batang,Litang and Kangding.The result shows that the helium isotope heat flux of the surface is in accord with the geophysical anomalies.The thermal contribution of crust and mantle in the Western Sichuan Plateau is relatively close:the average Qc/Q is 51.38℅in Batang,54.39℅in the Litang area,and 42.42℅in the Kangding area.The thermal structure pattern in this region tends to be a "warm crust-warm mantle"mode.The formation of high surface heat flux in the warm crust background need to be further studied.%川西高温水热活动区是地热学研究的典型地区之一,该区温泉、热泉点分布广泛,这些温泉大多沿北西-南东向的金沙江断裂、德格-乡城断裂、甘孜-理塘断裂、鲜水河三条主断裂呈条带状分布.对此高温水热活动区开展地热学分析,是研究青藏高原"东构造结"深部地球动力学过程、开发利用川西高原地热资源的重要基础,具有重要的科学研究意义和实际应用价值.本文针对实测氦同位素所发现的壳幔热结构比,依据重力、航空磁法、地震资料,采用反演莫霍面、居里面深度等方法,对川西高温水热活动区的深部热结构进行了分析.同时对中高温热泉密集出露的巴塘、理塘以及康定水热活动区的典型壳幔热流构与深部背景场进行了对比研究.结果表明,地表氦同位素热流估算的深部热流结构与地球物理异常结果保持较好的对应关系.川西地区热流结构中地壳、地幔的热贡献比较接近:巴塘地区平均Qc/Q=51.38℅、理塘地区平均Qc/Q=54.39℅、康定地区平均Qc/Q=42.42℅.川西地区表现出了"温壳温幔"型的深部热结构,但温地壳背景下形成较高地表热流的原因有待进一步研究.

摘要： 为了更好地使用岩石超声波速模型对地层地震成像过程进行约束,研究钻孔岩芯超声波速与地层地震波速的对应转换关系.提出一种根据岩体完整性状态,利用钻孔岩芯实测超声波速数据来获取地下介质地震波速分布的方法.首先根据岩性分类采用转换系数将岩块超声波速转换为岩块地震波速;再采用经验公式考虑水饱和度、温度和地压的影响,修正岩块波速;然后由钻孔岩芯图像确定岩体体积节理数工和岩体完整性系数Kv,将岩块地震波速转换为岩体地震波速.为了验证该转换方法的有效性,依托华阳川铀铌铅多金属矿床浅层勘探工程,选取5个钻孔开展对比研究.地层地震波速的空间分布采用背景噪声成像法反演得到,而通过超声波速转换得到的岩体波速数据,沿钻孔方向实施了空间平均化处理.两者对比结果表明,地震波速的空间趋势是一致的,数值上具有高度的相关性.

摘要： 南黄海盆地特殊的地震地质条件导致准确获得海相地层速度存在困难。基于本区的多道地震资料,联合使用多数据质控速度分析、初至波层析反演及层控网格层析反演3种方法,提取出了较为可靠的海相中-古生界速度信息。结合提取的地震速度及下扬子区已有钻井信息,总结了南黄海主要海相地层的速度分布。本区海相地层速度整体呈高低相间分布结构,存在速度的突变和倒转。海相下构造层速度高,不同地层间速度差异较小。由于海相下构造层无钻井资料约束,获得的速度信息仍存在不确定性。

摘要： 通过对本课题的研究,明确了控制爆破在泥岩砂岩地层和浅埋地层斜井施工技术,避免了爆破设计的不合理,根据实际爆破参数的验证,总结出了适合本工程泥岩砂岩地层的爆破参数计算方法及相关参数,对地震波允许参数的设定提供了参考.针对泥岩砂岩地层和浅埋地层开挖施工常见问题进行分析,提出合理的解决措施,帮助项目提前做好了合理施工方案,避免石油管道泄漏爆炸事故的发生,从而大大降低了工程施工成本和对当地环境不利的影响.%Based on the research of this subject, it is clear that the construction technology of blasting control in the mudstone sandstone strata and shallow buried strata avoids the unreasonable blasting design. According to the verification of actual blasting parameters, the calculation method of blasting parameters and related parameters of the mudstone sandstone strata is provided, which provides a reference for the setting of the allowable parameters of the seismic wave. This paper analyzes the common problems of excavation and construction of mudstone sandstone strata and shallow burial strata, and puts forward reasonable measures to help the project make a reasonable construction plan ahead of time, avoid the occurrence of oil pipeline leakage and explosion, thus greatly reducing the construction cost and adverse effects on the local environment.

摘要： 笔者建立了基于地震层析成像的热地幔对流模型,并分别计算了有、无地袁板块运动约束下的中国大陆岩石圈底部的地幔对流速度场和水平剪切应力场.结果表明,欧亚板块运动速度对中国大陆地幔浅部的对流速度场影响较大,而对对流应力场的影响较小.在岩石圈底部(约100 km深度),对流速度方向差异基本上达到了50％,幅度差异达到了80％以上,而地幔对流应力场的差异在方向上基本上小于10％、幅度上基本上小于20％.因此,如果仅研究岩石圈底部地幔对流应力场,则可以采用无板块速度约束的热对流模型.但如果要研究地幔浅部的对流状态和格局,则需要考虑板块运动速度.%In this paper,seismic tomography-based thermal mantle convection models are presented and then used to predict the mantle convective velocity and mantle convection-induced horizontal shear stress (HSS) at the base of the lithosphere below Mainland China,with or without constraint of surface plate motions.Our results show that motions of the Eurasian plate have larger effects on the convective velocity than on HSS in the uppermost mantle beneath Mainland China.At the base of the lithosphere (at depth ～ 100 kin),the percent difference between the directions of convective velocities predicted by the two models above is almost above 50％,and that between the amplitudes is almost larger than 80％;whereas for HSS,the percent differences are almost smaller than 10％ and 20％ in direction and amplitude,respectively.Therefore we conclude that,if it were only to predict the basal HSS of the lithosphere due to mantle convection,the thermal mantle convection model without constraint of the Eurasian plate motions could be adopted;however,if the convective state and pattern in the uppermost mantle were to be investigated,it is necessary to consider surface plate motions.

摘要： 岩石的地震波速是区域构造研究和浅部地震勘探的基础.哈密地区是我国重要的产油区域之一,但目前仍缺少地震波速方面的基础资料.这是在实验室条件下模拟了深层环境,分析了四组哈密地区砂岩样品的波速与压力之间的关系.结果表明:随着压力增加,砂岩样品的纵横波速度以对数关系增加;砂岩的孔隙度影响波速的变化速率;砂岩地震波速的滞后效应主要表现为:在相同压力下,降压过程中的波速大于升压过程中的波速;高压下的波速接近于用空间平均模型计算的理论波速.实验所得的波速数据与哈密地区的测井数据解释相吻合,其可以作为判断油层、水层和干层的辅助条件.%Seismic velocity of rocks have a wide range of applications in regional tectonic research and shallow seismic exploration.Hami region is one of the important production bases for oil and gas in China,but there is no data of seismic velocity.Deep environment is simulated in laboratory and the change rule of P-and S-wave velocity which come from four groups of sandstones in Hami is studied in this paper.The result shows that the P-wave and S-wave velocity increase (or decrease) with the pressure at the logarithmic rule.The porosity of sandstone affect the rate and size of increased wave velocity.The lag effect of sandstone's wave velocity mainly in the procession that wave velocity of decompression is greater than boost when they are under the same pressure.The wave velocity with high pressure is closed to the theoretical wave velocity of space averaging models.Wave velocity data measured by experiment is consistent with the log data in Hami.Therefore it can be used as the auxiliary condition to define oil layer,water layer and dry layer.

摘要： 通过一系列野外及实验室试验,对澳门地区土层地震波速进行了初步测量,从初步的现场试验结果来看,浅层的地震波速受到堆填物料的影响颇大,变化显得较为宽广,随着深度的增加,由于天然土层较为均一,地震波速的变化亦见收窄,现场剪切波速数据与标准贯入试验(SPT)数据可见一定吻合性.而通过等向及异向固结的实验室试验,发现澳门海泥的剪切波速在同一平均有效应力而不同应力路径下大致相同,此外,比较野外及实验室试验结果,亦显示澳门海泥的剪切波速受土体天然结构影响不大.

摘要： 利用中国地质块体结构和地震波速分布将岩石圈水平分块、垂向分层,反演出中国大陆岩石圈上地幔温度分布,并在中国大陆地区195个地表气象台站提供的地表地温等资料以及居里面深度和捕虏体测得的平衡温度作为深部温度等的约束下,根据三维稳态热传导方程,计算出了中国大陆及邻区岩石圈三维热结构.通过分析,青藏高原下方莫霍面处的典型温度范围为800-1000°C,而前寒武克拉通下方(如印度板块、四川盆地、华南、华北及塔里木)的典型温度范围则位于500-700°C之间.从深度超过60km处的热状态来看,中朝克拉通是破坏了的;青藏高原中部(特别是羌塘地区)的热结构支持下地壳或上地幔物质向东流动的看法.普遍而言,中国大陆内火山分布与地壳或上地幔内的高温区域是一致的.造山带区域存在很多明显的热转换区域,位于中国大陆地壳东西部之间的热转换区与南北地震带向吻合.

摘要： 井地地震CT测试技术是通过地面激振、井中接受的方式,利用地震波对地质体进行透视,通过地震波走时和衰减的观测和反演计算对地质体内部结构成像的物探新技术.由于波的传播路径相互交错,通过成像软件,可获得井周围一定深度范围内的地震波速变化的三维图像,从而查明地层中的岩溶的分布情况.

摘要： 为充分了解不同尺度地质系统的岩石动力学特征，更好利用不同规模勘探资料有效解释地质系统的特征，在岩石物理实验对岩石弹性模量、波速与应力的关系进行研究总结的基础上，探讨野外实际地震勘探中地震波速与地应力的关系。rn 研究表明，岩石、岩体或地质系统是多相(固液气)、多种矿物、各种构造的综合体，不论尺度如何均可以用多重介质孔隙裂隙地质模型描述，而波速是一定尺度不均匀地质体的等效波速，在一定地质条件下，各种尺度介质波速、岩石弹性模量与有效应力存在依存关系。坚硬的岩石通常很脆，如果岩石处于很大的应力之下，应力的释放将在岩石中产生微裂隙，微裂隙会降低地震速度。只要把岩石恢复到原来的应力状态就可能消除所引发的微裂隙。在岩体未达到破坏时，应力越大，波速愈高。压力释放引发的微裂隙有助于现场确定最大/最小水平应力方向，波速的椭圆长轴即为主应力方向。而上覆地层压力不变(或变化很小)，净储层压力的影响正好与储层压力对地震特性的影响相反。随着净储层压力的增加，所有岩石的地震特性都增高。这种增高的量级取决于若干其它的因素(孔隙形状、孔隙度、孔隙流体、岩性等等)。不论野外地震勘探还是实验研究，地震波速与地应力的一致性为地质系统的正反演解释提供了宏观控制依据：通过测量弄清岩石的应力史(它只能通过测量或应力反演方法来量化表示)有助于地震勘探的数据采集规划和解释，反过来，地震波速信息为恢复现代地应力场、储层识别的量化研究提供依据。

摘要： 尽管远离板块边缘,中国东北地区有着广泛的板内火山活动,长白山是其中最大的火山.关于长白山火山的成因,目前至少有两种解释模式:(1)俯冲板块的深部脱水作用(Zhao et al.,2009);(2)上涌热地幔物质的减压熔融作用(Tang et al.,2014).本研究中，从中国东北地区可用的流动和固定台站的波形数据中，采用相关的方法量取相对走时残差，用于远震走时层析成像。采用P和S波联合反演的方法，获得了中国东北地区深达800km的P和S波速度异常结构。最后提出了一个模型来解释长白山火山的成因，长白山火山热物质有两个来源：一部分源自地幔转换带中俯冲板块的脱水作用，另一部分源自下地幔热物质上涌所导致的减压熔融。这两部分物质的共同作用，使得长白山火山成为中国东北地区最大和最活跃的板内火山。

摘要： 东营凹陷北部陡坡带深层发育的巨厚致密砂砾岩体，纵向多期叠置，横向呈裙带扇形分布。扇体相带窄，物性变化快，为低孔低渗储层，其地球物理响应特征表现复杂。通过对十口取心井100余块超声波样品的选取、测定和实验数据的分析，认为影响砂砾岩体地震波传播速度的内因为岩性、孔隙度、孔隙结构和灰质含量；饱和不同流体对岩石纵波速度的影响较小，饱和油、水后弹性模量与干岩样有一定差别，而饱和油和饱和水时的岩石物理参数差异几乎不可分辨。测试分析结果为砂砾岩体储层地震敏感属性的选取奠定了基础。

摘要： 东营凹陷北部陡坡带深层发育的巨厚致密砂砾岩体，纵向多期叠置，横向呈裙带扇形分布。扇体相带窄，物性变化快，为低孔低渗储层，其地球物理响应特征表现复杂。通过对十口取心井100余块超声波样品的选取、测定和实验数据的分析，认为影响砂砾岩体地震波传播速度的内因为岩性、孔隙度、孔隙结构和灰质含量；饱和不同流体对岩石纵波速度的影响较小，饱和油、水后弹性模量与干岩样有一定差别，而饱和油和饱和水时的岩石物理参数差异几乎不可分辨。测试分析结果为砂砾岩体储层地震敏感属性的选取奠定了基础。

摘要： 东营凹陷北部陡坡带深层发育的巨厚致密砂砾岩体，纵向多期叠置，横向呈裙带扇形分布。扇体相带窄，物性变化快，为低孔低渗储层，其地球物理响应特征表现复杂。通过对十口取心井100余块超声波样品的选取、测定和实验数据的分析，认为影响砂砾岩体地震波传播速度的内因为岩性、孔隙度、孔隙结构和灰质含量；饱和不同流体对岩石纵波速度的影响较小，饱和油、水后弹性模量与干岩样有一定差别，而饱和油和饱和水时的岩石物理参数差异几乎不可分辨。测试分析结果为砂砾岩体储层地震敏感属性的选取奠定了基础。

摘要： 东营凹陷北部陡坡带深层发育的巨厚致密砂砾岩体，纵向多期叠置，横向呈裙带扇形分布。扇体相带窄，物性变化快，为低孔低渗储层，其地球物理响应特征表现复杂。通过对十口取心井100余块超声波样品的选取、测定和实验数据的分析，认为影响砂砾岩体地震波传播速度的内因为岩性、孔隙度、孔隙结构和灰质含量；饱和不同流体对岩石纵波速度的影响较小，饱和油、水后弹性模量与干岩样有一定差别，而饱和油和饱和水时的岩石物理参数差异几乎不可分辨。测试分析结果为砂砾岩体储层地震敏感属性的选取奠定了基础。

摘要： 一种基于地震波静力触探试验的复杂土层剪切波速计算方法，步骤如下：获取地震波静力触探试验中的修正锥尖阻力和侧摩阻力沿深度分布数据，并求得土性类型指数沿深度的分布值；计算各参数沿某一深度处的参数值修正为相邻的多个深度处参数值的平均值，获得参数沿深度的平滑分布曲线；综合应用土性类型指数和修正锥尖阻力进行土层划分，应用条件判断将土层分界线划分为四大类；依据具体的试验方法和遵循的规范，对四大类土层分界线进行合并和去除，并添加辅助分界线，完成用于剪切波速计算的土层划分；应用斯涅尔定律、实测的地震波传递时长和新的土体分层，应用正演模拟方法迭代求解每一土层的剪切波速，获得土体剪切波速沿深度的分布曲线。

摘要： 本发明公开了一种利用斜井直达地震波走时计算地层地震波速度的方法，该方法假设激发地震波的位置为炮点S，其二维坐标为(x

摘要： 本发明公开了一种地震岩石的地震波速度频散及衰减预测方法及装置，其中，该方法包括：获取地震岩石的井测数据；基于该井测数据，使用含裂隙的孔隙介质弹性波理论模型，确定地震岩石的饱和岩石体积模量；基于该饱和岩石体积模量，使用周期层状斑块饱和介质的斑块饱和模型，确定等效平面波模量；以及，基于该等效平面波模量确定地震波衰减和速度频散。通过本发明，实现了同时表征微观尺度的孔隙、裂隙结构与中观尺度的流体非均匀分布对地震波衰减的影响。

摘要： 本发明提供一种地震波速度频散及衰减的预测方法、装置、设备及系统，其中方法包括：构建用于计算孔隙岩石频变弹性参数的随机斑块饱和模型，基于储层物性参数确定随机斑块饱和介质对应的复频变纵波模量，构建饱含不同流体的裂缝介质模型，基于背景介质的岩石纵波模量确定垂直于裂缝方向的纵波模量，将随机斑块饱和模型中的随机斑块饱和介质替换裂缝介质模型的背景介质，建立岩石物理理论模型，根据岩石物理理论模型，预测岩石物性参数与流体类型对应的地震波速度频散和衰减响应。本发明实施例提供的预测方法、装置、设备及系统，能够更加准确地预测实际裂缝储层岩石中地震波的速度频散与衰减特性。

摘要： 本发明公开一种地震波速度的确定方法及装置，包括：采集多个地震波；对多个地震波进行频谱分析，得到对应的多个频谱；获取多个地震波对应的多个频散曲线，频散曲线指示频率和速度的关系；根据多个地震波对应的多个频散曲线和多个频谱确定基准频谱高点；根据基准频谱高点，将多个频谱分为正常频谱和非正常频谱；若地震波对应的频谱为正常频谱，获取地震波的频散曲线中第二高的速度值作为地震波的速度，若是地震波对应的频谱为非正常频谱，则获取地震波的频散曲线中最高的速度值作为地震波的速度。

摘要： 本发明公开了一种利用斜井直达地震波走时计算地层地震波速度的方法，该方法假设激发地震波的位置为炮点S，其二维坐标为(x

摘要： 本发明公开了一种地震岩石的地震波速度频散及衰减预测方法及装置，其中，该方法包括：获取地震岩石的井测数据；基于该井测数据，使用含裂隙的孔隙介质弹性波理论模型，确定地震岩石的饱和岩石体积模量；基于该饱和岩石体积模量，使用周期层状斑块饱和介质的斑块饱和模型，确定等效平面波模量；以及，基于该等效平面波模量确定地震波衰减和速度频散。通过本发明，实现了同时表征微观尺度的孔隙、裂隙结构与中观尺度的流体非均匀分布对地震波衰减的影响。

摘要： 本发明提供了一种基于深度学习的隧道地震波速反演方法及系统，本发明为深度学习地震反演方法添加隧道施工现场的环境信息、地质先验信息与隧道探测观测方式信息编码以帮助算法获得相应先验知识，并利用已有数据库对隧道地震波速反演网络进行初步的网络预训练，以辅助算法匹配当前隧道工程探测任务；同时建立网络更新同步隧道开挖的反演优化方式，在隧道边挖边探的同时，随着探测的进行扩充数据库并优化算法，训练方式逐步无监督化，获得更好的泛化性和更准确的反演结果。

摘要： 本发明提出一种基于观测数据自编码的多尺度无监督地震波速反演方法，本发明通过对观测数据自编码以提取数据中的大尺度信息，利用此信息引导反演网络完成对速度模型中不同尺度特征的恢复，降低反演的非线性程度，在此基础上将训练好的观测数据自编码器编码结构嵌入反演网络以完成反演网络前端对地震观测数据信息的有效提取，使得反演网络能更好的解析地震数据包含的信息，更好的建立地震数据与速度模型之间的映射关系，实现了反演方法完全无监督化，并为输入网络的观测数据添加位置编码以辅助网络感知观测系统的布设形式，便于实际工程应用。该方法能在无真实地质模型作网络训练标签的条件下取得地震速度模型较准确的反演效果。

摘要： 本发明提供了一种地震波速度自适应无网格场节点建立方法及系统，属于地震数值正演领域。该方法包括：输入模型上的各个坐标点处的速度值v、震源子波主频fp，该模型在测线方向的最大坐标值xmax和最小坐标值xmin，该模型在深度方向的最大坐标值zmax和最小坐标值zmin，最大场节点间距hmax、最小场节点间距hmin；从该模型的起始点出发，利用一个场节点的坐标点处的速度值获得横测线方向间距和深度方向场节点间距，利用该场节点的坐标值和所述横测线方向间距和深度方向场节点间距获得相邻场节点的坐标值，以此类推，最终获得该模型上的所有场节点的坐标值。利用本发明方法大幅减少了场节点数目，提高了地震波数值模拟计算效率。