摘要:
Method for determining best and worst cases for values of model parameters such as porosity and shale volume fraction generated by non-unique matrix inversion of physical data such as seismic reflection amplitudes. The matrix is diagonalized, and then orthonormal basis vectors associated with insignificant diagonal elements are used to generate upper and lower bounds on the solution. Best and worst case solutions are determined as linear combinations of the null basis vectors, where the expansion coefficients are determined by making a best fit to the upper and lower bounds.
摘要:
Method for constructing an integrated rock physics model that simulates both shale anisotropy and stress-induced anisotropy of clastic rocks. In the model, the total pore volume is divided into three parts according to the estimated shale volume and effective stress: (1) clay-related pores, (2) sand-related pores, and (3) microcracks (mainly in the sand component). The pore space is then partitioned into the clay-related and sand-related pores using a scheme first disclosed by Xu and White in 1995. The model simulates shale anisotropy via the preferred orientation of clay-related pores and stress-induced anisotropy via the preferred orientation of microcracks, which is controlled by the differential stresses. Laboratory measurements or well logs are needed to establish a relationship between crack density and the effective stress.
摘要:
Method for predicting lithology and porosity of subsurface rocks from seismic reflection data. The seismic data is inverted to yield elastic properties of the rocks such as the compressional and shear impedances. A rock physics model is built to relate porosity, the shale volume fraction, the fluid content of the rock and the elastic properties of the rock. The model is run backward in a second inversion process to solve for porosity and lithologic properties such as the shale volume fraction.
摘要:
Method for obtaining rock parameters such as porosity and vshale directly from inversion of seismic data corresponding to a single trace location. This method is distinguished from existing methods that obtain elastic properties from inversion of seismic data, then relate the elastic parameters to rock lithology parameters such as porosity or vshale because it is accomplished in one step, can incorporate anisotropy and does not require multiple trace locations for stability. The data are separated into partial stacks, and a wavelet is specified for each stack. A set of linearized equations are constructed relating seismic reflectivity to changes in elastic parameters, and another set of linearized equations is constructed relating the changes in elastic parameters to the lithologic parameters. The linearized reflectivity equations are combined with the linearized rock physics equations, convolved with the specified wavelets, and equated to the seismic data. The resulting matrix equations are then inverted and a solution is obtained for all offsets simultaneously.
摘要:
Method for obtaining rock parameters such as porosity and vshale directly from inversion of seismic data corresponding to a single trace location. This method is distinguished from existing methods that obtain elastic properties from inversion of seismic data, then relate the elastic parameters to rock lithology parameters such as porosity or vshale because it is accomplished in one step, can incorporate anisotropy and does not require multiple trace locations for stability. The data are separated into partial stacks, and a wavelet is specified for each stack. A set of linearized equations are constructed relating seismic reflectivity to changes in elastic parameters, and another set of linearized equations is constructed relating the changes in elastic parameters to the lithologic parameters. The linearized reflectivity equations are combined with the linearized rock physics equations, convolved with the specified wavelets, and equated to the seismic data. The resulting matrix equations are then inverted and a solution is obtained for all offsets simultaneously.
摘要:
Systems and methods which provide electromagnetic subsurface mapping to derive information with respect to subsurface features whose sizes are near to or below the resolution of electromagnetic data characterizing the subsurface are shown. Embodiments operate to identify a region of interest (203) in a resistivity image generated (202) using electromagnetic data (201). One or more scenarios may be identified for the areas of interest, wherein the various scenarios comprise representations of features whose sizes are near to or below the resolution of the electromagnetic data (204). According to embodiments, the scenarios are evaluated (205), such as using forward or inverse modeling, to determine each scenarios' fit to the available data and further to determine their geologic reasonableness (206). Resulting scenarios may be utilized in a number of ways, such as to be substituted in a resistivity image for a corresponding region of anomalous resistivity for enhancing the resistivity image (207).
摘要:
Method for estimating porosity and water saturation or other geological parameters of a subsurface region when the lithology of the region is unknown, requiring only geophysical data from remote surveys. The geophysical data are inverted (14-18) treating lithology as a third, and discrete, unknown model parameter (12) to be solved for in the inversion. A technique for solving mixed integer non-linear programs may be used. Suitable rock physics relationships (13) are used to relate the desired geological parameters to the geophysical model parameters required for simulating the geophysical data in the inversion process.
摘要:
Method for determining best and worst cases for values of model parameters such as porosity and shale volume fraction generated by non-unique matrix inversion of physical data such as seismic reflection amplitudes. The matrix is diagonalized, and then orthonormal basis vectors associated with insignificant diagonal elements are used to generate upper and lower bounds on the solution. Best and worst case solutions are determined as linear combinations of the null basis vectors, where the expansion coefficients are determined by making a best fit to the upper and lower bounds.
摘要:
Method for reducing a 3D joint inversion of at least two different types of geophysical data acquired by 3-D surveys (21) to an equivalent set of ID inversions. First, a 3D inversion is performed on each data type separately to the yield a 3-D model of a physical property corresponding to the data type (22). Next, a ID model of the physical property is extracted at selected (x,y) locations. A ID simulator (23) and the ID model of the physical property is then used at each of the selected locations to create a synthetic ID data set at each location (24). Finally, the ID synthetic data sets for each different type of geophysical data are jointly inverted at each of the selected locations, yielding improved values of the physical properties. Because the joint inversion is a ID inversion, the method is computationally advantageous, while recognizing the impact of 3-D effects.
摘要:
Systems and methods which provide electromagnetic subsurface mapping to derive information with respect to subsurface features whose sizes are near to or below the resolution of electromagnetic data characterizing the subsurface are shown. Embodiments operate to identify a region of interest (203) in a resistivity image generated (202) using electromagnetic data (201). One or more scenarios may be identified for the areas of interest, wherein the various scenarios comprise representations of features whose sizes are near to or below the resolution of the electromagnetic data (204). According to embodiments, the scenarios are evaluated (205), such as using forward or inverse modeling, to determine each scenarios' fit to the available data and further to determine their geologic reasonableness (206). Resulting scenarios may be utilized in a number of ways, such as to be substituted in a resistivity image for a corresponding region of anomalous resistivity for enhancing the resistivity image (207).