CN115629418B - 一种区域地震反应谱场构造方法 - Google Patents

一种区域地震反应谱场构造方法 Download PDF

Info

Publication number
CN115629418B
CN115629418B CN202211525226.4A CN202211525226A CN115629418B CN 115629418 B CN115629418 B CN 115629418B CN 202211525226 A CN202211525226 A CN 202211525226A CN 115629418 B CN115629418 B CN 115629418B
Authority
CN
China
Prior art keywords
point
spectrum
points
measuring
measuring points
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.)
Active
Application number
CN202211525226.4A
Other languages
English (en)
Other versions
CN115629418A (zh
Inventor
田源
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Science and Technology Beijing USTB
Original Assignee
University of Science and Technology Beijing USTB
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Science and Technology Beijing USTB filed Critical University of Science and Technology Beijing USTB
Priority to CN202211525226.4A priority Critical patent/CN115629418B/zh
Publication of CN115629418A publication Critical patent/CN115629418A/zh
Application granted granted Critical
Publication of CN115629418B publication Critical patent/CN115629418B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种区域地震反应谱场构造方法，包括：获取区域每一地震动测点的地震反应谱和位置；针对各区域内测点，分别基于其与各周边测点的距离及所有周边测点的地震反应谱，计算其所在位置在每一选定特征周期点处的试探反应谱；然后计算每一区域内测点的实际反应谱与试探反应谱在每个特征周期点处的比值；针对当前待构造反应谱的点，利用其与周围测点的位置关系及其周围测点对应的上述比值，确定各特征周期点处的反应谱加权系数；基于上述结果，针对当前待构造反应谱的点，计算其地震动反应谱；遍历区域内所有待构造反应谱的点，得到区域内地震反应谱场。本发明可直接给出地形相对平坦区域的地震反应谱场，且兼顾了局部小区域的场地特征。

Description

一种区域地震反应谱场构造方法
技术领域
本发明涉及土木工程结构抗震技术领域，特别涉及一种区域地震反应谱场构造方法。
背景技术
地震反应谱对于地震破坏力的初步评价具有重要价值。对于不同类型的结构，所关注的反应谱也存在一定区别。例如，研究人员一般通过地震加速度反应谱来判断地震对于低层结构的破坏力，通过地震速度和位移反应谱判断地震对于高层结构的破坏力。然而，地震发生后，快速的获取受灾区域的地震反应谱场面临诸多挑战。这主要是由于现实的强震台网密度不足以提供精细化的地震反应谱场，难以反映大面积范围内各处建筑所在区域的局部场地情况，导致评估结果出现偏差。因此，如何简单有效的构造发震区域的地震反应谱场对震害评估和震后应急具有重要意义。但目前尚缺乏相关的模型。
发明内容
本发明提供了一种区域地震反应谱场构造方法，以实现区域震后精细地震反应谱场的构造，进而为地震破坏力的快速评估提供重要数据，为震后应急提供更加精细、科学的依据。
为解决上述技术问题，本发明提供了如下技术方案：
一方面，本发明提供了一种区域地震反应谱场构造方法，包括：
获取待构造地震反应谱场区域的每一地震动测点的地震动记录和位置坐标，并基于测点的地震动记录，计算每一测点的地震反应谱；其中，所述地震动测点包括周边测点和区域内测点，所述周边测点指的是区域范围外轮廓附近的测点，所述区域内测点指的是被所述周边测点包围在内的测点；
基于测点的位置坐标，分别计算每一区域内测点与各周边测点的距离；
针对每一区域内测点，分别基于其与各周边测点的距离及所有周边测点的地震反应谱，计算其所在位置在每一选定特征周期点处的试探反应谱；然后计算每一区域内测点的实际反应谱与试探反应谱在每个特征周期点处的比值；
针对当前待构造反应谱的点，在其周围选定多个测点作为其周围测点，并利用其与周围测点的位置关系以及其周围测点的实际反应谱与试探反应谱在每个特征周期点处的比值，确定当前待构造反应谱的点在每个特征周期点处的反应谱加权系数；
针对当前待构造反应谱的点，基于其所在位置在每一选定特征周期点处的试探反应谱以及每个特征周期点处的反应谱加权系数，得到其地震动反应谱；
遍历区域内所有待构造反应谱的点，得到各点的地震动反应谱，以构造出区域内地震反应谱场。
进一步地，所述针对每一区域内测点，分别基于其与各周边测点的距离及所有周边测点的地震反应谱，计算其所在位置在每一选定特征周期点处的试探反应谱，包括：
通过以下公式得到每一区域内测点对应的试探反应谱：
Figure 284766DEST_PATH_IMAGE001
其中，
Figure 238946DEST_PATH_IMAGE002
表示第
Figure 281464DEST_PATH_IMAGE003
号测点所在位置在特征周期
Figure 459636DEST_PATH_IMAGE004
处的试探反应谱值；
Figure 410274DEST_PATH_IMAGE005
为 周边测点的总数；
Figure 964883DEST_PATH_IMAGE006
表示第
Figure 255050DEST_PATH_IMAGE007
号周边测点在特征周期
Figure 810315DEST_PATH_IMAGE008
处的实际反应谱值；
Figure 185933DEST_PATH_IMAGE009
表示 第
Figure 278654DEST_PATH_IMAGE010
号区域内测点与第
Figure 688907DEST_PATH_IMAGE007
号周边测点的距离；
Figure 474460DEST_PATH_IMAGE011
为预设指数。
进一步地，
Figure 337374DEST_PATH_IMAGE011
取值范围为1~2。
进一步地，所述针对当前待构造反应谱的点，在其周围选定多个测点作为其周围测点，并利用其与周围测点的位置关系以及其周围测点的实际反应谱与试探反应谱在每个特征周期点处的比值，确定当前待构造反应谱的点在每个特征周期点处的反应谱加权系数，包括：
针对当前待构造反应谱的点P，确定其周围最近的至少三个测点，作为其周围测点；
利用以下公式计算点P各周围测点在每个特征周期处的反应谱加权系数：
Figure 296103DEST_PATH_IMAGE012
其中，
Figure 560862DEST_PATH_IMAGE013
表示点P的第k个周围测点在特征周期
Figure 779966DEST_PATH_IMAGE004
处的反应谱加权系数；
Figure 130176DEST_PATH_IMAGE014
表示当点P周围第k个测点属于区域内测点时，第k个测点所在位置的实际反应 谱与试探反应谱在
Figure 564699DEST_PATH_IMAGE004
处的比值；
利用点P与周围测点的位置关系确定点P在每个特征周期点处的反应谱加权系数。
进一步地，所述利用点P与周围测点的位置关系确定点P在每个特征周期点处的反应谱加权系数，包括：
取所有周围测点的反应谱加权系数的几何平均值或所有周围测点的反应谱加权 系数的加权几何平均值，作为点P在特征周期
Figure 949544DEST_PATH_IMAGE008
处的反应谱加权系数。
进一步地，对于当前待构造反应谱的点，其周围测点形成的外包络能够包围当前待构造反应谱的点。
进一步地，针对当前待构造反应谱的点，基于其所在位置在每一选定特征周期点处的试探反应谱以及每个特征周期点处的反应谱加权系数，得到其地震动反应谱，包括：
对于当前待构造反应谱的点P，采用以下公式计算其地震动反应谱：
Figure 76900DEST_PATH_IMAGE015
其中，
Figure 179986DEST_PATH_IMAGE016
表示点P在特征周期
Figure 418200DEST_PATH_IMAGE008
处的地震动反应谱值；
Figure 655339DEST_PATH_IMAGE017
表示点P在 特征周期
Figure 536619DEST_PATH_IMAGE008
处的反应谱加权系数；N表示周边测点的总数；
Figure 119478DEST_PATH_IMAGE006
表示第J号周边测点在特 征周期
Figure 833487DEST_PATH_IMAGE008
处的实际反应谱值；
Figure 927345DEST_PATH_IMAGE018
表示点P与第J号周边测点的距离；
Figure 662083DEST_PATH_IMAGE019
为预设指数。
进一步地，遍历区域内所有待构造反应谱的点，得到各点的地震动反应谱，以构造出区域内地震反应谱场，包括：
对指定区域进行预设间距的离散；
遍历离散后的区域内所有的离散点，得到各点的地震动反应谱；利用所有点的地震动反应谱，构造出区域内地震反应谱场。
再一方面，本发明还提供了一种电子设备，其包括处理器和存储器；其中，存储器中存储有至少一条指令，所述指令由处理器加载并执行以实现上述方法。
又一方面，本发明还提供了一种计算机可读存储介质，所述存储介质中存储有至少一条指令，所述指令由处理器加载并执行以实现上述方法。
本发明提供的技术方案带来的有益效果至少包括：
本发明提供了一种区域地震反应谱场构造方法，利用现有的可获取的地震记录实现整片区域的震后精细地震反应谱场的构造，能够直接给出地形相对平坦区域的地震反应谱场，且兼顾了局部小区域的场地特征；本发明能够为地震破坏力的快速评估提供重要数据，为震后应急提供更加精细、科学的依据。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的区域地震反应谱场构造方法的执行流程示意图；
图2是测点分布示意图；
图3是四个周边测点的速度反应谱对比图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。
第一实施例
针对区域完整地震反应谱场构造的问题，本实施例提供了一种区域地震反应谱场构造方法。该方法的执行流程如图1所示，包括以下步骤：
S1，获取待构造地震反应谱场区域的每一地震动测点的地震动记录和位置坐标，并基于测点的地震动记录，计算每一测点的地震反应谱；
其中，本实施例将所有地震动测点分为两类：区域周边测点和区域内测点。其中，区域周边测点是指区域范围外轮廓附近的测点，这些测点处于所有测点的最外侧。区域内测点是指被周边测点包围在内的测点；
本实施例获取所有测点的分布情况并计算每个测点的地震动记录对应的反应谱，这里的反应谱可以是加速度反应谱、速度反应谱和位移反应谱中的任意一种，计算过程中所采用的阻尼比根据实际情况自定，整个过程中保持一致即可，本实施例并不限定具体的反应谱类型以及阻尼比的选取方案。
S2，基于测点的位置坐标，分别计算每一区域内测点与各周边测点的距离；
S3，针对每一区域内测点，分别基于其与各周边测点的距离及所有周边测点的地震反应谱，计算其所在位置在每一选定特征周期点处的试探反应谱；然后计算每一区域内测点的实际反应谱与试探反应谱在每个特征周期点处的比值；
具体地，在本实施例中，上述S3包括：
S31，根据分析需求，选定反应谱特征周期点；其中，所选定的反应谱特征周期点具 体数目根据需求自定，为方便说明，定义第s个特征周期点为
Figure 208602DEST_PATH_IMAGE008
；
S32，利用所有周边测点的地震反应谱，针对每个区域内测点所在位置的试探反应谱，采用如下加权函数计算其在每一选定特征周期点处的试探反应谱：
Figure 54198DEST_PATH_IMAGE001
其中，
Figure 999633DEST_PATH_IMAGE002
表示第
Figure 905272DEST_PATH_IMAGE003
号测点所在位置在特征周期
Figure 939087DEST_PATH_IMAGE004
处的试探反应谱值；N为周 边测点的总数；
Figure 322795DEST_PATH_IMAGE006
表示第J号周边测点在特征周期
Figure 922404DEST_PATH_IMAGE008
处的实际反应谱值；
Figure 671048DEST_PATH_IMAGE009
表示第
Figure 192159DEST_PATH_IMAGE003
号区域内测点与第J号周边测点的距离；
Figure 382488DEST_PATH_IMAGE011
为预设指数，一般取1~2。
S33，针对每个区域内测点，按以下形式计算其实际反应谱与试探反应谱在每个特征周期点的比值：
Figure 305444DEST_PATH_IMAGE020
其中，
Figure DEST_PATH_IMAGE022AA
表示第
Figure 716003DEST_PATH_IMAGE003
号测点在特征周期
Figure 84929DEST_PATH_IMAGE004
处的实际反应谱值；
Figure 751053DEST_PATH_IMAGE023
表示第
Figure 528516DEST_PATH_IMAGE010
号测点所在位置在特征周期
Figure 946859DEST_PATH_IMAGE004
处的试探反应谱值；
Figure 331310DEST_PATH_IMAGE024
表示第
Figure 80086DEST_PATH_IMAGE003
号测点处实际反应 谱与试探反应谱在特征周期
Figure 703266DEST_PATH_IMAGE004
处的比值。
S4，针对当前待构造反应谱的点，在其周围选定多个测点作为其周围测点，并利用其与周围测点的位置关系以及其周围测点的实际反应谱与试探反应谱在每个特征周期点处的比值，确定当前待构造反应谱的点在每个特征周期点处的反应谱加权系数；
具体地，在本实施例中，上述S4包括：
S41，针对当前待构造反应谱的点P，确定其周围最近的至少三个测点，作为其周围测点；其中，周围测点形成的外包络能够包含当前待构造反应谱的点；选定的这些周围测点可能包含两种类型，一种属于区域周边测点；另一种属于区域内测点。
S42，利用以下形式计算这些周围测点在每个特征周期处的加权系数：
Figure 292511DEST_PATH_IMAGE025
其中，
Figure 229505DEST_PATH_IMAGE013
表示点P的第k个周围测点在特征周期
Figure 293888DEST_PATH_IMAGE004
处的反应谱加权系数；
Figure 452468DEST_PATH_IMAGE014
表示当点P周围第k个测点属于区域内测点时，第k个测点所在位置的实际反应谱 与试探反应谱在
Figure 884718DEST_PATH_IMAGE008
处的比值；
S43，利用P与周围测点的位置关系确定每个特征周期点的反应谱加权系数。
为方便说明，这里定义点P处在特征周期
Figure 823855DEST_PATH_IMAGE008
处的反应谱加权系数为
Figure 917363DEST_PATH_IMAGE026
。这里 至少可以有以下两大类处理方式：
(1) 直接平均法：取所有周围测点加权系数的几何平均值，即：
Figure 727187DEST_PATH_IMAGE027
其中，
Figure 923813DEST_PATH_IMAGE028
表示点P周围选定的测点总数。
(2) 加权平均法：取所有周围测点加权系数的加权几何平均值，加权值可以是P点 至各测点的距离
Figure 488262DEST_PATH_IMAGE029
次方的倒数，即：
Figure 569482DEST_PATH_IMAGE030
其中，
Figure 30550DEST_PATH_IMAGE031
表示点P与周围第k个测点的距离；
Figure 866919DEST_PATH_IMAGE029
可根据实际情况选取，建议一般为 1或2。
除上述方法外，任何适用于小范围内标量插值的方法都可以尝试采用，这里不再赘述。
S5，针对当前待构造反应谱的点，基于其所在位置在每一选定特征周期点处的试探反应谱以及每个特征周期点的反应谱加权系数，得到其地震动反应谱；
具体地，在本实施例中，上述S5的实现过程为：针对区域内的任意待构造反应谱的点，将S3中指定的加权函数与每个特征周期点的反应谱加权系数进行组合，计算该点的地震动反应谱。对于区域内任意点P，可用以下公式计算其地震动反应谱：
Figure 46228DEST_PATH_IMAGE015
其中，
Figure 335999DEST_PATH_IMAGE016
表示通过本实施例方法构造的P点在特征周期
Figure 385995DEST_PATH_IMAGE008
处的地震动反应谱 值；
Figure 924424DEST_PATH_IMAGE017
表示点P处在特征周期
Figure 325449DEST_PATH_IMAGE008
处的反应谱加权系数；N表示周边测点的总数；
Figure 810788DEST_PATH_IMAGE006
表示第J号周边测点在特征周期
Figure 449711DEST_PATH_IMAGE008
处的实际反应谱值；
Figure 424621DEST_PATH_IMAGE018
表示点P与第J号周边测点的距 离；
Figure 247696DEST_PATH_IMAGE019
为预设指数，一般取1~2。
S6，遍历区域内所有待构造反应谱的点，得到各点的地震动反应谱，以构造出区域内地震反应谱场。
具体地，在本实施例中，上述S6包括：
S61，对指定区域进行适当间距的离散；
S62，遍历离散后的区域内所有的离散点，按照S4~S5得到各点的地震动反应谱；利用所有点的地震动反应谱，构造出区域内地震反应谱场。
本实施例的区域地震反应谱场构造方法，解决了区域完整地震反应谱场构造的问题，能够为地震破坏力的快速评估提供重要数据，为震后应急提供更加精细、科学的依据。下面结合具体的案例对本实施例的方法进行详细的阐述。
某区域测点分布及测点编号如图2所示。首先选定圆点处测点为区域周边测点，方点处测点为区域内测点，在本案例中，区域周边测点共有4个，区域内测点共有2个。另外设置一个测点用于后续说明本方法精度，用三角点标记。获取每个测点的地震动记录后计算地震反应谱，本案例以速度反应谱为例，4个周边测点与2个区域内测点的反应谱结果如图3所示。
针对每个区域内测点，分别计算其与所有周边测点的距离，得到表1结果：
表1 每个区域内测点与所有周边测点的距离
Figure 536726DEST_PATH_IMAGE032
为简单起见，本案例仅选定3个反应谱特征周期点，分别为1s、2s和4s。利用4个周边测点的地震反应谱根据S3中给出的公式，将r指数取为2，计算每个区域内测点位置的试探速度反应谱，结果如表2所示：
表2 每个区域内测点位置的试探速度反应谱
Figure 233417DEST_PATH_IMAGE033
针对每个区域内测点，计算实际反应谱与试探反应谱在每个特征周期点的比值，结果如表3所示：
表3 实际反应谱与试探反应谱在每个特征周期点的比值
Figure 379228DEST_PATH_IMAGE034
针对区域内的任意点，这里为说明方法精度，选取图2中的三角点P。选定的4个测点：边-2、边-3、内-1、内-2。分别计算P点与这4个点的距离为：0.8579、0.8218、0.5042、0.3874km。这4个点在3个特征周期点的加权系数如表4所示：
表4 加权系数
Figure 754845DEST_PATH_IMAGE035
这里分别采用直接平均法和加权平均法两种方法计算最终的加权系数，公式参见 S5，其中加权平均法的指数
Figure 319338DEST_PATH_IMAGE029
取为2。采用两种方法计算得到的最终加权系数结果如表5所 示：
表5 最终加权系数结果
Figure 870536DEST_PATH_IMAGE036
计算点P与4个周边测点（边-1~边-4）的距离，分别为：1.1869、0.8579、0.8218、 1.0877km。然后，采用S6中的公式计算，将
Figure 656089DEST_PATH_IMAGE011
指数取为2，最终在每个特征周期点的地震动反 应谱，计算结果与实际测点结果对比如表6所示：
表6 地震动反应谱计算结果与实际测点结果对比
Figure 659948DEST_PATH_IMAGE037
由上可见，本实施实例的方法具有较高的精度和可靠性。通过重复上述流程并遍历区域内所有点，最终可以得到区域内地震反应谱场，在此不再赘述。
第二实施例
本实施例提供一种电子设备，其包括处理器和存储器；其中，存储器中存储有至少一条指令，所述指令由处理器加载并执行，以实现第一实施例的方法。
该电子设备可因配置或性能不同而产生比较大的差异，可以包括一个或一个以上处理器（central processing units，CPU）和一个或一个以上的存储器，其中，存储器中存储有至少一条指令，所述指令由处理器加载并执行上述方法。
第三实施例
本实施例提供一种计算机可读存储介质，该存储介质中存储有至少一条指令，所述指令由处理器加载并执行，以实现上述第一实施例的方法。其中，该计算机可读存储介质可以是ROM、随机存取存储器、CD-ROM、磁带、软盘和光数据存储设备等。其内存储的指令可由终端中的处理器加载并执行上述方法。
此外，需要说明的是，本发明可提供为方法、装置或计算机程序产品。因此，本发明实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本发明实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质上实施的计算机程序产品的形式。
本发明实施例是参照根据本发明实施例的方法、终端设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、嵌入式处理机或其他可编程数据处理终端设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理终端设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理终端设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。这些计算机程序指令也可装载到计算机或其他可编程数据处理终端设备上，使得在计算机或其他可编程终端设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程终端设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
还需要说明的是，在本文中，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者终端设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者终端设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者终端设备中还存在另外的相同要素。
最后需要说明的是，以上所述是本发明优选实施方式，应当指出，尽管已描述了本发明优选实施例，但对于本技术领域的技术人员来说，一旦得知了本发明的基本创造性概念，在不脱离本发明所述原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。所以，所附权利要求意欲解释为包括优选实施例以及落入本发明实施例范围的所有变更和修改。

Claims (5)

1.一种区域地震反应谱场构造方法，其特征在于，包括：
获取待构造地震反应谱场区域的每一地震动测点的地震动记录和位置坐标，并基于测点的地震动记录，计算每一测点的地震反应谱；其中，所述地震动测点包括周边测点和区域内测点，所述周边测点指的是区域范围外轮廓附近的测点，所述区域内测点指的是被所述周边测点包围在内的测点；
基于测点的位置坐标，分别计算每一区域内测点与各周边测点的距离；
针对每一区域内测点，分别基于其与各周边测点的距离及所有周边测点的地震反应谱，计算其所在位置在每一选定特征周期点处的试探反应谱；然后计算每一区域内测点的实际反应谱与试探反应谱在每个特征周期点处的比值；
针对当前待构造反应谱的点，在其周围选定多个测点作为其周围测点，并利用其与周围测点的位置关系以及其周围测点的实际反应谱与试探反应谱在每个特征周期点处的比值，确定当前待构造反应谱的点在每个特征周期点处的反应谱加权系数；
针对当前待构造反应谱的点，基于其所在位置在每一选定特征周期点处的试探反应谱以及每个特征周期点处的反应谱加权系数，得到其地震动反应谱；
遍历区域内所有待构造反应谱的点，得到各点的地震动反应谱，以构造出区域内地震反应谱场；
所述针对每一区域内测点，分别基于其与各周边测点的距离及所有周边测点的地震反应谱，计算其所在位置在每一选定特征周期点处的试探反应谱，包括：
通过以下公式得到每一区域内测点对应的试探反应谱：
Figure QLYQS_1
其中，
Figure QLYQS_2
表示第i号测点所在位置在特征周期
Figure QLYQS_3
处的试探反应谱值；N为周边测点的总数；
Figure QLYQS_4
表示第J号周边测点在特征周期
Figure QLYQS_5
处的实际反应谱值；d i,J 表示第i号区域内测点与第J号周边测点的距离；r为预设指数；
所述针对当前待构造反应谱的点，在其周围选定多个测点作为其周围测点，并利用其与周围测点的位置关系以及其周围测点的实际反应谱与试探反应谱在每个特征周期点处的比值，确定当前待构造反应谱的点在每个特征周期点处的反应谱加权系数，包括：
针对当前待构造反应谱的点P，确定其周围最近的至少三个测点，作为其周围测点；
利用以下公式计算点P各周围测点在每个特征周期处的反应谱加权系数：
Figure QLYQS_6
其中，
Figure QLYQS_7
表示点P的第k个周围测点在特征周期
Figure QLYQS_8
处的反应谱加权系数；
Figure QLYQS_9
表示当点P周围第k个测点属于区域内测点时，第k个测点所在位置的实际反应谱与试探反应谱在
Figure QLYQS_10
处的比值；
利用点P与周围测点的位置关系确定点P在每个特征周期点处的反应谱加权系数；
针对当前待构造反应谱的点，基于其所在位置在每一选定特征周期点处的试探反应谱以及每个特征周期点处的反应谱加权系数，得到其地震动反应谱，包括：
对于当前待构造反应谱的点P，采用以下公式计算其地震动反应谱：
Figure QLYQS_11
其中，
Figure QLYQS_12
表示点P在特征周期
Figure QLYQS_13
处的地震动反应谱值；
Figure QLYQS_14
表示点P在特征周期
Figure QLYQS_15
处的反应谱加权系数；N表示周边测点的总数；
Figure QLYQS_16
表示第J号周边测点在特征周期
Figure QLYQS_17
处的实际反应谱值；d P,J 表示点P与第J号周边测点的距离；r为预设指数。
2.如权利要求1所述的区域地震反应谱场构造方法，其特征在于，r取值范围为1~2。
3.如权利要求1所述的区域地震反应谱场构造方法，其特征在于，所述利用点P与周围测点的位置关系确定点P在每个特征周期点处的反应谱加权系数，包括：
取所有周围测点的反应谱加权系数的几何平均值或所有周围测点的反应谱加权系数的加权几何平均值，作为点P在特征周期
Figure QLYQS_18
处的反应谱加权系数。
4.如权利要求1所述的区域地震反应谱场构造方法，其特征在于，对于当前待构造反应谱的点，其周围测点形成的外包络能够包围当前待构造反应谱的点。
5.如权利要求1所述的区域地震反应谱场构造方法，其特征在于，遍历区域内所有待构造反应谱的点，得到各点的地震动反应谱，以构造出区域内地震反应谱场，包括：
对指定区域进行预设间距的离散；
遍历离散后的区域内所有的离散点，得到各点的地震动反应谱；利用所有点的地震动反应谱，构造出区域内地震反应谱场。
CN202211525226.4A 2022-12-01 2022-12-01 一种区域地震反应谱场构造方法 Active CN115629418B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211525226.4A CN115629418B (zh) 2022-12-01 2022-12-01 一种区域地震反应谱场构造方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211525226.4A CN115629418B (zh) 2022-12-01 2022-12-01 一种区域地震反应谱场构造方法

Publications (2)

Publication Number Publication Date
CN115629418A CN115629418A (zh) 2023-01-20
CN115629418B true CN115629418B (zh) 2023-03-17

Family

ID=84910842

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211525226.4A Active CN115629418B (zh) 2022-12-01 2022-12-01 一种区域地震反应谱场构造方法

Country Status (1)

Country Link
CN (1) CN115629418B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102680193A (zh) * 2012-05-16 2012-09-19 东南大学 一种结构抗震分析的区间反应谱法
CN110927783A (zh) * 2019-11-14 2020-03-27 北京科技大学 一种二维场地二维地震反应谱绘制方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI660190B (zh) * 2017-03-15 2019-05-21 國立臺灣科技大學 震後診斷方法及震後診斷元件
JP2019109065A (ja) * 2017-12-15 2019-07-04 大成建設株式会社 設計用地震動の設定方法
CN109308373A (zh) * 2018-08-10 2019-02-05 清华大学 考虑场地-城市效应的区域建筑地震反应分析方法及装置
CN111999771B (zh) * 2020-08-27 2021-03-26 四川省地震局减灾救助研究所 一种实现地震烈度快速评估及自动修正的方法及系统
CN114442153B (zh) * 2022-01-30 2023-11-28 福州大学 一种近断层地震动拟合方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102680193A (zh) * 2012-05-16 2012-09-19 东南大学 一种结构抗震分析的区间反应谱法
CN110927783A (zh) * 2019-11-14 2020-03-27 北京科技大学 一种二维场地二维地震反应谱绘制方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
师黎静.区域场地近地表速度结构建模研究.2022,第44卷(第44期),360-367. *

Also Published As

Publication number Publication date
CN115629418A (zh) 2023-01-20

Similar Documents

Publication Publication Date Title
Rodríguez‐Iturbe et al. The design of rainfall networks in time and space
CN103336093B (zh) 一种区域空间质量分析方法
JP6001740B1 (ja) 構造物の伝達関数の高精度評価・地震応答予測・劣化診断システム及びその方法
CN108243425A (zh) 一种待建基站的位置确定方法及装置
CN104361157B (zh) 一种建筑物间风环境的评价方法
Schweitzer et al. Modeling of quality for engineering geodesy processes in civil engineering
Park et al. Model updating method for damage detection of building structures under ambient excitation using modal participation ratio
CN102012522A (zh) 地震采集数据的评价方法及装置
JP5281475B2 (ja) 常時微動計測に基づく建物の健全性診断法、診断装置及び診断プログラム
CN113340286B (zh) 一种土地规划勘察勘测项目测绘信息数据分析方法、设备及计算机存储介质
CN115629418B (zh) 一种区域地震反应谱场构造方法
CN109859323B (zh) 一种基于三角网模型林分空间格局加权的方法
JP2003296396A (ja) 建物の期待ライフサイクルコスト評価システムおよび期待ライフサイクルコスト評価プログラムを記録した記録媒体
JP3017154B2 (ja) 地震被害推定システム
CN114997028B (zh) 一种基于改进稳定图的空间网格结构省时随机子空间方法
JP2008309632A (ja) 浸水深現地調査補助システム
CN103376465A (zh) 建筑物楼层的地震即时分析系统及其方法与储存媒体
CN116522085A (zh) 一种全自动拉索频率提取、定阶和索力识别方法与应用
US20220173821A1 (en) Machine learning apparatus
JP2011501171A (ja) 放射線量の定量方法及び関連する等線量曲線を求める方法
CN109933921B (zh) 一种滚石灾害风险评估方法、装置、系统及存储介质
CN108304653B (zh) 晶体培植工艺环境振动影响定量快速评估技术
Weku et al. Optimal Bandwidth for Geographically Weighted Regression to Model the Spatial Dependency of Land Prices in Manado, North Sulawesi Province, Indonesia
Hellmanns et al. A planning tool for the automated quantification and visualization of blue space
JP2023065904A (ja) 海洋観測計画立案支援プログラム、海洋観測計画立案支援装置、海洋観測計画立案支援方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant