共查询到18条相似文献,搜索用时 187 毫秒
1.
2.
3.
4.
流动成像测井研究进展 总被引:2,自引:3,他引:2
流动成像测井属于对非均匀介质动态的非线性测量,通过对油气井内多相流体信息的投影扫描和反演处理,实时成像显示流体分布及流动状况,国外研制的电导法,电容法仪器采用电导探针或电容元件构成阵列测量探头,分别利用油气与水的导电特性和介电特性差异辨识井内流体,对流动截面的测量局限于个别点上,已在水平井中见到应用效果,但对物场信息投影测量的数据量和分辨率未能满足成像要求。我国提出的电磁法采用环状阵列电极构成测量探头,综合利用油气与水的导电和介电特性差异辨识井内流体,可以测量获得多相流体流动截面的清晰图像,为研究开发全新的流动成像测井系统奠定了良好基础。 相似文献
5.
6.
随钻地质导向钻井的关键在于对边界的探测,为提高仪器探测深度,时谐源激励的电磁波测井方法通常采用降低频率、增大源距的方式.瞬变电磁波测井信号源的突然关断产生会产生感应涡流,涡流随时间向地层深部扩散,与时谐源激励方式相比,其探测深度更大且测量过程不受信号源的干扰.因此,本文提出一种时间域瞬变电磁波测井边界远探测方法,采用余弦变换的数值滤波算法,模拟层状地层同轴发射接收线圈的瞬变电磁波测井响应,结果显示,地层电导率越大,电磁波传播速度越慢,测量晚期感应电动势与地层电导率线性相关;通过定义层状介质总场与线圈系所在当前层背景场的差值可方便提取界面信息,对界面的探测距离可达数十米;瞬变电磁波测井响应受源距的影响很小,为利用短源距实现远探测提供了可能.瞬变电磁波测井与时谐源电磁波测井相比优势明显,在电磁波测井领域中应用前景广阔. 相似文献
7.
电磁波测井(EPT测井)是美国八十年代推出的一种新型地球物理测井方法。由于它能够提供用于计算冲洗带地层含水饱和度的电磁波传播时间信息,因而越来越受到人们关注。本文详细介绍EPT测井的原理及它的TPO解释方法。EPT测井首先承认地层(被测介质)是导电介质,因为只有在传导电流存在的情况下,才可能测出电磁波幅度的衰减,从而反推地层的有关性质。在传导电流和自由电荷为零的理想介质中,是测不出电磁波的幅度变化的,这也即为什么EPT测井不能在油基泥浆井中进行的原因。其次,EPT测井认为对发射天线发射的电磁波可以当作均匀平面简谐波来处理,这样在保证一定测量精度的条件下,使TPO解释公式变得简单实用。 相似文献
8.
9.
实际井眼条件下过套管电阻率测井响应的传输线方程正演算法 总被引:3,自引:4,他引:3
过套管电阻率测井通过测量套管壁电势实现测量地层视电阻率,不同的地层模型对过套管电阻率测井会有不同的测井响应. 在传输线方程系数中本文考虑了径向含多个界面地层的影响,对轴向(井轴方向)呈层状、径向阶跃变化地层模型,给出了电势分布递推公式,计算了过套管电阻率测井响应,绘出了测井响应曲线. 算例表明,该方法有较强的轴向地层界面分辨能力;对径向均匀地层计算结果非常接近地层模型电导率,能较好地反映实际地层特征;对高电导地层,低电导水泥环有明显的测井响应;当径向存在多个地层界面时计算结果则为径向各层电导率的综合反映,不是某一径向层的电导率;算例显示该方法有较快(算例运行时间在1s内)的计算速度. 相似文献
10.
本文讨论了任意N层平面分层,每层又可具有任意柱面分层的轴对称二维非均匀有耗介质结构,分析电磁波在其中的传播、散射特性.由于采用电磁场的数值模式匹配理论,径向用离散数值本征模式替代连续模,纵向则用模式匹配方法导出场递推关系的解析表达式,大大减少了数值计算量,而计算精度不变.场的数值结果可很好的再现地下电参数剖面.相应的数值分析程序可用于从感应测井到电磁波测井的计算机辅助解释. 相似文献
11.
在分析地电场测量电极的极化电位产生机理的基础上,通过同类电极组合试验和不同类电极对比试验,对极化电位在地电场观测中的影响形成几点认识:①由于电极和土壤介质之间的接触电位和扩散电位的影响,极化电位的存在具有普遍性;②当工作状态稳定后,电极极化电位可能使得观测结果产生固定偏差,但在误差允许的范围内;③正常情况下,电极极化电位的存在不影响对地电场变化特性的观测。 相似文献
12.
L. ESKOLA 《Geophysical Prospecting》1979,27(3):616-627
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.
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.
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.
YUTAKA SASAKI 《Geophysical Prospecting》1989,37(4):395-406
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.
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. 相似文献