岩心核磁共振实验对低孔低渗碳酸盐岩储层的适应性研究

低孔低渗碳酸盐岩储层矿物成分复杂、岩石骨架参数难以确定、储集空间类型多样、孔隙结构及孔渗关系复杂、常规测井曲线响应特征不明显，使其测井评价极其困难．本文利用核磁共振测井定量评价低孔低渗碳酸盐岩储层岩心的孔隙结构、计算储层参数．利用T₂谱分布曲线分析孔隙结构、计算T₂截止值；在此基础上计算岩心总孔隙度、有效孔隙度、束缚水孔隙度、渗透率等储层参数，并与常规岩心实验结果进行对比分析；最后，总结出核磁共振测井在低孔低渗储层中的应用优势与局限，为核磁共振测井评价模型的建立提供基础数据．     

小白鼠肌肉组织的NMR质子自旋交换分析

本文在Zimmerman-Brittin两相质子交换核磁共振弛豫模型基础上,分析了NMR弛豫实验中检测信号与各相表现和本征弛豫多数的关系,编写了自动化处理实验数据的计算机程序,这一技术可用于复相系统中不同成分的NMR表现和本征弛豫特性研究中,本文中的样品是选用健康新鲜的小白鼠肌肉,没加任何处理,用h-h,s-h,s-s脉冲序列,反转恢复法(π-τ-π/2)在强场下(0.92T)做T₁、T₂测定实验,分析结果表明本征弛豫参数T₁=1050ms,T₂=4500μs的成分是由肌肉中的"自由水"引起的,其质子相对含量为69%;本征弛豫参数T₁=530ms,T₂=26μs的成分是由肌肉中的"束附水"引起的,其质子相对含量为9%,本征弛豫多数T₁=530ms,T₂=1250μs的成分是由肌肉中的各种大分子和有机物引起的,其质子相对含量为9%,本征弛豫参数T₁=470ms,T₂=1250μs的成分由样品中的脂肪引起的,其质子相对含量为13%,在肌肉组织中的质子与水中质子之间有强烈的交换作用,其交换率k=1000s^-1.在脂肪中的质子与其它成分之间没有交换作用。     

基于低场核磁共振的抚顺油页岩孔隙连通性演化研究

油页岩原位注热开采过程中，储层内部孔隙结构的连通性直接影响载热介质的流动行为和传热效率，同时对油气产物的扩散和流动行为起控制作用.本文利用低场核磁共振（LF NMR）技术，考察了不同热解终温（23~650℃）处理时，饱和水及束缚水状态下抚顺油页岩的T₂谱，分析了可动流体T₂截止值、束缚流体孔隙度、饱和流体孔隙度、渗透率等NMR孔隙参数，定量研究了随热解终温升高，抚顺油页岩孔隙结构的连通性演化规律.研究结果表明热解终温对抚顺油页岩孔隙连通性及渗透率的变化起控制作用，且可动流体孔隙度对总孔隙度的增加起主要促进作用，这说明热解终温升高加大了渗透率及油气产物的输运能力.本文为深入认识油页岩原位热解过程中孔隙结构的演化提供了依据.     

基于改进非线性拟合的核磁共振T2谱多指数反演

核磁共振T₂谱多指数反演算法是开展复杂体系样品核磁共振(NMR)弛豫研究最重要的数学工具. 常用的T₂谱多指数反演算法一般都是事先给出弛豫时间T₂分布的布点, 然后转化为线性拟合问题进行求解. 在求解的T₂谱较为分散的时候, 反演得到的T₂谱精确度不高, 分辨率较低. 非线性拟合是解决这个问题的有效办法. 本文针对分散T₂谱反演利用非线性拟合时遇到的初值依赖及运算复杂问题, 利用线性回归最小二乘方法, 改进了其中的带非负约束非线性优化模型, 将搜索的反演参数从T₂, f 减少为T₂, 加快了收敛速度, 减少了对初值的依赖, 提高了反演精度, 使算法更加稳健. 通过用改进的Levenberg-Marquardt算法和差分进化算法进行计算机模拟反演及实验数据反演, 验证了改进方法在核磁共振T₂ 谱反演中的有效性.     

玉米胚芽胚乳中油水动态行为的NMR研究

应用核磁共振时域二维弛豫相关谱分析的方法,对玉米胚芽的核磁共振弛豫特性进行了研究.对以下三种样品进行了测定:(1)胚芽不含油(用化学方法蒸馏出)但是含有不同百分含量的水;(2)胚芽不含水(真空蒸发)但是含有不同百分含量的油;(3)含有不同百分含量的水的胚乳水系统.研究表明随着水含量的增加水的纵向弛豫时间T₁减小,表明水与胚芽基质具有强烈的相互作用,而油的纵向弛豫时间T₁随着油含量的增加变化很小,表明油与胚芽基质之间没有大的相互作用,油几乎独立地处在胚芽中;对低水分(少于5%)玉米胚芽水系统,结果表明水是束缚于胚芽基质上,而油则是以油胞的形式存在于胚芽中.在玉米胚芽水系统中有两个T₁弛豫成分,长T₁弛豫成分从基质到水的交换率为57s^-1;从水到基质的交换率为113s^-1.短T₁成分从基质到水的交换率为245s^-1;从水到基质的交换率为821s^-1.在低水合状态下,胚乳水系统仅有一个T₁弛豫成分.胚乳到水的交换率高于37s^-1.     

表面活性剂CTAB水溶液中的T2弛豫分散研究

用CPMG脉冲序列测定了表面活性剂十六烷基三甲基溴化铵（CTAB）分子中的氮甲基(N-CH₃)质子的横向弛豫时间（T₂表观）,并发现测得的T₂表观\}与序列中的重聚脉冲间隔时间的一半τ _cp有关,说明存在横向弛豫分散现象. 当在τ_cp≤1 ms时,T₂表观与τ²_cp}呈线性关系;而当τ_cp≥4.6 ms时,T₂表观变得与τ_cp无关. 利用Luz-Meiboom两体化学交换模型计算了不同浓度的CTAB溶液中的N-CH₃质子的本征横向弛豫时间（T₂本征）和化学交换速率k_ex,发现k_ex与T₂本征和自扩散系数D一样,在临界胶束浓度（CMC）附近发生突变. 这个突变反映了CTAB分子在从单体到胶束的转变过程中其动力学特性发生了改变.      

油井样品NMR T2谱的影响因素及T2截止值的确定方法

T₂ 谱是核磁共振（NMR）测、录井技术应用与解释、评价的基础. 岩样T₂ 谱受仪器测量参数、样品性质（岩性、颗粒大小、样品粒度、样品干湿状态、孔隙流体含量及性质、磁化率、润湿性）及地层水矿化度等因素的影响. T₂ 截止值是T₂ 谱中最重要的参数之一,选取的科学性与准确性直接影响到核磁共振测量结果. 通过文献查询,对T2 谱的影响因素及T₂ 截止值的确定方法进行了分析.     

基于ADI-FDTD模拟的岩石核磁共振横向弛豫特征分析

将复杂的骨架-孔隙系统抽象成等效双孔介质,根据Bloch方程构建数学模型,用交替隐式时域有限差分(ADI-FDTD)和联合反演迭代法(SIRT)进行横向宏观磁化矢量的数值模拟与核磁共振T₂谱的反演,定量研究扩散系数、弛豫速率、孔隙组分比和孔隙宽度对核磁响应的影响.结果表明:横向宏观磁化矢量衰减速率与扩散系数和微孔隙分量成正比,与孔隙宽度成反比,与表面弛豫速率基本无关.当扩散系数较大、孔隙宽度较小时,核磁共振T₂谱难以直观反映孔隙组分及孔隙结构.应用核磁共振评价孔隙结构时需特别注意扩散系数和孔隙尺寸的影响.     

大分子中磁偶极弛豫的分离

提出采用仅与偶极弛豫有关的¹³C核同¹H核之间的核交叉弛豫速率来考察分子内部运动状况。给出了实际用于考察分子运动状况的关系式,从理论上阐明了本文方法比传统上直接采用常规测定的¹³C核纵向弛豫时间T₁、演向弛豫时间T₂及NOE增强因子η来研究分子运动状况的方法要合理。建立了既适用于小分子又适用于大分子自旋系统¹³C核总自旋弛豫中偶极弛豫分离的公式.     

嵌段共聚物溶致液晶相中水的2H-NMR动力学分析

采用D₂O 的²H-NMR线型和弛豫分析了PEO-PPO-PEO/D₂O/对二甲苯体系的层状和六角液晶相的动力学行为. 通过实验测得了两个不同体系的自旋 晶格弛豫时间T₁、自旋-自旋弛豫时间T₂和²H-NMR 谱. ²H-NMR 谱均为具有四极劈裂的粉末谱线型，且在谱图的中心，β_LD=54.7°时存在一个倒峰. 倒峰的出现直接表明引起体系中弛豫的主要动力学过程处于极窄化区域. 采用NMR弛豫模型，通过调节动力学参数，使理论模拟的²H-NMR谱、弛豫时间、倒峰的大小与实验的对应量相吻合，求得了体系的动力学参数.     

近红外光谱技术结合RCA和SPA方法检测土壤总氮研究

基于近红外光谱技术结合连续投影算法和回归系数分析对检测土壤总氮含量进行研究。采集农田土壤样本近红外光谱数据,土壤样本数量共394个。由于原始光谱数据量大,在500～2 500 nm光谱波长范围基础上,为简化模型,在原始光谱基础上采用连续投影算法和回归系数分析提取特征变量,以两种变量选择方法提取的特征变量作为输入,分别采用偏最小二乘回归(PLS)、 多元线性回归(MLR)和最小二乘支持向量机(LS-SVM)建模方法建立总氮预测模型,共建立了9个预测模型,最优预测集的决定系数为0.81,剩余预测偏差RPD为2.26。研究表明,基于连续投影算法和回归系数分析选择的特征波长可以应用于近红外光谱检测土壤总氮含量,同时可以大大简化模型,适合开发便携式土壤养分检测仪。     

Estimation of Permeability by Integrating Nuclear Magnetic Resonance (NMR) Logs with Mercury Injection Capillary Pressure (MICP) Data in Tight Gas Sands

It has been a great challenge to determine permeability in tight gas sands due to the generally poor correlation between porosity and permeability. The Schlumberger Doll Research (SDR) and Timur–Coates permeability models, which have been derived for use with nuclear magnetic resonance (NMR) data, also lose their roles. In this study, based on the analysis of the mercury injection experiment data for 20 core plugs, which were drilled from tight gas sands in the Xujiahe Formation of central Sichuan basin, Southwest China, two empirical correlations between the pore structure index ( $ \sqrt {{K \mathord{\left/ {\vphantom {K \varphi }} \right. \kern-0em} \varphi }} $ , defined by the square root of the ratio of rock permeability and porosity) and the R ₃₅ (the pore throat radius corresponding to 35.0 % of mercury injection saturation), the pore structure index and the Swanson parameter have been developed. To consecutively estimate permeability in field applications, based on the study of experimental NMR measurements for 36 core samples, two effective statistical models, which can be used to derive the Swanson parameter and R ₃₅ from the NMR T ₂ logarithmic mean value, have been established. These procedures carried out on the experimental data set can be extended to reservoir conditions to estimate consecutive formation permeability along the intervals with which NMR logs were acquired. The processing results of several field examples using the proposed technique show that the classification scale models are effective only in tight gas reservoirs, whereas the SDR and Timur–Coates models are inapplicable. The R ₃₅-based model is of significance in thin sands with high porosity and high permeability, but the predicted permeability curves in tight gas sands are slightly lower. In tight gas and thin sands, the Swanson parameter model is all credible.     

The influence of track modelling options on the simulation of rail vehicle dynamics

This paper investigates the effect of different models for track flexibility on the simulation of railway vehicle running dynamics on tangent and curved track. To this end, a multi-body model of the rail vehicle is defined including track flexibility effects on three levels of detail: a perfectly rigid pair of rails, a sectional track model and a three-dimensional finite element track model. The influence of the track model on the calculation of the nonlinear critical speed is pointed out and it is shown that neglecting the effect of track flexibility results in an overestimation of the critical speed by more than 10%. Vehicle response to stochastic excitation from track irregularity is also investigated, analysing the effect of track flexibility models on the vertical and lateral wheel–rail contact forces. Finally, the effect of the track model on the calculation of dynamic forces produced by wheel out-of-roundness is analysed, showing that peak dynamic loads are very sensitive to the track model used in the simulation.     

Streamer radius model and its assessment using two-dimensionalmodels

The variable radius method is proposed to approximate the radius of the ionization channel in the one and a half-dimensional (1.5-D) fluid models for studying streamer development, in which the unreasonable constant radius in the traditional 1.5-D fluid models is corrected. The streamer development and propagation between the 1.5-D fluid model with the variable radius method and the two-dimensional (2-D) fluid model using the same initial and external conditions are compared. The radius in each stage of streamer development from the 1.5-D fluid model with the variable radius method shows agreement to a certain degree with that of the 2-D fluid model. The purpose of this paper is not to negate the role of the 2-D fluid models, but to explore the potential of the 1.5-D fluid models and make them more useful and accurate as well as to understand the evolution of streamer radius. The streamer development from the 1.5-D fluid model with the variable radius method not only maintains simplicity of the 1.5-D fluid models, but also presents agreement with the 2-D fluid models for streamers     

Topographic Distribution of Photon Measurement Density Functions on the Brain Surface by Hybrid Radiosity-Diffusion Method

An accurate modelling of light propagation in the head is required to develop an algorithm to reconstruct the image of brain activity. Most previous studies have calculated the light propagation in two dimensional models because of their advantage in computation time and memory requirement over three dimensional models. However, in topographic imaging, the sensitivity distribution in the cross sections parallel to the brain surface which cannot be obtained from a two dimensional model is most important to reconstruct the image. In this study, the light propagation in three dimensional adult head models is calculated by finite element method and hybrid radiosity-diffusion method. The light propagation in the adult head is strongly affected by the non-scattering cerebrospinal fluid (CSF) layer surrounding the brain. The sensitive area is shifted toward the deeper region, and is spread around the CSF layer. The intensely sensitive region on the brain surface is broadly distributed between the source and detector. However, the sensitive region does not penetrate into the deeper part of the brain.     

基于CFD的射流抛光喷射距离的分析和优化

分析了喷射距离对射流抛光效果的影响,基于计算流体动力学进行了喷射距离的分析和优化.通过构建射流抛光不同喷射距离的物理模型,采用能更好地处理流线弯曲程度较大的流动的RNG k-ε紊流模型应用于射流抛光的数学建模,使用SIMPLEC算法对射流模型进行数值计算,得到了不同模型的射流抛光冲击射流流场及工件壁面上的冲击压力、紊动强度、壁流速度分布.根据射流抛光对冲击射流特性的要求,比较和分析不同喷距模型的数值仿真结果,结果显示,射流抛光最优化喷距范围为喷嘴口径的10倍至12倍之间.     

多孔介质中卡森流体的分形分析

研究了非牛顿流体中的卡森流体在多孔介质中的流动特性.基于服从分形分布的弯曲毛细管束模型,运用分形几何理论推导出了该流体在多孔介质中流动的流量、流速、启动压力梯度和有效渗透率的分形解析解.模型中的每一个参数都有明确的物理意义,它将卡森流体在多孔介质中的流动特性与多孔介质的微结构参数有机联系起来.文中给出了卡森流体的流速、启动压力梯度和有效渗透率随着各影响因素的变化趋势,并进行了讨论.所得分形模型可以更深刻地理解卡森流体在多孔介质中流动的内在物理机理. 关键词： 多孔介质 卡森流体 分形     

Estimation of Saturation Exponent from Nuclear Magnetic Resonance (NMR) Logs in Low Permeability Reservoirs

The resistivity experimental measurements of 36 core samples, which were drilled from low permeability reservoirs of southwest China, illustrate that the saturation exponents are not agminate, but vary from 1.627 to 3.48; this leads to a challenge for water saturation estimation in low permeability formations. Based on the analysis of resistivity experiments, laboratory nuclear magnetic resonance (NMR) measurements for all 36 core samples, and mercury injection measurements for 20 of them, it was observed that the saturation exponent is proportional to the proportion of small pore components and inversely proportional to the logarithmic mean of NMR T ₂ spectrum (T _2lm). For rocks with high proportion of small pore components and low T _2lm, there will be high saturation exponents, and vice versa. The proportion of small pore components is characterized by three different kinds of irreducible water saturations, which are estimated by defining 30, 40 and 50 ms as T ₂ cutoffs separately. By integrating these three different kinds of irreducible water saturations and using T _2lm, a technique of calculating the saturation exponent from NMR logs is proposed and the corresponding model is established. The credibility of this technique is confirmed by comparing the predicted saturation exponents with the results from the core analysis. For more than 85 % of core samples, the absolute errors between the predicted saturation exponents from NMR logs and the experimental results are lower than 0.25. Once this technique is extended to field application, the accuracy of water saturation estimation in low permeability reservoirs will be improved significantly.     

Perturbation theory for equilibrium properties of molecular fluids

Perturbation theory for the angular pair correlation function g(r₁₂ ω ₁ ω ₂), using a fluid with isotropic intermolecular forces as the reference system, is applied to the calculation of a variety of macroscopic properties. Comparisons with experiment are made for methane, oxygen and nitrogen (and carbon monoxide for infra-red and Raman band moments) in the dense fluid and liquid states. Theoretical expressions are given and calculations made for thermodynamic properties (isothermal compressibility, pressure, configurational energy, entropy and specific heat) both along and away from the vapour-liquid co-existence curve, for infra-red and Raman band moments, and for neutron scattering cross sections. Excellent agreement with experiment is obtained for all properties, except for the infra-red and Raman band moments; this latter comparison is inconclusive because of large experimental uncertainties. The anisotropic intermolecular forces are found to have very little effect on the liquid isothermal compressibility, in agreement with the first-order theory. Molecular anisotropy has a relatively small effect on the configurational energy and on the Helmholtz free energy, but the effect is large for pressure and specific heat. The pressure is more sensitive to short-range anisotropic forces than the other properties, whereas the specific heat is particularly sensitive to the long-range anisotropic forces. Mean squared torques (derived from infra-red and Raman band moments) are very sensitive to the strengths of the anisotropic forces, and are more sensitive to higher terms in the multipole series than are the other properties. The structure factors for oxygen and nitrogen are found to be little affected by the anisotropic forces.     

苹果树叶片叶绿素含量高光谱估测模型研究

叶片叶绿素含量是评估果树长势和产量的重要参数,实现快速、无损、精确的叶绿素含量估测具有重要意义。本研究以山东农业大学苹果园为试验区,采用高光谱分析技术探索苹果树叶片叶绿素含量的估测方法。通过分析叶片高光谱曲线特征,对原始光谱分别进行一阶微分、红边位置以及叶面叶绿素指数(LCI)变换,分别将其与叶绿素含量进行相关分析及回归分析,建立叶绿素含量估测模型并进行检验,从中筛选出精度最高的模型。结果显示,以LCI为变量的估测模型以及以一阶微分521和523nm组合为变量的估测模型拟合精度最高,决定系数R2分别为0.845和0.839,均方根误差RMSE分别为2.961和2.719,相对误差RE%分别为4.71%和4.70%。因此LCI及一阶微分是估测苹果树叶片叶绿素含量的重要指标。该模型对指导苹果树栽培生产具有积极意义。