位场异常分离的递减半径迭代法

本文提出了位场区域-剩余异常分离的空间域递减半径迭代法.由给定半径的圆周上八点位场值的算术平均导出一个新的八点圆周平均公式,它是一个由重磁数据计算区域异常的滤波器.该滤波器的传递函数有一主瓣和多个旁瓣,且半径越大,旁瓣越多,滤波器特性越差.由从大到小不同半径对应的传递函数的连乘积构造了递减半径迭代传递函数,它以大半径为参变量.递减半径迭代传递函数类似低通滤波器,其截止波数与大半径成反比.由递减半径迭代传递函数,给出空间域分离区域异常的递减半径线性迭代法,由重磁数据减去区域异常求得剩余异常.进一步,通过构造非线性修正系数,把空间域递减半径线性迭代法中线性迭代公式变成非线性迭代公式,得到空间域递减半径非线性迭代法.重磁理论模型及安徽泥河铁矿重磁资料的试验表明,空间域递减半径非线性迭代法可有效地压制假异常和高频干扰以及减少异常畸变,并能有效地分离出区域异常和剩余异常.     

连续胶结声速理论胶结半径表达式的一般性推广及应用

针对Dvorkin 和Nur （1996） 提出的连续胶结声速理论（CCT）中的胶结半径表达式仅适用于孔隙式胶结疏松砂岩的局限性,通过胶结物与颗粒之间的空间几何关系,在胶结物中心厚度不为零的前提下,推导了胶结物胶结半径的一般表达式,并考察了胶结物中心厚度对疏松砂岩声波速度的影响,探讨了该表达式在基底式胶结疏松砂岩纵横波速度预测中的应用.结果表明：随胶结物中心厚度变大,纵、横波速度均减小,且孔隙度越大,胶结物含量越少,速度下降越快;孔隙度或胶结物含量不变时,纵横波速比随胶结物中心厚度的增大而增大;纵横波速比随孔隙度的变化规律与胶结物中心厚度有关,当胶结物中心厚度较小时,孔隙度增大,纵横波速比略微减小,而当胶结物中心厚度较大时,孔隙度的增大会引起纵横波速比略微增大;对于基底式胶结疏松砂岩,使用推广后表达式修正的CCT模型相对原CCT模型能够更好地预测其声波速度,同时发现,利用其纵横波速比无法预测由孔隙度或胶结物含量变化引起的疏松砂岩软硬的变化,但是在胶结物含量或孔隙度不变的情况下,纵横波速比可指示由胶结物中心厚度变化引起的疏松砂岩软硬的变化.     

Uniform head in horizontal and vertical wells

The steady-state head within a fully penetrating well may be estimated by evaluating the Thiem equation at the radius of the well. A method is presented here to extend results from the Thiem equation to horizontal wells and to partially penetrating wells. The particular model used in this investigation is based upon the analytic element method; it accurately reproduces a boundary condition of uniform head along the cylindrical surface at the perforated face of the well. This model is exercised over a representative range of parameters including the well's length, radius, and pumping rate, and the aquifer's hydraulic conductivity and thickness. Results are presented in a set of figures and tables that compare the well's drawdown to the drawdown that would have been obtained using the Thiem solution with the same pumping rate and radius. A methodology is presented to estimate the head within a horizontal or partially penetrating well by adding a correction term to results that can be readily obtained from computer models of vertical fully penetrating wells. This approach may also be used to contrast the differences in head between horizontal and vertical wells of various lengths, radii, and placement elevations.     

Determination of the formation temperature from shut-in logs: Estimation of the radius of thermal influence

As a consequence of drilling and completion operations most wells have reduced permeability near the borehole (skin zone). This may result in a corresponding change of thermal conductivity of formations near the borehole. At the same time the mathematical models, which describe the process of thermal recovery of formations are based on the assumption that the thickness of the radius of heated zone (during circulation of the drilling mud) is significantly larger than the well radius (or skin zone). It is shown below that the radius of the heated zone (radius of investigation or radius of thermal influence) can be approximated by simple formulas. It is shown that for six field bottom-hole temperature surveys the radius of thermal influence (RTI) is too small, and for this reason, the duration of the drilling mud circulation period should increased.     

填土地基强夯加固效果研究

针对回填土地基,通过量纲分析得到强夯加固深度和累积夯沉量的计算公式,结合已有工程资料对所提出的计算公式进行验证。结果表明:该公式能够较为精确地预测填土地基的加固深度与累积夯沉量;分析表明:累积夯沉量与加固深度比值随着诸多要素而变化,仅在确定的落距和夯锤半径范围内的才可以作为常数进行推广;最后分析夯锤半径及厚度对强夯效果的影响。     

Radial thickness variation in impact crater ejecta: implications for lunar basin deposits

A review of cratering data and available semi-empirical calculations suggests that the variation of ejecta thickness,t, with increasing range from lunar craters may be approximately modelled by the expression: t=0.14R^0.74(r/R^?3.0 wherer is range from the center of the crater andR, the crater radius, all in meters. This equation has been used to estimate the thickness of ejecta deposits at each of the Apollo sites contributed from the large multi-ringed frontside lunar basins. Predicted average thickness of Imbrium ejecta at Apollo 15 is 812 m; at Apollo 14, 130 m; at Apollo 17, 102 m; and at Apollo 16, 50 m. Since the sequence of formation of these basins is known, the stratigraphic column resulting from superimposed ejecta blankets can be calculated. Results suggest that pre-Nubium crustal material at upland Apollo sites lies at depths greater than 280 (Apollo 14) to 1940 m (Apollo 17). Predicted stratigraphic sections for the Apollo sites are tabulated.     

A gradient‐ratio method for the semi‐automatic interpretation of gravity map data sets

The application of semi‐automatic interpretation techniques to potential field data can be of significant assistance to a geophysicist. This paper generalizes the magnetic vertical contact model tilt‐depth method to gravity data using a vertical cylinder and buried sphere models. The method computes the ratio of the vertical to the total horizontal derivative of data and then identifies circular contours within it. Given the radius of the contour and the contour value itself, the depth to the source can be determined. The method is applied both to synthetic and gravity data from South Africa. The Matlab source code can be obtained from the author upon request.     

EFFECTS OF WELL CASING ON POTENTIAL FIELD MEASUREMENTS USING DOWNHOLE CURRENT SOURCES1

A mathematical formulation for the electric potential from point current-sources coaxial with a metal casing has been obtained. The excitation caused by the axial point-sources will produce currents in the pipe. By assuming that the pipe can be divided into many cylindrical ring segments with constant axially-directed current, the solution of the fields inside and outside the pipe can be formulated in an integral form. The integral equation applied to the segmented pipe yields a set of simultaneous linear equations which are solved for the currents in the pipe; these are then used to calculate the potentials anywhere outside the pipe in the medium. This solution has been used to study the distribution of the potentials in a half-space for a single current-source at and beyond the bottom of a finite length of casing. For a casing 0.1 m in radius and 0.006 m in wall thickness with a conductivity of 10⁶ S/m, in a half-space of 10^-2 S/m, it was found that only in a region very near the pipe does the pipe exert substantial influence on the fields of a point-source 100 casing diameters beyond the end of the pipe. It appears that cross-hole resistivity surveys can be implemented without corrections for the casing if the source is located at least 50–100 casing diameters beyond the end of the casing. Hole-to-surface surveys are much more affected by the pipe. For a pipe-source separation of 100 casing diameters, the surface measurements must not be closer than a half pipe length for a 5% or less field distortion.     

平面P波入射下深埋圆形复合式衬砌隧道抗减震机理研究

基于Fourier-Bessel级数展开法,研究深埋圆形三层复合式衬砌洞室在平面P波入射下的动应力集中问题,并给出三层衬砌洞室动应力集中系数级数解析解;依托某IX度地震区管道隧道实际工程,分析不同衬砌刚度组合和厚度组合对洞室动应力集中系数的影响。研究表明:注浆加固洞室围岩和设置减震层都可以降低二次衬砌动应力集中系数;增大围岩注浆区弹性模量和厚度,有利于减小衬砌动应力集中系数,最优围岩注浆区厚度为1倍洞室净空半径;减震层弹性模量降低,减震层厚度增大,二次衬砌动应力集中系数变小,减震层弹性模量宜低于围岩弹模1/20,最优减震层厚度宜取1/50的洞室内净空半径。最后针对实际管道隧道抗减震技术,考虑围岩稳定性,提出"围岩-加固圈-减震层-衬砌"新型减震结构,分析结果表明:对比其他三种抗减震措施,新型减震结构的减震效果最好。     

Radial groundwater flow to a finite diameter well in a leaky confined aquifer with a finite‐thickness skin

A mathematical model that describes the drawdown due to constant pumpage from a finite radius well in a two‐zone leaky confined aquifer system is presented. The aquifer system is overlain by an aquitard and underlain by an impermeable formation. A skin zone of constant thickness exists around the wellbore. A general solution to a two‐zone leaky confined aquifer system in Laplace domain is developed and inverted numerically to the time‐domain solution using the modified Crump (1976) algorithm. The results show that the drawdown distribution is significantly influenced by the properties and thickness of the skin zone and aquitard. The sensitivity analyses of parameters of the aquifer and aquitard are performed to illustrate their effects on drawdowns in a two‐zone leaky confined aquifer system. For the negative‐skin case, the drawdown is very sensitive to the relative change in the formation transmissivity. For the positive‐skin case, the drawdown is also sensitive to the relative changes in the skin thickness, and both the skin and formation transmissivities over the entire pumping period and the well radius and formation storage coefficient at early pumping time. Copyright © 2009 John Wiley & Sons, Ltd.     

The impact of reservoir scale on amplitude variation with offset

To investigate the vertical and horizontal impact of reservoir scale on the amplitude‐versus‐offset characteristics, we conduct seismic numerical simulations on models containing spatially confined lithologic units with different scales. We find that the reservoir scale has a nonlinear effect on the amplitude‐versus‐offset intercepts and gradients. As the reservoir width increases, amplitude‐versus‐offset intercept and gradient both first increase, then decrease, and finally remain stable. The amplitude‐versus‐offset intercept is maximum when the reservoir width is 80% larger than the Fresnel zone radius, whereas the amplitude‐versus‐offset gradient peaks at 1.5 times the Fresnel zone radius. Both amplitude‐versus‐offset intercept and gradient are approximately proportional to the reservoir width prior to reaching their maxima. When the lateral extent of the reservoir is more than three times the Fresnel zone radius, the amplitude‐versus‐offset attributes are constant. Modelling the reservoir thickness shows that intercept and gradient behave in a manner similar to that of tuning of thin beds. Both the amplitude‐versus‐offset intercept and gradient first increase and then decrease with the thickness, peaking at the tuning thickness. The thickness contribution to amplitude‐versus‐offset variations is negligible when the thickness is larger than 1.6 times of the tuning thickness. Considering the magnitude of the changes in amplitude‐versus‐offset intercept and gradient caused by reservoir scale, the width causes a maximum 433% intercept increase and a 344% gradient increase, whereas the thickness causes a maximum 100% intercept increase and a 73% gradient increase. Cross‐plotting the amplitude‐versus‐offset intercept and gradient shows the reservoir scale change gives rise to an anti‐clockwise spiraling effect. In conclusion, the lateral and vertical extents of the reservoir both play an important role in amplitude variation with offset. Our analysis shows that the lateral reservoir extent has a larger impact on the amplitude variation with offset than the vertical tuning effect.     

Quantifying the spatial variability of melting seasonal ground ice and its influence on potential evapotranspiration spatial variability in a boreal peatland

Melting seasonal ground ice (SGI) in western Boreal Plains (WBP) peatlands can reduce the available energy at the surface by reducing potential evapotranspiration (PET). PET often exceeds annual precipitation in the WBP. Including this effect in hydrological models may be important in assessing water deficits. However, SGI melt and the timing of ice-free conditions vary spatially, which suggests PET spatial variability could be influenced by SGI. Understanding this potential linkage can help improve site scale PET in peatland hydrological models. The objectives of this paper were (a) to quantify the effect of ice thickness and melt rate on peatland PET; (b) quantify the spatial variability of SGI thickness and melt rate across spatial scales; and (c) assess how/if spatial variability in SGI thickness/melt rate affects site scale PET. Results from the sensitivity analysis indicated that SGI thickness had a bigger impact on reducing PET compared with the melt rate. Two SGI thickness values were used that were observed on site: 0.32 m, which was measured in a more treed area, and 0.18 m, which was in a more open area. The 0.32 m had an average PET reduction of 14 mm (±0.7), over the month of May, compared with 9 mm (±1 mm) when there was 0.18 m of SGI, which are 13.7 and 8.8% reductions, respectively. SGI thickness and melt rate, both exhibited large- and small-scale spatial variability. At the large scale, spatial patterns in SGI thickness appeared to be influenced by extensive shading from the adjacent hillslopes. Small scale, SGI thickness may be a function of tree proximity and the snowpack. Finally, net radiation, rather than SGI, appeared to be the main driver behind PET spatial variability. This work enhances our conceptual understanding of the role of SGI in WBP peatlands. Future work can use the findings to better inform peatland hydrological models, allowing for better representation of peatlands in regional-scale models.     

Vertical vibration of ring foundations

A comprehensive analysis is made of the harmonic response of vertically excited, massless, rigid ring foundations supported at the surface of an homogeneous elastic halfspace. The parameters considered include the thickness to radius ratio for the ring, the frequency of the exciting force and of the resulting steady-state response, and Poisson's ratio for the supporting medium. The response quantities examined include the stiffness and damping coefficients of the system in an equivalent spring–dashpot representation of the supporting medium, the displacements of the foundation and of points on the ground surface away from the foundation, and the normal pressure at the foundation–medium interface. The results in each case are compared with those obtained for a solid disk having the same radius as the outer radius of the ring, and a simple approximate model is used to interpret the results and to provide insight into the action of the system. The reported data are computed by a method of analysis that takes due account of the mixed boundary conditions at the surface of the halfspace, and are believed to be of high accuracy.     

Parameter Estimation/Sensitivity Analysis for an Aquifer Test with Skin Effect

Aquifer information carried by aquifer test data may be affected by the presence of a finite thickness skin around the wellbore. The mathematical treatment for an aquifer accounting for the skin zone can be characterized by five parameters, that is, the outer radius of the skin zone and the transmissivity and storativity for each of the skin and aquifer zones. Sensitivity analysis was performed to examine the ground water flow behavior in the skin and aquifer zones in terms of the constant-head test (CHT) data. The simulated annealing procedure was applied to simultaneously determine the skin and aquifer parameters from the analysis of CHT data. Toward the previously mentioned goals, four suites of CHT data were analyzed in this article. The analyses of wellbore flow rate at the test well and the specific drawdown at the observation well gave accurate estimates for the skin and aquifer parameters, respectively. Only the skin thickness and both the skin and the aquifer diffusivities could be accurately estimated from the analysis of drawdown data in the observation well. The estimates for all skin and aquifer parameters from the composite analysis of flow rate and drawdown data were the most accurate. The results of sensitivity analyses and parameter estimations provide instructive references in the analysis of the skin-affected CHT data.     

井中磁源瞬变电磁响应特征研究

井中瞬变电磁波勘探是一个全空间地球物理场问题.采用Gaver-Stehfest逆拉氏变换方法,正演计算了瞬变信号激励下接收线圈上的电磁场响应.分析了包含井眼泥浆、套管、水泥环和地层的轴对称多层介质模型的电磁场响应特征,考察了各层介质参数对井中瞬变电磁响应的影响.不同电导率井眼泥浆的电磁场响应衰减曲线表明,井眼泥浆电导率...     

Slug Tests in Wells Screened Across the Water Table: Some Additional Considerations

The majority of slug tests done at sites of shallow groundwater contamination are performed in wells screened across the water table and are affected by mechanisms beyond those considered in the standard slug‐test models. These additional mechanisms give rise to a number of practical issues that are yet to be fully resolved; four of these are addressed here. The wells in which slug tests are performed were rarely installed for that purpose, so the well design can result in problematic (small signal to noise ratio) test data. The suitability of a particular well design should thus always be assessed prior to field testing. In slug tests of short duration, it can be difficult to identify which portion of the test represents filter‐pack drainage and which represents formation response; application of a mass balance can help confirm that test phases have been correctly identified. A key parameter required for all slug test models is the casing radius. However, in this setting, the effective casing radius (borehole radius corrected for filter‐pack porosity), not the nominal well radius, is required; this effective radius is best estimated directly from test data. Finally, although conventional slug‐test models do not consider filter‐pack drainage, these models will yield reasonable hydraulic conductivity estimates when applied to the formation‐response phase of a test from an appropriately developed well.     

地球外核的电荷密度

为计算地球磁极处的磁感应强度,建立地球的磁场是由带电的地球外核的旋转产生的模型.先根据毕奥-萨伐尔定律计算球形模型绕自转轴旋转时在自转轴直径上产生的磁感应强度;再利用已知的地球外核的内外半径及地球半径和磁极处的磁感应强度值,计算出地球外核的电荷体密度及面密度.结果表明:若外核的电荷呈均匀的体密度分布,则其电荷体密度为3.5507 C/m³;若外核的电荷均匀分布在外核的外表面,则其面密度为2.4581×10⁶ C/m².通过地球表面的磁感应强度信息利用物理规律和地球物理数据推测地球内部难以直接进行探测的相关信息,具有实际意义.根据地震学方法对地球外核厚度、转向等变化的最新研究数据按该文模型可推测地球磁场强度、极性等的变化.而地球磁场的变化对地球上的人类生活颇有影响.     

Estimating the thickness of a faulted mechanical layer based on fault spacing: An example from Miyako‐jima Island,southern Ryukyu Arc,Japan

The thickness of the mechanical layer that hosts a group of faults can be estimated from the spacing of saturated faults (i.e. the constant spacing between faults when the fault system is fully developed and has attained its final fault density). We measured fault spacing for a group of saturated active normal faults on Miyako‐jima Island (southern Ryukyu Arc, Japan) and estimated the thickness of the faulted mechanical layer. The measured fault spacing is 1.30 ±0.14 km, and the thickness of the mechanical layer is < 3 km, which is the upper limit of the seismogenic zone. This faulted mechanical layer corresponds to a sedimentary layer in which earthquakes cannot occur. Results indicate that the shallow (< 3 km depth) normal faults on Miyako‐jima Island do not have the potential to cause medium‐size earthquakes as individual faults. The origin of the shallow normal faults might be related to the presence of a larger‐scale, deeper fault. The results indicate that fault spacing provides important information on the potential magnitude of earthquakes associated with active faults.     

Inverse problem for the current loop model: Possibilities and restrictions

The possibilities of determining arbitrary current loop parameters based on the spatial structures of the magnetic field components generated by this loop on a sphere with a specified radius have been considered with the use of models. The model parameters were selected such that anomalies created by current loops on a sphere with a radius of 6378 km would be comparable in value with the different-scale anomalies of the observed main geomagnetic field (MGF). The least squares method was used to solve the inverse problem. Estimates close to the specified values were obtained for all current loop parameters except the current strength and radius. The radius determination error can reach ±120 km; at the same time, the magnetic moment value is determined with an accuracy of ±1%. The resolvability of the current force and radius can to a certain degree be improved by decreasing the observation sphere radius such that the ratio of the source distance to the current loop radius would be at least smaller than eight, which can be difficult to reach when modeling MGF.     

Oblique flexural-gravity wave scattering by a submerged semi-circular ridge

This study examines the flexural-gravity wave scattering by a semi-circular ridge submerged in the ocean covered by a thin ice sheet. The ice sheet is treated as an elastic plate and its deflection is modelled using the Euler-Bernoulli beam equation. The fluid motion is described by linear potential theory. The series solutions of velocity potentials for obliquely and normally incident waves are both developed using multipole expansions. The unknowns in the series solutions are determined in terms of the impermeable condition on the ridge surface. Accurate results for the reflection and transmission coefficients of ridge and the deflection of ice sheet are presented. The effects of ridge radius, ice thickness, wave frequency and wave incident angle on the hydrodynamic quantities are shown.