共查询到19条相似文献，搜索用时 187 毫秒
1.
伴随着复电阻率法的广泛应用,发展精确和快速的正演和反演算法成为复电阻率法研究的重点.本文采用基于三角单元剖分的有限单元法进行了复电阻率二维数值模拟研究.为了提高计算速度,对无穷远边界进行了近似处理.整个正演计算过程分为两步,首先采用有限单元法计算四个不同频率的视复电阻率数据,然后对前一步得到的视复电阻率数据采用递推算法计算视Cole-Cole参数.采用这种正演算法与一维正演的结果进行了对比,验证了本文方法的正确性.设计了两个二维极化模型,数值模拟结果表明视复电阻率和Cole-Cole视参数等值线断面图对于异常目标体都有比较明显的反映. 相似文献
2.
渤海唐山地区的地下热结构及地表热流分布的数值模拟 总被引：1，自引：0，他引：1
本文采用空间有限元和时间有限差分法模拟了渤中盆地区域的演化历史和地下热结构,并分别计算了渤中上地幔的热隆起与浅部的沉积层对渤中盆地内部温度场的影响。文中得到的计算结果较好地解释了渤中盆地地区热流正异常的观测事实,据此,对渤中盆地在形成过程中产生的热应力与地震孕震力源的关系也进行了一些探讨 相似文献
3.
4.
孔隙压力对含随机缺陷岩石破坏过程及全部变形特征的影响 总被引：2，自引：0，他引：2
对于平面应变压缩条件下含有随机缺陷的岩样,利用FLAC研究了孔隙压力对岩样破坏过程、全部变形特征及前兆的影响。以前编写的若干FISH函数,被用于生成缺陷和计算轴向、侧向、体积应变及侧向应变与轴向应变比值的负值（计算得到的泊松比）。在峰值应力之后,密实岩石单元服从线性应变软化行为及随后的理想塑性行为,而材料缺陷呈现理想塑性行为。当孔隙压力较高时,应力-侧向应变曲线具有一平台;破坏的前兆更明显;变形后岩样的体积总是大于原始体积;在初始加栽阶段、均匀变形阶段及峰后变形阶段,由于明显的侧向膨胀,计算得到的泊松比远大于0．5。当孔隙压力较低时,在峰值应力之前,变形后岩样的体积小于原始体积,体积扩容出现于峰值应力之后,引起了负的体积应变。利用广义虎克定律,解释了平面应变弹性状态下数值结果的合理性。对岩样进行带状区域扫描后,确认随机缺陷的初始分布与岩样的最终破坏形态紧密相关。 相似文献
5.
固体传压高温高压实验中,岩样环境温度是重要参数。本文用有限元法进行了理论估算,并用热电偶实测了岩样区的温度场分布。给出了场的特征分析。 相似文献
6.
为了计算带任意地形的各向异性介质中二维大地电磁响应,本文在非结构化网格的基础上,采用有限体积法,开发了二维大地电磁各向异性正演模拟的新算法.首先,从Maxwell方程出发,推导二维各向异性介质中大地电磁场的边值问题;然后,采用三角网格自动生成技术对求解区域进行非结构化网格剖分,进而构建节点中心控制体积单元,利用有限体积方法,得到求解边值问题的大型稀疏线性方程组;最后,利用Pardiso精确地计算了大地电磁响应值.三个各向异性模型的计算结果表明,本文开发的有限体积算法,不仅能够高精度求解带任意地形的大地电磁电导率各向异性问题,而且对于同一模型,该方法的计算消耗和精度都与有限单元法相当.因此,有限体积法是处理电磁法各向异性问题的一种有效方法. 相似文献
7.
计算火灾下空心混凝土楼板内温度场的关键是计算空腔内的热量传递。假设火灾过程中空腔内气体压力为常数,空气流动引起的热量损失忽略不计,空腔内的热量传递通过壁面辐射和空气传导实现。基于奥本海姆网络法,本文提出了计算空心混凝土楼板空腔内辐射热流的算法,推导了矩阵形式的空腔内有效辐射方程,考虑到系数矩阵的对称性,采用LDLT分解法求解有效辐射方程。利用空腔边界上的有效辐射热流,计算空腔边界上单元的温度节点荷载。这种方法已经嵌入作者开发的温度场分析软件TFIELD之中,用TFIELD对底面按ISO834曲线进行加热的空心混凝土楼板进行计算,结果表明,本文算法十分有效,对研究空心混凝土楼板的抗火性能具有积极意义。 相似文献
8.
电阻率法有限单元正演模拟中,采用第三类边界条件时为保证精度仍要求较大范围的计算域.无单元法为地球物理领域的新兴正演模拟方法,其计算效率低,但其中采用的移动最小二乘(MLS)形函数相比于有限单元法形函数具有良好的连续性,模拟精度高.本文将MLS形函数应用于电阻率法有限单元2.5维正演的第三类边界条件处理,提出电阻率法有限单元-移动最小二乘(FEM-MLS)耦合正演方法.通过不同正演方法的模型算例模拟结果对比,验证了本文算法的有效性,并讨论了各个参数选择对模拟结果的影响.本文数值模拟结果表明采用第三类边界条件时,在同等计算精度前提下,FEM-MLS耦合法相比于有限单元法可进一步缩小计算域并提高了计算效率,相比于采用较大计算域满足边界条件的有限单元法计算效率提高了约一倍,相比于采用相同小范围计算域的有限单元法平均精度提高了约一倍. 相似文献
9.
结构阻尼与材料阻尼的关系 总被引：5，自引：1，他引：4
针对固体材料在弹性阶段的内阻尼,采用材料粘弹性本构关系,通过有限单元法建立了考虑剪切变形影响的杆单元阻尼矩阵,得到了材料拉压粘滞系数和剪切粘滞系数之间的关系。我们可以通过实验测试材料的损耗因子或粘滞阻尼系数(材料阻尼),像结构质量矩阵、刚度矩阵一样通过比较确切的计算得到弹性阶段结构内阻尼的阻尼矩阵(结构阻尼),而且此阻尼矩阵既基于材料阻尼实验测试,又便于数学处理,且物理意义明确。其次,本文分析了材料阻尼对结构阻尼比影响,得到了材料损耗因子与结构模态阻尼比间的关系,并通过10层钢筋混凝土剪切型框架结构进一步给出了具体的数值结果。 相似文献
10.
为了提高海洋可控源电磁的地形响应模拟精度,探讨了三维交错网格剖分有限体积法求解频率域电磁场解的数值方法.在似稳场近似条件下,推导了有限差分法和有限体积法低频电磁场控制方程的离散表达式.通过对离散表达式的分析比较表明,有限体积法较之有限差分具有更高的计算精度,同时也具有与有限差分相当的计算效率.对含地形界面网格单元电导率采用加权平均处理,通过与二维有限元程序的计算结果对比,验证了该方法可以有效提高有限体积法对地形变化的模拟精度.应用有限体积法计算了一个三维带地形储层模型的可控源电磁响应,分析表明地形变化对电场分量影响明显,磁场分量对地形变化不敏感. 相似文献
11.
电导率是表征岩石电学性质的重要物理参数,在地质资源勘查和测井解释等领域发挥着巨大作用.快速、准确地确定岩石电导率具有重要的理论和实践意义.作为近年来发展的一种岩石物理数值模拟工具,数字岩心技术在定量计算电导率等物性参数方面应用广泛.三维微观结构的准确获取是数字岩心技术计算岩石电导率的关键,但传统获取岩石三维微观结构的方法较为复杂费时.为了方便快速地通过数字岩心技术计算岩石的电导率,本文研究了岩石二维与三维数字岩心的电导率联系.我们基于微米级X射线CT扫描得到的三个砂岩样品的微观结构信息建立了三维数字岩心,并通过有限元法计算的三维数字岩心电导率与实验数据的对比验证数值计算方法的有效性.随后我们数字地扩展了岩石的孔隙,产生了较大孔隙度的三维数字岩心样本,在此基础上,计算了三维数字岩心和相应二维数字岩心的电导率,并通过Archie公式分别拟合了电导率与孔隙度之间的关系,得到了相应的胶结系数.结果表明,三维数字岩心的胶结系数小于二维数字岩心的胶结系数,且二者的比值与岩石实测孔隙度呈线性负相关关系.以该联系为纽带,通过二维图像快速计算得到的电导率与孔隙度关系,确定了三维数字岩心的电导率与孔隙度关系,并进一步通过三维数字岩心的孔隙度计算其电导率.该方法计算得到的人工砂岩样品的电导率与其三维数字岩心电导率相关系数高于96%,验证了基于二维图像的数字岩心电导率计算方法的有效性.本文的研究结果为快速、准确地计算岩石电导率提供了新的思路,在油气勘探开发中有广阔的应用前景. 相似文献
12.
以８００ＭＰａ高温高压三轴室为例，介绍了简体热应力计算方法和结果。分析讨论了热应力对筒体安全性的影响，指出高温差低载荷是最危险的情况，并提出预防与解决的方法。 相似文献
13.
某高温气冷堆核电厂结构地震反应分析 总被引：4，自引：3，他引：1
以某高温气冷堆核电厂结构为原型,利用有限元软件建立三维结构实体模型,开展了模态分析和弹性动力时程分析的数值计算,以探讨某高温气冷堆核电厂的结构特性和抗震性能,并重点分析了在三向地震动作用下高温气冷堆核电厂的加速度、位移反应时程和楼层反应谱.总体上看,高温气冷堆核电厂在两个水平向的刚度比较均匀,楼层反应接近,布局较为合理;在三向地震动作用下,顶层中心点的竖向楼层反应均明显大于两水平方向楼层反应,因此在高温气冷堆核电厂结构设计中,应关注竖向地震动对核电厂地震反应带来的不利影响. 相似文献
14.
EXPERIMENTAL STUDY OF TEMPERATURE EVOLUTION AND IDENTIFICATION OF INSTABILITY POSITION OF PLANAR STRIKE-SLIP FAULT DURING PROCESS OF STICK-SLIP 
下载免费PDF全文
下载免费PDF全文 Stick-slip of fault in laboratory accompanies change of temperature. Temperature change is not only concerned with sliding friction, but also with the stress state of the sample. In this article, we use infra-red thermal imaging system as wide-range observation means to study the temperature variation of different stages during the deformation of sample. The rock sample for the experiment is made of granodiorite from Fangshan County with a size of 300mm×300mm×50mm. It is cut obliquely at an angle of 45°, forming a planar fault. Two-direction servo-control system was used to apply load on the sample. The load in both directions was forced to 5MPa and maintained constant (5MPa) in the X direction, then the load in the Y direction was applied by a displacement rate of 0.5μm/s, 0.1μm/s and 0.05μm/s successively. The left and below lateral of the sample were fixed, and the right and top lateral of the sample were slidable when loaded. The experiment results show not only the temperature change from increase to decrease caused by conversion of stress accumulation to relaxation before and after the peak stress, but also opposite variation of temperature increase on fault and temperature decrease in rock during instability stage. Most important of all, we have found the temperature precursor identifying the position of instability through the temperature variation with time along the fault. It shows that rate of temperature increase of instability position keeps relative high value since the stage of strongly off-linear stage, and accelerates in stage of meta-instability. After separating the effect of friction and stress, we found that temperature increase occurs in the rock near the fault instead of on the fault, which means the mechanism of temperature increase is stress accumulation. Temperature of fault at the instability position does not increase, which means the position is locked. We speculate that the position of locked area on fault with high stress accumulation near the fault may be the future instability position. It is of significance of studying temperature variation during stick-slip to the monitoring of earthquake precursors. Heat caused by friction of earthquake needs long time to transfer to the surface and could not be detected as a precursor. While the stress of surface rock near the fault would change as the stress of interior rock changes, which could cause detectable temperature variations. The research purpose of this article is to find special change positions before instability. As the temperature variations are caused by stress and slip of fault, the results are also meaningful to analysis of stress and displacement data related to earthquake precursors. 相似文献
15.
高围压条件下岩石力学性质的温度效应 总被引：10，自引：0，他引：10
以往的研究表明，深部高温高压条件下岩石力学性质与地面所看到的有很大差别。深部岩石力学与深部的高温高压环境有关。本文介绍几个典型实验实例，叙述了的岩石样品在不同的温度条件下力学性质发生的变化。 相似文献
16.
数字岩石物理利用三维成像技术和数学方法,建立数字岩心,开展多场物理响应模拟,计算岩石的等效弹性参数,为岩石物理学研究开拓了新的领域.本文发展了一种基于高阶有限差分的岩石模量数值计算的新方法,该方法便于理解,容易实现,占用内存相对较少,计算效率高,结果合理,弥补了常规岩石物理实验周期长,成本高,误差大等不足.该方法将三维数字岩心样本嵌入一个具有与岩心样本骨架颗粒相同弹性性质的区域中扩充成一个新模型,测量由数字岩石样本非均匀结构带来的纵波或横波峰值振幅的时间差,利用该时间差(相对参考模型)估算纵横波的等效速度,进而求取等效弹性模量.理论模型和实际岩心的计算结果表明,数值模拟结果与实验结果有较高的吻合度,验证了该方法的合理性. 相似文献
17.
18.
本文从总电位法和异常电位法两个方面,分别就直流电法有限元中所采用的三类边界条件对计算结果的影响进行了分析,为采用合适的边界条件及选取最佳边界距离来提高正演精度提供了基础.结果表明,混合边界条件精度最高,可大大缩小求解区域而不影响计算精度,其次是Dirichlet边界条件精度较高,但测点越靠边界误差会越大,必须取足够大的边界区域,齐次边界条件的误差比较大,但如果采用非二极装置,通过电位差计算得到的视电阻率,由于无穷远边界对电位差的影响基本消除,视电阻率计算误差与混合边界条件下的接近,在反演中,为了节省计算时间,经常使用齐次边界条件进行有限元正演. 相似文献
19.
间断Galerkin有限元法（DG-FEM）作为一种有效的高阶有限元法受到了国内外学者的广泛关注.本文基于任意高阶间断Galerkin有限元法对弹性波方程进行空间离散,并将离散后所得的非齐次线性常微分方程系统齐次化,最后结合针对齐次问题的强稳定性保持龙格库塔（SSP Runge-Kutta）算法,将DG-FEM推广至时间任意高阶精度.另外,借鉴近最佳匹配层（NPML）的思想,基于复频移（CFS）拉伸坐标变换推导了一种新的PML吸收边界条件（简称为CFS-NPML）,该CFS-NPML能够与DG-FEM算法很好地结合,形成有效的起伏地表地震波传播数值模拟技术.数值试验结果表明,DG-FEM具有高阶精度,可以适应任意复杂起伏地表和复杂构造情况下的弹性波传播数值模拟.同时,CFS-NPML对包括面波等震相的人为边界反射都具有良好的吸收效果. 相似文献