共查询到20条相似文献，搜索用时 140 毫秒
1.
本文针对三维复杂介质的两点射线追踪的适用性和有效性,提出了网格模型和层状模型相结合的建模方式,文中利用规则网格化的层状模型对复杂地质单元进行数值描述,一方面避免了用纯粹复杂地质网格体描述时需要的大量内存,另一方面结合运用规则化的层状模型,必有利于提高计算效率,在计算单元自动生成的射线追踪方法基础上,本文结合波前法提出一种快速、有效的两点射线追踪方法,利用反射波波前在地表形成的射线方位角和倾角坐标体系,能快速、准确地收敛到所有从源点到接收点的射线路径,亦可有效地解决射线路径的多值解难题,从而可适用于复杂地质模型的特点. 相似文献
2.
《地球物理学进展》2020,(3)
球形磁性地质体是地质勘探中广泛遇到的基本地质体之一,选择合适的磁异常分析方法对其磁场进行准确正演计算具有重要的意义.为了准确快速地求解球形磁性地质体磁场特征,本文提出了一种基于磁偶极子构造原理求解球形磁性地质体磁场的数值计算方法,对球形磁性地质体按半径由小增大进行多层单元分割,再对每层单元进行块单元分割,将每个块单元视为"磁偶极子",利用基本磁偶极子公式计算了所有块单元在探测平面(线)的磁场大小,同时分析了球体磁场的剖面特征,并将数值解与解析解结果进行了对比和误差分析,最后建立双球组合模型结合Comsol多物理场仿真工具对其仿真对比,验证了磁偶极子构造法求解球体磁场的准确性.研究结果表明:数值解法结果与解析解结果磁异常波动趋势完全一致;不同测点处绝对误差有所差别,但磁异常值均在1 nT以下.本文提出的数值解法无需复杂的数学推导,计算结果稳定可靠. 相似文献
3.
《地球物理学进展》2017,(3)
为了解决复杂速度模型中的走时正、反演问题,例如:含不规则起伏地表、不规则地下波阻抗界面、以及不规则速度异常体的复杂地学模型,本文采用三角网格单元模型参数化下的分区多步改进型最短路径算法,实现了多震相地震射线的追踪计算,结合共轭梯度法求解带约束的阻尼最小二乘反演问题,实现了多震相走时联合同时反演成像的方法技术.当界面起伏较大时出现散射,从而造成散射点所在区域射线密度过密,导致该区域内速度和界面的过度更新.为了克服上述问题,我们在同时反演中引入了射线密度的概念,从而有效地解决了上述过度更新问题.数值模拟实验表明:采用三角网格单元进行模型参数化,可保证在复杂模型中的正演计算具有较高的计算精度;同时反演中可以准确地刻画不规则异常体和不规则反射界面.因此,本文提出的走时成像方法技术具有较广的实用价值. 相似文献
4.
地震直达波走时层析成像可归结为求解一个大型的、稀疏的、常常是病态的线性方程组.求解方程组常用的迭代法,需要一个比较合理的初始猜测解,也即是初始速度模型.初始模型关系到反演的效率甚至成像的正确性.本文在前人研究基础上提出一种生成模型网格节点初始速度方法,假定震源到检波点路径为直线,记录每条射线穿过的单元和统计每个网格单元穿过的射线数目、自动拾取网格节点所在单元的数目等.实例中,由原始数学模型的正演旅行时资料生成节点初始速度模型,效果可以.最后,分别采用均匀模型和本文方法生成的初始模型进行迭代反演,通过比较,证实该自动生成节点初始模型的可行性和可靠性,并对存在的问题进行讨论和解释. 相似文献
5.
工程电法勘探中，常需要探测埋深不大、规模较小、分布复杂的目标体，对探测方法的精度要求较高.由于三维高密度电阻率法数据采集密度大，可对目标体进行多方位观测，是目前工程探测中常用的方法之一.本文对三维高密度电阻率E SCAN法的视电阻率异常特征进行研究.设计几种典型地质体模型，利用有限单元法进行正演计算，并分析其异常特征和分布规律.算例结果表明，E SCAN法对低阻体的分辨能力强于高阻体.无需对所有电极进行扫描供电观测，即可分辨地质体电性特征及水平位置，但观测结果难以分辨地质体纵向延伸. 相似文献
6.
工程电法勘探中，常需要探测埋深不大、规模较小、分布复杂的目标体，对探测方法的精度要求较高.由于三维高密度电阻率法数据采集密度大，可对目标体进行多方位观测，是目前工程探测中常用的方法之一.本文对三维高密度电阻率ESCAN法的视电阻率异常特征进行研究.设计几种典型地质体模型，利用有限单元法进行正演计算，并分析其异常特征和分布规律.算例结果表明，ESCAN法对低阻体的分辨能力强于高阻体.无需对所有电极进行扫描供电观测，即可分辨地质体电性特征及水平位置，但观测结果难以分辨地质体纵向延伸. 相似文献
7.
本文进一步扩展了合成地震图的广义射线方法.由于算法中将地球模型划分成震源区、接收区和地幔区,因此可用于计算震源区与接收区速度结构不同时的体波合成地震图.新算法把广义射线理论与传播矩阵理论结合起来,并采用了分层Q值模型,从而可以自动包括指定介质层内所有多次波的混响.为简便数据准备,设计了射线编码自动生成算法. 本文算法与Kind扩展的反射率法用相同的地球模型进行了数值结果对比,结果表明,矩阵-射线方法用于计算体波地震图时,不但有足够好的精度,而且还具有速度快并适于单个震相研究的优点.在计算由10个接收点组成的地震剖面时,其计算时间仅为反射率法的7％左右. 相似文献
8.
本文进一步扩展了合成地震图的广义射线方法.由于算法中将地球模型划分成震源区、接收区和地幔区,因此可用于计算震源区与接收区速度结构不同时的体波合成地震图.新算法把广义射线理论与传播矩阵理论结合起来,并采用了分层Q值模型,从而可以自动包括指定介质层内所有多次波的混响.为简便数据准备,设计了射线编码自动生成算法. 本文算法与Kind扩展的反射率法用相同的地球模型进行了数值结果对比,结果表明,矩阵-射线方法用于计算体波地震图时,不但有足够好的精度,而且还具有速度快并适于单个震相研究的优点.在计算由10个接收点组成的地震剖面时,其计算时间仅为反射率法的7％左右. 相似文献
9.
直流电阻率数值模拟对于实测电阻率数据反演、解释地下真实电阻率结构具有重要意义.为解决复杂地电模型剖分以及处理截断边界时面临的困难,本文提出一种新的三维自然单元-无限元耦合法.该算法在自然单元区域灵活布设自然节点以填充复杂模型,并引入无限单元替代截断边界条件以促使刚度矩阵与电源位置无关.首先推导了稳定电流场基本方程;然后介绍了二、三维自然单元法及其形函数的构建过程;再引入三维多向映射无限单元并实现自然单元-无限元的有效耦合;最后通过4个数值模型,结合传统有限单元法、有限元-无限元耦合法以及COMSOL Multiphysics软件的计算结果,验证了算法的正确性及有效性.数值结果表明,本文所提算法满足精度要求,能有效减小计算区域,无需二次更新刚度矩阵,对复杂模型适应性强. 相似文献
10.
本文使用改进后的最短路径算法（MSPM）结合分区多步计算技术实现了三维复杂层状起伏介质中的多次透射、反射及转换波波前传播的数值模拟,以及相应走时和射线路径的跟踪计算.其原理是将三维复杂层状模型按速度界面分成若干个独立的计算区域,采用分步计算技术进行多次波的跟踪计算.基于多次波是通过速度界面简单的入射、透射、反射及转换波按一定规律及原理的不同组合,因此可实施分区多步计算技术.数值模拟实例及误差分析表明分区多步计算技术具有单步最短路径算法中的诸多优点:算法简单、数值计算稳健、计算精度高、速度快及全球解等,因此是解决多次波跟踪计算行之有效的方法. 相似文献
11.
针对2008年1月至2011年12月发生在新疆北天山西段的地震事件,分析新源、温泉和巴仑台3个台站记录到的波形资料,采取尾波归一化方法,计算出P、S波在多个频率点f的衰减系数Q-1P、Q-1S,利用最小二乘法拟合得到Q-1=Q0fη,得到该区域P波衰减系数和S波衰减系数。结果表明,新源研究区关于P波、S波的衰减指数为-0.982 1、-0.978 1;温泉研究区关于P波、S波的衰减指数为-1.049 7、-0.773 8;巴仑台研究区关于P波、S波的衰减指数为-0.870 0、-0.956 8。与中国云南地区关于P波、S波衰减结果相比,新源研究区域与云南的禄劝研究区地下介质衰减特性比较一致。 相似文献
12.
It is demonstrated that the blind deconvolution method is fully capable of recovering the unknown Greens function and of estimating the source time functions from observed seismic data of small earthquakes. Based on the assumption of the Gaussian-mixture model of the Greens function, the newly-formulated algorithm is evaluated using synthetic seismic data along with those of the May 8, 1996 Mexico earthquake (Mc = 4.6). Since the estimated results closely match the theoretical input very well, the method is then employed to analyze the source time functions of the July 7, 1995 Pu-Li, Taiwan earthquake (ML = 5.3). The stations triggered by this event were azimuthally well covered. Using the estimated source time functions, information pertaining to the directivity effect is readily obtained, and the actual fault plane of this event is identified, thus clearly indicating that this method provides a most efficient way to estimate the source time function of a small earthquake.Acknowledgment The authors would like to express their thanks to two anonymous reviewers for their valuable suggestions and Dr. I. Santamarias courteous assistance. They also appreciate the efforts of Drs. H.C. Chiu and R.J. Rau, who provided the seismic data and the fault plane solutions. The National Science Council, Taiwan, has supported this research (NSC 91-2119-M-194-011). 相似文献
13.
14.
The Resonance Capacity Method is proposed for the earthquake response analysis of hysteretic structures. Resonance Capacity is a physical quantity of structures which is related to the hysteretic energy absorbed by structures in one cycle and is equated to the acceleration, velocity and displacement amplitudes α0, d0 and d0 of earthquake ground motions at resonance.1 According to the idealized trapezoidal approximation of earthquake ground motions in the logarithmic period–velocity plane as proposed by Veletsos and Newmark,8 the Resonance Capacity property applies in each period range, short, medium and long, where α0, v0 and d0 respectively are approximately constant. In the medium range of periods, the energy dissipated in hysteretic loops and the deformation amplitudes of a single-degree system with elasto–plastic force–deformation relationships are calculated for the case of El Centro 1940, 18 May earthquake, by this Resonance Capacity Method. The result is compared with results from conventional numerical response analyses obtained by Berg and Thomaides,14 Kato and Akiyama12 and Veletsos and Newmark,8 and the general agreement is seen to be good. Therefore, it may be possible to apply this Resonance Capacity Method over the entire range of periods. By means of this method the earthquake response analysis of hysteretic systems can be performed easily, and the hysteretic energy and fatigue characteristics of structures may be taken into account directly, up to the point of fracture. 相似文献
15.
利用华北地区1970年～2009年的ML3.0级以上地震资料,使用相关分析方法,研究了大同及邻区的中小地震空间分布模式,及其在该区中强地震前的异常情况。研究工作首先是确定空间范围、时间窗和时间滑动步长,然后将地震的空间分布进行网格化再转化为序列;对一个时段和它的稍长时段的地震空间分布进行序列化,计算两个序列的相关系数;并分析相关系数在中强地震前是否明显变化。计算结果表明：当空间取以1989年大同-阳高5.9级地震震中为中心,经向长度和纬向长度均为1.6度的矩形区域,时段长度取3600天,稍长时段取3700天,滑动步长取100天的情况下,1980年～2009年的相关系数时间序列数值整体上稳定,在0.94～1之间波动,但在1989年山西大同-阳高5.9级地震、1991年山西大同5.8级地震与1999年山西浑源5.6级地震前2年均出现过相关系数小于0.94的异常现象,表明地震的空间分布上出现明显变化。并讨论了相关系数异常对于中强地震的前兆意义,可以得到以下3点认识：（1）大同及周边地区的中强地震前,中小地震空间分布模式发生显著变化,可作为2年尺度的震兆异常,为该区中期地震预测提供时间上的参考;（2）考虑到研究区的限定范围为震中周边0.8度以内,地震时空分布的相关分析结果对于地震地点的预测也具有指示意义;（3）相关分析可以为识别此类地震空间分布模式的异常变化提供帮助。 相似文献
16.
对SQIP地震学中期预测方法在华北地区的中期预测效果进行了重新评价,回溯性全时空扫描检验报准率R值约为0.67~0.68。做了方法中短期和短期预测指标的提取尝试。结果表明综合参量P值异常区由最大出现收缩后,平均约7个月左右发生未来主震,其中6个月内的震例约占有“收缩”震例的53%~54%,在3个月内的约占33%~35%;部分震例发生在异常区消失后7~8个月内。表明P值异常区“收缩”(或消失)一定程度上可以作为中、强地震的中短期预测标志。 相似文献
17.
18.
Lava flows from Mauna Loa and Huallai volcanoes are a major volcanic hazard that could impact the western portion of the island of Hawaii (e.g., Kona). The most recent eruptions of these two volcanoes to affect Kona occurred in a.d. 1950 and ca. 1800, respectively. In contrast, in eastern Hawaii, eruptions of neighboring Klauea volcano have occurred frequently since 1955, and therefore have been the focus for hazard mitigation. Official preparedness and response measures are therefore modeled on typical eruptions of Klauea.The combinations of short-lived precursory activity (e.g., volcanic tremor) at Mauna Loa, the potential for fast-moving lava flows, and the proximity of Kona communities to potential vents represent significant emergency management concerns in Kona. Less is known about past eruptions of Huallai, but similar concerns exist. Future lava flows present an increased threat to personal safety because of the short times that may be available for responding.Mitigation must address not only the specific characteristics of volcanic hazards in Kona, but also the manner in which the hazards relate to the communities likely to be affected. This paper describes the first steps in developing effective mitigation plans: measuring the current state of peoples knowledge of eruption parameters and the implications for their safety. We present results of a questionnaire survey administered to 462 high school students and adults in Kona. The rationale for this study was the long lapsed time since the last Kona eruption, and the high population growth and expansion of infrastructure over this time interval. Anticipated future growth in social and economic infrastructure in this area provides additional justification for this work.The residents of Kona have received little or no specific information about how to react to future volcanic eruptions or warnings, and short-term preparedness levels are low. Respondents appear uncertain about how to respond to threatening lava flows and overestimate the minimum time available to react, suggesting that personal risk levels are unnecessarily high. A successful volcanic warning plan in Kona must be tailored to meet the unique situation there. 相似文献
19.
A numerical evaluation of the complete Navier-Stokes equations of motion for steady-state, incompressible flow past an infinite circular cylinder is given in terms of the stream function, vorticity, and pressure distribution past such bodies. A method is described which allows use of these flow characteristics: (1) to approximate the characteristics of air flow past hexagonal columnar ice crystals falling under gravity at terminal velocity in air, (2) to compute the trajectory of supercooled cloud drops relative to such ice crystals, and (3) to determine the efficiency with which short columnar ice crystals and needle shaped ice crystals collide with supercooled cloud drops. It is found that for all columnar type ice crystals riming is negligible if the cloud drop size is less than 5 m, and that for riming to commence short columnar crystals must have diameters larger than 50 m, while needle crystals must have diameters larger than 40 m. It is further shown that the collision efficiency cut-offs at the small drop radius and at the large drop radius end of the collision efficiency diagram can be explained on the basis of the cloud drop trajectories for these drop size ranges. 相似文献
20.
A unified approach for the design of high ductility steel frames with concentric braces in the framework of Eurocode 8 
下载免费PDF全文
下载免费PDF全文 Edoardo M. Marino 《地震工程与结构动力学》2014,43(1):97-118
Eurocode 8 (EC8) stipulates design methods for frames with diagonal braces and for chevron braced frames, which differ as regards the numerical model adopted, the value of the behavior factor q and the estimation of the lateral strength provided by braces. Instead, in this paper, the use of the same design method is suggested for both types of concentrically braced frames. The design method is a generalization of the one proposed for chevron braced frames in a previous study. A numerical investigation is conducted to assess the reliability of this design method. A set of concentrically braced frames is designed according to the EC8 and proposed design methods. The seismic response of these frames is determined by nonlinear dynamic analysis. Finally, it is demonstrated that the proposed design method is equivalent to those provided by EC8, because it can ensure the same level of structural safety which would be expected when using EC8. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献