三维波动方程的数值频散关系及其叠前和叠后数值模拟

采用适于三维各向同性介质正演的有限差分格式 ,对声波方程进行了差分离散 ,导出了差分格式的数值频散关系。研究表明 ,三维正演的数值频散主要受差分精度及网格空间步长的影响 ,尤以后者影响更甚。在相同的条件下 ,剖分步长越大及差分精度越低 ,数值频散现象越明显。研究还表明 ,剖分网格有低通滤波作用 ,因而在正演时应使激发源的主要频率远离网格的截止频率 ,以保证精度。在此基础上对某地质模型进行了叠前和叠后数值模拟。计算结果表明 ,叠前模拟能够比较真实地反映地震波在地下介质中的传播规律 ,但计算量过大 ;叠后模拟时 ,爆炸界面法能较好地反映地下反射界面形状 ,但不宜计算较大的模型 ;平面波照射法可用于计算较大的模型 ,且能较好地反映波的动力学特征 ,但不能精确地反映界面形态     

任意偶数阶精度交错网格声波方程数值模拟

以波动理论为理论基础,利用泰勒级数展开推导出了一阶应力-速度声波方程组的空间任意偶数阶精度交错网格差分格式。选用衰减边界条件,进行了边界效果对比及差分精度选取测试。结果表明差分精度越高,频散越弱,数值模拟的效果越好。应用空间八阶、时间二阶精度差分格式,实现了各向同性介质模型的声波方程数值模拟。正演模拟结果可以清晰地反映出声波波场中波传播的运动学和动力学特征,可推广应用于三维声波高精度波场正演模拟中。     

基于紧致交错差分格式的二维声波及黏滞声波方程数值模拟

针对常规网格差分难以适用于地震波场数值模拟中复杂介质的问题,首次将紧致交错有限差分格式应用于黏滞声波方程的数值模拟研究并同声波方程的数值模拟进行了模拟精度、频散关系和稳定性分析等方面的比较。理论研究结果表明：当差分精度相同时,紧致交错网格所需节点数要少于常规的中心差分和交错差分格式,计算效率更高;同常规的交错差分与中心差分格式相比,紧致差分的截断误差更小,数值频散也更低,能够适用于粗网格计算;差分精度相同情况下时进行数值模拟,紧致交错格式所需要的时间网格更小,稳定性条件也更为严格;紧致交错差分格式在完全匹配层(perfectly matched layer, PML)条件下,能够对边界反射进行有效吸收。最后,对均匀、水平层状介质以及Marmousi模型进行了黏滞声波方程的数值模拟和波场特征分析,实验结果证明了该方法对于复杂介质的数值模拟的适应性和有效性,并具有较高的模拟精度及计算效率。     

优化的时空域等效交错网格有限差分正演模拟

地震波有限差分数值模拟中,使用优化的差分系数可以有效地减少数值频散(相速度误差)。为了同时压制空间频散与时间频散,提出一种时空优化的等效交错网格有效差分正演模拟方法。采用以减少给定波数范围内所有传播方向的相对频散误差为目标函数,并借助牛顿法快速获取优化的差分系数。频散分析表明等效交错网格正演在低速介质中容易出现空间频散,在高速介质中高波数分量容易出现时间频散,数值模拟显示该方法既能有效地压制空间频散,也可以减少时间频散。Marmousi模型测试表明该方法适用于复杂地质模型的地震波场正演模拟。相比传统差分方法,该时空优化差分方法在不增加计算量的前提下,可以有效地提高数值模拟精度。     

基于起伏地形的瑞利面波有限差分数值模拟

利用瑞利面波的频散特性及其不同频率瑞利面波能量随探测深度衰减变化规律,能够用来进行浅地表的地质勘探和地球内部构造的研究。利用有限差分数值求解的算法,模拟起伏条件下的地下介质中的瑞利面波;建立速度-应力弹性波动方程,利用交错网格的差分方式,引入完全匹配层(Perfectly matched layer,简称PML)的人工边界条件,对几种简单的起伏地表模型进行正演模拟。结果表明交错网格有限差分方法有比较好的精度,PML边界条件能够对来自边界的虚反射有较好的吸收处理。     

裂缝介质中瑞雷面波传播的渐变非均匀交错网格数值模拟

提出了渐变非均匀交错网格高阶差分方法,模拟裂缝介质中瑞雷面波的传播过程,并引入具有二阶时间精度和六阶空间精度的附加差分公式解决非均匀交错网格的不对称问题。这种差分方法可描述为：在裂缝处采用细网格,在其他地方采用粗网格,细网格到粗网格之间在空间上采用渐变步长,达到既细致地描述了局部细微结构、又保持了计算量适当的目的。在变步长的选取中提出了采用汉宁窗函数,这样既可以较快地从粗网格向细网格过渡,又避免了粗细网格过渡时发生频散,算例表明该方法完全适用于模拟裂缝介质中瑞雷面波的传播,计算精度较高,且能节约时间和内存,数值模拟结果非常直观地揭示了瑞雷面波在裂缝介质中的波场传播特征。     

弹性波传播的低数值频散波场模拟方法

压制数值频散以提高计算精度是检验地震波数值模拟方法的一个重要标准。基于弹性波传播方程,建立了低数值频散波场模拟的八阶FNRK方法。该方法以Runge-Kutta方法对时间导数进行三阶离散,以近似解析离散算子替代差分算子对空间偏导数进行八阶离散,结合通量校正传输技术消除离散后的数值频散。弹性波场模拟结果表明,与高阶有限差分方法相比,该方法能在压制数值频散方面具有明显的优势,计算精度提高,且适应于地震波在大规模复杂介质中传播的波场模拟。     

瑞雷波数值模拟中的数值频散校正策略及实例分析

为寻求快速、高精度的瑞雷波数值模拟方法,将常用于流体动力学中的通量校正传输技术(FCT)与交错网格有限差分方法相结合,对浅层各向同性弹性介质进行包括瑞雷面波和体波在内的全波场模拟。模拟结果表明:该方法有效地压制了粗网格条件下的数值频散现象,并保留了真实的波场振荡,特别是层状介质中瑞雷波的频散特征被凸显;在同一精度下,该方法与细化网格方法相比,计算效率提高1.3倍。     

弹性波方程的紧致差分方法

在对弹性波方程进行数值模拟时 ,低阶差分格式往往产生严重的数值频散 ,高阶显示差分格式需要用较多的网格点 ,不利于边界的处理。而紧致差分格式吸收了它们的优点 ,弥补了它们的不足。为此该文应用紧致差分格式的思想 ,发展了二维情况下弹性波方程初值问题的紧致差分方法 ,研究了它的稳定性 ,并用 Fourier方法分析了显示差分格式和紧致差分格式的相速度误差 ,最后利用紧致差分方法在粗网格条件下对地震波传播进行了数值模拟 ,并同五点四阶中心差分方法的计算结果进行了对比。结果表明 ,求解弹性波方程的紧致差分方法有效 ,且具有比同网格点差分格式更高的计算精度和较小的数值频散。     

三维各向异性介质中弹性波方程交错网格高阶有限差分法数值模拟

三维数值模拟是研究各向异性介质中复杂弹性波波场特征和传播规律的重要手段。根据Taylor展开式,推导出了交错网格一阶空间导数的任意偶数阶精度有限差分近似式和相应的差分系数计算式,给出了三维各向异性介质中弹性波一阶双曲型应力速度方程交错网格任意偶数阶精度差分格式和稳定性条件,并推导出了三维各向异性介质PML法吸收层边界条件公式和相应的交错网格差分格式。对方位各向异性介质模型和正交各向异性介质模型中弹性波的传播进行的三维数值模拟结果表明,弹性波在三维各向异性介质中传播时存在拟P波、拟SV波和拟SH波,并出现横波分裂、横波分裂盲区、波面三分叉等特殊现象,另外,弹性波场在空间是变化的,其拟P波、拟SV波和拟SH波的耦合关系比较复杂。     

共轭梯度法应用于解电磁场的散射和辐射问题

本文将共轭梯度法(CGM)应用于电大尺度(几个波长以))的散射和辐射问题的矩量法中.在二维、TM 情况下,对于用列表法有效地实现节约内贮作了比较详细的讨论;通过多个典型例子加以说明,并和传统方法进行了比较,归纳出应用CGM 的长处和存在的问题.在此基础上,将CGM 应用于解决散射体内部谐振对外影响的边界扩展中,并通过实例加以说明.     

再生式空气预热器的温度场及传热的计算与分析

通过对我国大型锅炉中常用的再生式空气预热器中传热元件的传热过程的分析，建立了它的传热数学模型，通过计算机求解，以《锅炉机组热力计算标准方法》￣［１］中回转式空气预热器的设计计算为例，进行了核算，平均数值符合很好。     

拉性载荷作用下煤和瓦斯突出的失稳机理及数值模拟

本文从变形力学系统平衡的稳定性观点对拉性载荷作用下煤和瓦斯突出的机理及判据进行了研究和讨论。图６、参２。     

反倾层状边坡倾倒变形破坏模式的岩层等厚度特性研究

为检验等厚层状特性这一假定条件在反倾岩质边坡研究中的适用性,以重庆市硝洞槽-郑家大沟反倾层状岩质岸坡为实例进行数值模拟验证。首先建立二维简化工程地质模型;并利用离散元软件(UDEC)对反倾边坡在相邻岩层不同厚度比值(1.0∶0.1)~(1.0∶1.0)条件下分别进行数值模拟计算。然后对比分析不同岩层厚度比值条件下边坡弯曲倾倒变形特征,得出了反倾层状边坡倾倒变形随相邻岩层厚度比值变化规律。最后通过工程实例对研究结做进一步验证。主要得出以下结论:1相邻岩层厚度越相近边坡越易发生弯曲倾倒变形;2仅当相邻岩层厚度比值处于(1.0∶0.8)~(1.0∶1.0)范围内时,反倾边坡才能可视为等厚层状分布。     

FMS物流系统的动态仿真及评价(Ⅰ)——动态仿真及“物流场”概念

首次提出了“物流场”的概念,对动态的物流系统进行计算机离散仿真,探讨了FMS物流系统的研究理论和方法,研究了加工系统的随机性。本文的结果为物流系统的研究提供了一个理论依据和实用综合软件。     

基于模拟数据分析地形校正模型效果及检验

对目前已有的多种地形校正模型(余弦校正模型、C校正模型、Minnaert模型、统计校正模型、SCS校正模型以及SCS C校正模型)的优缺点进行了详细的评价.分析了各种地形校正模型的起源和物理意义.采用数值模拟(Li-Strahler几何光学模型(GOMS))技术生成不同地形条件下的冠层反射率,并采用上述多种模型进行地形影响的消除,从理论上研究分析各种地形校正模型的特点.结果表明统计模型和C校正模型是这些地形校正模型中较好的模型方法.利用结果最优的模型方法对中国江西兴国县的地形进行校正试验,结果表明引入的C校正模型可以较优地进行地形影响的消除.     

成套电器面板计算机辅助设计与制造系统

成套电器面板计算机辅助设计与制造系统ＪＤＰ－ＣＡＤＡＭ是一个ＣＡＤ与ＣＡＭ集成数控自动编程系统，在ＣＡＤ部分，利用人工神经网络，输入细表，可以实现面板的自动布置，并自动形成开孔图，在ＣＡＭ部分，读取面板开孔图后，系统自动选择加工刀具和确定加工路线，自动生成ＮＣ程序并用动画模拟的方法校核其正确性，文章在总体概括ＪＤＰ－ＣＡＤＡＭ系统的基础上，重点介绍了数控自动编程部分的功能及实现方法。     

具有常数移民和急慢性阶段的SIS模型的研究

研究了具有常数移民以及具有急性和慢性两个阶段的SIS传染病模型.针对p=0和0相似文献   

浅埋地铁隧道CRD法与三台阶法施工可行性比选

针对兰州地铁一号线东岗站后配线区间浅埋隧道,对其施工开挖分别采用三台阶法与CRD法进行数值模拟研究。计算结果表明,当采用CRD法施工时,目标断面所在地表最大横向沉降量以及沿隧道轴线方向拱顶最大下沉量分别较三台阶法施工时小7 mm和6.4 mm;衬砌在各选定断面的最小安全系数,除入口断面外,其余均满足规范要求,而三台阶法施工时衬砌只在中间断面能满足要求;两种工法在隧道拱脚处均出现塑性区,但采用CRD法施工时的塑性区面积明显较小。据此建议采用CRD法施工并为同类工程施工提供参考。