Integrating seismic, CSEM, and well-log data for reservoir characterization

Size: px
Start display at page:

Download "Integrating seismic, CSEM, and well-log data for reservoir characterization"

Transcription

1 H o n o r a r y L e c t u r e Integrating seismic, CSEM, and well-log data for reservoir characterization Lucy MacGregor, RSI S urely seismic tells you everything you need to know about the Earth? This is a widely held view and there is no doubt that seismic methods are extremely powerful, generally applicable, and can be used to generate high-resolution images of structure and stratigraphy which can guide exploration, appraisal and development projects. However, there are some situations where seismic data fail to answer the geophysical question of interest. It is these situations where complementary sources of information must be used to obtain the information required. Well-log data provide a useful starting point to understand the requirement for, and power of, data integration. A variety of measurements are routinely made when a well is logged. The example in Figure 1 shows resistivity, VP, VS, acoustic impedance, density, and Poisson s ratio. A petrophysicist analyzing this well will integrate these measurements to provide an interpretation of the lithology and fluid properties. In this example a clean sand reservoir (shown in yellow in the second track) is overlain by a shale layer (grey), and is charged with oil in its upper part. Looking first at the measurements of seismic properties, the largest variations in the curves occur at lithology boundaries, for example between the shales and the underlying sand reservoir. In contrast, the effect on the seismic properties of a change in fluid is small. The black curve in the VP, VS, acoustic impedance and Poisson s ratio plots represents the in-situ measurement; the blue curve represents the equivalent response when the oil in the reservoir is replaced by water. The difference between the curves is small, illustrating the extremely poor sensitivity of seismic properties to changes in the fluid type. That is not to say that it is impossible to extract fluid information from seismic data; sometimes this is an extremely successful strategy. However, in many situations it is difficult or impossible to do so with certainty. The resistivity measurement, in contrast, shows a high degree of sensitivity to fluid content, but a relatively low sensitivity to changes in lithology. In the oil-saturated region of the reservoir, the resistivity is 1 2 two orders of magnitude higher than in the surrounding water-saturated sediments, where the variation in resistivity is small. Thus, the resistivity log provides a key measurement, used by petrophysicists in determining water saturation and fluid properties in a well. Often the properties of a reservoir are required away from the well bore, and to do this remote sensing methods are required. The lower left plot of Figure 1 shows an example of an acoustic impedance volume derived from inversion of a 3D seismic volume. Areas of low impedance can be correlated with known hydrocarbon accumulations at well locations, and by looking for areas with similar seismic properties other hydrocarbon-saturated zones can perhaps be predicted. However, the seismic method is sensitive to impedance and impedance contrasts, which the fluid-substitution study above 268 The Leading Edge March 2012 illustrates are inherently low in sensitivity to saturation effects (red circled areas in Figure 1). Electromagnetic methods allow remote measurement of resistivity which is extremely sensitive to the properties and distribution of fluids in a structure. The lower right panel in Figure 1 shows a resistivity section derived from inversion of marine CSEM data acquired over a gas reservoir (shown in red). In the same way that multiple data types are used in an integrated interpretation of well-log data, such resistivity information from electromagnetic surveys is complementary to seismic data, and can improve the constraint on fluid properties when used in an integrated geophysical interpretation. Of course no method is perfect. There are potential pitfalls in each when applied to resolve a given geophysical requirement. In the case of CSEM data analysis, it is an obvious (but surprisingly often forgotten) fact that the method measures resistivity and not hydrocarbons directly. Whereas zones of high resistivity can be caused by the presence of hydrocarbon fluids, there are many other causes, such as the presence of cemented sandstone, limestones, volcanics or coals. An interpretation of a CSEM data set to provide a profile, section, or volume of the Earth (itself a nonunique process) must be followed by an equally careful interpretation of the resulting resistivity data in terms of rock and fluid properties. The same is true of seismic data that provide measurements of reflectivity, or impedance (or a host of derived attributes and properties). These seismic properties must also Figure 1. (top) Example well-log suite. Tracks show (left to right) water saturation, lithology, resistivity, VP, VS, acoustic impedance, and Poisson s ratio. (lower left) Impedance volume derived from inversion of a 3D seismic volume. (lower right) Resistivity section derived from inversion of marine CSEM data.

2 H o n o r a r y L e c t u r e Figure 2. (left) Well-log suite from a North Sea gas reservoir. The effect of fluid substitution on VP, VS, and density is shown in brown (for 80% gas saturation) and red (for 20% gas saturation). (right top panel) The effect of gas saturation plotted in the P-impedance Poisson s ratio domain. Blue shows the response when the gas is water wet, red the response for a fizz gas (20% saturation), and brown for a commercial 80% saturation. (right bottom) AVA response for the same three cases. be carefully interpreted in terms of the underlying rock and fluid properties that give rise to them, and once again this interpretation has uncertainty associated with it. For example AVO and amplitude anomalies may be caused by lithology or thin-bed effects rather than fluids, and saturation is often difficult to interpret on the basis of seismic data alone. For any given geophysical question, the most robust answer will be obtained by using the tool, or combination of tools best suited to the task, and integrating the resulting data within a rock physics framework, to provide a shared Earth model that is geologically reasonable, and consistent with each of the geophysical data types available. This idea will be illustrated first using a synthetic case study, before moving on to a real example from the North Sea. Synthetic case study The synthetic example is based on a gas reservoir in the North Sea (Figure 2) in which a clean sand reservoir saturated with gas lies between 2100 and 2250 m. The question we wish to answer using geophysical methods is: What is the gas saturation of the reservoir away from the well location? The left side of Figure 2 shows the effect on the seismic response of changing the gas saturation, plotted in the Pimpedance Poisson s ratio domain, for a water-saturated reservoir (blue), a reservoir containing a low saturation of gas (20% = red) and an 80% saturation (brown). It is clear that the seismic response can detect the presence of the gas; there is good separation between the responses in the water-saturated and gas-saturated cases. However, this does not answer the original question because the responses of the reservoirs saturated at 80% and 20% are almost identical. Although seismic data can be used to detect the gas, it cannot usefully constrain the saturation. The conclusions are the same when the AVA response is considered. Whereas there is a large AVA effect caused by the presence of the gas, the difference between the responses at low and high saturations is small, and would certainly be hard to identify in real data (the well known fizz-gas problem). Turning now to resistivity, there are a number of models describing the electrical properties of rocks saturated with fluids (e.g., Han et al., 2011; Gelius and Wang, 2008; Carcione et al., 2007; Waxman and Smits, 1968; Hashin and Shtrikman, 1962). In this case, we use the simplest of these, Archie s law (Archie, 1942), an empirical equation describing the relationship between bulk resistivity of a two phase medium, RT, and its porosity φ, water saturation Sw and the resistivity of the interstitial fluids, Rw. The constants a, m, and n (respectively, the tortuosity, cementation exponent, and saturation exponent) are empirical, and must be calibrated on a case-bycase basis from well-log data. March 2012 The Leading Edge 269

3 Figure 3. Variation of resistivity with hydrocarbon saturation for a range of porosities, calculated from Archie s law using typical parameters of m = n = 2 and a = 1. Figure 4. Sensitivity of the CSEM method to a change in gas saturation for the reservoir shown in Figure 2. Sensitivity is calculated as the percentage difference between the inline electric field response for the gas-saturated reservoir, compared to the water-wet case, and is contoured as a function of source-receiver separation and signal frequency. The white solid contour shows a typical noise floor for the water depth considered here (1 km). Areas of the plot above and to the right of this line are below the noise floor and would not be usable in survey data. (left) 20% gas saturation. (right) 80% gas saturation. Bulk resistivity, R T, derived from Archie s law is plotted in Figure 3 against hydrocarbon saturation for a range of porosities. Taking the 25% porosity curve as an example, Figure 3 shows that the resistivity varies between about 1 Ωm for a saturation of 20%, up to about 25 Ωm when the gas saturation in the reservoir is 80%. This large difference in physical properties clearly would distinguish between low and high saturations if it could be measured remotely. However, note that there is a trade-off between saturation and porosity. A low-saturation, low-porosity medium will have the same resistivity as one with a much higher porosity and higher saturation. Measurements of resistivity can be made using a variety of geophysical methods. Here we concentrate on the marine controlled-source electromagnetic (CSEM) method (Constable and Srnka, 2007, provide a review of the method). The marine CSEM method uses a high-powered source to transmit a low-frequency signal through the sea floor to an array of receivers. Signal frequencies in the range Hz are typically employed. By studying the variation in the received fields as a function of source-receiver separation, geometry, and signal frequency, the resistivity of the subsurface can be determined to depths of typically about 3 km below mudline (depending on the structure and resistivity of the overburden). Figure 4 shows the sensitivity of a CSEM data set to changes in gas saturation for the same reservoir as above. The sensitivity is shown by plotting the difference between the inline electric field amplitude calculated for the gas-saturated reservoir compared to that calculated when the reservoir is water wet. This is plotted as a function of source-receiver separation (between 500 m and 20 km) and signal frequency (between 0.01 Hz and 5 Hz), covering the range of acquisition parameters typically used in the field. Looking first at the plot for the low gas-saturation (20%) case, it is clear that the effect of the gas on the CSEM response is small (at most 3%). Such a small difference would not be resolved in survey data, and so a fizz-gas-saturated reservoir would be indistinguishable from a waterwet reservoir using CSEM data alone. In contrast the high gas saturation (80%) has a large effect on the CSEM response (>15%). Returning to the geophysical question we are trying to answer (What is the gas saturation?), it is clear that whereas the seismic method can identify the presence of gas (but not its saturation), the CSEM method should be sensitive to the difference between low and high saturations. However, although sensitivity to the property of interest is a necessary condition for successful application of the CSEM method, it is not sufficient. It is also important that this property can be recovered from a survey data set which is finite, noisy, and subject to error. Conversion of a processed CSEM data set into interpretable resistivity sections or volumes relies inversion approaches that are themselves nonunique and subject to uncertainty, and must be carefully constrained to ensure a meaningful answer. This will be illustrated using a synthetic data set generated from a simplified 1D model of the gas reservoir shown in Figure 2, which has been block averaged to provide a background resistivity structure. Embedded in this is a 25% porosity reservoir, which has a resistivity of 25 Ωm when saturated with 80% gas, and 2 Ωm when water-saturated. From this simplified model, synthetic data at 0.1 Hz, 0.3 Hz and 0.7 Hz are generated and contaminated with realistic survey noise. These data can then be inverted to establish how accurately, and under what conditions the saturation in the reservoir can be determined. In this case, the Occam inversion algorithm described by Constable et al. (1987) is used. The left panels of Figure 6 show an unconstrained inversion of the synthetic data shown in Figure 5, along with the 270 The Leading Edge March 2012

4 H o n o r a r y L e c t u r e Figure 5. (left) Simplified 1D model representing the gas reservoir and surrounding structure. The reservoir has a resistivity of 25 Ωm when saturated with 80% gas, and 2 Ωm when water-saturated. (right) Synthetic CSEM amplitude and phase data generated from this model at three frequencies, and contaminated with realistic survey noise. Figure 6. Inversion of synthetic data shown in Figure 5. In each case, the inversion result is shown in red and the true model in dashed gray. (top) Gas-charged reservoir (80% saturation). (bottom) Watersaturated background model. (left) Unconstrained inversion results. (center) The top of the reservoir is assumed known from seismic and is included as a smoothness break. (right) Both top and base reservoir interval are assumed constrained by the seismic and are included. corresponding inversion of data generated from a background model in which the reservoir was water-wet for comparison. The inversion does a good job of recovering the smoothly varying background resistivity trend. However, in the gassaturated case the inversion returns a broad area of elevated resistivity, approximately centered on the known position of the reservoir, but with a peak resistivity considerably lower than in the input model. Because this feature is not present in the background inversion, we can have confidence in attributing it to the high resistivity associated with the gas. However, if the background model were not known (in each case the odd harmonics of the fundamental were also available for processing), interpretation of the result would be much less clear. Is the elevated resistivity shown in the upper left panel of Figure 7 the result of a localized hydrocarbon accumulation, or is it perhaps the result of a more dispersed variation in porosity or lithology? Without further information, the interpretation risk of an unconstrained inversion result such as this would be high. In most situations where CSEM data are acquired, seismic data are also available, and can be used to constrain structure in the CSEM inversion. Whereas an unconstrained inversion can only answer the question What is the resistivity of the Earth?, a poorly posed problem, inclusion of structural information from seismic modifies this to the better posed question What is the resistivity within a seismically defined structure? Constrained inversions for both the gas-charged and background cases are also shown in Figure 6. The center panels of Figure 6 show the result when the depth of top reservoir is assumed constrained by seismic, and is therefore included as a break in the smoothness constraint in the inversion. This constraint has little effect on the inversion for the background structure (lower panels) indicating that such constraints will not create structure where there is none. In the gas-saturated case, the introduction of structural constraints localizes the recovered resistivity to the depth of the known structure. Such a result would provide a clearer indication of high resistivity within this structure. The right panels of Figure 6 show the results when the depths of both the top and bottom of the reservoir are constrained by seismic. For the gas-saturated case, the presence, depth, and resistivity of the gas-charged zone are now clearly resolved. When the structure is known, providing quantitative constraints on the resistivity properties within that structure becomes possible. The models shown in the top row of Figure 6 are mathematically equivalent and based solely on the CSEM data and inversion; we have no reason to prefer one over the other. However, the geological interpretation resulting from each model is likely to be significantly different. So we must question what is actually constrained in the inversion process? The answer is that the parameter best constrained is the transverse resistance, defined as the vertically integrated resistivity across an interval (Harris et al., 2009). After calculating the transverse resistance across the interval of elevated resistivity in the unconstrained inversion result and comparing it to the transverse resistance calculated over the same interval for the true model, the inversion gives a result accurate to approximately 12%. Doing the same comparison for the fully constrained model gives a result accurate to 5 6%. In both cases, the transverse resistance is well constrained; however, with good structural constraint from seismic the accuracy of the recovered transverse resistance, and hence any rock or fluid properties derived from it, improves considerably. This will be illustrated next using a case study from the North Sea. March 2012 The Leading Edge 271

5 H o n o r a r y L e c t u r e Figure 7. (left) Well-log data from the chalk layer. Tracks show resistivity, density, porosity, lithology, and water saturation. (right) Histograms showing the distribution of resistivity within porous oilcharged chalk, porous water-wet chalk, and tight chalk. Determining the fluid in a North Sea chalk layer The second example comes from a chalk reservoir in the North Sea. The area of porous chalk forming the reservoir interval lies in the top m of the chalk layer, which itself varies between 100 and 500 m in thickness and lies approximately 2.5 km below sea floor. Although seismic data can be used to determine the structure and porosity of the chalk, they are extremely insensitive to the fluid type and distribution within it, because the impedance of porous water-saturated chalk is almost the same as that of porous oil-saturated chalk. Figure 7 shows a well log through the chalk layer, penetrating areas of porous oil-saturated, porous water-saturated, and tight chalk. Histograms illustrating the distribution of resistivity in each of these three zones are also shown. The average bulk resistivity of the porous chalk is considerably higher in the presence of hydrocarbons. However, resistivity alone remains ambiguous, because the distribution of resistivity within a tight chalk is extremely similar to that in a porous hydrocarbon-bearing chalk. There is a trade-off between porosity and saturation that cannot be resolved on the basis of resistivity measurements taken in isolation. Given both impedance and resistivity measurements (provided respectively by seismic and CSEM data), this ambiguity can be resolved. The seismic data can be used to identify porous areas of the chalk. Once the porosity is constrained, the CSEM data can be used to establish whether porous areas identified correspond to a high resistivity (potentially indicating hydrocarbon charge) or low resistivity. The CSEM survey layout is shown in Figure 8. The analysis presented here concentrates on line 2, a line trending NESW that passes from an area characterized by low-to-intermediate porosity chalk to an area of higher-porosity chalk at its northern end. Two source tows along this line were performed, the first at a fundamental frequency of 0.05 Hz and 272 The Leading Edge March 2012 Figure 8. Map showing the layout of the CSEM survey. Black circles represent receiver locations, and the source tow line is shown by the continuous black line, which also extended 10 km outside the receiver array (not shown). The colors represent the average porosity in the upper 24 m of the chalk layer, derived from calibrated inversion of seismic data. Note the salt diapir in the southwest corner of the survey area. the second at a fundamental frequency of 0.2 Hz (within each case, the odd harmonics of the fundamental also available for processing). The water depth in the survey area was m. The interpretation of CSEM data proceeds in stages, starting with simple 1D approaches and gradually adding complexity as required to explain the data. For this survey, the results of the 1D reconnaissance modeling stages indicated that the overburden structure (from sea floor down to the top of the chalk) was laterally continuous along the line. Within the chalk, however, the resistivity increases gradually to the north. The next stage in the interpretation is to apply a 2.5D inversion. Seismic structural constraints on the depths of the top and base of chalk were included in the CSEM inversion as breaks in smoothness. The resulting resistivity section is shown in Figure 9. There is a clear increase in resistivity toward the northern end of the chalk layer, confirming the results of the 1D modeling. The transverse resistance between the top and base chalk is also shown in Figure 9, and again increases to the northern end of the line. However, the cause of the increase in resistivity remains ambiguous when only the CSEM data are considered. It could be the result of a decrease in porosity of the chalk, a change in the chalk thickness or a change in the fluid content of the chalk from water to more resistive hydrocarbons. By incorporating the seismic data in the analysis, this ambiguity can be resolved. Seismic data can be used to constrain the porosity of the chalk layer. The left panel of Figure 10 shows the result of inversion of the seismic data along the CSEM line calibrated to porosity. The high-porosity zone in the upper few tens of meters of the chalk forms the reservoir interval, underlain by a much lower-porosity section. Applying Archie s law and assuming 100% water saturation allows us to convert the seismically derived porosity into an equivalent seismically derived resistivity (Figure 10).

6 March 2012 The Leading Edge 273

7 Figure 9. (bottom) 2D resistivity section resulting from inversion of 0.05 Hz, 0.25 Hz and 0.6 Hz amplitude and phase data from CSEM line 2. Structural constraints from seismic data are included as breaks in the smoothness constraint at the top and base of the chalk. The resistivity within the chalk increases toward the northern end (right side in the figure) of the line. (top) Transverse resistance between the top and base chalk calculated from the CSEM inversion result. This parameter is well constrained by CSEM data. Figure 10. (left) Porosity section within the chalk along line 2, derived from inversion of seismic data. (middle and right) The low-baseline and high-baseline (respectively) seismically derived resistivity, calculated from the porosity section using Archie s law calibrated at wells in the area, and assuming 100% saturation with water. Low-baseline case assumes Rw = 0.04 Ωm and the high-baseline case assumes Rw = 0.08 Ωm. The empirical Archie parameters are calibrated using well-log data from the survey area to give a value of the cementation parameter, m=1.6. The interstitial water resistivity, Rw, varies across the survey from values below 0.04 Ωm in areas close to the salt diapir where pore fluids are likely to be more saline, to values above 0.08 Ωm. Given the significant variation in Rw, low- and high-baseline resistivity sections were calculated (Figure 10, middle and right). The seismically derived transverse resistance in the chalk can then be calculated as before (Figure 11, blue curves). This reflects the effect of variations in porosity and chalk thickness. Any deviation from this curve is therefore indicative of a change in pore-fluid properties or type. This is also illustrated in Figure 11, which shows the equivalent transverse resistance when a synthetic hydrocarbon reservoir, consistent with horizontal logs through the known hydrocarbon accumulation in the area, is included in the model (red and pink curves for the low- and high-baseline curves, respectively). By comparing the CSEM-derived transverse resistance with the seismically derived values, the effects of fluids in the chalk can be isolated. This comparison is shown by black in Figure 11. At the southern end of the line, outside the area known to be charged with hydrocarbons, the CSEM result follows the low-baseline seismically derived curve closely. This area also corresponds to the portion of the line closest to the salt diapir, so a low Rw is plausible if the pore fluids have a higher salinity than elsewhere as a result. The well log from which the synthetic reservoir was constructed runs between 14 and 19 km along the line. In this area, the CSEM result follows the high-baseline seismically derived curve closely. Outside this area (from 19 to 25 km along the line), the large difference between the CSEM and seismically derived transverse resistance suggests an extension to the north of the hydrocarbon accumulation, either beneath or adjacent to the line. There are often ambiguities in the interpretation of a single geophysical data type. Here we have illustrated this with a case study from a North Sea chalk layer. Seismic data can be used to identify porous areas in the chalk, but cannot distinguish between water-wet and hydrocarbon-charged chalk. In contrast CSEM data can be used to identify areas of high resistivity; however, these may be caused either by low-porosity zones, or by porous zones that are charged with resistive hydrocarbons. By integrating the analysis of the two data types, the ambiguities in each can be mitigated, providing in this case an effective indication of the fluid content of the chalk. The future The ultimate goal of a geophysical analysis is to find a constrained model of geology, lithology, and fluid properties on which commercial decisions about reservoir exploration, development, and management can be based. To achieve this, the Earth can be interrogated with a number of tools. In this paper we have considered seismic, CSEM, and well-log data. However, gravity data, magnetic data, and magnetotellurics (for example) also have an important role to play (De Stefano et al., 2011; Jegen et al., 2009). Each data type must 274 The Leading Edge March 2012

8 March 2012 The Leading Edge 275

9 Figure 11. Comparison between the CSEM and seismically derived transverse resistance curves. Light blue and dark blue curves = high- and low-baseline seismically derived transverse resistance assuming 100% water saturation. Pink (high-baseline) and red (low-baseline) curves = same but with a synthetic hydrocarbon reservoir based on horizontal well-log information included. Black curve = CSEM-derived transverse resistance. In the southern (left) and central part of the line, the CSEM result agrees well with the seismic and well-log predictions. At the northern end of the line, where no well log is available, the large difference between CSEM and seismically derived transverse resistance suggests an extension of the resistive hydrocarbon layer to the north. be interpreted within an integrated framework so that the resulting Earth model is consistent with all the data used in its construction. Within this integrated framework, CSEM can provide valuable complementary information. However, there are a few considerations in applying this method to ensure that it delivers its full potential. The first, and perhaps most important, point is to ensure the geophysical question to which the CSEM method is applied is well posed and suitable for the technique. If the question is poorly posed, the answer will be equally poorly constrained, and interpretation risk of the result will be high. Good survey design is essential. Acquisition parameters such as required offsets and transmission frequencies are critical to the success of a CSEM survey. Survey parameters must be carefully targeted to answer the question posed above. If they are not, then the survey may fail to be sensitive to the structures of interest, leading to a disappointing result. In order to design the survey correctly, a good background model is required. CSEM survey parameters depend critically on factors such as overburden structure, resistivity, and anisotropy, and therefore constraint on these is critical to a well designed survey. The applicability of CSEM methods must be considered with this in mind. CSEM methods can in principle be applied at any stage during the life cycle of an oil field, although the question being addressed varies: Frontier exploration: What is the resistivity structure of the Earth? Exploration: What is the resistivity within the seismically defined structure? Appraisal: Given the seismic structure and the properties of the reservoir at a well, how do these properties vary across the field? Monitoring: Given the seismic structure, and properties in the well, how do these properties vary across the field and through time? The interpretation uncertainty and risk associated with CSEM declines as more information is available to construct an integrated interpretation. Therefore, of these possible applications, the last three, where both well and seismic data are available, will produce results with the lowest uncertainties. To date, the CSEM industry has been focused on exploration, with only a few surveys targeted at appraisal, and none directly at monitoring, for which the technology may be much better suited. Perhaps in the future, the most widely applied use of CSEM methods will shift toward these applications. Indeed several recent studies (for example Orange et al., 2009; Andreis et al., 2011; Liang et al., 2011) have demonstrated the potential power of CSEM methods as a complement to seismic in monitoring applications. Conclusions There are many geophysical questions that cannot be answered unambiguously (or at all) by a single data type. The examples presented here illustrate situations where integration of seismic, well and CSEM data, exploiting the strengths of each, is required. Many other examples, incorporating a variety of data types, exist. The field of data integration and shared Earth modeling is large, and rapidly expanding with technology moving from qualitative to quantitative approaches. There are many challenges, for example combining 276 The Leading Edge March 2012

10 different physical properties, using technologies that sense the Earth at a range of different scales is not straightforward. However, there are also many opportunities, if these challenges can be overcome, to improve the quality and constraints of geophysical information on which commercial decisions are based. References Andreis, D. and L. MacGregor, 2011, Using CSEM to monitor production from a complex 3D gas reservoir: a synthetic case study: The Leading Edge, 30, no. 9, , doi: / Archie, G. E., 1942, The electrical resistivity log as an aid in determining some reservoir characteristics: Transactions of the American Institute of Mining Engineers, 146, Carcione, J. M., B. Ursin, and J. I. Nordskag, 2007, Cross-property relations between electrical conductivity and the seismic velocity of rocks: Geophysics, 72, no. 5, E193 E204, doi: / Constable, S. and L. J. Srnka, 2007, An introduction to marine controlled-source electromagnetic methods for hydrocarbon exploration: Geophysics, 72, no. 2, WA3 WA12, doi: / Constable, S. C., R. L. Parker, and C. G. Constable, 1987, Occam s inversion: A practical algorithm for generating smooth models from EM sounding data: Geophysics, 52, no. 3, , doi: / De Stefano, M., F. G. Andreasi, S. Re, M. Virgilio, and F. Snyder, 2011, Multiple domain simultaneous joint inversion of geophysical data with application to subsalt imaging: Geophysics, 76, no. 3, R69 R80, doi: / Gelius, L. J. and Z. Wang, 2008, Modeling production caused changes in conductivity for a siliciclastic reservoir: A differential effective medium approach: Geophysical Prospecting, 56, no. 5, , doi: /j x. Han, T. C., A. I. Best, L. M. MacGregor, J. Sothcott, and T. A. Minshull, 2011, Joint elastic electrical effective medium models of reservoir sandstones: Geophysical Prospecting, 59, no. 4, , doi: /j x. Harris, P., Z. Du, L. MacGregor, W. Olsen, R. Shu, and R. Cooper, 2009, Joint interpretation of seismic and CSEM data using well log constraints: an example from the Luva Field: First Break, 27, Hashin, Z. and S. Shtrikman, 1962, A variational approach to the theory of the effective magnetic permeability of multiphase materials: Journal of Applied Physics, 33, no. 10, , doi: / Jegen, M. D., R. W. Hobbs, P. Tarits, and A. Chave, 2009, Joint inversion of marine magnetotelluric and gravity data incorporating seismic constraints. Preliminary results of sub-basalt imaging off the Faroe Shelf: Earth and Planetary Science Letters, 282, no. 1 4, 47 55, doi: /j.epsl Liang, L., A. Abubakar, and T. Habashy, 2011, Production monitoring using joint inversion of marine CSEM data and production data: 81st Annual International Meeting, SEG, Expanded Abstracts, Orange, A., K. Key, and S. Constable, 2009, The feasibility of reservoir monitoring using time-lapse CSEM: Geophysics, 74, no. 2, F21 F29, doi: / Waxman, M. H. and L. J. M. Smits, 1968, Electrical conductivities in oil-bearing shaly sands: Society of Petroleum Engineers Journal, 8, no. 2, , doi: /1863-a. especially David Andreis, Paola Vera de Newton, Uwe Strecker, and Anyela Morecote who all contributed to the work presented. The North Sea case study data were acquired by OHM Ltd., and thanks to Maersk, Chevron, and Shell for permission to publish the case study. Many of the workflows described were developed under the WISE consortium project, sponsored by Maersk, Chevron, Det Norske, Total, DONG, and the UK s Department for Energy and Climate Change (DECC). The time-lapse modeling included in the lecture was performed as part of a research project in collaboration with BP, and funded by the UK Technology Strategy Board. Finally, thanks to SEG for organizing the Honorary Lecture program, and to Shell for its generous sponsorship. Corresponding author: lucy.macgregor@rocksolidimages.com Acknowledgments: I thank my colleagues at RSI for their help, March 2012 The Leading Edge 277

2011 SEG SEG San Antonio 2011 Annual Meeting 771. Summary. Method

2011 SEG SEG San Antonio 2011 Annual Meeting 771. Summary. Method Geological Parameters Effecting Controlled-Source Electromagnetic Feasibility: A North Sea Sand Reservoir Example Michelle Ellis and Robert Keirstead, RSI Summary Seismic and electromagnetic data measure

More information

Rock physics integration of CSEM and seismic data: a case study based on the Luva gas field.

Rock physics integration of CSEM and seismic data: a case study based on the Luva gas field. Rock physics integration of CSEM and seismic data: a case study based on the Luva gas field. Peter Harris*, Zhijun Du, Harald H. Soleng, Lucy M. MacGregor, Wiebke Olsen, OHM-Rock Solid Images Summary It

More information

Reducing Uncertainty through Multi-Measurement Integration: from Regional to Reservoir scale

Reducing Uncertainty through Multi-Measurement Integration: from Regional to Reservoir scale Reducing Uncertainty through Multi-Measurement Integration: from Regional to Reservoir scale Efthimios Tartaras Data Processing & Modeling Manager, Integrated Electromagnetics CoE, Schlumberger Geosolutions

More information

Tu B3 15 Multi-physics Characterisation of Reservoir Prospects in the Hoop Area of the Barents Sea

Tu B3 15 Multi-physics Characterisation of Reservoir Prospects in the Hoop Area of the Barents Sea Tu B3 15 Multi-physics Characterisation of Reservoir Prospects in the Hoop Area of the Barents Sea P. Alvarez (RSI), F. Marcy (ENGIE E&P), M. Vrijlandt (ENGIE E&P), K. Nichols (RSI), F. Bolivar (RSI),

More information

Summary. Introduction

Summary. Introduction Noel Black*, Glenn A. Wilson, TechnoImaging, Alexander V. Gribenko and Michael S. Zhdanov, TechnoImaging and The University of Utah Summary Recent studies have inferred the feasibility of time-lapse controlled-source

More information

Vertical and horizontal resolution considerations for a joint 3D CSEM and MT inversion

Vertical and horizontal resolution considerations for a joint 3D CSEM and MT inversion Antony PRICE*, Total E&P and Don WATTS, WesternGeco Electromagnetics Summary To further explore the potential data content and inherent limitations of a detailed 3D Controlled Source ElectroMagnetic and

More information

Determination of reservoir properties from the integration of CSEM and seismic data

Determination of reservoir properties from the integration of CSEM and seismic data Determination of reservoir properties from the integration of CSEM and seismic data Peter Harris, 1 Rock Solid Images, and Lucy MacGregor, 2 Offshore Hydrocarbons Mapping, discuss the advantages in reservoir

More information

Using seismic guided EM inversion to explore a complex geological area: An application to the Kraken and Bressay heavy oil discoveries, North Sea

Using seismic guided EM inversion to explore a complex geological area: An application to the Kraken and Bressay heavy oil discoveries, North Sea Using seismic guided EM inversion to explore a complex geological area: An application to the Kraken and Bressay heavy oil discoveries, North Sea Zhijun Du*, PGS Summary The integrated analysis of controlled

More information

Improved Exploration, Appraisal and Production Monitoring with Multi-Transient EM Solutions

Improved Exploration, Appraisal and Production Monitoring with Multi-Transient EM Solutions Improved Exploration, Appraisal and Production Monitoring with Multi-Transient EM Solutions Folke Engelmark* PGS Multi-Transient EM, Asia-Pacific, Singapore folke.engelmark@pgs.com Summary Successful as

More information

Anisotropic 2.5D Inversion of Towed Streamer EM Data from Three North Sea Fields Using Parallel Adaptive Finite Elements

Anisotropic 2.5D Inversion of Towed Streamer EM Data from Three North Sea Fields Using Parallel Adaptive Finite Elements Anisotropic 2.5D Inversion of Towed Streamer EM Data from Three North Sea Fields Using Parallel Adaptive Finite Elements K. Key (Scripps Institution of Oceanography), Z. Du* (PGS), J. Mattsson (PGS), A.

More information

Downloaded 08/29/13 to Redistribution subject to SEG license or copyright; see Terms of Use at

Downloaded 08/29/13 to Redistribution subject to SEG license or copyright; see Terms of Use at New approach to 3D inversion of MCSEM and MMT data using multinary model transform Alexander V. Gribenko and Michael S. Zhdanov, University of Utah and TechnoImaging SUMMARY Marine controlled-source electromagnetic

More information

A Broadband marine CSEM demonstration survey to map the Uranus salt structure

A Broadband marine CSEM demonstration survey to map the Uranus salt structure A Broadband marine CSEM demonstration survey to map the Uranus salt structure M. Vöge ( 1 ), A. A. Pfaffhuber ( 1 ), K. Hokstad ( 2 ), B. Fotland( 2 ) ( 1 ) Norwegian Geotechnical Institute, Sognsveien

More information

The elastic properties such as velocity, density, impedance,

The elastic properties such as velocity, density, impedance, SPECIAL SECTION: Rr ock Physics physics Lithology and fluid differentiation using rock physics template XIN-GANG CHI AND DE-HUA HAN, University of Houston The elastic properties such as velocity, density,

More information

Uncertainty analysis for the integration of seismic and CSEM data Myoung Jae Kwon & Roel Snieder, Center for Wave Phenomena, Colorado School of Mines

Uncertainty analysis for the integration of seismic and CSEM data Myoung Jae Kwon & Roel Snieder, Center for Wave Phenomena, Colorado School of Mines Myoung Jae Kwon & Roel Snieder, Center for Wave Phenomena, Colorado School of Mines Summary Geophysical inverse problems consist of three stages: the forward problem, optimization, and appraisal. We study

More information

Towed Streamer EM data from Barents Sea, Norway

Towed Streamer EM data from Barents Sea, Norway Towed Streamer EM data from Barents Sea, Norway Anwar Bhuiyan*, Eivind Vesterås and Allan McKay, PGS Summary The measured Towed Streamer EM data from a survey in the Barents Sea, undertaken in the Norwegian

More information

We Simultaneous Joint Inversion of Electromagnetic and Seismic Full-waveform Data - A Sensitivity Analysis to Biot Parameter

We Simultaneous Joint Inversion of Electromagnetic and Seismic Full-waveform Data - A Sensitivity Analysis to Biot Parameter We-09-04 Simultaneous Joint Inversion of Electromagnetic and Seismic Full-waveform Data - A Sensitivity Analysis to Biot Parameter J. Giraud* (WesternGeco Geosolutions), M. De Stefano (WesternGeco Geosolutions)

More information

Downloaded 10/29/15 to Redistribution subject to SEG license or copyright; see Terms of Use at

Downloaded 10/29/15 to Redistribution subject to SEG license or copyright; see Terms of Use at Electrical anisotropy drivers in the Snøhvit region of the Barents Sea Michelle Ellis*, Lucy MacGregor, Rolf Ackermann, Paola Newton, Robert Keirstead, Alberto Rusic, Slim Bouchrara, Amanda Geck Alvarez,

More information

Downloaded 10/29/15 to Redistribution subject to SEG license or copyright; see Terms of Use at

Downloaded 10/29/15 to Redistribution subject to SEG license or copyright; see Terms of Use at CSEM based Anisotropy trends in the Barents Sea. Slim Bouchrara *, Lucy MacGregor, Amanda Alvarez, Michelle Ellis, Rolf Ackermann, Paola Newton, Robert Keirstead, Alburto Rusic,Yijie Zhou and Hung-Wen

More information

Estimating the hydrocarbon volume from elastic and resistivity data: A concept

Estimating the hydrocarbon volume from elastic and resistivity data: A concept INTERPRETER S CORNER Coordinated by Rebecca B. Latimer Estimating the hydrocarbon volume from elastic and resistivity data: A concept CARMEN T. GOMEZ, JACK DVORKIN, and GARY MAVKO, Stanford University,

More information

23855 Rock Physics Constraints on Seismic Inversion

23855 Rock Physics Constraints on Seismic Inversion 23855 Rock Physics Constraints on Seismic Inversion M. Sams* (Ikon Science Ltd) & D. Saussus (Ikon Science) SUMMARY Seismic data are bandlimited, offset limited and noisy. Consequently interpretation of

More information

SEG Houston 2009 International Exposition and Annual Meeting

SEG Houston 2009 International Exposition and Annual Meeting The role of EM rock physics and seismic data in integrated 3D CSEM data analysis I. Brevik*, StatoilHydro, Pål T. Gabrielsen, Vestfonna and Jan Petter Morten, EMGS Summary An extensive 3D CSEM dataset

More information

Geophysical model response in a shale gas

Geophysical model response in a shale gas Geophysical model response in a shale gas Dhananjay Kumar and G. Michael Hoversten Chevron USA Inc. Abstract Shale gas is an important asset now. The production from unconventional reservoir like shale

More information

G008 Advancing Marine Controlled Source Electromagnetics in the Santos Basin, Brazil

G008 Advancing Marine Controlled Source Electromagnetics in the Santos Basin, Brazil G008 Advancing Marine Controlled Source Electromagnetics in the Santos Basin, Brazil M.P. Buonora* (Petrobras), A. Zerilli (WesternGeco Electromagnetics), T. Labruzzo (WesternGeco Electromagnetics) & L.F.

More information

Estimating vertical and horizontal resistivity of the overburden and the reservoir for the Alvheim Boa field. Folke Engelmark* and Johan Mattsson, PGS

Estimating vertical and horizontal resistivity of the overburden and the reservoir for the Alvheim Boa field. Folke Engelmark* and Johan Mattsson, PGS Estimating vertical and horizontal resistivity of the overburden and the reservoir for the Alvheim Boa field. Folke Engelmark* and Johan Mattsson, PGS Summary Towed streamer EM data was acquired in October

More information

Downloaded 03/06/15 to Redistribution subject to SEG license or copyright; see Terms of Use at

Downloaded 03/06/15 to Redistribution subject to SEG license or copyright; see Terms of Use at Integration of seismic, CSEM and well data in the Perdido basin, Gulf of Mexico Philippe Nivlet, Luis Sánchez Pérez, EMGS Services México, José Antonio Escalera Alcocer, José Trinidad Martínez Vázquez,

More information

Downloaded 09/09/15 to Redistribution subject to SEG license or copyright; see Terms of Use at

Downloaded 09/09/15 to Redistribution subject to SEG license or copyright; see Terms of Use at Reservoir properties estimation from marine broadband seismic without a-priori well information: A powerful de-risking workflow Cyrille Reiser*, Matt Whaley and Tim Bird, PGS Reservoir Limited Summary

More information

The prediction of reservoir

The prediction of reservoir Risk Reduction in Gas Reservoir Exploration Using Joint Seismic-EM Inversion GAS EXPLORATION By Yoram Rubin, G. Michael Hoversten, Zhangshuan Hou and Jinsong Chen, University of California, Berkeley A

More information

Use of Seismic and EM Data for Exploration, Appraisal and Reservoir Characterization

Use of Seismic and EM Data for Exploration, Appraisal and Reservoir Characterization Use of Seismic and EM Data for Exploration, Appraisal and Reservoir Characterization Anton Ziolkowski and Folke Engelmark Petroleum Geo-Services CSEG, Calgary, 6 May 2009 Outline Exploration, appraisal,

More information

Interpretation and Reservoir Properties Estimation Using Dual-Sensor Streamer Seismic Without the Use of Well

Interpretation and Reservoir Properties Estimation Using Dual-Sensor Streamer Seismic Without the Use of Well Interpretation and Reservoir Properties Estimation Using Dual-Sensor Streamer Seismic Without the Use of Well C. Reiser (Petroleum Geo-Services), T. Bird* (Petroleum Geo-Services) & M. Whaley (Petroleum

More information

Integrating rock physics and full elastic modeling for reservoir characterization Mosab Nasser and John B. Sinton*, Maersk Oil Houston Inc.

Integrating rock physics and full elastic modeling for reservoir characterization Mosab Nasser and John B. Sinton*, Maersk Oil Houston Inc. Integrating rock physics and full elastic modeling for reservoir characterization Mosab Nasser and John B. Sinton*, Maersk Oil Houston Inc. Summary Rock physics establishes the link between reservoir properties,

More information

Modeling and interpretation of CSEM data from Bressay, Bentley and Kraken area of East Shetland Platform, North Sea

Modeling and interpretation of CSEM data from Bressay, Bentley and Kraken area of East Shetland Platform, North Sea from Bressay, Bentley and Kraken area of East Shetland Platform, North Sea Anwar Bhuiyan*, Rune Sakariassen, Øystein Hallanger and Allan McKay, PGS Summary The heavy oil reservoirs of the Bentley, Bressay

More information

Feasibility and design study of a multicomponent seismic survey: Upper Assam Basin

Feasibility and design study of a multicomponent seismic survey: Upper Assam Basin P-276 Summary Feasibility and design study of a multicomponent seismic survey: Upper Assam Basin K.L.Mandal*, R.K.Srivastava, S.Saha, Oil India Limited M.K.Sukla, Indian Institute of Technology, Kharagpur

More information

Shaly Sand Rock Physics Analysis and Seismic Inversion Implication

Shaly Sand Rock Physics Analysis and Seismic Inversion Implication Shaly Sand Rock Physics Analysis and Seismic Inversion Implication Adi Widyantoro (IkonScience), Matthew Saul (IkonScience/UWA) Rock physics analysis of reservoir elastic properties often assumes homogeneity

More information

We apply a rock physics analysis to well log data from the North-East Gulf of Mexico

We apply a rock physics analysis to well log data from the North-East Gulf of Mexico Rock Physics for Fluid and Porosity Mapping in NE GoM JACK DVORKIN, Stanford University and Rock Solid Images TIM FASNACHT, Anadarko Petroleum Corporation RICHARD UDEN, MAGGIE SMITH, NAUM DERZHI, AND JOEL

More information

Reservoir properties inversion from AVO attributes

Reservoir properties inversion from AVO attributes Reservoir properties inversion from AVO attributes Xin-gang Chi* and De-hua Han, University of Houston Summary A new rock physics model based inversion method is put forward where the shaly-sand mixture

More information

Downloaded 10/02/18 to Redistribution subject to SEG license or copyright; see Terms of Use at

Downloaded 10/02/18 to Redistribution subject to SEG license or copyright; see Terms of Use at Multi-scenario, multi-realization seismic inversion for probabilistic seismic reservoir characterization Kester Waters* and Michael Kemper, Ikon Science Ltd. Summary We propose a two tiered inversion strategy

More information

Detection, Delineation and Characterization of Shallow Anomalies Using Dual Sensor Seismic and Towed Streamer EM data

Detection, Delineation and Characterization of Shallow Anomalies Using Dual Sensor Seismic and Towed Streamer EM data Detection, Delineation and Characterization of Shallow Anomalies Using Dual Sensor Seismic and Towed Streamer EM data A.J. McKay* (Petroleum Geo-Services ASA), M. Widmaier (Petroleum Geo- Services ASA),

More information

Towed Streamer EM Integrated interpretation for accurate characterization of the sub-surface. PETEX, Tuesday 15th of November 2016

Towed Streamer EM Integrated interpretation for accurate characterization of the sub-surface. PETEX, Tuesday 15th of November 2016 Towed Streamer EM Integrated interpretation for accurate characterization of the sub-surface. PETEX, Tuesday 15th of November 2016 Joshua May Sales and Marketing Manager Agenda Introduction to Towed Streamer

More information

B033 Improving Subsalt Imaging by Incorporating MT Data in a 3D Earth Model Building Workflow - A Case Study in Gulf of Mexico

B033 Improving Subsalt Imaging by Incorporating MT Data in a 3D Earth Model Building Workflow - A Case Study in Gulf of Mexico B033 Improving Subsalt Imaging by Incorporating MT Data in a 3D Earth Model Building Workflow - A Case Study in Gulf of Mexico E. Medina* (WesternGeco), A. Lovatini (WesternGeco), F. Golfré Andreasi (WesternGeco),

More information

A014 Uncertainty Analysis of Velocity to Resistivity Transforms for Near Field Exploration

A014 Uncertainty Analysis of Velocity to Resistivity Transforms for Near Field Exploration A014 Uncertainty Analysis of Velocity to Resistiity Transforms for Near Field Exploration D. Werthmüller* (Uniersity of Edinburgh), A.M. Ziolkowski (Uniersity of Edinburgh) & D.A. Wright (Uniersity of

More information

Rock physics and AVO applications in gas hydrate exploration

Rock physics and AVO applications in gas hydrate exploration Rock physics and AVO applications in gas hydrate exploration ABSTRACT Yong Xu*, Satinder Chopra Core Lab Reservoir Technologies Division, 301,400-3rd Ave SW, Calgary, AB, T2P 4H2 yxu@corelab.ca Summary

More information

Evaluation of Rock Properties from Logs Affected by Deep Invasion A Case Study

Evaluation of Rock Properties from Logs Affected by Deep Invasion A Case Study Evaluation of Rock Properties from Logs Affected by Deep Invasion A Case Study Jahan Zeb a, Reece Murrell b a CGG Australia, 1 Ord Street, West Perth, WA 6005 Contact email: Jahan.Zeb@cgg.com b Esso Australia,

More information

Integration of seismic and fluid-flow data: a two-way road linked by rock physics

Integration of seismic and fluid-flow data: a two-way road linked by rock physics Integration of seismic and fluid-flow data: a two-way road linked by rock physics Abstract Yunyue (Elita) Li, Yi Shen, and Peter K. Kang Geologic model building of the subsurface is a complicated and lengthy

More information

Quantitative Interpretation

Quantitative Interpretation Quantitative Interpretation The aim of quantitative interpretation (QI) is, through the use of amplitude analysis, to predict lithology and fluid content away from the well bore. This process should make

More information

Exploring with Controlled Source Electro- Magnetic (CSEM) methods: from 2D profiling to 3D multi-azimuth surveying

Exploring with Controlled Source Electro- Magnetic (CSEM) methods: from 2D profiling to 3D multi-azimuth surveying Exploring with Controlled Source Electro- Magnetic (CSEM) methods: from 2D profiling to 3D multi-azimuth surveying M. Darnet, P. Van Der Sman, R.E. Plessix, J.L. Johnson, M. Rosenquist ( 1 ) ( 1 ) Shell

More information

Reservoir properties prediction using CSEM, pre-stack seismic and well log data: Case Study in the Hoop Area, Barents Sea, Norway

Reservoir properties prediction using CSEM, pre-stack seismic and well log data: Case Study in the Hoop Area, Barents Sea, Norway Reservoir properties prediction using CSEM, pre-stack seismic and well log data: Case Study in the Hoop Area, Barents Sea, Norway ABSTRACT We present an example from the Hoop area of the Barents Sea showing

More information

SEG Houston 2009 International Exposition and Annual Meeting

SEG Houston 2009 International Exposition and Annual Meeting : The Santos Basin Project - Brazil Andrea Zerilli * and Tiziano Labruzzo, WesternGeco EM, Marco Polo Buonora, Paulo de Tarso Luiz Menezes and Luiz Felipe Rodrigues, Petrobras E&P/GEOF/MP Summary The Santos

More information

Statistical Rock Physics

Statistical Rock Physics Statistical - Introduction Book review 3.1-3.3 Min Sun March. 13, 2009 Outline. What is Statistical. Why we need Statistical. How Statistical works Statistical Rock physics Information theory Statistics

More information

An empirical method for estimation of anisotropic parameters in clastic rocks

An empirical method for estimation of anisotropic parameters in clastic rocks An empirical method for estimation of anisotropic parameters in clastic rocks YONGYI LI, Paradigm Geophysical, Calgary, Alberta, Canada Clastic sediments, particularly shale, exhibit transverse isotropic

More information

Downloaded 11/02/16 to Redistribution subject to SEG license or copyright; see Terms of Use at Summary.

Downloaded 11/02/16 to Redistribution subject to SEG license or copyright; see Terms of Use at   Summary. in thin sand reservoirs William Marin* and Paola Vera de Newton, Rock Solid Images, and Mario Di Luca, Pacific Exploración y Producción. Summary Rock Physics Templates (RPTs) are useful tools for well

More information

Summary. Introduction

Summary. Introduction Effect of over- and under-burden on time-lapse CSEM monitoring capabilities Arash JafarGandomi* and Andrew Curtis, School of GeoSciences, The University of Edinburgh, Kings Buildings, Edinburgh EH9 3JW,

More information

Geological Classification of Seismic-Inversion Data in the Doba Basin of Chad*

Geological Classification of Seismic-Inversion Data in the Doba Basin of Chad* Geological Classification of Seismic-Inversion Data in the Doba Basin of Chad* Carl Reine 1, Chris Szelewski 2, and Chaminda Sandanayake 3 Search and Discovery Article #41899 (2016)** Posted September

More information

Downloaded 09/16/16 to Redistribution subject to SEG license or copyright; see Terms of Use at

Downloaded 09/16/16 to Redistribution subject to SEG license or copyright; see Terms of Use at Ehsan Zabihi Naeini*, Ikon Science & Russell Exley, Summit Exploration & Production Ltd Summary Quantitative interpretation (QI) is an important part of successful Central North Sea exploration, appraisal

More information

RC 1.3. SEG/Houston 2005 Annual Meeting 1307

RC 1.3. SEG/Houston 2005 Annual Meeting 1307 from seismic AVO Xin-Gong Li,University of Houston and IntSeis Inc, De-Hua Han, and Jiajin Liu, University of Houston Donn McGuire, Anadarko Petroleum Corp Summary A new inversion method is tested to directly

More information

SEAM2: 3D NON-SEISMIC MODELING OF A COMPLEX MIDDLE EAST O&G PROSPECT

SEAM2: 3D NON-SEISMIC MODELING OF A COMPLEX MIDDLE EAST O&G PROSPECT 6th International Symposium on Three-Dimensional Electromagnetics Berkeley, California, USA, March 28-3, 217 SEAM2: 3D NON-SEISMIC MODELING OF A COMPLEX MIDDLE EAST O&G PROSPECT M. D. Watts 1, R. L. Mackie

More information

Keywords. CSEM, Inversion, Resistivity, Kutei Basin, Makassar Strait

Keywords. CSEM, Inversion, Resistivity, Kutei Basin, Makassar Strait Noor Jehan Ashaari Muhamad* (EMGS Asia Pacific), Ritesh Mohan Joshi (EMGS Asia Pacific), Muhamad Afifie Chan Mahadie Chan (EMGS Asia Pacific) mmuhamad@emgs.com Keywords CSEM, Inversion, Resistivity, Kutei

More information

6298 Stress induced azimuthally anisotropic reservoir - AVO modeling

6298 Stress induced azimuthally anisotropic reservoir - AVO modeling 6298 Stress induced azimuthally anisotropic reservoir - AVO modeling M. Brajanovski* (Curtin University of Technology), B. Gurevich (Curtin University of Technology), D. Nadri (CSIRO) & M. Urosevic (Curtin

More information

A11 Planning Time-lapse CSEM-surveys for Joint Seismic-EM Monitoring of Geological Carbon Dioxide Injection

A11 Planning Time-lapse CSEM-surveys for Joint Seismic-EM Monitoring of Geological Carbon Dioxide Injection A11 Planning Time-lapse CSEM-surveys for Joint Seismic-EM Monitoring of Geological Carbon Dioxide Injection T. Norman* (NTNU Trondheim and EMGS ASA), H. Alnes (StatoilHydro Research Trondheim), O. Christensen

More information

ROCK PHYSICS MODELING FOR LITHOLOGY PREDICTION USING HERTZ- MINDLIN THEORY

ROCK PHYSICS MODELING FOR LITHOLOGY PREDICTION USING HERTZ- MINDLIN THEORY ROCK PHYSICS MODELING FOR LITHOLOGY PREDICTION USING HERTZ- MINDLIN THEORY Ida Ayu PURNAMASARI*, Hilfan KHAIRY, Abdelazis Lotfy ABDELDAYEM Geoscience and Petroleum Engineering Department Universiti Teknologi

More information

water L v i Chapter 4 Saturation

water L v i Chapter 4 Saturation 4. Resistivity The presence of hydrocarbons is identified by the electrical resistance of the formation. These electrical properties of rocks depend on the pore geometry and fluid distribution. That is,

More information

Sub-basalt exploration in the Kutch-Saurashtra basin using EM

Sub-basalt exploration in the Kutch-Saurashtra basin using EM Sub-basalt exploration in the Kutch-Saurashtra basin using EM Paper ID: 2003081 Deepankar Borgohain 1, Krishna Kumar 2, U G Marathe 3, Pradipta Mishra 3, Deepak Kumar 1 1.EMGS Asia Pacific Sdn Bhd, 1009

More information

Enhancing the resolution of CSEM inversion using seismic constraints Peter Harris*, Rock Solid Images AS Lucy MacGregor, OHM Surveys Ltd

Enhancing the resolution of CSEM inversion using seismic constraints Peter Harris*, Rock Solid Images AS Lucy MacGregor, OHM Surveys Ltd Enhancing the resolution of CSEM inversion using seismic constraints Peter Harris*, Rock Solid Images AS Lucy MacGregor, OHM Surveys Ltd Summary The vertical resolution of inverted CSEM data remains an

More information

Rock physics and AVO analysis for lithofacies and pore fluid prediction in a North Sea oil field

Rock physics and AVO analysis for lithofacies and pore fluid prediction in a North Sea oil field Rock physics and AVO analysis for lithofacies and pore fluid prediction in a North Sea oil field Downloaded 09/12/14 to 84.215.159.82. Redistribution subject to SEG license or copyright; see Terms of Use

More information

QUANTITATIVE INTERPRETATION

QUANTITATIVE INTERPRETATION QUANTITATIVE INTERPRETATION THE AIM OF QUANTITATIVE INTERPRETATION (QI) IS, THROUGH THE USE OF AMPLITUDE ANALYSIS, TO PREDICT LITHOLOGY AND FLUID CONTENT AWAY FROM THE WELL BORE This process should make

More information

Keywords. PMR, Reservoir Characterization, EEI, LR

Keywords. PMR, Reservoir Characterization, EEI, LR Enhancing the Reservoir Characterization Experience through Post Migration Reprocessed (PMR) Data A case study Indrajit Das*, Ashish Kumar Singh, Shakuntala Mangal, Reliance Industries Limited, Mumbai

More information

Determination of the Laminar, Structural and Disperse Shale Volumes Using a Joint Inversion of Conventional Logs*

Determination of the Laminar, Structural and Disperse Shale Volumes Using a Joint Inversion of Conventional Logs* Determination of the Laminar, Structural and Disperse Shale Volumes Using a Joint Inversion of Conventional Logs* Ambrosio Aquino-López 1, Aleksandr Mousatov 1, Mikhail Markov 1, and Elena Kazatchenko

More information

New Frontier Advanced Multiclient Data Offshore Uruguay. Advanced data interpretation to empower your decision making in the upcoming bid round

New Frontier Advanced Multiclient Data Offshore Uruguay. Advanced data interpretation to empower your decision making in the upcoming bid round New Frontier Advanced Multiclient Data Offshore Uruguay Advanced data interpretation to empower your decision making in the upcoming bid round Multiclient data interpretation provides key deliverables

More information

B-31 Combining 4D seismic and reservoir

B-31 Combining 4D seismic and reservoir B-31 Combining 4D seismic and reservoir simulation: key to effective reservoir management ANDREY BAKULIN 1, NICK DRINKWATER 1, CLAUDE SIGNER 1, SARAH RYAN 2 & ANDY O DONOVAN 3 1 Schlumberger Cambridge

More information

The GIG consortium Geophysical Inversion to Geology Per Røe, Ragnar Hauge, Petter Abrahamsen FORCE, Stavanger

The GIG consortium Geophysical Inversion to Geology Per Røe, Ragnar Hauge, Petter Abrahamsen FORCE, Stavanger www.nr.no The GIG consortium Geophysical Inversion to Geology Per Røe, Ragnar Hauge, Petter Abrahamsen FORCE, Stavanger 17. November 2016 Consortium goals Better estimation of reservoir parameters from

More information

Downloaded 11/20/12 to Redistribution subject to SEG license or copyright; see Terms of Use at

Downloaded 11/20/12 to Redistribution subject to SEG license or copyright; see Terms of Use at AVO crossplot analysis in unconsolidated sediments containing gas hydrate and free gas: Green Canyon 955, Gulf of Mexico Zijian Zhang* 1, Daniel R. McConnell 1, De-hua Han 2 1 Fugro GeoConsulting, Inc.,

More information

Novel approach to joint 3D inversion of EM and potential field data using Gramian constraints

Novel approach to joint 3D inversion of EM and potential field data using Gramian constraints first break volume 34, April 2016 special topic Novel approach to joint 3D inversion of EM and potential field data using Gramian constraints Michael S. Zhdanov 1,2, Yue Zhu 2, Masashi Endo 1 and Yuri

More information

ERTH2020 Introduction to Geophysics The Seismic Method. 1. Basic Concepts in Seismology. 1.1 Seismic Wave Types

ERTH2020 Introduction to Geophysics The Seismic Method. 1. Basic Concepts in Seismology. 1.1 Seismic Wave Types ERTH2020 Introduction to Geophysics The Seismic Method 1. Basic Concepts in Seismology 1.1 Seismic Wave Types Existence of different wave types The existence of different seismic wave types can be understood

More information

Multiple horizons mapping: A better approach for maximizing the value of seismic data

Multiple horizons mapping: A better approach for maximizing the value of seismic data Multiple horizons mapping: A better approach for maximizing the value of seismic data Das Ujjal Kumar *, SG(S) ONGC Ltd., New Delhi, Deputed in Ministry of Petroleum and Natural Gas, Govt. of India Email:

More information

National Geophysical Survey Science Opportunities

National Geophysical Survey Science Opportunities National Geophysical Survey Science Opportunities Anton Ziolkowski Professor of Petroleum Geoscience University of Edinburgh BGS Keyworth 4 April 2016 Outline Brief biography Science opportunities High

More information

PETROLEUM GEOSCIENCES GEOLOGY OR GEOPHYSICS MAJOR

PETROLEUM GEOSCIENCES GEOLOGY OR GEOPHYSICS MAJOR PETROLEUM GEOSCIENCES GEOLOGY OR GEOPHYSICS MAJOR APPLIED GRADUATE STUDIES Geology Geophysics GEO1 Introduction to the petroleum geosciences GEO2 Seismic methods GEO3 Multi-scale geological analysis GEO4

More information

The role of seismic modeling in Reservoir characterization: A case study from Crestal part of South Mumbai High field

The role of seismic modeling in Reservoir characterization: A case study from Crestal part of South Mumbai High field P-305 The role of seismic modeling in Reservoir characterization: A case study from Crestal part of South Mumbai High field Summary V B Singh*, Mahendra Pratap, ONGC The objective of the modeling was to

More information

TIV Contrast Source Inversion of mcsem data

TIV Contrast Source Inversion of mcsem data TIV Contrast ource Inversion of mcem data T. Wiik ( ), L. O. Løseth ( ), B. Ursin ( ), K. Hokstad (, ) ( )epartment for Petroleum Technology and Applied Geophysics, Norwegian University of cience and Technology

More information

D034 On the Effects of Anisotropy in Marine CSEM

D034 On the Effects of Anisotropy in Marine CSEM D034 On the Effects of Anisotropy in Marine CSEM L.O. Løseth* (Statoil ASA, formerly Dept. of Physics, NTNU), B. Ursin (Dept. of Appl. Geophys., NTNU) & L. Amundsen (Statoil ASA) SUMMARY In marine CSEM

More information

Joint inversion of borehole electromagnetic and sonic measurements G. Gao, A. Abubakar, T. M. Habashy, Schlumberger-Doll Research

Joint inversion of borehole electromagnetic and sonic measurements G. Gao, A. Abubakar, T. M. Habashy, Schlumberger-Doll Research Joint inversion of borehole electromagnetic and sonic measurements G. Gao, A. Abubakar, T. M. Habashy, Schlumberger-Doll Research SUMMARY New-generation electromagnetic (EM) and sonic logging tools are

More information

Seismic Guided Drilling: Near Real Time 3D Updating of Subsurface Images and Pore Pressure Model

Seismic Guided Drilling: Near Real Time 3D Updating of Subsurface Images and Pore Pressure Model IPTC 16575 Seismic Guided Drilling: Near Real Time 3D Updating of Subsurface Images and Pore Pressure Model Chuck Peng, John Dai and Sherman Yang, Schlumberger WesternGeco Copyright 2013, International

More information

The SPE Foundation through member donations and a contribution from Offshore Europe

The SPE Foundation through member donations and a contribution from Offshore Europe Primary funding is provided by The SPE Foundation through member donations and a contribution from Offshore Europe The Society is grateful to those companies that allow their professionals to serve as

More information

Simultaneous Inversion of Clastic Zubair Reservoir: Case Study from Sabiriyah Field, North Kuwait

Simultaneous Inversion of Clastic Zubair Reservoir: Case Study from Sabiriyah Field, North Kuwait Simultaneous Inversion of Clastic Zubair Reservoir: Case Study from Sabiriyah Field, North Kuwait Osman Khaled, Yousef Al-Zuabi, Hameed Shereef Summary The zone under study is Zubair formation of Cretaceous

More information

Uncertainties in rock pore compressibility and effects on time lapse seismic modeling An application to Norne field

Uncertainties in rock pore compressibility and effects on time lapse seismic modeling An application to Norne field Uncertainties in rock pore compressibility and effects on time lapse seismic modeling An application to Norne field Amit Suman and Tapan Mukerji Department of Energy Resources Engineering Stanford University

More information

Cross-well seismic modelling for coal seam delineation

Cross-well seismic modelling for coal seam delineation P-134 Sanjeev Rajput, CSIRO, P. Prasada Rao*, N. K. Thakur, NGRI Summary Finite-difference analyses is attempted to simulate a multi layered complex coal seam model in order to differentiate top and bottom

More information

Downloaded 09/16/16 to Redistribution subject to SEG license or copyright; see Terms of Use at

Downloaded 09/16/16 to Redistribution subject to SEG license or copyright; see Terms of Use at Data Using a Facies Based Bayesian Seismic Inversion, Forties Field, UKCS Kester Waters* (Ikon Science Ltd), Ana Somoza (Ikon Science Ltd), Grant Byerley (Apache Corp), Phil Rose (Apache UK) Summary The

More information

P1488 DECC Relinquishment Report OMV (U.K.) Ltd.

P1488 DECC Relinquishment Report OMV (U.K.) Ltd. 1. Synopsis: P1488 DECC Relinquishment Report OMV (U.K.) Ltd. UK offshore Licence P1488 was awarded to OMV, Dana Petroleum and Rocksource in the 24 th Licence Round on the 24 th of January 2007 as a frontier

More information

SEG Houston 2009 International Exposition and Annual Meeting SUMMARY

SEG Houston 2009 International Exposition and Annual Meeting SUMMARY Simultaneous joint inversion of MT and CSEM data using a multiplicative cost function Aria Abubakar, Maokun Li, Jianguo Liu and Tarek M. Habashy, Schlumberger-Doll Research, Cambridge, MA, USA SUMMARY

More information

Recent advances in application of AVO to carbonate reservoirs: case histories

Recent advances in application of AVO to carbonate reservoirs: case histories Recent advances in application of AVO to reservoirs: case histories Yongyi Li, Bill Goodway*, and Jonathan Downton Core Lab Reservoir Technologies Division *EnCana Corporation Summary The application of

More information

Some consideration about fluid substitution without shear wave velocity Fuyong Yan*, De-Hua Han, Rock Physics Lab, University of Houston

Some consideration about fluid substitution without shear wave velocity Fuyong Yan*, De-Hua Han, Rock Physics Lab, University of Houston ain enu Some consideration about fluid substitution without shear wave velocity Fuyong Yan*, De-Hua Han, Rock Physics Lab, University of Houston Summary When S-wave velocity is absent, approximate Gassmann

More information

IPTC PP Challenges in Shallow Water CSEM Surveying: A Case History from Southeast Asia

IPTC PP Challenges in Shallow Water CSEM Surveying: A Case History from Southeast Asia IPTC-11511-PP Challenges in Shallow Water CSEM Surveying: A Case History from Southeast Asia Sandeep K. Chandola*, Rashidah Karim, Amy Mawarni, Russikin Ismail, Noreehan Shahud, Ramlee Rahman, Paul Bernabe

More information

Petroleum Exploration

Petroleum Exploration Petroleum Exploration Upstream Petroleum Exploration The role of exploration is to provide the information required to exploit the best opportunities presented in the choice of areas, and to manage research

More information

Modeling of 3D MCSEM and Sensitivity Analysis

Modeling of 3D MCSEM and Sensitivity Analysis PIERS ONLINE, VOL. 3, NO. 5, 2007 641 Modeling of 3D MCSEM and Sensitivity Analysis Zhanxiang He 1, 2, Zhigang Wang 2, Gang Yu 3, Kurt Strack 3, and Haiying Liu 2 1 Department of Geosciences, University

More information

Verification of Archie Constants Using Special Core Analysis and Resistivity Porosity Cross Plot Using Picket Plot Method

Verification of Archie Constants Using Special Core Analysis and Resistivity Porosity Cross Plot Using Picket Plot Method Int'l Journal of Computing, Communications & Instrumentation Engg. (IJCCIE) Vol. 4, Issue (207) ISSN 2349-469 EISSN 2349-477 Verification of Archie Constants Using Special Core Analysis and Resistivity

More information

Fr Reservoir Monitoring in Oil Sands Using a Permanent Cross-well System: Status and Results after 18 Months of Production

Fr Reservoir Monitoring in Oil Sands Using a Permanent Cross-well System: Status and Results after 18 Months of Production Fr-01-03 Reservoir Monitoring in Oil Sands Using a Permanent Cross-well System: Status and Results after 18 Months of Production R. Tondel* (Statoil ASA), S. Dümmong (Statoil ASA), H. Schütt (Statoil ASA),

More information

Summary. Seismic Field Example

Summary. Seismic Field Example Water-saturation estimation from seismic and rock-property trends Zhengyun Zhou*, Fred J. Hilterman, Haitao Ren, Center for Applied Geosciences and Energy, University of Houston, Mritunjay Kumar, Dept.

More information

Seismic reservoir and source-rock analysis using inverse rock-physics modeling: A Norwegian Sea demonstration

Seismic reservoir and source-rock analysis using inverse rock-physics modeling: A Norwegian Sea demonstration 66 Seismic reservoir and source-rock analysis using inverse rock-physics modeling: A Norwegian Sea demonstration Kenneth Bredesen 1, Erling Hugo Jensen 1, 2, Tor Arne Johansen 1, 2, and Per Avseth 3, 4

More information

Site Characterization & Hydrogeophysics

Site Characterization & Hydrogeophysics Site Characterization & Hydrogeophysics (Source: Matthew Becker, California State University) Site Characterization Definition: quantitative description of the hydraulic, geologic, and chemical properties

More information

Applications of finite-difference modelling to coalscale seismic exploration

Applications of finite-difference modelling to coalscale seismic exploration Applications of finite-difference modelling to coalscale seismic exploration Shaun Strong 1,2, Troy Peters 1 1. Velseis Pty Ltd 2 University of Queensland Introduction Geological environments with significant

More information

Rock Physics and Quantitative Wavelet Estimation. for Seismic Interpretation: Tertiary North Sea. R.W.Simm 1, S.Xu 2 and R.E.

Rock Physics and Quantitative Wavelet Estimation. for Seismic Interpretation: Tertiary North Sea. R.W.Simm 1, S.Xu 2 and R.E. Rock Physics and Quantitative Wavelet Estimation for Seismic Interpretation: Tertiary North Sea R.W.Simm 1, S.Xu 2 and R.E.White 2 1. Enterprise Oil plc, Grand Buildings, Trafalgar Square, London WC2N

More information

Relinquishment Report for Licence Number P1356, Block 48/8c March 2008

Relinquishment Report for Licence Number P1356, Block 48/8c March 2008 Relinquishment Report for Licence Number P1356, Block 48/8c March 2008 Table of Contents 1.0 Header...3 2.0 Synopsis...3 3.0 Exploration Activities...4 4.0 Prospectivity Analysis...11 5.0 Conclusions...12

More information