基于螺旋线上谱因式分解的地震波场隐式辛算法   总被引：14，自引：8，他引：14       下载免费PDF全文

罗明秋  刘洪  李幼铭 《地球物理学报》2001,44(3):379-388

均匀介质、复杂各向同性介质和各向异性介质中的地震波传播过程,可用统一形式的标量声波方程描述.考虑到在无损耗条件下,地震波方程描述了地震波场这一个无穷维的哈密顿体系随时间的演化过程,该过程为一个单参数连续辛变换,因而可以在其哈密顿形式表述下导出其辛格式.与显式辛算法相比,隐式辛格式对应的隐式辛几何算法具有无条件稳定的特点,可以允许较大的计算步长.但是由于隐式算法不可避免地面临高阶矩阵的求逆,其每一步的计算速度较慢.为实现矩阵快速求逆,文中采用了螺旋边界条件下谱因式分解的方法.在螺旋边界条件下,需要求逆的矩阵化为带状矩阵,而且其各列非零元素的位置和大小具有非常好的相似性,因而可以采用谱因式分解的方法实现快速LU分解.文中采用二阶精度的隐式蛙跳辛格式和谱因式分解方法,计算了常速度、层状介质和Marmousi模型中的波场.计算表明,隐式辛算法不失为波场计算的一种好方法.  相似文献   

各向异性弹性波有限元叠前逆时偏移   总被引：26，自引：10，他引：26       下载免费PDF全文

张美根  王妙月 《地球物理学报》2001,44(5):711-719

利用有限元法和最小走时射线追踪的界面点法，实现了各向异性弹性波的叠前逆时偏移.理论模型资料的偏移结果清晰准确，证实了该系统的可靠性.通过对各向异性模型资料的各向同性偏移处理，发现常规偏移剖面存在较大误差，地质体的垂向深度和横向位置与实际模型有偏离.  相似文献   

广东地区现代构造应力场的模拟计算   总被引：2，自引：0，他引：2  

闻则刚  王正尚  杨马陵 《华南地震》2001,21(3):29-36

运用已知的断裂边界条件、区内震源机制解、小震应力降和环境应力场等资料作为基本数据,建立地质模型和计算模型,用有限元方法模拟计算了华南地区的现代地壳构造应力场.初步得到了本区的外力作用方式、断裂活动情况、应力集中状况、区应力场分布状况,为划分地震危险区段提供了依据.  相似文献   

扁椭球形雨滴和冰粒后向散射实验与理论研究   总被引：2，自引：0，他引：2  

王庆安  刘黎平  巢增明  张尧培 《大气科学》1996,20(6):735-744

本文介绍了对扁椭球形雨滴和冰（雹）粒散射能力的实验模拟和测量技术，报告了一个系列扁水(冰)椭球后向散射截面的实验测量结果；并用扩张边界条件法(EBCM)计算出相应的理论值，用Mie公式计算出等容积球形后向散射截面。将实验结果与理论计算作了综合对比分析，得到了若干有用的结论。  相似文献   

参数化的多次散射雷达方程和它的反演理论 II:反演理论     

邱金桓  吕达仁 《大气科学》1992,16(6):725-736

本文根据参数化的多次散射雷达方程,分析了激光多次散射对激光回波信号和一次散射雷达方程反演结果的效应,提出了迭代前向积分法和迭代后向积分法以求解参数化的多次散射雷达方程,并根据多次散射信号对云和气溶胶的前向散射相函数的高度敏感性提出了一个从双接收视场的多次散射激光回波信号,同时确定云和气溶胶消光系数分布和前向散射相函数的方法.此外,本文还根据数值试验分析了这些反演方法在地基激光雷达和空间站激光雷达两种情形下的适用性.  相似文献   

参数化的多次散射雷达方程和它的反演理论 Ⅱ:反演理论     

邱金桓  吕达仁 《大气科学》1992,(6)

本文根据参数化的多次散射雷达方程,分析了激光多次散射对激光回波信号和一次散射雷达方程反演结果的效应,提出了迭代前向积分法和迭代后向积分法以求解参数化的多次散射雷达方程,并根据多次散射信号对云和气溶胶的前向散射相函数的高度敏感性提出了一个从双接收视场的多次散射激光回波信号,同时确定云和气溶胶消光系数分布和前向散射相函数的方法.此外,本文还根据数值试验分析了这些反演方法在地基激光雷达和空间站激光雷达两种情形下的适用性.  相似文献   

浅谈地下水流模拟中河流的处理方法   总被引：3，自引：0，他引：3  

谭杰  王杰  杨操静 《地下水》2006,28(4):28-30

本文分析了河水与地下水不同作用形式下的补排关系,给出了河流在地下水流模拟中的两种处理方法.根据河流与地下水的相互作用规律将其概化为四种边界条件,并给出了在不同边界条件中河水与地下水交换量的计算公式,初步建立了河水与地下水的解析耦合模型.  相似文献   

一个可供现代数学分析研究的气候动力学模型   总被引：6，自引：0，他引：6  

曾庆存 《大气科学》1998,22(4):408-417

针对以往大气动力学方程和气候动力学方程数学研究所使用模式的不足之处,本文作了改进,使数学模型更接近于实际情况。这些包括,（1）取大气上界气压为零,而非某个正的（小）常数;（2）消除了整层无辐散近似,从而使地面气压为直接的预报量,但在地面气压预报方程中加入了水平湍流扩散项（平滑）;（3）对大气运动的内热源的辐射传输过程和水汽相变过程作了合理的处理,而非看作已知的外强迫,用适当的方法可以用解析函数较好地逼近这些过程。另外还提出合理的边界条件。我们的目的是要证明这样提出的初—边值问题的解的存在性,以后还可研究其基本性质和长期行为。  相似文献   

Tests and applications of an approach to absorbing reflected waves towards incident boundary     

张洪生  王炎  许春辉  商辉  于小伟 《中国海洋工程》2013,27(6):703-718

If the upstream boundary conditions are prescribed based on the incident wave only, the time-dependent numerical models cannot effectively simulate the wave field when the physical or spurious reflected waves become significant. This paper describes carefully an approach to specifying the incident wave boundary conditions combined with a set sponge layer to absorb the reflected waves towards the incident boundary. Incorporated into a time-dependent numerical model, whose governing equations are the Boussinesq-type ones, the effectiveness of the approach is studied in detail. The general boundary conditions, describing the down-wave boundary conditions are also generalized to the case of random waves. The numerical model is in detail examined. The test cases include both the normal one-dimensional incident regular or random waves and the two-dimensional oblique incident regular waves. The calculated results show that the present approach is effective on damping the reflected waves towards the incident wave boundary.  相似文献   

基于多尺度有限元的海上升压站脚靴式桩套筒连接结构分析研究     

张伟  张坤  姜贞强  孙震洲  王永发  吕国儿  陈杰峰 《海洋工程》2022,40(2):55-66

脚靴式海上升压站导管架的桩套筒与主导管的连接结构是关键荷载传递部位。整体梁系有限元模型无法反应局部板型结构的失效状态;截断板系有限元模型则受制于边界条件的合理性,结果不够精确。文中基于某工程海上升压站,建立不同精细化范围的多尺度有限元模型及不同边界条件的截断模型,对结构屈服失效及屈曲失效复核的同时论证了局部网格对计算精细程度的影响及截断模型对边界条件的敏感性;此外,研究了轭板板厚及应力比对屈曲失效的敏感程度,为今后海上升压站脚靴结构不同阶段的设计、优化提供参考和借鉴。  相似文献