Inversion of the power spectrum from magnetic anomalies

Geophysics ◽  
1994 ◽  
Vol 59 (3) ◽  
pp. 391-401 ◽  
Author(s):  
Juan Garca‐Abdeslem ◽  
Gordon E. Ness

We develop methods for the modeling and inversion of the power‐density spectrum from magnetic anomaly data assuming that the crustal magnetic field is caused by an ensemble of vertical‐sided and uniformly magnetized prisms. The solution of the forward problem is achieved in the wavenumber domain, where a synthetic spectrum is given by the product of the mathematical expectations of single‐valued functions that describe depth, thickness, and horizontal dimensions of prisms in the ensemble. We use Gaussian and uniform distributions to describe the ensemble and provide a variety of functions from which different statistical models can be obtained. The solution of the inverse problem is achieved iteratively, starting from an initial set of model parameters. It is based on the ridge‐regression algorithm, and its usefulness is assessed in a number of examples with numeric, synthetic and real data spectra. The methods are first tested on the spectrum obtained from a simple artificial magnetic anomaly and on the artificial spectrum caused by an ensemble of source bodies and are found to be capable of recovering the source parameters. Next, the methods are applied to marine magnetic data from a survey offshore of the Yucatán Peninsula, Mexico. The results of this last application are consistent with the crustal structure observed at Chicxulub hole.

Geophysics ◽  
2007 ◽  
Vol 72 (3) ◽  
pp. L21-L30 ◽  
Author(s):  
Soraya Lozada Tuma ◽  
Carlos Alberto Mendonça

We present a three-step magnetic inversion procedure in which invariant quantities with respect to source parameters are inverted sequentially to give (1) shape cross section, (2) magnetization intensity, and (3) magnetization direction for a 2D (elongated) magnetic source. The quantity first inverted (called here the shape function) is obtained from the ratio of the gradient intensity of the total-field anomaly to the intensity of the anomalous vector field. For homogenous sources, the shape function is invariant with source magnetization and allows reconstruction of the source geometry by attributing an arbitrary magnetization to trial solutions. Once determined, the source shape is fixed and magnetization intensity is estimated by fitting the total gradient of the total-field anomaly (equivalent to the amplitude of the analytic signal of magnetic anomaly). Finally, the source shape and magnetization intensity are fixed and the magnetization direction is determined by fitting the magnetic anomaly. As suggested by numerical modeling and real data application, stepped inversion allows checking whether causative sources are homogeneous. This is possible because the shape function from inhomogeneous sources can be fitted by homogeneous models, but a model obtained in this way fits neither the total gradient of the magnetic anomaly nor the magnetic anomaly itself. Such a criterion seems effective in recognizing strongly inhomogeneous sources. Stepped inversion is tested with numerical experiments, and is used to model a magnetic anomaly from intrusive basic rocks from the Paraná Basin, Brazil.


Geophysics ◽  
2010 ◽  
Vol 75 (4) ◽  
pp. L79-L90 ◽  
Author(s):  
Daniela Gerovska ◽  
Marcos J. Araúzo-Bravo ◽  
Kathryn Whaler ◽  
Petar Stavrev ◽  
Alan Reid

We present an automatic procedure for interpretation of magnetic or gravity gridded anomalies based on the finite-difference similarity transform (FDST). It is called MaGSoundFDST (magnetic and gravity sounding based on the finite-difference similarity transform) and uses a “focusing” principle in contrast to deriving multiple clusters of many solutions as in the widely used Euler deconvolution method. The source parameters are characterized by isolated solutions, and the interpreter obtains parallel images showing the horizontal position, depth, and structural index [Formula: see text] value. The underlying principle is that the FDST of a potential field anomaly becomes zero or linear at all observation points when the central point of similarity (CPS) of the transform coincides with a source field’s singular point and a correct [Formula: see text] value is used. The procedure involves calculating a 3D function that evaluates the linearity of the FDST for a series of [Formula: see text] values, using a moving window and sounding the subsurface along a verticalline under each window center. We then combine the 3D results for different [Formula: see text] values into a single map whose minima determine the horizontal position of the sources. The [Formula: see text] value and the CPS depth associated with each minimum determine the [Formula: see text] value and depth of the corresponding source. Only one estimate characterizes a simple source, which is a major advantage over other window-based procedures. MaGSoundFDST uses only the measured anomalous field and its upward continuation, thus avoiding the direct use of field derivatives. It is independent of the magnetization-vector direction in the magnetic data case. The procedure accounts for a linear background of local gravity or magnetic anomalies and has been applied effectively to several cases of synthetic and real data. MaGSoundFDST shares common features with the magnetic and gravity sounding based on the differential similarity transform (MaGSoundDST) but is more stable in estimating depth and structural index in the presence of random noise.


Author(s):  
M M Serim ◽  
Ö C Özüdoğru ◽  
Ç K Dönmez ◽  
Ş Şahiner ◽  
D Serim ◽  
...  

Abstract We investigate timing and spectral characteristics of the transient X-ray pulsar 2S 1417−624 during its 2018 outburst with NICER follow up observations. We describe the spectra with high-energy cut-off and partial covering fraction absorption (PCFA) model and present flux-dependent spectral changes of the source during the 2018 outburst. Utilizing the correlation-mode switching of the spectral model parameters, we confirm the previously reported sub-critical to critical regime transitions and we argue that secondary transition from the gas-dominated to the radiation pressure-dominated disc do not lead to significant spectral changes below 12 keV. Using the existing accretion theories, we model the spin frequency evolution of 2S 1417−624 and investigate the noise processes of a transient X-ray pulsar for the first time using both polynomial and luminosity-dependent models for the spin frequency evolution. For the first model, the power density spectrum of the torque fluctuations indicate that the source exhibits red noise component (Γ ∼ −2) within the timescales of outburst duration which is typical for disc-fed systems. On the other hand, the noise spectrum tends to be white on longer timescales with high timing noise level that indicates an ongoing accretion process in between outburst episodes. For the second model, most of the red noise component is eliminated and the noise spectrum is found to be consistent with a white noise structure observed in wind-fed systems.


2017 ◽  
Author(s):  
Khalid S. Essa ◽  
Mahmoud El-Hussein

Abstract. A robust Particle Swarm Optimization (PSO) investigation for magnetic data by a 2D dipping dike has been presented. The interpretive model parameters are: the amplitude coefficient (K), the depth to the top of the dipping dike (z), exact origin of the dipping dike (x0), and the width of dipping dike (w). The inversion procedure is actualized to gauge the parameters of a 2D dipping dike structures where it has been confirmed first on synthetic models without and with different level of random noise. The results of the inversion demonstrate that the parameters derived from the inversion concur well with the true ones. The root mean square (RMS) is figured by the strategy which is considered as the misfit between the measured and computed anomalies. The technique has been warily and effectively applied to real data examples from China and UK with the presence of ore bodies. The present technique can be applicable for mineral exploration and ore bodies of dike-like structure embedded in the shallow and deeper subsurface.


2019 ◽  
Vol XVI (2) ◽  
pp. 1-11
Author(s):  
Farrukh Jamal ◽  
Hesham Mohammed Reyad ◽  
Soha Othman Ahmed ◽  
Muhammad Akbar Ali Shah ◽  
Emrah Altun

A new three-parameter continuous model called the exponentiated half-logistic Lomax distribution is introduced in this paper. Basic mathematical properties for the proposed model were investigated which include raw and incomplete moments, skewness, kurtosis, generating functions, Rényi entropy, Lorenz, Bonferroni and Zenga curves, probability weighted moment, stress strength model, order statistics, and record statistics. The model parameters were estimated by using the maximum likelihood criterion and the behaviours of these estimates were examined by conducting a simulation study. The applicability of the new model is illustrated by applying it on a real data set.


Galaxies ◽  
2021 ◽  
Vol 9 (2) ◽  
pp. 25
Author(s):  
Debjit Chatterjee ◽  
Arghajit Jana ◽  
Kaushik Chatterjee ◽  
Riya Bhowmick ◽  
Sujoy Kumar Nath ◽  
...  

We study the properties of the faint X-ray activity of Galactic transient black hole candidate XTE J1908+094 during its 2019 outburst. Here, we report the results of detailed spectral and temporal analysis during this outburst using observations from Nuclear Spectroscopic Telescope Array (NuSTAR). We have not observed any quasi-periodic-oscillations (QPOs) in the power density spectrum (PDS). The spectral study suggests that the source remained in the softer (more precisely, in the soft–intermediate) spectral state during this short period of X-ray activity. We notice a faint but broad Fe Kα emission line at around 6.5 keV. We also estimate the probable mass of the black hole to be 6.5−0.7+0.5M⊙, with 90% confidence.


Mathematics ◽  
2021 ◽  
Vol 9 (16) ◽  
pp. 1850
Author(s):  
Rashad A. R. Bantan ◽  
Farrukh Jamal ◽  
Christophe Chesneau ◽  
Mohammed Elgarhy

Unit distributions are commonly used in probability and statistics to describe useful quantities with values between 0 and 1, such as proportions, probabilities, and percentages. Some unit distributions are defined in a natural analytical manner, and the others are derived through the transformation of an existing distribution defined in a greater domain. In this article, we introduce the unit gamma/Gompertz distribution, founded on the inverse-exponential scheme and the gamma/Gompertz distribution. The gamma/Gompertz distribution is known to be a very flexible three-parameter lifetime distribution, and we aim to transpose this flexibility to the unit interval. First, we check this aspect with the analytical behavior of the primary functions. It is shown that the probability density function can be increasing, decreasing, “increasing-decreasing” and “decreasing-increasing”, with pliant asymmetric properties. On the other hand, the hazard rate function has monotonically increasing, decreasing, or constant shapes. We complete the theoretical part with some propositions on stochastic ordering, moments, quantiles, and the reliability coefficient. Practically, to estimate the model parameters from unit data, the maximum likelihood method is used. We present some simulation results to evaluate this method. Two applications using real data sets, one on trade shares and the other on flood levels, demonstrate the importance of the new model when compared to other unit models.


2021 ◽  
Vol 11 (15) ◽  
pp. 6998
Author(s):  
Qiuying Li ◽  
Hoang Pham

Many NHPP software reliability growth models (SRGMs) have been proposed to assess software reliability during the past 40 years, but most of them have focused on modeling the fault detection process (FDP) in two ways: one is to ignore the fault correction process (FCP), i.e., faults are assumed to be instantaneously removed after the failure caused by the faults is detected. However, in real software development, it is not always reliable as fault removal usually needs time, i.e., the faults causing failures cannot always be removed at once and the detected failures will become more and more difficult to correct as testing progresses. Another way to model the fault correction process is to consider the time delay between the fault detection and fault correction. The time delay has been assumed to be constant and function dependent on time or random variables following some kind of distribution. In this paper, some useful approaches to the modeling of dual fault detection and correction processes are discussed. The dependencies between fault amounts of dual processes are considered instead of fault correction time-delay. A model aiming to integrate fault-detection processes and fault-correction processes, along with the incorporation of a fault introduction rate and testing coverage rate into the software reliability evaluation is proposed. The model parameters are estimated using the Least Squares Estimation (LSE) method. The descriptive and predictive performance of this proposed model and other existing NHPP SRGMs are investigated by using three real data-sets based on four criteria, respectively. The results show that the new model can be significantly effective in yielding better reliability estimation and prediction.


2013 ◽  
Vol 2013 ◽  
pp. 1-9 ◽  
Author(s):  
Ali Kargarnejad ◽  
Mohsen Taherbaneh ◽  
Amir Hosein Kashefi

Tracking maximum power point of a solar panel is of interest in most of photovoltaic applications. Solar panel modeling is also very interesting exclusively based on manufacturers data. Knowing that the manufacturers generally give the electrical specifications of their products at one operating condition, there are so many cases in which the specifications in other conditions are of interest. In this research, a comprehensive one-diode model for a solar panel with maximum obtainable accuracy is fully developed only based on datasheet values. The model parameters dependencies on environmental conditions are taken into consideration as much as possible. Comparison between real data and simulations results shows that the proposed model has maximum obtainable accuracy. Then a new fuzzy-based controller to track the maximum power point of the solar panel is also proposed which has better response from speed, accuracy and stability point of view respect to the previous common developed one.


Geophysics ◽  
2016 ◽  
Vol 81 (4) ◽  
pp. U25-U38 ◽  
Author(s):  
Nuno V. da Silva ◽  
Andrew Ratcliffe ◽  
Vetle Vinje ◽  
Graham Conroy

Parameterization lies at the center of anisotropic full-waveform inversion (FWI) with multiparameter updates. This is because FWI aims to update the long and short wavelengths of the perturbations. Thus, it is important that the parameterization accommodates this. Recently, there has been an intensive effort to determine the optimal parameterization, centering the fundamental discussion mainly on the analysis of radiation patterns for each one of these parameterizations, and aiming to determine which is best suited for multiparameter inversion. We have developed a new parameterization in the scope of FWI, based on the concept of kinematically equivalent media, as originally proposed in other areas of seismic data analysis. Our analysis is also based on radiation patterns, as well as the relation between the perturbation of this set of parameters and perturbation in traveltime. The radiation pattern reveals that this parameterization combines some of the characteristics of parameterizations with one velocity and two Thomsen’s parameters and parameterizations using two velocities and one Thomsen’s parameter. The study of perturbation of traveltime with perturbation of model parameters shows that the new parameterization is less ambiguous when relating these quantities in comparison with other more commonly used parameterizations. We have concluded that our new parameterization is well-suited for inverting diving waves, which are of paramount importance to carry out practical FWI successfully. We have demonstrated that the new parameterization produces good inversion results with synthetic and real data examples. In the latter case of the real data example from the Central North Sea, the inverted models show good agreement with the geologic structures, leading to an improvement of the seismic image and flatness of the common image gathers.


Sign in / Sign up

Export Citation Format

Share Document