Abstract: Methods and systems for optimizing drilling operations in a wellbore using a drill string are provided. The methods and systems include measuring a first formation characteristic with at least one sensor, measuring a second formation characteristic by means of a hydraulic test, the at least one second formation characteristic being different from the at least one first formation characteristic, generating a model to represent a formation around the wellbore, the model incorporating the first formation characteristic and the second formation characteristic, and performing a drilling operation based on the generated model.

Abstract: A method for partial differential equation inversion, the method including receiving measured data do; selecting an objective function having first and second measures N1 and N2, wherein the objective function depends on three independent variables V, u, and f, V being a perturbation of a wave equation operator L from a background operator L0, u being a wavefield that satisfies the wave equation operator L, and f being a source function that describes the source of the waves; optimizing with a processor the objective function by finding a minimum or a maximum using the inversion; calculating with the processor solutions V*, u*, and f* of the three independent variables V, u, and f; and generating with the processor an image of an object based on the solutions V*, u*, and f*.

Abstract: In an example implementation, first seismic energy is generated using first seismic sources positioned on an earth's surface. First data including measurements of the first seismic energy is obtained from first geophones positioned at a first depth below the earth's surface. Second data including measurements of the first seismic energy is obtained from second geophones positioned on the earth's surface. Second seismic energy is generated using second seismic sources positioned on an earth's surface and proximal to the second geophones. Third data including measurements of the second seismic energy is obtained from third geophones positioned at the first depth below the earth's surface. A propagation of the first seismic energy along a first path is estimated based on the first, second and third data. One or more characteristics of the target are determined based on the estimate.

Abstract: An interface autonomously collects input data to control the operation or travel path of an object, includes a plurality of independently operating sensors that detect the occurrence of an event. A master fractal controller initiates a space-variant fractalization protocol in order to collect the input data related to the detected event. In response to the input, the master fractal controller sends control signals to a servient controller to control the operation or travel path of the object. The master fractal controller initiates the space-variant fractalization protocol by identifying or defining a fractal space, and by further initiating sub-servient fractals in sub-servient spaces, within which an intersection of a graphical function layer that includes an executable grid with a graphical data layer causes a command or function be executed within the sub-servient fractal space.

Abstract: Systems and methods are provided for reducing error in data compression and decompression when data is transmitted over low bandwidth communication links, such as satellite links. Embodiments of the present disclosure provide systems and methods for variable block size compression for gridded data, efficiently storing null values in gridded data, and eliminating growth of error in compressed time series data.

Abstract: A method and apparatus for forecasting output power of wind turbine in a wind farm. The present invention provides a method for forecasting output power of a wind turbine in a wind farm, including: generating a corrected data set based on environmental data collected from at least one sensor in the wind farm; correcting a weather forecasting model by using the corrected data set; obtaining a forecast value of wind information at the wind turbine based on the corrected weather forecasting model; and forecasting the output power of the wind turbine based on the forecast value and a power forecasting model.

Abstract: Strata surfaces can be identified in well logs and correlated across the well logs taking into account manual corrections. For example, a computing device can receive well logs. The computing device can determine multiple stratum-surface correlations based on the well logs. Then, the computing device can receive user input indicating a correction to a particular stratum-surface correlation. Based on the correction to the particular stratum-surface correlation, the computing device can update some or all of the other stratum-surface correlations.

Abstract: Disclosed are method and apparatus for obtaining residual gravity anomaly. The method comprises obtaining Bouguer gravity anomaly of a target region, determining a first pre-set range corresponding to each sampling point in the target region, obtaining a first regional field value of each sampling point in the first pre-set range through a surface fitting method based on the coordinate and field value of the sampling point, traversing the target region to obtain first regional gravity anomaly of the target region according to the first regional field values of the sampling points in the target region, and obtaining first residual gravity anomaly of the target region according to the Bouguer gravity anomaly and the first regional gravity anomaly. The residual gravity anomaly thus obtained is more accurate, thus enabling an accurate prediction of an underground geological body.

Abstract: A method and apparatus for machine learning for use with automated seismic interpretation include: obtaining input data; extracting patches from a pre-extraction dataset based on the input data; transforming data of a pre-transformation dataset based on the input data and geologic domain knowledge and/or geophysical domain knowledge; and generating augmented data from the extracted patches and the transformed data. A method and apparatus for machine learning for use with automated seismic interpretation include: a data input module configured to obtain input data; a patch extraction module configured to extract patches from a pre-extraction dataset that is based on the input data; a data transformation module configured to transform data from a pre-transformation dataset that is based on the input data and geologic domain knowledge and/or geophysical domain knowledge; and a data augmentation module configured to augment data from the extracted patches and the transformed data.

Abstract: Systems and methods are provided to correct seismic data for the undesired effects caused by near surface velocity variations. In one embodiment, a method includes receiving travel time data for a near surface region and estimating an initial velocity model for the near surface region using the travel time data. The method can include updating the velocity model by performing an inversion iteration of including inversion of travel times to estimate slowness. The process can also include calculating at least one long-wavelength static for the near surface region. The long-wavelength statics may be used to correct for undesired effects caused by near surface velocity variations.

Abstract: Methods for full waveform inversion (FWI) in the midpoint-offset domain include using a computer system to sort seismic traces into common midpoint-offset bins (XYO bins). For each XYO bin, a linear moveout correction is applied to a collection of seismic traces within the XYO bin. The collection of seismic traces is stacked to form a pilot trace. The computer system determines a surface-consistent residual static correction for each seismic trace. The computer system determines that the surface-consistent residual static correction for each seismic trace is less than a threshold. Responsive to the determining that the surface-consistent residual static correction is less than the threshold, the computer system stacks the collection of seismic traces to provide the pilot trace. The computer system groups the pilot traces for the XYO bins into a set of virtual shot gathers. The computer system performs one-dimensional FWI based on each virtual shot gather.

Abstract: Systems and methods develop a three-dimensional model of a subterranean formation based on vertical seismic profiles at a plurality of well locations.

Abstract: A new method for iteratively picking the seismic first breaks and conducting imaging of the near-surface velocity structures in an iterative fashion is provided that the first-break picks of the input seismic data are applied to image the near-surface velocity structures and the calculated travel times associated with the updated velocity structures are applied to help refine the first-break picks in the first break picking process until first-break picks satisfy a number of quality control criteria, statics solutions are optimized, and the near surface imaging reaches an acceptable data misfit. This invention produces a velocity model that can be used for near surface statics corrections or for the prestack depth migration.

Abstract: Systems and methods for editing multi-Z horizons interpreted from seismic data are provided. A multi-Z horizon having a plurality of surfaces is visualized within a two-dimensional (2D) representation of seismic data displayed via a graphical user interface (GUI) of an application executable at a computing device of a user. Input is received via the GUI from the user for editing one or more of the plurality of surfaces of the multi-Z horizon within a current view of the displayed 2D representation of the seismic data. A location of the received input relative to each of the plurality of surfaces within the current view is determined. The one or more surfaces of the multi-Z horizon are modified based on the location of the received input within the current view. The visualization of the multi-Z horizon within the GUI is updated, based on the modified one or more surfaces.

Abstract: Provided in this disclosure are systems and methods for selection of petrophysical techniques to model reservoirs. Reservoir properties were calculated using two techniques—a deterministic technique and an optimizing petrophysics technique, and simulation models were developed. The technique that yielded a simulation model that aligned more accurately and consistently with the reservoir data was selected to further develop the reservoir model.

Abstract: The invention provides a method of detecting a seismic event, which comprises acquiring (110) a digital signal x characteristic of a signal measured by at least one seismic sensor, and calculating (130) a time-frequency distribution for at least one section of a given duration of said signal, in a given frequency band. For each frequency of said frequency band, the calculated time-frequency distribution is normalized. The method also comprises calculating (150) the moving average of the normalized time-frequency distribution ZD, in said frequency band and in a time window, given reference L, centered on the time n; and detecting (160) a seismic event when the average exceeds a predefined threshold value. The invention also provides a corresponding detection system.

Abstract: A method can include receiving data acquired via survey equipment disposed in an environment; analyzing at least a portion of the data via an iterative adaptive approach that implements matching pursuit to generate analysis results; and identifying at least one feature in the environment based at least in part on a portion of the analysis results.

Abstract: The invention notably relates to a computer-implemented method for processing a 4D seismic signal relative to a subsoil, the subsoil including a zone subject to extraction and/or injection, the method comprising: providing the 4D seismic signal; identifying a part of the 4D seismic signal corresponding to a zone of the subsoil distinct from the zone subject to extraction and/or injection; determining a noise model of the 4D seismic signal based on the identified part of the 4D seismic signal; and processing the 4D seismic signal based on the noise model. This improves the field of 4D seismic data processing.

Abstract: A method is described for assessing subsurface structure uncertainty based on at least one subsurface horizon. The method calculates seismic continuity attributes to determine a mappability of the subsurface horizon(s); determines horizontal uncertainty for each fault in vertical uncertainty for each horizon; generates probabilistic scenarios for a subsurface geometry for at least one conceptual model; and generates a map of geological model uncertainty based on the probabilistic scenarios. In some embodiments, the probabilistic scenarios are stochastic simulations. In some embodiments, generating a map of geological model uncertainty is based on information entropy. The method may be executed by a computer system.

Abstract: A method for directional Q compensation of seismic data may comprise calculating angle-dependent subsurface travel times; applying directional Q compensation to the prestack seismic data to obtain Q-compensated data in time-space domain, wherein the directional Q compensation is based on the angle-dependent subsurface travel times; and using the Q-compensated data to generate an image of the subsurface. Directional Q compensation may comprise determining an angle-dependent forward E operator and an angle-dependent adjoint E* operator using the angle-dependent subsurface travel times; and applying a sparse inversion algorithm using the angle-dependent operators to obtain a model of Q-compensated data.

Abstract: A system and method of scanning an environment and acquiring an image is provided. The system includes a mobile device having a camera and a first position indicator. A scanner having a light emitter and a light receiver is provided. The scanner determining coordinates of surfaces in an environment in response to emitting light with the light emitter and receiving light with the light receiver, the scanner having a second position indicator. One or more processors are provided that determine the position of the mobile computing device and transmits the data between the scanner in response to the first position indicator engaging the second position indicator.

Abstract: An information processing apparatus that is communicable with an in-vehicle device includes a detection section to detect a change of a specific feature at every predetermined timing based on probe information, the specific feature being a feature on a road and used in control of automated driving. The detection section includes a first calculation section to calculate a first index based on the probe information; a second calculation section to calculate a second index based on the probe information; and a comparison section to detect a change of the specific feature by comparing the first index and the second index.

Abstract: Systems and methods for training a model that refines estimated parameter values include computer processors and non-transitory electronic storage that stores subsurface map data sets that correspond to different subsurface volumes of interest, the system configured to obtain training data including unrefined subsurface map data sets specifying estimated parameter values of a first parameter as a function of position within corresponding subsurface volumes of interest, obtain an initial seismic mapping model, generate a conditioned seismic mapping model, and store the conditioned seismic mapping model in the electronic storage.

Abstract: Changes to interpreted horizons in a subsurface representation are propagated to synthetic horizons in the subsurface representation by maintaining relative distances between the synthetic horizons and interpreted horizons. Distance ratios of synthetic horizons between interpreted horizons are applied to changed interpreted horizons to derive new locations of synthetic horizons.

Abstract: A method is disclosed and includes receiving a seismic cube. The seismic cube includes a three-dimensional image of a portion of a subsurface area. The method further includes providing the seismic cube to a machine learning process. The machine learning process includes one or more neural networks used for predicting a location of a subsurface seismic layer in the received seismic cube. The method also includes receiving, from the machine learning process, the prediction of the location of the subsurface seismic layer in the seismic cube.

Abstract: Various embodiments described herein provide methods of hydrocarbon management and associated systems and/or computer readable media including executable instructions. Such methods (and by extension their associated systems and/or computer readable media for implementing such methods) may include obtaining geophysical data (e.g., seismic or other geophysical data) from a prospective subsurface formation (that is, a potential formation or other subsurface region of interest for any of various reasons, but in particular due to potential for production of hydrocarbons) and using a trained machine learning (ML) system for direct hydrocarbon indicators (DHI) analysis of the obtained geophysical data. Hydrocarbon management decisions may be guided by the DHI analysis.

Abstract: A pile head analysis system captures a plurality of images of a construction site including a pile using a camera mounted on an unmanned aerial vehicle (UAV), acquires the images, generates a three-dimensional model of the construction site, detects a pile head from the three-dimensional model, and determines consistency between the detected pile head and preliminarily acquired design information.

Abstract: A method of at least partially deghosting recorded seismic s-waves, wherein recorded seismic data is provided, wherein the recorded seismic data has been recorded at a receiver located beneath the Earth's surface, and wherein the recorded seismic data includes s-wave data. The method may include the steps of finding a model of the Earth's crust for use in deghosting the recorded seismic data using the s-wave data, wherein the model includes at least one region and wherein the model includes the Earth's surface and the location of the receiver, using the model to find a deghosting operator that, when applied to the s-wave data, at least partially deghosts the s-wave data, and applying the deghosting operator to the s-wave data to at least partially deghost the s-wave data.

Abstract: Identifying a low permeable conglomerate diagenetic trap can be implemented according to a method that comprises: determining a first relation curve between a depth and a critical physical property of a known diagenetic trap in a target work area, and a second relation curve between a reservoir physical property of the known diagenetic trap and a designated seismic attribute; determining a third relation curve between the depth and the critical physical property in the target work area and the designated seismic attribute according to the first relation curve and the second relation curve; and performing a diagenetic trap identification of the target work area according to the third relation curve. Identification accuracy of a low permeable conglomerate diagenetic trap can thereby be improved.

Abstract: The method for producing a geological vector model (GVM) based on seismic data includes the step of forming a Model-Grid, which includes creating a network of small units, called patches, to which a relative geological age is assigned, a set of patches with the same relative geological age corresponding to a geological layer, called the geological horizon. The method includes the step of sampling the Model-Grid in two directions perpendicular to each other, enabling the Model-Grid to be sampled in a plurality of vertical planes and the step of forming two-dimensional geological vector models (2DGVM). The step of forming includes forming two-dimensional horizons (Hb) with distinct relative geological ages using the patches belonging to each sampled plane, each two-dimensional geological vector model (2DGVM) corresponding to a vertical plane originating from the sampling of the Model-Grid.

Abstract: System and method for scalable and reliable scheduling of iterative seismic full wavefield inversion algorithms with alternating parallel and serial stages of computation on massively parallel computing systems. The workers are independent, initiating actions and unaware of each other, and given limited information. This enables application of optimal scheduling, load-balancing, and reliability techniques specific to seismic inversion problems. The central dispatcher specifies the structure of the inversion, including task dependency, and keeps track of progress of work. Management tools enable the user to make performance and reliability improvements during the execution of the seismic inversion.

Abstract: The embodiments of the present application include acquiring a monitoring region and each observation point therein; partitioning the monitoring region into N layers of grids according to a seismic source positioning accuracy, wherein a side length of a grid cell of an i-th layer of grid is D/2i-1, i=1, . . . N, and D is an initial side length of the grid cell and not more than a double of a distance between the respective observation points; searching all nodes in a first layer of grid to acquire a node satisfying a preset condition therefrom; from i=2, determining and searching nodes satisfying a first preset requirement in the i-th layer of grid, to acquire a node satisfying the preset condition therefrom, until a search in an N-th layer of grid is completed, wherein a node satisfying the preset condition acquired in the N-th layer of grid is a seismic source point location.

Abstract: Automatic calibration for modeling a field includes performing model calibration iterations based on measurements of the field. Each model iteration includes obtaining, during the field operation, a current measurement of the field, generating likelihood values corresponding to model realizations of the field, where each likelihood value is generated by at least comparing the current measurement of the field to a modeling result of a corresponding model realization, selecting, based on the likelihood values, at least one selected model realization from the model realizations, generating, by at least adjusting a first input value of the at least one selected model realization, an adjusted model realization of the field based on the at least one selected model realization, and adding the adjusted model realization to the model realizations. Automatic calibration for modeling a field further includes generating a calibrated modeling result of the field based on the adjusted model realization.

Abstract: A subsoil region is represented using a grid having columns of cells arranged in an array of horizontal positions. The cells of a column are delineated by geological layering of the subsoil region. The method, used for estimating elastic parameters in the subsoil region, comprises obtaining seismic traces and performing seismic inversion for the columns of cells. The seismic inversion for a column comprises selecting a set of elastic parameter values that minimizes a cost function regardless of statistical distribution of the elastic parameter values. The set of elastic parameter values is selected among a plurality of candidate sets of elastic parameter values each having values of elastic parameters for the cells of the column.

Abstract: A method for generating a plurality of structural models for a geological setting involves identifying a poorly imaged portion of an input seismic data set. A set of geologically valid kinematic base models are selected to define a kinematic evolution scenario. A correspond set of kinematic base model parameters is defined for the models. A structural model is generated for both a well-imaged portion of the seismic data and the poorly imaged portion. A misfit between the input horizons and the and the modeled horizons is calculated and steps are repeated for a predetermined number of iterations to produce a best-fit model. The steps are repeated then to produce a plurality of best-fit models of geologically plausible solutions for the geological setting. The method is particularly suitable for complex geological settings.

Abstract: An inspection specifying unit specifies a member position that is a point and a type of damage that is a soundness determination target on the basis of building structure information. An inspection data analysis unit acquires the image data from an inspection data acquisition unit, analyzes the image data, and determines from the image data, whether or not there is damage of the type of damage specified as the inspection target, and the degree of damage in a case where there is the damage. A damage influence degree calculation unit calculates the degree of damage influence for each damage. A soundness determination unit determines the soundness of the entire building on the basis of the degree of damage influence for each damage. A soundness output unit outputs the soundness determined in to a display, a printer, or the like.

Abstract: The present invention discloses a method and system for processing gravity and magnetic data in geological resource exploration. The method includes: acquiring first (i) potential field data and (ii) gradient data of an observation surface, performing upward continuation of the acquired data using a wave-number domain conversion method to obtain second and third gradient data and second potential field data, and determining third potential field data using a fourth-order explicit scheme Milne method according to the first, second, and third gradient data, and the second potential field data; calculating fourth gradient data using an ISVD method according to the third potential field data; and correcting the third potential field data using a fourth-order implicit scheme Simpson method according to the fourth gradient data, the first potential field data, and the first and second gradient data to obtain corrected third potential field data.

Abstract: A system and method for attenuating surface waves in common shot gathers of seismic data recorded by a set of geophones by: iteratively executing a genetic algorithm over a plurality of generations to generate an optimal one-dimensional (1D) Earth model based on the common shot gather data by, successively refining a pool of candidate Earth models to better fit the common shot gather data, until optimal Earth models in sequential generations converge; generating synthetic surface wave data based on the optimal Earth model and canceling the synthetic surface wave data from the common shot gather data to generate new common shot gather data that reduces the noise due to surface waves; and iteratively executing the genetic algorithm over each new common shot gather data until optimal Earth models generated in sequential iterations of the genetic algorithm converge.

Abstract: This disclosure describes processes and systems for generating a seismic image of a subterranean formation from recorded seismic data gathers obtained in a marine seismic survey of the subterranean formation. The seismic data comprises recorded pressure and vertical velocity wavefields that are used to separate the recorded pressure wavefield into upgoing and downgoing pressure wavefields. A seismic image is computed from the subterranean formation based on a product of the downgoing pressure wavefield and a migration operator applied to the upgoing pressure wavefield. The downgoing pressure wavefield is a boundary source wavefield and the upgoing pressure wavefield is boundary receiver wavefield of the migration operator. The seismic image is iteratively updated by computing a residual seismic image based on the upgoing pressure wavefield and adding the residual seismic image to the seismic image.

Abstract: In an example, the autonomous vehicle (“AV”) can be configured among the other vehicles and railway to communicate with a rider on a peer to peer basis to pick up the rider on demand from a location on a track, like a railway, tram or other track, rather than the rider being held hostage to -a fixed, railway schedule. The rider can have an application on his/her cell phone, which tracks each of the AVs. and contact them using the application on the cell phone. In an example, the AV is configured for both on-track and off track operation with different operating parameters for on-track and off track, including speed, degree of autonomy, sensors used etc.

Abstract: Methods of identifying a flying object using digital imaging that may include: obtaining data of a propagating wavefield through a propagating volume that includes a volume above the earth's surface; obtaining a reference digital image of the propagating volume; selecting a holographic computational method of wavefield imaging; selecting a wavefield based on one or more parameters; calculating a sampling ratio by dividing a number of data samples in the data subset by a number of image samples in the data subset; decimating the data subset; generating a new digital image based on the selected holographic computational method of imaging, the decimated data subset, and parameters corresponding to the data subset; and determining a quantitative difference measure between the reference digital image and the new digital image, and image quality.

Abstract: Automatic propagation of real-world parent seismic images to efficiently generate a collection of realistic synthetic child training images to train a model for accurate automatic seismic interpretation. A 3D structural model in a present-day geological space (e.g., GB) depicting subsurface locations of particles (e.g., in region B) may be transformed by a 3D coordinate space transformation (e.g., uvtB) to a depositional space (e.g., G*B) depicting past depositional locations of those particles (e.g., corresponding depositional region B). A real-world parent image depicting subsurface locations of particles (e.g., in region A) may be transformed, via a forward transformation (e.g., uvtA), to a depositional seismic image in the depositional space of the three-dimensional structural model (e.g., G*A=G*B). A reverse transformation (e.g., uutB?1) may transform the depositional seismic image from the depositional space into synthetic child training images in the present-day geological space (e.g.

Abstract: A state determination device determines a state of a paint device including an air motor that rotates by a supply of a gas, a rotational shaft that serves as a rotational center of the air motor, and a rotary atomizer head that is connected to the rotational shaft and sprays paint. The paint device further includes a vibration detection unit that detects a vibration of the paint device, and the state determination device includes a frequency component acquisition unit and an abnormality determination unit. The frequency component acquisition unit acquires a frequency component caused by a contact between the rotational shaft and the air motor from a detection result by the vibration detection unit. The abnormality determination unit determines whether an abnormality of the air motor exists or not based on an acquisition result by the frequency component acquisition unit.

Abstract: A coarse error correction system for detecting, predicting, and correcting errors in neural networks is provided. The coarse error correction system receives a first set of statistics that are computed from values collected from a neural network during a training phase of the neural network. The coarse error correction system computes a second set of statistics based on values collected from the neural network during a run-time phase of the neural network. The coarse error correction system detects an error in the neural network during the run-time phase of the neural network by comparing the first set of statistics with the second set of statistics. The coarse error correction system increases a voltage setting to the neural network based on the detected error.

Abstract: A method includes receiving desired locations of nodes for deployment on a seabed of a seismic survey where each of the nodes includes a sealed housing and, within the sealed housing, at least one battery and spaced seismic sensors electrically powered by the at least one battery; determining locations of the nodes as deployed on the seabed where at least some of the determined locations differ from their corresponding desired locations; acquiring seismic data sensed by the spaced seismic sensors of the nodes where the acquired seismic data corresponds to the determined locations; and, based at least in part on the acquired seismic data, a spacing of the spaced seismic sensors and the desired locations, generating seismic data for the desired locations.

Abstract: A computer executable algorithm adapted to propagate a boundary surface of a seed that is placed within a region of interest of a visual representation of a 3D seismic data so as to follow a natural contour of said region of interest, wherein said algorithm is executable to: (i) generate at least one attribute volume comprising at least on attribute derivable from said 3D seismic data set; (ii) generate at least one characteristic parameter for a plurality of candidate events of said 3D seismic data within a predefined gate region located forward of said propagating boundary surface; (iii) generate and assign a probability characteristic for said plurality of candidate events based on said at least one attribute volume and said at least one characteristic parameter; and propagate said boundary surface towards and incorporating any one of said plurality of candidate events that fulfils an acceptance criteria of said probability characteristic so as to generate a surface along the natural contour of said region

Abstract: Inversion for marine seismic imaging of a full reflected wavefield can include generating an image of a subsurface formation by full wavefield migrating a recorded seismic wavefield and generating numerically modeled data using the image. A mismatch between the numerically modeled data and the seismic wavefield can be determined. Responsive to determining that the mismatch exceeds an inversion match threshold, the image can be updated using the mismatch between the numerically modeled data and the seismic wavefield.

Abstract: This disclosure is directed to processes and systems that generate enhanced-resolution seismic images by interpolating sparsely recorded seismic data. Structured dictionary learning is employed to train a set of basis vectors, called “atoms,” and corresponding sparse coefficients on patches of the recorded seismic data. The atoms are constrained to represent the geometric structure of reflection events in the recorded seismic data gather. Linear combinations of the atoms are used to compute interpolated patches over a finer receiver-coordinate grid. The interpolated patches replace the original patches in the recorded seismic data to obtain interpolated seismic data that can be used to generate an image of the subterranean formation.

Abstract: A seismic rock physics inversion method based on a large area tight reservoir includes steps: building a multi-scale rock physics model; analyzing fluid sensitivities of rock physics parameters in two scales of acoustic logging and ultrasonic wave, and sifting the rock physics parameters that are most sensitive to a porosity and a gas saturation in a plurality of observation scales; building a single-well rock physics template, preferably a standard template; considering lateral variations and heterogeneity of reservoir geological features, fine-tuning input parameters of the rock physics template according to gas testing situations of all wells in a large work area, optimizing the whole work area and building a three-dimensional work area rock physics template data volume, and combining the data volume with pre-stack seismic inversion to calculate a porosity and a saturation of a target layer; and smoothing a result and finally outputting a reservoir parameter inversion data volume.

Abstract: The present disclosure provides a technique for marine seismic imaging that processes data acquired from two or more different seismic surveys in a combined manner to advantageous effect. The different seismic surveys may use seabed sensors at same positions on the seabed, but they may have different shot locations and may be performed at different times. In one use case, the technique may be used to image a subsurface location that is difficult to image using either survey alone. In another use case, the technique may be used to image a subsurface location under an obstruction. The technique may also be utilized to efficiently monitor a reservoir over time.