Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Data virtualization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 937 records · Page 52

Radiographic Imaging and Tomography

Radiographic imaging and tomography (RadIT) come in many flavors such as X-ray imaging and tomography (IT), proton IT, neutron IT, muon IT, neutrino IT, etc. We identify five RadIT themes: Physics, Sources, Detectors, Methods, and Data Science, which are integral parts of image interpretation and three-dimensional (3D) tomographic reconstruction. Traditionally, RadIT have been driven by medicine, non-destructive testing, material sciences, and security applications. The latest thrusts of growth come from automation, machine vision, additive manufacturing and virtual reality (‘metaverse’). The five RadIT themes parallel their counterparts in optical IT. Synergies among different forms of RadIT and with optical IT motivate further advances towards multi-modal IT and quantum IT.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

EPCAPE-PT-LANL Measurements: Photoacoustic Extinctometer

Coastal cities offer a unique environment for studying aerosol-cloud interactions and the effects of urban emissions on cloud properties. As part of the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE), the Partitioning Thrust by Los Alamos National Laboratory (EPCAPE-PT-LANL) was conducted. Our campaign focused on measuring the optical and chemical properties of aerosols and their interactions within marine stratocumulus clouds in La Jolla, California. EPCAPE-PT-LANL enhances the primary goals of EPCAPE through innovative observations of vapor-phase transitions between aerosols and cloud droplets, the impact of black carbon on aerosol-cloud dynamics, and the effects of cloud processing on aerosol optical properties. Instrument: Photoacoustic Extinctiometer at 870 nm (Droplet Measurement Technology) Header: - AbsorptionNoise_870nm[1/Mm]: Noise level associated with the absorption measurement at 870 nm, expressed in inverse megameters (Mm⁻¹). - Absorption_870nm[1/Mm]: Absorption coefficient of aerosols measured at 870 nm, expressed in inverse megameters (Mm⁻¹). - Extinction_870nm[1/Mm]: Extinction coefficient of aerosols measured at 870 nm, expressing the sum of scattering and absorption by aerosols, in inverse megameters (Mm⁻¹). - Scattering_870nm[1/Mm]: Scattering coefficient of aerosols measured at 870 nm, expressed in inverse megameters (Mm⁻¹). - BackgroundAbsorption_870nm[1/Mm]: Background absorption measurement at 870 nm, used to correct the primary absorption data, expressed in inverse megameters (Mm⁻¹). - BackgroundScattering_870nm[1/Mm]: Background scattering measurement at 870 nm, used to correct the primary scattering data, expressed in inverse megameters (Mm⁻¹). - Temperature[C]: The ambient temperature at the time of the measurement, expressed in degrees Celsius. - RelativeHumidity[%]: The relative humidity at the time of the measurement, expressed as a percentage. - CVI_Flag[bool]: A boolean flag indicating whether the Counterflow Virtual Impactor (CVI) was active (true) or inactive (false) during the measurement

54 ENVIRONMENTAL SCIENCES↗

A comparison of cosmological models using strong gravitational lensing galaxies

Strongly gravitationally lensed quasar-galaxy systems allow us to compare competing cosmologies as long as one can be reasonably sure of the mass distribution within the intervening lens. In this paper, we assemble a catalog of 69 such systems from the Sloan Lens ACS and Lens Structure and Dynamics surveys suitable for this analysis, and carry out a one-on-one comparison between the standard model, ΛCDM, and the R{sub h}=ct universe, which has thus far been favored by the application of model selection tools to other kinds of data. We find that both models account for the lens observations quite well, though the precision of these measurements does not appear to be good enough to favor one model over the other. Part of the reason is the so-called bulge-halo conspiracy that, on average, results in a baryonic velocity dispersion within a fraction of the optical effective radius virtually identical to that expected for the whole luminous-dark matter distribution modeled as a singular isothermal ellipsoid, though with some scatter among individual sources. Future work can greatly improve the precision of these measurements by focusing on lensing systems with galaxies as close as possible to the background sources. Given the limitations of doing precision cosmological testing using the current sample, we also carry out Monte Carlo simulations based on the current lens measurements to estimate how large the source catalog would have to be in order to rule out either model at a ∼99.7% confidence level. We find that if the real cosmology is ΛCDM, a sample of ∼200 strong gravitational lenses would be sufficient to rule out R{sub h}=ct at this level of accuracy, while ∼300 strong gravitational lenses would be required to rule out ΛCDM if the real universe were instead R{sub h}=ct. The difference in required sample size reflects the greater number of free parameters available to fit the data with ΛCDM. We point out that, should the R{sub h}=ct universe eventually emerge as the correct cosmology, its lack of any free parameters for this kind of work will provide a remarkably powerful probe of the mass structure in lensing galaxies, and a means of better understanding the origin of the bulge-halo conspiracy.

79 ASTRONOMY AND ASTROPHYSICS↗

X-ray fluorescence microscopy methods for biological tissues

Abstract Synchrotron-based X-ray fluorescence microscopy is a flexible tool for identifying the distribution of trace elements in biological specimens across a broad range of sample sizes. The technique is not particularly limited by sample type and can be performed on ancient fossils, fixed or fresh tissue specimens, and in some cases even live tissue and live cells can be studied. The technique can also be expanded to provide chemical specificity to elemental maps, either at individual points of interest in a map or across a large field of view. While virtually any sample type can be characterized with X-ray fluorescence microscopy, common biological sample preparation methods (often borrowed from other fields, such as histology) can lead to unforeseen pitfalls, resulting in altered element distributions and concentrations. A general overview of sample preparation and data-acquisition methods for X-ray fluorescence microscopy is presented, along with outlining the general approach for applying this technique to a new field of investigation for prospective new users. Considerations for improving data acquisition and quality are reviewed as well as the effects of sample preparation, with a particular focus on soft tissues. The effects of common sample pretreatment steps as well as the underlying factors that govern which, and to what extent, specific elements are likely to be altered are reviewed along with common artifacts observed in X-ray fluorescence microscopy data.

Pushie, M. Jake (ORCID:0000000174945427)↗

A Grid Intensity Map-Based Source Model for Monte Carlo Modeling of IMRT Field Irradiation - Paper 5

An accurate source model is important for TPS (Treatment Planning System) dose calculation and verification, especially for IMRT (Intensity Modulated Radiotherapy) irregular and inhomogeneous fields. Monte Carlo simulation is the most accurate way to calculate dose distributions in the heterogeneous media like human bodies. A full radiotherapy modeling generally includes Linac modeling and phantom dose calculation. Full Linac modeling is very time-consuming and requires highly detailed knowledge of the accelerator composition and structure, however, which is very difficult to obtain for commercial interests. Often all relevant parameters are not known with sufficient accuracy, and thus a time-consuming tuning of e.g. the electron beam spot size and energy must be performed. A multiple source model (MSM) regards particles coming from different parts of an accelerator as if they are coming from different virtual sources. MSM uses a set of virtual sources representing the treatment beam (e.g. primary photons, head scatter, electron contamination) by a set of analytic functions or histograms with parameters derived from dosimetric measurements or a full MC simulation. Tillikainen et al built a three sub-source model by matching with measurement data using an automatic process to optimize the free model parameters for a Varian 2100 C/D accelerator. Lin et al have built a convenient five-photon-sub-source model for Varian 2300C 6 MV-X beam based on Monte Carlo simulated Linac for ARTS (Accurate Radiotherapy System). The particle generation of MSM is very fast compared to full Linac simulation, and its sampling efficiency is not limited by the size of the phase space file (PSF). However, generally MSM is generated based on regular field PSF and needs an irregular component to form the beam shape meanwhile apt to bringing error Linac modeling and MSM producing from PSF analysis actually both are multi-freedom optimization process by adjusting several parameters (e.g. electron energy and angular divergence). Linac modeling adjusts its parameters from the beginning of electron hitting target through all components to fit the simulation to the measurement. MSM tunes its PSF-derived parameters by optimizing to minimize the deviations between the simulation and the measurement. Since the full Linac modeling and the MSM model both need to fine tune their parameters to match the simulation to the measurement, and sometimes maybe offsets the real accelerator for its multi-parameter optimization process. So, why not fine tune directly all source model parameters just behind all Linac components to fit the simulation to the measurement? A grid intensity-based surface source model was suggested here, which assumed the particle emitting from the bottom surface of Linac last beam trimmer MLC (Multi Leaf Collimator). The particle emitting position is sampled on the grid intensity map optimized by TPS, and the particle weight is adjusted by its original grid intensity. This source model not only takes advantage of the Monte Carlo method in the heterogeneous media, but also avoids effectively the transportation speed shortcoming of Monte Carlo simulation through Linac components. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

NDNSEG (Named Data Network Scalable Environment Generator) [SWR-21-91]

NDNSEG (Named Data Network Scalable Environment Generator) automates the deployment of NDN environments in NREL's Cyber-Energy Emulation Platform (CEEP). Specifying the number of clients, servers, gateways, and pv_inverters and running NDNSEG will automate the generation of required files to deploy the number of specified virtual machines, the underlying networking infrastructure, and configuration for GitLab CI/CD.

Peterson, Jordan↗

Influences of lidar scanning parameters on wind turbine wake retrievals in complex terrain

Abstract. Scanning lidars enable the collection of spatially distributed measurements of turbine wakes and the estimation of wake properties such as magnitude, extent, and trajectory. Lidar-based characterizations, however, may be subject to distortions due to the observational system. Distortions can arise from the resolution of the measurement points across the wake, the projection of the winds onto the beam, averaging along the beam probe volume, and intervening evolution of the flow over the scan duration. Using a large-eddy simulation and simulated measurements with a virtual lidar model, we assess how scanning lidar systems may influence the properties of the retrieved wake using a case study from the Perdigão campaign. We consider three lidars performing range-height indicator sweeps in complex terrain, based on the deployments of lidars from the Danish Technical University (DTU) and German Aerospace Center (DLR) at the Perdigão site. The unwaked flow, measured by the DTU lidar, is well-captured by the lidar, even without combining data into a multi-lidar retrieval. The two DLR lidars measure a waked transect from different downwind vantage points. In the region of the wake, the observation system reacts to the smaller spatial and temporal variations of the winds, allowing more significant observation distortions to arise. While the measurements largely capture the wake structure and trajectory over its 4–5 D extent, limited spatial resolution of measurement points and volume averaging lead to a quicker loss of the two lobes in the near wake, smearing of the vertical bounds of the wake (< 30 m), wake center displacements up to 10 m, and dampening of the maximum velocity deficit by up to a third. The virtual lidar tool, coupled with simulations, provides a means for assessing measurement capabilities in advance of measurement campaigns.

17 WIND ENERGY↗

HDG-1 Experiment Irradiation Monitoring Data Qualification Final Report

SUMMARY The U.S. Department of Energy (DOE) Advanced Reactor Technologies (ART) Graphite Research and Development (GRD) Program is conducting a series of six experiments to quantify the effects of irradiation on nuclear-grade graphite. This report documents the qualification of irradiation monitoring data for the fifth experiment, High Dose Graphite-1 (HDG-1). Qualified monitoring data are required by the ART program to support the design and licensing of the first high-temperature reactor (HTR) nuclear plant. Data are classified as Qualified if they meet the usage requirements described in the experiment planning and quality assurance (QA) documents, Failed if they do not meet those requirements and provide no usable information, or Trend if they do not fully meet all requirements but still provide useful information subject to an assessment of how any deficiencies may affect a particular use of the data. HDG-1 irradiation began with Advanced Test Reactor (ATR) Cycle 168B on August 24, 2020, and concluded after Cycle 173C on January 27, 2025. The HDG-1 capsule was removed from the reactor core twice—during core internal change (CIC) Cycle 170A and powered axial locator mechanism (PALM) Cycle 172A—to prevent overheating of the graphite specimens during high-power PALM cycles. The capsule was therefore irradiated during a total of seven normal ATR cycles: 168B, 169A, 171A, 171B, 173A, 173B, and 173C. Irradiation monitoring data evaluated in this report include thermocouple (TC) temperature, gas flow rate, gas moisture, gas pressure, specimen load, and graphite stack displacement. Temperature. A total of 14,508,065 TC temperature records were captured. Of these, 13,901,785 (95.8%) are Qualified and 606,280 (4.2%) are Failed. The principal source of failed temperature data was the instrument failure of TC-9 (Zone 2) on June 24, 2024, and TC-10 (Zone 1) on July 5, 2024, near the end of Cycle 173A, which resulted in 595,554 Failed readings. An additional 379 missing values and 10,347 slightly negative values from TC-13 during ATR outages are also Failed. Neither TC-9 nor TC-10 was used as a temperature-control TC, and their failures did not compromise capsule condition monitoring. Correlation analysis of all 13 TCs found no evidence of virtual junction formation. Control chart analysis revealed clear downward drift of approximately 80°C for TC-6 (Zone 3) relative to other stable TCs, and possible downward drift of approximately 60°C for TC-13 relative to the Zone 5 control TC (TC-1), though TC-13 remained consistent with the Zone 2 control TC (TC-12). Gas flow. A total of 20,088,090 gas flow rate records were captured. Of these, 19,941,463 (99.3%) are Qualified and 146,627 (0.7%) are Failed due to missing values. All argon, helium, and total gas flow data were within expected ranges throughout the irradiation. Gas moisture. A total of 1,116,005 outlet gas moisture values were captured. Of these, 1,101,421 (98.7%) are Qualified and 14,584 (1.3%) are Failed, comprising 14,556 out-of-range values and 28 missing values. The out-of-range moisture values exceeded 22,000 ppmv for approximately 1 week at the beginning of Cycle 173A, when accumulated moisture evaporated after the capsule was retrieved from water storage during PALM Cycle 172A and reinserted into the east flux trap. Moisture levels returned to below 25 ppmv for the remaining three cycles, and the transient high-moisture event did not affect the integrity of specimen irradiation. Gas pressure. A total of 7,812,035 gas pressure values were captured. Of these, 6,642,048 (85.0%) are Qualified and 1,169,987 (15.0%) outlet pressure values are Failed, comprising 718,537 zero outlet pressure values due to sensor failure from Cycle 168B through Cycle 171B, 54,550 missing values, and 396,900 too-low outlet pressure values, ranging from 1.1 to 1.6 psia after sensor replacement during Cycle 173A. Load. A total of 6,696,030 load values were captured. Of these, 6,694,580 (99.98%) are Qualified and 1,450 (0.02%) are Failed due to missing values. Applied loads to the six specimen stacks were stable throughout the irradiation. Stack displacement. A total of 6,696,030 displacement values were captured. Of these, 5,713,297 (85.32%) are Qualified and 3,781 (0.06%) are Failed due to missing values. Stack displacement increased consistently throughout the irradiation, reaching approximately 3.08 in. for Channels 5 and 6 by the end of irradiation. 978,952 (14.62%) substantially elevated displacements observed for Channel 6 beginning in Cycle 171A and for Channel 5 beginning in Cycle 173A are assigned Trend status. Raising pressure. A total of 1,115,999 raising pressure values were captured. Of these, 1,115,430 (99.95%) are Qualified and 569 (0.05%) are Failed due to missing values. Ram pressure. A total of 6,696,030 ram pressure values were captured. Of these, 6,692,249 (99.95%) are Qualified and 3,484 (0.05%) are Failed due to missing values. Stack raising was perf

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Towards time-resolved MicroED grid preparation using mix-and-inject gas dynamic virtual nozzles

Recent progress in gas dynamic virtual nozzle (GDVN) technologies in combination with high-brilliance synchrotron and X-ray free-electron lasers (XFELs) has allowed the visualization of protein dynamics in crystallo by mixing macromolecular protein crystals with a substrate using tunable mixing times on the order of milliseconds to seconds prior to serial X-ray diffraction data collection. This has become the method of choice for high-resolution structure determination of intermediate states. However, such experiments require large counts of crystals of proper sizes for high-resolution data collection, and premium beam times for screening efforts. Cryogenic microcrystal electron diffraction (MicroED) represents a complementary technique that may be a more accessible avenue for time-resolved nanocrystallography compared with serial X-ray diffraction experiments. MicroED can produce full diffraction datasets from just a few submicrometre-thick crystals, and the approach is more readily accessible, requiring standard cryogenic transmission electron microscopy (TEM) equipment available at many universities and institutes. Cryogenic MicroED, like other forms of cryo-EM, begins with rapidly freezing biological material on electron microscopy grids. In the case of MicroED, micro- to nano-crystals (<500 nm thick) are deposited onto electron microscopy grids and plunge-frozen for subsequent electron diffraction data collection. Here, we have incorporated GDVN technology developed originally for XFEL experiments into the freezing process as a first step towards time-resolved studies. We describe the limited deposition efficiency of the model MicroED protein proteinase K on TEM grids using GDVNs, preceding sample vitrification and successful MicroED data collection. We discuss both the initial results from such experiments and the methodological challenges in developing this approach into a reliable workflow for millisecond-to-second time-resolved structural studies of macromolecules. Our results promise a strategy to deposit crystals on grids using GDVNs and determine high-resolution structures by MicroED, constituting a first step towards development of time-resolved MicroED experiments.

MicroED↗

A Study on the Use of XCT and FEA to Predict the Elastic Behavior of Additively Manufactured Parts of Cylindrical Geometry

Defining general criteria for the acceptability of defects within industrial components is often complicated, since the specific load conditions and the criticality of the given application should be considered individually. In order to minimize the risk of failure, high safety factors are commonly adopted during quality control. However this practice is likely to cause the rejection of components whose defects would be instead acceptable if a more sound knowledge of the component behaviour were achieved. Parts produced by additive manufacturing (AM) may suffer from various defects, including micro- or macro-holes, delamination and microstructural discontinuities. Such processes, which are specially suitable for one-off components, require robust and reliable inspection before a part is accepted or rejected, since the refusal of even a single part at the end of the production process represents a significant loss. For this reason, it would be very useful to simulate in a reliable way whether a certain defect is truly detrimental to the proper working of the part during operation or whether the component can still be used, despite the presence of a defect. To this purpose, the paper highlights the benefits of a synergistic interaction between Industrial X-ray computed tomography (XCT) and finite element analysis (FEA). Internal defects of additively manufactured parts can be identified in a non-destructive way by means of XCT. Then FEA can be performed on the XCT-based virtual model of the real component, rather than on the ideal CAD geometry. A proof of concept of this approach is proposed here for a reference construct produced in an Aluminium alloy by AM. Numerical results of the proposed combined XCT–FEA procedure are contrasted with experimental data from tensile tests. The findings sustain the reliability of the method and allow to assess its full provisional accuracy for parts of cylindrical geometry designed to operate in the elastic field. The paper moves a step beyond the present application limits of tomography as it is currently employed for AM parts and it evidences instead the possibility of extending the usage of tomography to acceptance testing and prediction of operative behaviour.

42 ENGINEERING↗

Redshift remapping and cosmic acceleration in dark-matter-dominated cosmological models

The standard relation between the cosmological redshift and cosmic scale factor underlies cosmological inference from virtually all kinds of cosmological observations, leading to the emergence of the LambdaCDM cosmological model. This relation is not a fundamental theory and thus observational determination of this function (redshift remapping) should be regarded as an insightful alternative to holding its standard form in analyses of cosmological data. We present non-parametric reconstructions of redshift remapping in dark-matter-dominated models and constraints on cosmological parameters from a joint analysis of all primary cosmological probes including the local measurement of the Hubble constant, Type Ia supernovae, baryonic acoustic oscillations (BAO), Planck observations of the cosmic microwave background (CMB) radiation (temperature power spectrum) and cosmic chronometers. The reconstructed redshift remapping points to an additional boost of redshift operating in late epoch of cosmic evolution, but affecting both low-redshift observations and the CMB. The model then predicts a significant difference between the actual Hubble constant, h=0.48±0.02, and its local determination, h obs =0.73±0.02. The ratio of these two values coincides closely with the maximum expansion rate inside voids formed in the corresponding open cosmological model with Ω m =0.87±0.03, whereas the actual value of the Hubble constant implies the age of the Universe that is compatible with the Planck LambdaCDM cosmology. The new dark-matter-dominated model with redshift remapping provides excellent fits to all data and eliminates recently reported tensions between the Planck LambdaCDM cosmology, the local determination of the Hubble constant and the BAO measurements from the Ly α forest of high-redshift quasars.

79 ASTRONOMY AND ASTROPHYSICS↗

Beam-spin Asymmetries from Semi-inclusive Pion Electroproduction

We have measured the moment A{sup sin{phi}}{sub LU} corresponding to the polarized electron beam-spin asymmetry in SIDIS. A{sup sin{phi}}{sub LU} is a twist-3 quantity providing information about quark-gluon correlations. Data were taken with the CLAS Spectrometer at Jefferson Lab using a 5.498 GeV longitudinally polarized electron beam and an unpolarized liquid hydrogen target. All three pion channels (pi{sup +}, pi{sup 0} and pi{sup -}) were measured simultaneously over a large range of kinematics within the virtuality range Q{sup 2} ~ 1.0-4.5 GeV{sup 2}. The observable was measured with good statistical precision over a large range of z, P{sub T}, x{sub B}, and Q{sup 2}, which permits comparison with several reaction models. The discussed measurements provide an upgrade in statistics over previous measurements, and serve as the first evidence for the negative sign of the {pi}{sup -} sin{phi} Moment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

VERAView

In order to support engineering analysis of Virtual Environment for Reactor Analysis (VERA) model results, the Consortium for Advanced Simulation of Light Water Reactors (CASL) needs a tool that provides visualizations of HDF5 files that adhere to the VERAOUT specification. VERAView provides an interactive graphical interface for the visualization and engineering analyses of output data from VERA. The Python-based software provides instantaneous 2D and 3D images, 1D plots, and alphanumeric data from VERA multi-physics simulations.

Lee, RonaldW.↗

Secure Enclaves: An Isolation-centric Approach for Creating Secure High Performance Computing Environments

High performance computing environments are often used for a wide variety of workloads ranging from simulation, data transformation and analysis, and complex workflows to name just a few. These systems may process data at various security levels but in so doing are often enclaved at the highest security posture. This approach places significant restrictions on the users of the system even when processing data at a lower security level and exposes data at higher levels of confidentiality to a much broader population than otherwise necessary. The traditional approach of isolation, while effective in establishing security enclaves poses significant challenges for the use of shared infrastructure in HPC environments. This report details current state-of-the-art in virtualization, reconfigurable network enclaving via Software Defined Networking (SDN), and storage architectures and bridging techniques for creating secure enclaves in HPC environments.

97 MATHEMATICS AND COMPUTING↗

VERAView

In order to support engineering analysis of Virtual Environment for Reactor Analysis (VERA) model results, the Consortium for Advanced Simulation of Light Water Reactors (CASL) needs a tool that provides visualizations of HDF5 files that adhere to the VERAOUT specification. VERAView provides an interactive graphical interface for the visualization and engineering analyses of output data from VERA. The Python-based software provides instantaneous 2D and 3D images, 1D plots, and alphanumeric data from VERA multi-physics simulations.

Lee, Ronald W.↗

On the effective spectral emissivity of clear skies and the radiative cooling potential of selectively designed materials

Thermophotonic devices are optically designed to be spectrally selective in order to reject heat to outer space through atmospheric windows of low thermal absorption. The determination of thermal equilibrium temperatures for thermophotonic devices requires the knowledge of the effective spectral emissivity of the sky. In this work, individual contributions of participating gases and aerosols to the spectral values of the sky emissivity are analyzed in the entire infrared spectrum as well as in seven distinct bands for which water vapor either dominates or is virtually transparent to infrared radiation. We also propose high-fidelity correlations for the effective sky emissivity as functions of the normalized ambient partial pressure of water vapor () for both broadband and for the seven spectral bands. Here, the correlations are derived using a combination of ground experimental data, high resolution spectral data for the main atmospheric constituents and spectral models to reconstruct the spectral distribution of infrared thermal radiation from the atmosphere to the ground. These results enable direct calculation of the equilibrium temperature and cooling efficiency of radiative cooling devices in terms of meteorological conditions observed at the surface level. For hot and dry conditions, the passive radiative coolers have a cooling potential of 150.8 W m–2 while for humid conditions, the coolers are mostly ineffective.

14 SOLAR ENERGY↗

Shear reflectivity compensation in full-waveform inversion using least-squares reverse-time migration

SUMMARY The computational cost of elastic-waveform inversion is too high for inverting PP reflections, while using acoustic full-waveform inversion (FWI) is inaccurate because it does not depend on the shear modulus/velocity/impedance that affects elastic PP wavefield amplitudes. To solve this problem, we develop a waveform inversion method that uses acoustic least-squares reverse-time migration (LSRTM) to compensate the shear reflectivity for acoustic FWI. Our method is based on the quasi-elastic-wave equation developed by Chapman et al. (2014). The quasi-elastic-wave equation uses a linearized acoustic-wave equation with shear modulus μ as a virtual source to correct the acoustic PP wavefield amplitudes toward elastic ones. Our waveform inversion method inverts for elastic parameters by minimizing the L2 norm of the difference between recorded and predicted PP reflections modelled using the quasi-elastic-wave equation. Numerical tests on synthetic and field data show that our method can properly handle the amplitudes of elastic PP reflections and provides an accurate estimate of the P- and S-wave velocities/impedances and, in some cases, the density. The method does not need the computationally expensive numerical solution to the elastic-wave equation. It also gives a better estimate of elastic parameters than a pure LSRTM method for elastic PP reflections.

Feng, Zongcai↗