Revisiting the Lowrie–Fuller test: alternating field demagnetization characteristics of single-domain through multidomain glass–ceramic magnetite

This paper reports the alternating field demagnetization characteristics of glass–ceramic magnetite assemblages carrying weak-field thermoremanent magnetization (TRM), weak-field anhysteretic remanent magnetization (ARM), and saturation remanence (J_rs). Average grain sizes vary from less than 0.1 μm to approximately 100 μm, and hysteresis parameters indicate that these assemblages encompass single-domain (SD) through truly multidomain (MD) behavior. In all assemblages, weak-field TRM and weak-field ARM are more stable to alternating field demagnetization than is (J_rs). This response is especially remarkable in the 100 μm assemblage, which otherwise displays truly MD behavior. Although the SD samples pass the Lowrie–Fuller test for SD behavior, calculations presented here show that populations of noninteracting, uniaxial SD grains should behave in just the opposite sense to that reported originally by Lowrie and Fuller. This discrepancy could indicate that SD, glass–ceramic magnetite populations are more affected by magnetic interactions than would be expected for magnetite crystals that nucleated individually from a silicate matrix. This interpretation is supported by the SD assemblages failing the ‘Cisowski' test: that is, the curves for acquisition and AF demagnetization of (J_rs) intersect well below the 50% mark. However, a second and intriguing explanation of the SD-like results obtained from all samples is that alternating field demagnetization characteristics reflect a strong dependence of local energy minimum domain state, and its associated stability, on the state of magnetization.     

改进的地壳速度结构遗传反演方法

遗传算法的收敛速度与模型参数初始搜索范围的划分精度有关.本文提出了在低精度的基础上开始迭代,在选代过程中逐步缩小搜索范围的改进措施.这种改进使得遗传算法的收敛速度和反演精度同时得到提高,并且又不增加模型空间的大小.一维地壳速度结构的反演计算结果表明方法是有效的.最后给出了实际应用中应采用的反演策略.     

三维密度界面的正反演研究和应用

重力位场的界面反演是位场处理解释中的重要问题.本文将基于快速傅里叶变换的频率域界面反演方法Parker-Oldenburg公式推广到物性可随深度变化的三维情况,得出了密度可以横向、纵向任意变化的重力界面正反演公式.该方法在计算时可以合理地选取地面下某一深度作为基准面以减小界面起伏,使迭代易于收敛.理论模型试验表明该方法反演精度高,收敛速度快,在密度界面反演中具有广泛的实用价值.最后利用该方法反演华北地区莫霍面的深度,反演结果得到了地震测深数据的验证.     

拟合加密控制点反应谱的地震动时程合成

在传统地震动反应谱频域拟合方法的基础上,提出了拟合加密控制点目标反应谱的地震动合成方法,并在此基础上利用实际算例与传统方法进行了对比.本文方法在每次迭代调整过程中均考虑了对地震动加速度强度包络函数的模拟,并且对调整后的地震动进行了基线校正,从而确保最终合成的地震动时程既能以一定精度拟合加密目标反应谱,又能满足目标包络函数的要求,同时保证积分所得位移曲线不产生基线漂移.数值算例分析结果表明,本文方法与传统非加密控制点目标谱拟合方法相比,在地震动时程曲线合成、反应谱拟合精度以及收敛速度等方面均具有一定的优势.      

重磁位场转换计算中迭代法的综合分析与研究

处理转换计算在重磁资料解释中发挥着重要的作用,但一些计算如向下延拓、化极等有时是很不稳定的,在频率域中则表现为其转换因子具有明显的放大作用,所以其FFT理论计算结果是不稳定的.因此,很多研究工作都是围绕增加计算的稳定性、提高计算效果进行的,其中迭代法是近来在研究中受到普遍重视的方法技术,并取得了较好的成果.但也存在对迭代法研究还不够深入,对其存在的缺点认识不够充分、客观等问题,例如,迭代法进行延拓及化极等计算时,对一些具体应用虽能在一定程度上获得较好的计算结果,但却存在计算结果并不会随着迭代次数的增加而得到持续改善的问题,对于原本不稳定的计算,迭代法在迭代次数比较大时,所得的结果依然是不稳定的.为此,本文在对迭代法进行分析研究的基础上,进一步推导了迭代法的通式,并分析了对迭代法收敛性影响的各种因素.分析结果表明:迭代法收敛到FFT理论直接计算结果的决定因素是计算过程中如何选择原始数据到目标数据的映射函数;在选择了合适的映射函数的情况下,迭代次数不仅仅是决定计算成本,而是决定结果好坏的关键因素;增加迭代次数虽然能够使计算收敛到FFT直接计算理论结果,但如果该理论结果本身就是不稳定的,则迭代法计算如果收敛,也是收敛到一个不稳定的结果.所以针对位场处理转换中一些不稳定计算采用迭代法,并没有从根本上解决计算的不稳定性问题.     

积分迭代法的正则性分析及其最优步长的选择

位场积分迭代法是一种解决位场大数据量、大深度向下延拓的有效方法.本文基于Kirsch正则化子理论,推导了积分迭代法对应的正则化滤子函数,并证明积分迭代法为一种求解位场向下延拓不适定反问题的正则化方法.针对积分迭代法迭代步长固定、迭代次数较多影响收敛速度的问题,提出该迭代法最优迭代步长的选择原理.理论模型和实测数据对比分...     

结合欧拉矢量的反演算法构建青藏高原东北缘地壳运动速度场模型

探讨地壳运动速度场模型的构建方法,提出结合欧拉矢量的维多样性动态权重粒子群算法构建地壳运动速度场模型。通过模拟算例验证该算法的稳定性和有效性,建立的速度场模型与线性权重粒子群算法和非线性权重粒子群算法的计算结果相比具有较高的精度,且收敛速度较快。利用青藏高原东北缘1999—2013年中国地壳运动观测网络观测到的GPS水平速率结果,在块体划分和模型辨识的基础上,建立青藏高原东北缘地壳运动速度场模型,并将其与最小二乘配置法的计算结果进行比较,结果表明改进的粒子群算法建立的地壳运动速度场模型具有较高的精度。     

地层衰减对地震波速度逆散射反演的影响研究

在利用地震波数据进行地球物理反演时,地层对地震波的吸收衰减效应会对地层物性参数的准确反演产生较大的影响,因此利用黏弹性声波方程进行反演更符合实际情形.本文在考虑地层衰减效应进行频率空间域正演模拟的基础上,提出基于黏弹性声波方程的频率域逆散射反演算法并对地震波传播速度进行反演重建,在反演过程中分别用地震波传播复速度和实速度来表征是否考虑地层吸收衰减效应.基于反演参数总变差的正则化处理使反演更加稳定,在反演中将低频反演速度模型作为高频反演的背景模型进行逐频反演,由于单频反演过程中背景模型保持不变,故该方法不需要在每次迭代中重新构造正演算子,具有较高的反演效率;此外本文在反演过程中采用了基于MPI的并行计算策略,进一步提高了反演计算的效率.在二维算例中分别对是否考虑地层吸收衰减效应进行了地震波速度反演,反演结果表明考虑衰减效应可以得到与真实模型更加接近的速度分布结果,相反则无法得到正确的地震波速度重建结果.本文算法对复杂地质模型中浅层可以反演得到分辨率较高的速度模型,为其他地震数据处理提供比较准确的速度信息,在地层深部由于地震波能量衰减导致反演分辨率不太理想.     

Equivalent mechanical models of sloshing fluid in arbitrary‐section aqueducts

This paper reports on a semi‐analytical/numerical method to model sloshing water in an arbitrarily shaped aqueduct. The water motion is assumed to be inviscid, compressible, and linear (small displacement). The transverse sloshing fluid in an aqueduct is equivalently simplified as a fixed rigid mass M₀ and a mass–spring system (M₁, K₁). According to a rule that the actual fluid (computed with finite element model) and its equivalent mechanical model have the same first sloshing frequency and acting effects on the aqueduct, the analytical solutions of the fixed (impulsive) mass M₀, sloshing (convective) massM₁, spring stiffness K₁, and their locations in the aqueduct body are acquired by the least squares (curve fitting) algorithm. Applying this equivalent principle, the equivalent mechanical models are respectively obtained for the sloshing water in rectangular, semicircular, U‐shaped, and trapezoid aqueducts. The equivalent principle and fluid models are validated through comparison investigations involving rectangular and U‐shaped aqueducts. The dynamic properties and seismic responses of the original and equivalent systems are simulated, compared, and discussed for a U‐shaped aqueduct bridge. The main purpose of this paper is to provide a simplified model of sloshing fluid for the seismic/wind‐resistant computation of the support structures of the aqueduct bridge. Copyright © 2011 John Wiley & Sons, Ltd.     

Decomposition methods for water resources optimization models with fixed costs

This paper illustrates the application of two decomposition algorithms, generalized Benders decomposition (GBD) and outer approximation (OA), to water resources problems involving cost functions with both discrete and nonlinear terms. Each algorithm involves the solution of an alternating finite sequence of nonlinear programming subproblems and relaxed versions of a mixed-integer linear programming master problem. Three example models, involving capacity expansion of a conjunctively managed surface and groundwater system, are formulated and solved to demonstrate the performance of the algorithms. The results show that OA obtains solutions in far fewer iterations than GBD, but OA requires more computational resources per iteration. As a result, depending on the mixed-integer programming and nonlinear programming solvers available, GBD may be better suited for solving larger planning problems.     

Lp模反褶积及其在地震勘探中的应用

在文献[1,2]的基础上,本文提出一种矩阵逆迭代和增强强反射的自适应算法,方法简明,使用方便。理论模拟和实际资料处理结果说明该法是很有效的。     

反馈迭代法压制表面多次波效果分析

地震数据处理通常假定反射数据仅由一次波组成,地震勘探资料中的表面多次波通常被视为相干噪音予以压制.基于波动方程预测的反馈迭代法为数据驱动的方法,可有效压制复杂地下介质的表面多次波,通常可由级数展开法和迭代法两种方法实现.文中分别从理论阐述以及应用效果分析方面对基于波动方程预测的反馈迭代法的两种实现途径进行研究,并采用GPU/CPU协同并行加速计算预测表面多次波,改善以往CPU计算效率低的状况.文中对含表面多次波的炮数据和复杂的SMAART模型进行了计算,并将利用级数展开法与迭代法压制表面多次波的效果进行对比分析,结果表明,迭代形式反馈迭代法压制表面多次波的方法效果更佳.由于实现迭代法和级数展开法需要完成的自适应匹配相减的次数不同,迭代法的计算成本略高于级数展开法.     

利用微分进化法确定海洋磁场向下延拓中的最优参数

位场的解析延拓是实现不同高度海洋地磁场相互转换的主要途径,是构建海洋三维磁空间背景场模型的关键技术.针对位场向下延拓迭代法中最优正则化参数及最佳迭代次数难以确定问题,尝试引入微分进化法,以正则化参数及迭代次数为种群变量,以延拓结果的熵值为目标函数,以目标函数最小化为搜索准则,实现两种参数的并行全局寻优.采用实测数据对微分进化法在几种常用的迭代法中最优正则化参数及最佳迭代次数的确定进行了分析,与传统L-曲线准则确定的最优正则化参数及多次试验确定的最佳迭代次数进行对比,结果表明：微分进化法确定的最优参数能使三种迭代法取得最佳迭代效果,延拓结果与真实地磁场最为接近,并且该法计算稳定、自适应强,建议在海洋磁场数据向下延拓中应用.     

弹性波场数值模拟的隐式差分多重网格算法

为了得到稳定的弹 性波数值模拟,而不得不选择隐式差分方程;为了提高解的精度,又不得不增加节点数目, 但同时也降低了隐式迭代求解的收敛速度. 为此,本文使用隐式差分的多重网格算法进行弹 性波数值模拟,多重网格算法通过粗网格收敛较快的迭代过程求出近似解,以近似解为初值 使用细网格进行精确的迭代求解,从而加速了隐式迭代求解的过程,能够以较高计算速度、 精度、稳定性完成弹性波传播过程的数值模拟.     

Development of the universal and simplified soil model coupling heat and water transport

It is very important to develop a universal soil model with higher simplicity and more accuracy, which can be widely applied to very general cases such as wet or dry soil, frozen or unfrozen soil and homogeneous or heterogeneous soil. Firstly in this study, based on analysis of both magnitude order and the numerical simulation results, the universal and simplified soil model (USSM) coupling heat and mass transport processes is developed. Secondly, in order to avoid the greater uncertainty caused by the phase change term in numerical iteration process for the model solution obtaining, new version of the universal simplified soil model (NUSSM) is further derived through variables transformation, and accordingly a more efficient numerical scheme for the new version is designed well. The simulation results from the NUSSM agree with the results from more complicated and accurate soil model very well, also reasonably reproduce the observed data under widely real conditions. The new version model, because of its simplicity, will match for the development of land surface model.     

电离层层析重构的一种新算法

自适应联合迭代重构算法是电离层层析成像中一种收敛速度较快的反演算法.然而,在电离层电子密度重构过程中,该方法对迭代初值的精度要求较高.针对上述问题,本文提出了一种约束自适应联合迭代重构算法,该方法通过附加合理的平滑约束,减弱了没有观测信息的格网对迭代初值精度的依赖性,有效地提高了反演结果精度.数值模拟实验和实测数据的反演结果证实了该算法的可行性和在重构精度上的优越性.     

钢筋混凝土结构抗震能力评估的简化能力谱方法

本文通过对现有能力谱法的研究,在吸收前人研究成果的基础上,探讨了一种简化的能力谱方法。该方法不是根据需求谱与能力谱有无交点(性能点)为评估依据,而是以能力谱为根据求出不同延性状态下既有建筑物相应的抗震能力,并与需求谱相应的谱加速度比较,判断结构的抗震能力是否满足要求。该方法不需要复杂的迭代计算来求结构的性能点,计算过程简单;同时,可以考虑既有建筑物的老化及损伤。最后,本文应用自行编制的计算程序,通过一个工程实例说明了该方法的应用及其特点。     

双侧向电阻率测井反演算法研究

双侧向测井仪因其激励源的复杂性给其测量数据的反演带来了极大的困难．本文利用变形玻恩（Born）迭代法对轴对称二维非均匀介质中双侧向测井仪的位场数据进行了严格反演．在每一次迭代过程中采用了高效数值模式匹配的正演计算方法,它可半解析地求出反演中所需计算的格林（Green）函数及其偏导数,同时利用所推导出的位场非线性积分方程中积分运算的半解析形式,大大提高了计算效率和反演质量．通过反演实例表明,反演结果与电导率原状剖面吻合较好,并且受噪声影响较小．     

改进的高精度匹配追踪方法研究及应用

本文在常规算法的基础上,提出了一种改进的匹配追踪方法.该方法通过引入新的完备库构建策略,并以雷克子波作为原子,利用原始地震数据最大相关性估计完备库的原子的位置和能量,使得分解精度进一步提高,极大地提高了算法的适应性;为避免原子间隔过小问题,引进了最小原子间距,使得分解效率和分解质量进一步提高.模型试算和实际资料应用表明,本文算法不仅提高了信号稀疏表示的质量,加快了收敛速度,而且算法的适应性和分解精度也得到提高;该方法能够较好地挖掘地震有效信息,提高地震解释精度.     

Approximate traveltime inversion in downhole microseismic monitoring

In downhole microseismic monitoring, accurate event location relies on the accuracy of the velocity model. The model can be estimated along with event locations. Anisotropic models are important to get accurate event locations. Taking anisotropy into account makes it possible to use additional data – two S-wave arrivals generated due to shear-wave splitting. However, anisotropic ray tracing requires iterative procedures for computing group velocities, which may become unstable around caustics. As a result, anisotropic kinematic inversion may become time consuming. In this paper, we explore the idea of using simplified ray tracing to locate events and estimate medium parameters. In the simplified ray-tracing algorithm, the group velocity is assumed to be equal to phase velocity in both magnitude and direction. This assumption makes the ray-tracing algorithm five times faster compared to ray tracing based on exact equations. We present a set of tests showing that given perforation-shot data, one can use inversion based on simplified ray-tracing even for moderate-to-strong anisotropic models. When there are no perforation shots, event-location errors may become too large for moderately anisotropic media.