首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到18条相似文献,搜索用时 187 毫秒
1.
油井多相流电磁成像测量敏感场仿真   总被引:1,自引:1,他引:1       下载免费PDF全文
为了考察油井多相流电磁成像测量的敏感场分布,本文根据阵列电极结构和测量电磁场特性,应用有限元求解流体截面的电磁场问题,对于不同介质分布模型,计算电势分布,进而模拟测量敏感场.仿真结果表明,测量敏感场呈马鞍面状分布,在靠近发射电极和测量电极的区域敏感性较强,在弧形区域的中间敏感性较小.  相似文献   

2.
油井流动成像束状探测电磁场   总被引:2,自引:0,他引:2       下载免费PDF全文
为了克服油井流动成像电磁探测的"软场"效应,提高成像质量,本文从测量物理模型出发,采用在轴向上屏蔽和径向上聚焦的方法构建束状探测场,通过有限元仿真计算电势分布以及测量敏感场,发现轴向上屏蔽后的电势分布均匀平坦,从而探测电磁场可以简化为二维场;径向上聚焦后敏感区域主要集中在发射电极和接收电极之间的束状区域内,敏感场呈马鞍状,其中间敏感性明显加强.模拟流动实验测量结果表明,束状探测电磁场的测量信号较强,成像质量明显提高.  相似文献   

3.
为了探索电磁流动成像测井仪器激发的相位场信号的传播特性,本文根据实际电磁流动成像传感器的阵列电极结构,应用有限元方法进行正演仿真模拟,得到在不同持水率的层流环境下各电极所测量到的电位相位值,以及6种典型的相位敏感场分布.正演仿真结果表明,相位测量值随着持水率的增加而减小,相位敏感场只对靠近接收电极区域的敏感性较强.不同流型下的相位特征区别明显,表明了利用相位场数据进行电磁流动成像是可行的.  相似文献   

4.
流动成像测井研究进展   总被引:2,自引:3,他引:2  
流动成像测井属于对非均匀介质动态的非线性测量,通过对油气井内多相流体信息的投影扫描和反演处理,实时成像显示流体分布及流动状况,国外研制的电导法,电容法仪器采用电导探针或电容元件构成阵列测量探头,分别利用油气与水的导电特性和介电特性差异辨识井内流体,对流动截面的测量局限于个别点上,已在水平井中见到应用效果,但对物场信息投影测量的数据量和分辨率未能满足成像要求。我国提出的电磁法采用环状阵列电极构成测量探头,综合利用油气与水的导电和介电特性差异辨识井内流体,可以测量获得多相流体流动截面的清晰图像,为研究开发全新的流动成像测井系统奠定了良好基础。  相似文献   

5.
电磁全息测量数据包含两种模态,因此重建流动图像也需采用"双模"融合的敏感场。首先,通过电磁全息探测物理场分析,结合层析成像数学理论Radon反变换,从定解问题推导了全息测量敏感场函数表达式。其次,将有限元仿真计算得到的全息测量敏感场应用于模拟流动试验和全息成像。结果表明,基于复电位φ关于极径r的偏导数的全息测量敏感场契合了Radon反变换的数学表达,且充分体现了幅度、相位测量敏感性;反演所得流动图像更加吻合实际流型。全息测量敏感场构建有效克服了传统单模敏感场在计算精度和计算效率等方面的局限性。  相似文献   

6.
瞬变电磁波测井边界远探测方法研究   总被引:1,自引:0,他引:1       下载免费PDF全文
随钻地质导向钻井的关键在于对边界的探测,为提高仪器探测深度,时谐源激励的电磁波测井方法通常采用降低频率、增大源距的方式.瞬变电磁波测井信号源的突然关断产生会产生感应涡流,涡流随时间向地层深部扩散,与时谐源激励方式相比,其探测深度更大且测量过程不受信号源的干扰.因此,本文提出一种时间域瞬变电磁波测井边界远探测方法,采用余弦变换的数值滤波算法,模拟层状地层同轴发射接收线圈的瞬变电磁波测井响应,结果显示,地层电导率越大,电磁波传播速度越慢,测量晚期感应电动势与地层电导率线性相关;通过定义层状介质总场与线圈系所在当前层背景场的差值可方便提取界面信息,对界面的探测距离可达数十米;瞬变电磁波测井响应受源距的影响很小,为利用短源距实现远探测提供了可能.瞬变电磁波测井与时谐源电磁波测井相比优势明显,在电磁波测井领域中应用前景广阔.  相似文献   

7.
电磁波测井(EPT测井)是美国八十年代推出的一种新型地球物理测井方法。由于它能够提供用于计算冲洗带地层含水饱和度的电磁波传播时间信息,因而越来越受到人们关注。本文详细介绍EPT测井的原理及它的TPO解释方法。EPT测井首先承认地层(被测介质)是导电介质,因为只有在传导电流存在的情况下,才可能测出电磁波幅度的衰减,从而反推地层的有关性质。在传导电流和自由电荷为零的理想介质中,是测不出电磁波的幅度变化的,这也即为什么EPT测井不能在油基泥浆井中进行的原因。其次,EPT测井认为对发射天线发射的电磁波可以当作均匀平面简谐波来处理,这样在保证一定测量精度的条件下,使TPO解释公式变得简单实用。  相似文献   

8.
电磁成像正演仿真与传感器参数优化设计   总被引:1,自引:0,他引:1       下载免费PDF全文
应用有限元方法对电磁层析成像测量进行正演模拟,计算不同测量组合下各种典型流型的电势实部和虚部分布,考察电极个数、电极张角以及电极厚度对模拟测量响应的影响,并采用神经网络方法对传感器各结构参数进行优化设计.研究结果表明,介质分布对电势影响起主要作用,其中虚部分布特征更加明显;传感器各参数中电极数目对测量值影响较大,电极张角和电极厚度有一定影响.  相似文献   

9.
过套管电阻率测井通过测量套管壁电势实现测量地层视电阻率,不同的地层模型对过套管电阻率测井会有不同的测井响应. 在传输线方程系数中本文考虑了径向含多个界面地层的影响,对轴向(井轴方向)呈层状、径向阶跃变化地层模型,给出了电势分布递推公式,计算了过套管电阻率测井响应,绘出了测井响应曲线. 算例表明,该方法有较强的轴向地层界面分辨能力;对径向均匀地层计算结果非常接近地层模型电导率,能较好地反映实际地层特征;对高电导地层,低电导水泥环有明显的测井响应;当径向存在多个地层界面时计算结果则为径向各层电导率的综合反映,不是某一径向层的电导率;算例显示该方法有较快(算例运行时间在1s内)的计算速度.  相似文献   

10.
电磁波对轴对称二维层状介质的散射   总被引:41,自引:4,他引:41       下载免费PDF全文
本文讨论了任意N层平面分层,每层又可具有任意柱面分层的轴对称二维非均匀有耗介质结构,分析电磁波在其中的传播、散射特性.由于采用电磁场的数值模式匹配理论,径向用离散数值本征模式替代连续模,纵向则用模式匹配方法导出场递推关系的解析表达式,大大减少了数值计算量,而计算精度不变.场的数值结果可很好的再现地下电参数剖面.相应的数值分析程序可用于从感应测井到电磁波测井的计算机辅助解释.  相似文献   

11.
在分析地电场测量电极的极化电位产生机理的基础上,通过同类电极组合试验和不同类电极对比试验,对极化电位在地电场观测中的影响形成几点认识:①由于电极和土壤介质之间的接触电位和扩散电位的影响,极化电位的存在具有普遍性;②当工作状态稳定后,电极极化电位可能使得观测结果产生固定偏差,但在误差允许的范围内;③正常情况下,电极极化电位的存在不影响对地电场变化特性的观测。  相似文献   

12.
A numerical method is proposed for solving the problem of steady current flow. The electrodynamic model is replaced by the equivalent stationary charge distribution obtained by Poisson's analysis, in which the surface integral equation for field intensity is reduced to a set of simultaneous linear algebraic equations by means of the method of sub-areas. The solution of the set allows the calculation of an approximation for the charge density distribution on the discontinuity surfaces of conductivity. The method is valid for complex conductivities, whereby the apparent phase shift of IP can be calculated from the complex potential or field intensity. The phase shift anomaly calculated as an application is very similar to the corresponding frequency effect anomaly. The method allows the calculation of the mise-à-la-masse effect as a solution to a potential problem, in which the primary current electrode is located within the body to be surveyed.  相似文献   

13.
3D resistivity inversion using 2D measurements of the electric field   总被引:3,自引:0,他引:3  
Field and 'noisy' synthetic measurements of electric-field components have been inverted into 3D resistivities by smoothness-constrained inversion. Values of electrical field can incorporate changes in polarity of the measured potential differences seen when 2D electrode arrays are used with heterogeneous 'geology', without utilizing negative apparent resistivities or singular geometrical factors. Using both the X - and Y -components of the electric field as measurements resulted in faster convergence of the smoothness-constrained inversion compared with using one component alone. Geological structure and resistivity were reconstructed as well as, or better than, comparable published examples based on traditional measurement types. A 2D electrode grid (20 × 10), incorporating 12 current-source electrodes, was used for both the practical and numerical experiments; this resulted in 366 measurements being made for each current-electrode configuration. Consequently, when using this array for practical field surveys, 366 measurements could be acquired simultaneously, making the upper limit on the speed of acquisition an order of magnitude faster than a comparable conventional pole–dipole survey. Other practical advantages accrue from the closely spaced potential dipoles being insensitive to common-mode noise (e.g. telluric) and only 7% of the electrodes (i.e. those used as current sources) being susceptible to recently reported electrode charge-up effects.  相似文献   

14.
XIE Tao  LU Jun 《地震地质》2015,37(4):1125-1135
We calculate three-dimensional sensitivity coefficients distribution of apparent resistivity observation when Schlumberger array is used by using finite element method. Analysis results suggest that for the situation of one-dimensional positive or minus coefficient of surface medium, three-dimensional sensitivity coefficients distribution at surface shows similar patterns, and sensitivity coefficients distributions of different layered electric structures are also similar. There are two approximate ellipses at the two-dimensional surface plane between current electrodes and potential electrodes, where sensitivity coefficients are minus, and sensitivity coefficients at other areas are positive. Sensitivity coefficients at two approximate ellipses between current electrodes and potential electrodes are minus at the vertical section along monitoring line, while others are positive. From the three-dimensional view, minus sensitivity coefficients are at the two approximate half ellipsoids between current electrodes and potential electrodes when arrays are applied at surface. And coefficients near the electrodes are much greater than other areas. When resistivity of local areas at surface changes, we can qualitatively analyze the disturbing effects caused by the areas using three dimensional sensitivity coefficients distribution, and the analysis result can serve as reference for further experiment and numerical model quantitative analysis.  相似文献   

15.
The relation between the horizontal profiles of the subsurface resistivity and surface magnetotelluric data can be described by the input and output of a moving-average filter. The impulse response of this spatial filter, which characterizes the averaging process of the magnetotelluric measurements, is given by the sensitivity profile. Thus, the sensitivity analysis can provide insight into the characteristics of the measurements and hence the mechanism of the static effects. The sensitivity analysis presented here consists of constructing the vertical section of the sensitivity distribution using the finite-element method and then Fourier transforming the selected horizontal profiles. When the dipole is assumed for measuring the electric field, the static effects can be explained by the high-pass filter characteristics for the near-surface. When the electrode separation is taken into account, the sensitivity can be obtained by averaging the sensitivities for the dipoles over the horizontal distance equal to the electrode separation. Therefore, the higher-frequency components at each depth decrease with increasing electrode separation. Thus, although the static effects can be reduced simply by increasing the electrode separation, information on the resistivity variation at depth is also lost. However, such an adverse effect can be reduced by making the EMAP-type measurements followed by the spatial filtering of the profile data using the tapered weighting function.  相似文献   

16.
River bifurcations are key nodes within braided river systems controlling the flow and sediment partitioning and therefore the dynamics of the river braiding process. Recent research has shown that certain geometrical configurations induce instabilities that lead to downstream mid‐channel bar formation and the formation of bifurcations. However, we currently have a poor understanding of the flow division process within bifurcations and the flow dynamics in the downstream bifurcates, both of which are needed to understand bifurcation stability. This paper presents results of a numerical sensitivity experiment undertaken using computational fluid dynamics (CFD) with the purpose of understanding the flow dynamics of a series of idealized bifurcations. A geometric sensitivity analysis is undertaken for a range of channel slopes (0.005 to 0.03), bifurcation angles (22° to 42°) and a restricted set of inflow conditions based upon simulating flow through meander bends with different curvature on the flow field dynamics through the bifurcation. The results demonstrate that the overall slope of the bifurcation affects the velocity of flow through the bifurcation and when slope asymmetry is introduced, the flow structures in the bifurcation are modified. In terms of bifurcation evolution the most important observation appears to be that once slope asymmetry is greater than 0.2 the flow within the steep bifurcate shows potential instability and the potential for alternate channel bar formation. Bifurcation angle also defines the flow structures within the bifurcation with an increase in bifurcation angle increasing the flow velocity down both bifurcates. However, redistributive effects of secondary circulation caused by upstream curvature can very easily counter the effects of local bifurcation characteristics. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

17.
Aquifer vulnerability to pesticide migration through till aquitards   总被引:1,自引:0,他引:1  
This study investigates the influence of key factors-mainly recharge rate and degradation half-life--on downward migration of the widely used pesticide mecoprop (MCPP) through a typical clayey till aquitard. The study uses the numerical model FRAC3Dvs, which is a three-dimensional discrete fracture/matrix diffusion (DFMD) numerical transport model. The model was calibrated with laboratory and field data from a site near Havdrup, Denmark, but the overall findings are expected to be relevant to many other sites in similar settings. Fracture flow and MCPP transport parameters for the model were obtained through calibration using well-characterized laboratory experiments with large (0.5 m diameter by 0.5 m high) undisturbed columns of the fractured till and a field experiment. A second level of upscaling and sensitivity analysis was then carried out using data on hydraulic head, fracture spacing, and water budget from the field site. The simulations of downward migration of MCPP show that MCPP concentration and mass flux into the underlying aquifer, and hence the aquifer vulnerability to this pesticide compound, is mainly dependent on the degradation rate of the pesticide, the overall aquitard water budget, and the ground water recharge rate into the aquifer. The influence of flow rate, matrix diffusion, and degradation rate are intertwined. This results in one to four orders of magnitude higher MCPP flux into the aquifer from aquifer recharge rates of 20 and 120 mm/yr, respectively, for no degradation and MCPP half-life of 0.5 yr. From a sensitivity analysis it was found that the range of MCPP flux into the aquifer varied less than one order of magnitude due to (1) changing fracture spacing from 1 to 10 m, or (2) preferential flow along inclined thin sand layers, which represent common conditions for the current and other settings of clayey till in Denmark and other glaciated areas in Europe and North America. The results indicate that for aquifers overlain by fractured clayey tills, the vulnerability to contamination with pesticides (pesticide flux into the aquifer) and other widespread agricultural contaminants is going to vary strongly in the watershed as a function of the distribution of aquitard water budget (flow rate) and aquitard redox environment (controlling contaminant degradation rates), even if the thickness of the till is relatively constant. DFMD modeling of cause-effect relationships within such systems has great potential to support decisions in planning, regulation, and contaminant remediation.  相似文献   

18.
Three choices of control variables for meteorological variational analysis (3DVAR or 4DVAR) are associated with horizontal wind: (1) streamfunction and velocity potential, (2) eastward and northward velocity, and (3) vorticity and divergence. This study shows theoretical and numerical differences of these variables in practical 3DVAR data assimilation through statistical analysis and numerical experiments. This paper demonstrates that (a) streamfunction and velocity potential could potentially introduce analysis errors; (b) A 3DVAR using velocity or vorticity and divergence provides a natural scale dependent influence radius in addition to the covariance; (c) for a regional analysis, streamfunction and velocity potential are retrieved from the background velocity field with Neumann boundary condition. Improper boundary conditions could result in further analysis errors; (d) a variational data assimilation or an inverse problem using derivatives as control variables yields smoother analyses, for example, a 3DVAR using vorticity and divergence as controls yields smoother wind analyses than those analyses obtained by a 3DVAR using either velocity or streamfunction/velocity potential as control variables; and (e) statistical errors of higher order derivatives of variables are more independent, e.g., the statistical correlation between U and V is smaller than the one between streamfunction and velocity potential, and thus the variables in higher derivatives are more appropriate for a variational system when a cross-correlation between variables is neglected for efficiency or other reasons. In summary, eastward and northward velocity, or vorticity and divergence are preferable control variables for variational systems and the former is more attractive because of its numerical efficiency. Numerical experiments are presented using analytic functions and real atmospheric observations.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号