Apparent-density mapping using entropic regularization

Geophysics ◽  
2007 ◽  
Vol 72 (4) ◽  
pp. I51-I60 ◽  
Author(s):  
João B. C. Silva ◽  
Francisco S. Oliveira ◽  
Valéria C. F. Barbosa ◽  
Haroldo F. Campos Velho

We present a new apparent-density mapping method on the horizontal plane that combines the minimization of the first-order entropy with the maximization of the zeroth-order entropy of the estimated density contrasts. The interpretation model consists of a grid of vertical, juxtaposed prisms in both horizontal directions. We assume that the top and the bottom of the gravity sources are flat and horizontal and estimate the prisms’ density contrasts. The minimization of the first-order entropy favors solutions presenting sharp borders, and the maximization of the zeroth-order entropy prevents the tendency of the source estimate to become a single prism. Thus, a judicious combination of both constraints may lead to solutions characterized by regions with virtually constant estimated density contrasts separated by sharp discontinuities. We apply our method to synthetic data from simulated intrusive bodies in sediments that present flat and horizontal tops. By comparing our results with those obtained with the smoothness constraint, we show that both methods produce good and equivalent locations of the sources’ central positions. However, the entropic regularization delineates the boundaries of the bodies with greater resolution, even in the case of 100-m-wide bodies separated by a distance as small as [Formula: see text]. Both the proposed and the global smoothness constraints are applied to real anomalies from the eastern Alps and from the Matsitama intrusive complex, northeastern Botswana. In the first case, the entropic regularization delineates two sources, with a horizontal and nearly flat top being consistent with the known geologic information. In the second case, both constraints produce virtually the same estimate, indicating, in agreement with results of synthetic tests, that the tops of the sources are neither flat nor horizontal.

Geophysics ◽  
2010 ◽  
Vol 75 (2) ◽  
pp. L39-L50 ◽  
Author(s):  
João B. Silva ◽  
Suzan S. Vasconcelos ◽  
Valeria C. Barbosa

A new apparent-magnetization mapping method on the horizontal plane combines minimization of first-order entropy with maximization of zeroth-order entropy of the estimated magnetization. The interpretation model is a grid of vertical, juxtaposed prisms in both horizontal directions. To estimate the magnetization of the prisms, assume that the top and bottom of the magnetic sources are horizontal. Minimization of the first-order entropy favors solutions with sharp borders, and the maximization of zeroth-order entropy prevents the tendency of the estimated source to become a single prism with large magnetization. Thus, a judicious combination of both constraints can lead to solutions characterized by regions with virtually constant magnetizations separated by sharp discontinuities. This is applied to synthetic data from simulated intrusive bodies in sediments that have horizontal tops. By comparing the results with those obtained with the common Tikhonov regularization (smoothness constraint) method, it is shown that both methods produce good and equivalent locations of the central positions of the sources. However, entropic regularization delineates the boundaries of the bodies with greater detail. Both the proposed and the smoothness constraints are applied to real anomaly data over a magnetic skarn in Butte Valley, Nevada, U.S.A. Entropic regularization produced an estimated magnetization distribution with sharper boundaries, smaller volume, and higher apparent magnetization as compared with results produced by incorporating the smoothness constraint.


Geophysics ◽  
2010 ◽  
Vol 75 (3) ◽  
pp. I29-I35 ◽  
Author(s):  
João B. Silva ◽  
Alexandre S. Oliveira ◽  
Valéria C. Barbosa

We have developed a gravity interpretation method for estimating the discontinuous basement relief of a sedimentary basin. The density contrast between the basement and the sediments is assumed to be known, and it could be either constant or vary monotonically with depth. The interpretation model consists of a set of vertical, juxtaposed prisms, whose thicknesses are the parameters to be estimated. We used the entropic regularization that combines the minimization of the first-order entropy measure with the maximization of the zeroth-order entropy measure of the solution vector. We validated the method by applying it to synthetic data produced by a simulated basin bordered by high-angle step faults; we obtained a good definition of the relief, particularly of the discontinuities. We also applied the method to a profile across the Büyük Menderes Valley in West Turkey and obtained a solution exhibiting a gravity fault with large slip on the northern border of the valley. When applied to the interpretation of a discontinuous basement relief, the method has a better performance than the global smoothness method. It is comparable to the weighted smoothness method, but it does not require the a priori knowledge about the maximum basin depth.


Author(s):  
Amarjot Singh Bhullar ◽  
Gospel Ezekiel Stewart ◽  
Robert W. Zimmerman

Abstract Most analyses of fluid flow in porous media are conducted under the assumption that the permeability is constant. In some “stress-sensitive” rock formations, however, the variation of permeability with pore fluid pressure is sufficiently large that it needs to be accounted for in the analysis. Accounting for the variation of permeability with pore pressure renders the pressure diffusion equation nonlinear and not amenable to exact analytical solutions. In this paper, the regular perturbation approach is used to develop an approximate solution to the problem of flow to a linear constant-pressure boundary, in a formation whose permeability varies exponentially with pore pressure. The perturbation parameter αD is defined to be the natural logarithm of the ratio of the initial permeability to the permeability at the outflow boundary. The zeroth-order and first-order perturbation solutions are computed, from which the flux at the outflow boundary is found. An effective permeability is then determined such that, when inserted into the analytical solution for the mathematically linear problem, it yields a flux that is exact to at least first order in αD. When compared to numerical solutions of the problem, the result has 5% accuracy out to values of αD of about 2—a much larger range of accuracy than is usually achieved in similar problems. Finally, an explanation is given of why the change of variables proposed by Kikani and Pedrosa, which leads to highly accurate zeroth-order perturbation solutions in radial flow problems, does not yield an accurate result for one-dimensional flow. Article Highlights Approximate solution for flow to a constant-pressure boundary in a porous medium whose permeability varies exponentially with pressure. The predicted flowrate is accurate to within 5% for a wide range of permeability variations. If permeability at boundary is 30% less than initial permeability, flowrate will be 10% less than predicted by constant-permeability model.


1998 ◽  
Vol 13 (39) ◽  
pp. 3169-3177 ◽  
Author(s):  
IOANNIS GIANNAKIS ◽  
K. KLEIDIS ◽  
A. KUIROUKIDIS ◽  
D. PAPADOPOULOS

We study string propagation in an anisotropic, cosmological background. We solve the equations of motion and the constraints by performing a perturbative expansion of the string coordinates in powers if c2 — the worldsheet speed of light. To zeroth order the string is approximated by a tensionless string (since c is proportional to the string tension T). We obtain exact, analytical expressions for the zeroth- and first-order solutions and we discuss some cosmological implications.


1988 ◽  
Vol 53 (2) ◽  
pp. 554-570 ◽  
Author(s):  
Kosta Došen ◽  
Peter Schroeder-Heister

This paper is meant to be a comment on Beth's definability theorem. In it we shall make the following points.Implicit definability as mentioned in Beth's theorem for first-order logic is a special case of a more general notion of uniqueness. If α is a nonlogical constant, Tα a set of sentences, α* an additional constant of the same syntactical category as α and Tα, a copy of Tα with α* instead of α, then for implicit definability of α in Tα one has, in the case of predicate constants, to derive α(x1,…,xn) ↔ α*(x1,…,xn) from Tα ∪ Tα*, and similarly for constants of other syntactical categories. For uniqueness one considers sets of schemata Sα and derivability from instances of Sα ∪ Sα* in the language with both α and α*, thus allowing mixing of α and α* not only in logical axioms and rules, but also in nonlogical assumptions. In the first case, but not necessarily in the second one, explicit definability follows. It is crucial for Beth's theorem that mixing of α and α* is allowed only inside logic, not outside. This topic will be treated in §1.Let the structural part of logic be understood roughly in the sense of Gentzen-style proof theory, i.e. as comprising only those rules which do not specifically involve logical constants. If we restrict mixing of α and α* to the structural part of logic which we shall specify precisely, we obtain a different notion of implicit definability for which we can demonstrate a general definability theorem, where a is not confined to the syntactical categories of nonlogical expressions of first-order logic. This definability theorem is a consequence of an equally general interpolation theorem. This topic will be treated in §§2, 3, and 4.


1983 ◽  
Vol 27 (01) ◽  
pp. 13-33
Author(s):  
Francis Noblesse

A new slender-ship theory of wave resistance is presented. Specifically, a sequence of explicit slender-ship wave-resistance approximations is obtained. These approximations are associated with successive approximations in a slender-ship iterative procedure for solving a new (nonlinear integro-differential) equation for the velocity potential of the flow caused by the ship. The zeroth, first, and second-order slender-ship approximations are given explicitly and examined in some detail. The zeroth-order slender-ship wave-resistance approximation, r(0) is obtained by simply taking the (disturbance) potential, ϕ, as the trivial zeroth-order slender-ship approximation ϕ(0) = 0 in the expression for the Kochin free-wave amplitude function; the classical wave-resistance formulas of Michell [1]2 and Hogner [2] correspond to particular cases of this simple approximation. The low-speed wave-resistance formulas proposed by Guevel [3], Baba [4], Maruo [5], and Kayo [6] are essentially equivalent (for most practical purposes) to the first-order slender-ship low-Froude-number approximation, rlF(1), which is a particular case of the first-order slender-ship approximation r(1): specifically, the first-order slender-ship wave-resistance approximation r(1) is obtained by approximating the potential ϕ in the expression for the Kochin function by the first-order slender-ship potential ϕ1 whereas the low-Froude-number approximation rlF(1) is associated with the zero-Froude-number limit ϕ0(1) of the potentialϕ(1). A major difference between the first-order slender-ship potential ϕ(1) and its zero-Froude-number limit ϕ0(1) resides in the waves that are included in the potential ϕ(1) but are ignored in the zero-Froude-number potential ϕ0(1). Results of calculations by C. Y. Chen for the Wigley hull show that the waves in the potential ϕ(1) have a remarkable effect upon the wave resistance, in particular causing a large phase shift of the wave-resistance curve toward higher values of the Froude number. As a result, the first-order slender-ship wave-resistance approximation in significantly better agreement with experimental data than the low-Froude-number approximation rlF(1) and the approximations r(0) and rM.


2019 ◽  
Vol 3 (3) ◽  
Author(s):  
Jin-Yi Li ◽  
Xiang-Wei Lv ◽  
Guo-Qiang Zhong ◽  
Hong-Hong Ke

Abstract Background Micro-reentry tachycardia usually emerges in scar tissues related to post-atrial fibrillation ablation and cardiomyopathy. It is difficult to identify the micro-reentry circuit accurately by conventional mapping method. Case summary A 74-year-old man presented with paroxysmal atrial tachycardia (AT) presenting as palpitations. He was evaluated by an electrophysiological examination using a high-density CARTO mapping system. The mapping results showed the AT with a cycle length of 184 ms was focused on his right atrial fossa ovalis (FO). In this small area, the high-density mapping demonstrated a significant micro-reentrant tachycardia. Radiofrequency ablation at the centre of the micro-reentrant circuit successfully terminated the AT. No recurrences were observed during a 12-month follow-up. Discussion This case demonstrated a micro-reentrant AT originates from the FO without cardiomyopathy or previous ablation with specific loops. This is an unusual location for AT though and can cause difficulty for operators if it terminates or is non-sustained. High-density mapping using a PentaRay catheter can effectively characterize micro-reentrant circuits and determine the real target for ablation therapy.


2012 ◽  
Vol 490-495 ◽  
pp. 3733-3737
Author(s):  
Shu Hong Jing ◽  
Shou Zhi Pu ◽  
Shi Qiang Cui

A new photochromic diarylethene compound, 1-(2,4-dimethoxy-5-pyrimidine)-2-[2-methyl-5-(3-pyridine)-3-thienyl]perfluorocyclopentene(1a), was synthesized, and its photochromic reactivity, fluorescent and electrochemical property were also investigated. Diarylethene 1a changed the color from colorless to pink upon irradiation with UV light, in which absorption maxima were observed at 520 and 519 nm in hexane and PMMA film, respectively. The the photochromic reaction kinetics indicated that the cyclization processes of 1 belong to the zeroth order reaction and the cycloreversion process belong to the first order reaction. This new photochromic system also exhibited remarkable fluorescence switching in hexane solution and this new photochromic system also exhibited remarkable optical storage character.


2019 ◽  
Vol 2019 ◽  
pp. 1-11 ◽  
Author(s):  
Cheng-Biao Fu ◽  
Hei-Gang Xiong ◽  
An-Hong Tian

Discussion on the application of fractional derivative algorithm in monitoring organic matter content in field soil is scarce. This study is aimed at improving the accuracy of soil organic matter (SOM) content estimation in arid region, and the undesirable model precision caused by the missing information associated with the larger discrepancy between conventional integer-order, i.e., first order and second order, derivative, and raw spectral data. We utilized fractional derivative (of zeroth order to second order in 0.2-order interval) processing on the field spectral reflectance (R) of the salinized soil sample from Fukang, Xinjiang, and its square root-transformed (R), log-transformed (lgR), inverse-transformed (1/R), and inverse log-transformed (1/lgR) values. The correlation coefficient of each fractional derivative of transformed value with SOM content was calculated. The simulation showed the derivative reflectance value approximates zero. When increasing from zeroth order to first order, the derivative curve gradually aligns to the first-order curve, and the destination alignment was also seen while increasing from first order to second order. The significance test of 0.05 showed initial increase and later decay of bands in the five spectral transformations as the order increases. For specific bands, the derivative algorithm clearly justifies the correlation between soil spectra and organic matter content, and all of the absolute highest correlation coefficient values were obtained at fractional orders. When compared with integer-order derivative, fractional derivative is significantly better in improving correlation, showing overall superiority. The result supports the application of fractional derivative in the hyperspectral remote monitor of SOM in arid zone, which may in turn realize the timely and accurate SOM monitor in arid zone, and provides the basis for ecological restoration.


2003 ◽  
Vol 12 (06) ◽  
pp. 767-779 ◽  
Author(s):  
Jörg Sawollek

It is an open question whether there are Vassiliev invariants that can distinguish an oriented knot from its inverse, i.e., the knot with the opposite orientation. In this article, an example is given for a first order Vassiliev invariant that takes different values on a virtual knot and its inverse. The Vassiliev invariant is derived from the Conway polynomial for virtual knots. Furthermore, it is shown that the zeroth order Vassiliev invariant coming from the Conway polynomial cannot distinguish a virtual link from its inverse and that it vanishes for virtual knots.


Sign in / Sign up

Export Citation Format

Share Document