Source author record

Michael E. Glinsky

Michael E. Glinsky appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

6works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

6 published item(s)

preprint2022arXiv

Physics-assisted Generative Adversarial Network for X-Ray Tomography

X-ray tomography is capable of imaging the interior of objects in three dimensions non-invasively, with applications in biomedical imaging, materials science, electronic inspection, and other fields. The reconstruction process can be an ill-conditioned inverse problem, requiring regularization to obtain satisfactory results. Recently, deep learning has been adopted for tomographic reconstruction. Unlike iterative algorithms which require a distribution that is known a priori, deep reconstruction networks can learn a prior distribution through sampling the training distributions. In this work, we develop a Physics-assisted Generative Adversarial Network (PGAN), a two-step algorithm for tomographic reconstruction. In contrast to previous efforts, our PGAN utilizes maximum-likelihood estimates derived from the measurements to regularize the reconstruction with both known physics and the learned prior. Compared with methods with less physics assisting in training, PGAN can reduce the photon requirement with limited projection angles to achieve a given error rate. The advantages of using a physics-assisted learned prior in X-ray tomography may further enable low-photon nanoscale imaging.

preprint2016arXiv

A novel workflow for seismic net pay estimation with uncertainty

This paper presents a novel workflow for seismic net pay estimation with uncertainty. It is demonstrated on the Cassra/Iris Field. The theory for the stochastic wavelet derivation (which estimates the seismic noise level along with the wavelet, time-to-depth mapping, and their uncertainties), the stochastic sparse spike inversion, and the net pay estimation (using secant areas) along with its uncertainty; will be outlined. This includes benchmarking of this methodology on a synthetic model. A critical part of this process is the calibration of the secant areas. This is done in a two step process. First, a preliminary calibration is done with the stochastic reflection response modeling using rock physics relationships derived from the well logs. Second, a refinement is made to the calibration to account for the encountered net pay at the wells. Finally, a variogram structure is estimated from the extracted secant area map, then used to build in the lateral correlation to the ensemble of net pay maps while matching the well results to within the nugget of the variogram. These net pay maps are then integrated, over the area of full saturation gas, to give the GIIP distribution (Gaussian distributions for the porosity, gas expansion factor, and gas saturation for the sand end member are assumed and incorporated in the estimate of GIIP). The method is demonstrated on the Iris (UP5 turbidite) interval. The net pay is corrected for reduction in the amplitudes over part of the area due to shallow gas. The sensitivity of the GIIP to the independent stochastic variables is estimated (determining the value of information) so that business decisions can be made that maximize the value of the field.

preprint2014arXiv

A new perspective on renormalization: the scattering transformation

The scattering transformation developed by Mallat is put into the perspective of field theory. It is shown to be a simultaneous transformation of the field and the "time" parameter explicitly used in the definition of path integrals central to the Lagrangian approach and in the definition of the "time" ordered products used in the Hamiltonian (or canonical) approach. This transformation preserves the form of the S-matrix as "time" ordered products. The transformed "time" coordinate is the inverse "time" scale. This is traditionally the UV cutoff or renormalization parameter in standard approaches to renormalization. The critical calculation will be the determination of the form of the effective action in this transformed coordinate system. This action will now be expressed as an integral of a Lagrangian density that is a function of the renormalization parameter or transformed "time". Other symmetries of the action can be explicitly built into the transformation. It will be demonstrated on a simple 1D $ϕ^4$ field theory. This non-perturbative approach has great potential in possibly being used to renormalize quantum gravity and obtaining expressions for the strongly coupled limit of QCD.

preprint2014arXiv

Focal mechanism estimation by classification

A classification technique for identifying focal mechanism type and fault plane orientation based on the polarity of P-wave "first motion" data is derived. A support vector machine is used to classify the polarity data in the space of spherical harmonic functions. The classification is non-parametric in the sense that there is no requirement to make a priori assumptions source mechanism. A metric of similarity potentially able to distinguish shear versus tensile dislocation without requiring estimation of the fault plane orientation is a natural consequence of this procedure. Going further, correlation functions between template source mechanism is derived, gives an estimate of fault plane orientation assuming a particular source mechanism.

preprint2014arXiv

Rock physics and geophysics for unconventional resource, multi-component seismic, quantitative interpretation

An extension of a previously developed rock physics model is made that quantifies the relationship between the ductile fraction of a brittle/ductile binary mixture and the isotropic seismic reflection response. By making a weak scattering (Born) approximation and plane wave (eikonal) approximation, with a subsequent ordering according to the angle of incidence, singular value decomposition analysis are done to understand the stack weightings, number of stacks, and the type of stacks that will optimally estimate the two fundamental rock physics parameters. Through this angle ordering, it is found that effective wavelets can be used for the stacks up to second order. Finally, it is concluded that the full PP stack and the "full" PS stack are the two optimal stacks needed to estimate the two rock physics parameters. They dominate over both the second order AVO "gradient" stack and the higher order (4th order) PP stack (even at large angles of incidence). Using this result and model based Bayesian inversion, the detectability of the ductile fraction (shown by others to be the important quantity for the geomechanical response of unconventional reservoir fracking) is demonstrated on a model characteristic of the Marcellus shale play.

preprint2011arXiv

Turbidity current flow over an obstacle and phases of sediment wave generation

We study the flow of particle-laden turbidity currents down a slope and over an obstacle. A high-resolution 2D computer simulation model is used, based on the Navier-Stokes equations. It includes poly-disperse particle grain sizes in the current and substrate. Particular attention is paid to the erosion and deposition of the substrate particles, including application of an active layer model. Multiple flows are modeled from a lock release that can show the development of sediment waves (SW). These are stream-wise waves that are triggered by the increasing slope on the downstream side of the obstacle. The initial obstacle is completely erased by the resuspension after a few flows leading to self consistent and self generated SW that are weakly dependant on the initial obstacle. The growth of these waves is directly related to the turbidity current being self sustaining, that is, the net erosion is more than the net deposition. Four system parameters are found to influence the SW growth: (1) slope, (2) current lock height, (3) grain lock concentration, and (4) particle diameters. Three phases are discovered for the system: (1) "no SW", (2) "SW buildup", and (3) "SW growth". The second phase consists of a soliton-like SW structure with a preserved shape. The phase diagram of the system is defined by isolating regions divided by critical slope angles as functions of current lock height, grain lock concentration, and particle diameters.