Geolog6.6.1 Multimin Tutorial
Geolog6.6.1 Multimin Tutorial
1
Multimin Tutorial
Copyright © 1998-2006, Paradigm Geotechnology B.V. All Rights Reserved.
The information in this document is subject to change without notice and should not be construed as a commitment by Paradigm
Geotechnology B.V. or any of its affiliates (collectively, "Paradigm"). Paradigm assumes no responsibility for any errors that may
appear in this document.
The Copyright Act of the United States, Title 17 of the United States Code, Section 501 prohibits the reproduction or transmission
of Paradigm’s copyrighted material in any form or by any means, electronic or mechanical, including photocopying and recording,
or by any information storage and retrieval system without permission in writing from Paradigm. Violators of this statute will be
subject to civil and possible criminal liability. The infringing activity will be enjoined and the infringing articles will be impounded.
Violators will be personally liable for Paradigm’s actual damages and any additional profits of the infringer, or statutory damages
in the amount of up to $150,000 per infringement. Paradigm will also seek all costs and attorney fees. In addition, any person
who infringes this copyright willfully and for the purpose of commercial advantage or private financial gain, or by the reproduction
or distribution of one or more copies of a copyrighted work with a total retail value of over $1,000 shall be punished under the
criminal laws of the United States of America, including fines and possible imprisonment.
The following are trademarks or registered trademarks of Paradigm Geotechnology B.V. or any of its affiliates (collectively,
"Paradigm") in the United States or in other countries: Paradigm™, Paradigm logo, Artist™, Cement-IT™, EarthStudy 2D™, e-
Geoscience™, Epos™, Explorer™, EZdip™, Focus™, GeoDepth™, GeoLAN™, Geolog®, GeoScene™, GeoSec™, iMap™,
IntegralPlus™, LogLAN™, MapLAN™, NexModel™, OpenGeo Tools™, PlaNet™, Probe™, Pump-IT™, Reservoir Navigator™,
ResSim™, Rock & Fluid Canvas™, SeisEarth™, SeisFacies™, SeisX™, Slider™, SolidGeo™, Stratimagic™, Surveyor™,
Sysdrill™, Sysdrill Director™, Sysdrill DirectorGeo™, Sysdrill IDEAS™, Sysdrill SaG™, Vanguard™, Velocity Navigator™, VG
Tools™, VolLAN™, VoxelGeo™, VXplot™, and/or other Paradigm products referenced herein..
The following are trademarks or registered trademarks of their respective companies or organizations:
Adobe Reader - Adobe Systems, Inc.
corEVAL is a trademark of ACS Laboratories Pty Ltd.
FLEXlm - Macrovision Corporation
Exceed and Exceed 3D - Hummingbird Communications
GOCAD - Earth Decision Sciences
GeoFrame, GeoQuest, GeoNet, GeoShare - Schlumberger Ltd.
HP, HPGL, HPGL/2 and PaintJet - Hewlett-Packard Company
IBM AIX - IBM Corporation
IRIS, IRIX, Onyx, Origin - Silicon Graphics, Inc.
JustCGM and Justplot - Justcroft International Ltd.
Linux - Linus Torvalds
Microsoft, MS, MS-DOS, Windows, Windows NT, Visual C++-Microsoft Corporation
Motif and OSF/Motif - Open Software Foundation, Inc.
Netscape - Netscape Communications Corp.
NuTCRACKER - Mortice Kern Systems, Inc.and MKS Software Inc.
NVIDIA, Quadro, GeForce - NVIDIA Corporation
OpenWindows - Sun Microsystems, Inc.
OpenWorks - Landmark Graphics Corporation
Operton, Athlon - AMD
Pentium, Xeon - Intel Corporation
Red Hat - Red Hat, Inc.
RP66/A - Schlumberger Ltd.
SDI and PrintMaster - System Development, Inc.
SunOS, Solaris - Sun Microsystems, Inc.
SPARCstation, UltraSPARC - SPARC International, Inc.
UNIX - The Open Group
WOW - Landmark Graphics Corporation
XVision Eclipse - The Santa Cruz Operation, Inc.
X Window System - Massachusetts Institute of Technology
This tutorial is designed, with a hands-on, practical approach, to assist you in learning how to
perform sophisticated formation evaluation using Multimin, the advanced petrophysics
application of Paradigm’s Geolog system.
The tutorial is taught by professional trainers to small groups of interested people from the
petroleum industry. The material is presented by following worked examples on a workstation.
Prerequisites
To gain the most from this tutorial, the following prerequisites are advised to avoid
unnecessary frustration:
• A working knowledge of the Geolog system and its conventions, the methods to
import data into the software, and the way data is stored.
• A working knowledge of the methods used in Geolog to perform basic
preparation tasks for petrophysical data including:
— depth-matching, editing and despiking for petrophysical logs;
— environmental corrections;
— data displays;
— crossplotting.
Document Conventions
In this document, all INPUT to the computer is in Bold Courier New, while all OUTPUT from
the computer is in Courier New, but not bold.
Tutorial Data
The following additional files (files not supplied with software) are used in this tutorial:
Multimin is a software application that provides advanced formation evaluation answers for
petroleum geologists, engineers, and petrophysicists. The application is based in the
probabilistic, or optimizing, approach to modelling wireline and rock data. The name
"Multimin" is an abbreviation of "Multiple Mineral", revealing the application’s strong
geological heritage and nature. However, current petrophysical theory, rapid numerical
programming, and an intuitive user interface underpin this geological nature. Multimin users
are able to rapidly produce detailed petrophysical analyses of formations using a wide range
of wireline and rock data.
The following steps will guide you through your first Multimin analysis, and then explore,
review, challenge and change this analysis. After completing these steps, you will have a
broad understanding of how to accomplish a Multimin analysis.
The example you will follow is from a well that was drilled in the Alberta Basin of Canada during
the late eighties. The data you will analyze are the original core and wireline logging
measurements that were available at the time the completion decision was made on this well.
The well was drilled to test the Upper Triassic Halfway Formation, and the zone to be analyzed
is the only possible candidate for a gas or oil completion. Based on your petrophysical analysis
and recommendation, this expensive exploration well would either be abandoned, or
additional expense equal to half the original well cost will be incurred to complete the zone you
have predicted to be hydrocarbon productive.
Procedure
Before performing an analysis, you need to gather the required information (provided below
for this tutorial) and then you will:
Required Information
Before analyzing the example well, you need to gather all the information and data
that will be required. This includes loading the "hard" data such as wireline logs
and core analyses. It also includes reviewing all fluid analyses that have been
performed for the formation, including formation water, oil and gas analyses. The
"soft" data would include knowledge of the geology of the zone for analysis.
Initially, all geological and engineering information from adjacent wells should be
reviewed before a Multimin analysis. Once a few wells in the area have been
analyzed, you will develop particular models that can be applied to analysis
intervals. In this context, model means an appropriate combination of:
• Choice of minerals for available measurements.
• Choice of fluids for available measurements.
• Response equation parameters for all minerals and fluids on each measurement.
• Fluid analysis parameters.
• Electrical properties such as cementation factor.
• Mineral properties such as cation-exchange capacity.
Exercise 1
Note also the "standard" names for the regular core analysis results—
PHIT_CR, KA_CR, RHOG_CR and SWT_CR have been applied on
loading.
These are the names used in the Paradigm-supplied Multimin layouts,
crossplots, and models. You could choose to use alternative names, but you
would then need to edit the layouts, crossplots and models, or include the
alternative names in the alias.alias file. See "Appendix A - Multimin
Aliases’ on Page 78 for a review of the implementation of aliases in relation
to Multimin.
Exercise 2
Depth matching and curve editing are processes that are common in
petrophysical analyses. They are a normal part of a structured approach to
any petrophysical analysis, including one that includes a Multimin analysis.
The curve editing process may include despiking, splicing, and even
wholesale graphical and text replacement of clearly spurious data.
Various Paradigm-supplied layouts (e.g., composite) can be used to view the
data to verify its quality, or you can build your own layout. The zone of interest
is from 4820-5085 feet (1470-1550 meters).
For the sake of this tutorial, let us consider this data set to have no need for
either depth matching or editing. If you disagree, please feel free to perform
these tasks now using the normal Geolog processing options.
Note the value of DFD is high due to the presence of barites in the mud.
Environmental Corrections
The next phase of the structured processing chain would normally be the application
of appropriate environmental corrections. These would be applied to the edited and
shifted curves to produce a new group of curves that had the suffix _COR.
For this example, and at this point, we are not going to perform any environmental
corrections. This is a matter of expedient choice, and could be argued on the basis
of the data to hand. The hole appears to be mainly in gauge, has an optimal
calibration diameter (200mm), and shows only moderate invasion. Even the heavy
barites mud would not provide a significant neutron correction.
Procedure
In this step, the Multimin analysis is performed. All ancillary data which is, in this example, rock
descriptions, fluid analyses, and analysis intervals should be reviewed. You then:
• Determine the minerals and fluids to be used, and take into consideration what
unusual effects whole mud invasion will have on log responses.
• Build a Multimin model for the HALFWAY_TEST well.
• Run an analysis.
• Review the results of the analysis.
• Make a completion decision based on the review.
Rock Descriptions
Multimin is a petrophysical package that calculates values for porosity, and the
relative volumes of water and hydrocarbons contained in the porosity. It is also a
geological package, calculating the volumes of each of the minerals that you
choose to include in the analysis.
The appropriate choice of which minerals to model is often very easy. However,
sometimes the accuracy and time to perform an analysis is critically dependent on
the choice of the right minerals, and it is worthwhile to use all sources of
knowledge to make this choice:
• The most definitive data is obtained from core descriptions and analyses—X-Ray
Diffraction (XRD), Scanning Electron Microscopy (SEM), particularly for clay
types, and thin section petrography.
• Extremely useful insights are also gained from reviewing binocular microscope
descriptions of drill cuttings.
• In difficult zones, SEM viewing of cuttings can provide vital clues. In using the
SEM, Energy Dispersive X-Ray Spectra (EDX) can be measured on particular
mineral grains. This not only identifies a mineral, but also gives the composition,
enabling a better choice of response parameters for wireline logging equations.
In new exploration areas, or in areas of restricted information, there will be times
when the only data available to determine the correct mineralogy and fluids will be
the wireline logs themselves. In these cases, your Multimin analysis will be very
much from a log analyst's perspective!
4822.83 - 4831.03 1470.0 - 1472.5 Tight siltstone with major clay, 40% dolomite, no
fluorescence.
4831.03 - 4839.24 1472.5 - 1475.0 Very fine grained marl with 10% dolomite, minor clay,
6% intergranular porosity and fair permeability.
Moderately well sorted quartz grains. Poor
fluorescence. Coring recommended.
4839.24 - 4884.19 1475.0 - 1488.7 Cored interval. 100% recovery. Separate core
description. Mud was oil flecked during coring
operation.
4884.19 - 4896.65 1488.7 - 1492.5 Fine grained marl with well sorted quartz grains.
10 - 20% dolomite. Porosity from 6 - 15% decreasing
with degree of dolomitization. Good permeability.
Poor to moderate fluorescence.
4896.65 - 4954.07 1492.5 - 1510.0 Fine grained marl. Very well sorted quartz grains.
10 - 20% dolomite. Porosity from 6 - 9% decreasing
with degree of dolomitization. Fair permeability and
poor fluorescence.
4954.07 - 5019.69 1510.0 - 1530.0 Similar to 1492.5 - 1510.0 meters with decreasing
fluorescence with depth.
5019.69 - 5085.30 1530.0 - 1550.0 Similar to 1492.5 -1510.0 meters with no
fluorescence.
5085.30 - 5150.92 1550.0 - 1570.0 Tight siltstone with major clay, 40% dolomite, no
fluorescence.
4839.24 - 4840.88 1475.0 - 1475.5 Upper very fine grained, moderately well sorted
quartzose sandstone variably (5-15%) dolomitized.
Minor illitic clay. Fair permeability with 5%
inter-granular porosity. Poor oil stain
4840.88 - 4854.00 1475.5 - 1479.5 Tight marl with abundant unidentified clay. Matrix is
composed of roughly equal quartz and dolomite. No
oil stain.
4854.00 - 4884.19 1479.5 - 1488.7 Upper fine grained quartzose sandstone. Minor illitic
clay. Very well sorted. Dolomitization varies from
5 - 25%. Dolomitization reduces porosity, which
varies from 9 - 15%. Good permeability with whole
barites mud invasion seen in open pores. Poor oil
stain.
Exercise 1
3. The other information that piques our interest is the invasion of the pore
system by the whole mud, presumably carrying barites into the formation. In
setting up the model, the following question should be considered:
Fluid Analyses
We have already chosen gas as the likely hydrocarbon present in this well, and it
is desirable to know the gas specific gravity so that appropriate logging responses
to gas can be estimated.
If the zone was subjected to a production test, a drill stem test, or a wireline
formation fluid test, a sample of the hydrocarbon might have been recovered.
However, for the purposes of this example, no such recovery was available at the
time of performing the Multimin analysis. There are 3 possible sources of
information:
• From the mud logging data.
• Using data from nearby exploration wells that tested gas, or from wells within
the same gas field if this were a development well.
• By inference from reviewing data from elsewhere in the basin.
Nearby producing wells have gas analyses with a specific gravity of 0.65.
No formation water was recovered by testing. Fortunately, data from offset wells is
present, and a comprehensive formation water catalog of the formation suggests a
water resistivity of 0.06 ohmm at 77 degF (25 degC).
Resistivities of mud, mud filtrate, and mud cake were measured by the wireline
logging contractor.
Analysis Intervals
Multimin analyses are generally performed on an interval basis, although a simple
analysis could be run over a defined range. For the HALFWAY_TEST well, only
one interval will be analyzed, but it is common practice to have many analysis
intervals in a single well. Interval subdivision could be a result of lithological or fluid
changes, or possibly as a result of hole conditions, etc. Interval subdivision would
be performed by reviewing the wireline logs, mud logs, core descriptions, test
results, and supported by crossplot analysis.
The HALFWAY interval has been preset in the INTERVAL.INTERVAL log which
was loaded from the Geolog ASCII file.
1. Select Petrophysics > Multimin > Model Maintenance from the Well menu
to display the Multimin Model Editor window (see Figure 2).
The Model Maintenance option can be used to build new models or change
existing models. Models are well-independent, and so can be used in any well
where that particular lithology/fluid model is applicable.
2. Select Model > New (Old Defaults) to open the default model. The initial
display is the Unknowns+ > Responses screen (see Figure 3).
To start, you can choose either of 2 default models. The difference between
these 2 default models lies in the wireline log responses of the different
minerals (see the "Note" on Page 20 for further information).
) When you select any option from the Unknowns+ menu, the Unknowns (top)
section is always displayed and you can therefore select/change the
Unknowns options from any of these combination screens.
Note user defined special minerals and fluids are possible. More on this in
"Complex Example" on Page 45 and " Verification Example" on Page 59.
The Responses section allows you to assign response values to each of the
selected Unknowns (volumes) for each of the chosen logs. We will do this in
the next section.
5. Logs are turned on or off under the "Selected" column in the Wireline
Equations section. Ensure the following equations are on (set to YES):
— Formation density
— Neutron
— Sonic transit time
— Photoelectric absorption
— Total gamma
— Unflushed conductivity
— Flushed conductivity
6. We will accept the default parameters for our chosen volumes and equations.
The green cells in the "Method" column can be toggled to allow the choice of
different response equations. There are options for Neutron, Sonic transit
time, Unflushed conductivity, and Flushed conductivity.
7. At this stage, leave the Neutron at LINEAR, the Sonic transit time at WYLLIE
LINEAR and select the DUAL-WATER LINEAR for both conductivity equations.
Each of the available options will be explained by your instructor.
8. In the Unc. Method (Uncertainty Method) column, leave the defaults as shown
(Value for all logs except CT and CXO, which are Interval).
The uncertainty method can be toggled between value and log for all wireline
equations. There is a third (Interval) option for CT and CXO. Uncertainty
values are presented for each response equation. The role of uncertainties is
to provide a relative weighting for the effect of each of the equations in the
Multimin analysis. A higher uncertainty will cause the equation to have a
lesser effect on the solution. See "Multimin Uncertainties" in the Multimin
online help documentation for further information.
9. The seven equations chosen represent a typical suite of wireline logging data.
However, other tools are represented by equations here, and you should
scroll through the complete list of equations. Below the partial list shown are
many more types, including user defined response equations.
The default values you see are dependent upon the initial default model chosen. The "Old
Defaults" (as we have chosen for this model) are based on various published sources. Three
references are suggested to represent (but not to limit to) the myriad of data sources:
"Element Mineral Rock Catalog", O. Serra. Schlumberger Manual SMP- 7040, 1990.
"Mineral Logging Parameters: Nuclear and Acoustic", A. Ellis et al. The Technical Review,
Schlumberger, v. 36, no. 1, p. 38-52, 1988.
The "Herron Defaults" are based on more recent work described in the following paper.
Which defaults you use will depend on local experience. After choosing the correct minerals
and fluids, setting these responses is the most important factor in obtaining a good Multimin
analysis. As a Multimin user you will probably build up a small library of references containing
tool responses to various fluids and minerals. Multimin is a modeling program where a
modeled or predicted set of logs is compared to a measured set of logs (and we include core
measurements as part of these available logs). The closer the match, the better the analysis.
The match can be improved by varying the responses shown here. However, only reasonable
changes, based on real evidence, should be made. The change in responses should be within
the range of known values for a mineral or fluid.
It is suggested that the Herron defaults would be more applicable with authigenic clays while the
old defaults might in fact be composite values (possibly including small percentages of heavy
radioactive minerals and other impurities). For these reasons, you might find that the Herron
defaults work well in carbonates while the old defaults work well in shaley sands. There is no
hard and fast rule. You will need to see which ones work best in different areas or even fields,
and it might also depend on the number of logs you have available—with more logs you can
work with more minerals and so use more pure minerals rather than pseudos or composites.
Clay bound water responses are set equal to the responses of the flushed zone free water
(except for conductivity), and the responses for clays (illite, in this example) are those for dry
clays. When the Dual Water saturation equation is used, Multimin internally computes the
volume of clay bound water associated with the clay as a function of CEC, temperature and
alphax/alphau. Response parameters for wet clays are then calculated. See the "Multimin
Technical Reference" documentation in Geolog's online help for details. When the analysis is
executed, Multimin will initially solve for volumes of wet clays. In storing answers to the well
database, it saves the volumes of each dry clay, and the combined volume of clay bound
water. Whether bound water volumes, clay or shale volumes are computed depends on the
saturation equation chosen.
10. We will not use any core data in the analysis at this stage, so we will not
select any equations from the Core, XRD-IR, or Petrography screens.
However, you should look at these screens (under the Equations menu) to
see what other types of data can be used to help in a Multimin analysis.
11. Select Method > Parameters to display the Dual Water Saturation
Parameters screen.
12. In the Method row in the Dual Water Saturation Parameters section, ensure
the Saturation factor (W) is set to Value.
Notice the values for Cementation exponent (M) and Saturation exponent (N)
are grayed out—WHY?
Setting Saturation factor (W) to Log would give you the option to account for a
variable cementation exponent, but leave as Value for now.
At this stage, ensure "Allow expansion of Clay Bound Water" and "Allow
application of Beta-dilute" are toggled ON—WHY?
Note these are mineral properties, and not wireline equation responses. The
grain densities listed here have different values from the formation density
wireline responses—WHY?
The Cation Exchange Capacity (CEC) values listed here for various clays
should be treated as starting points only. Clays are chemically extremely
variable, and their CECs are also extremely variable. It is a good strategy to
vary these values to assist in the matching of measured to predicted
conductivity curves. For instance, if the predicted unflushed zone and flushed
zone conductivity curves were consistently less than the measured data,
consider increasing the CEC of the dominant clay. This will put in place more
clay bound water, increasing the predicted conductivities.
The program constraints (which are not user editable) are displayed.
Note the "Prog Water Mud" would be replaced by "Prog Oil Mud" if Mud Type
under Verify > Edit Parameters had been set to Oil instead of Water.
Look at the Prog X Bndwat and Prog U Bndwat constraints. What do the
numbers mean?
User constraints are defined by the user, and are shown on the menu to be
available for user editing. They must be toggled to "Yes" to be made active.
When this is done, the entire row on the screen is highlighted in green. See
"Constraints" in the Multimin online help documentation for information on the
application of User Constraints.
17. Bounds are placed on the upper and lower values the chosen volumes can
have. Leave the default values.
19. Leave "Calculate fluid responses" set to YES. Note the color of the
corresponding Units field is gray.
22. Select the WIRE set and click the Log Headers tab.
23. Locate the mid-zone values for Formation temperature (FTEMP) and
Formation pressure (FPRESS), and enter these values.
24. Copy both logs to the Interval set, and from there enter the values into the
Model FTEMP and FPRESS parameters.
27. Check the Responses. Note that some of the fluid fields are colored gray.
The gray Response fields and the uncertainties (Unc. Value/Log column) for
CT and CXO (since the method is set to "Interval") change value. They have
been calculated from the parameters entered in the Edit Parameters screen
using the fluid properties package within Multimin. When "Calculate fluid
responses" is set to YES, none of the gray Response fields are user-editable.
Do not be alarmed by this mass of numbers! They are explained fully in the
"Model Quality" section of the Multimin Technical Reference and Multimin
online help documentation.
) The major feature is the interval condition number, reported as 3.25033. The
condition number is a method of describing how well the measured curves
are capable of discriminating the unknown components. One can have the
best quality measurements, but if they are not suited to discriminating the
chosen unknowns, then this number will tell you. The smaller the number
the better.
Multimin describes this model as having a good interval condition number so
we can feel confident in running this analysis. However, we should still run
some models with bad condition numbers, even up to values approaching 10,
the limit at which Multimin will cut out for fully linear models (you can set this to
a higher value in the Verify > Edit Parameters screen—experienced users
only).
) Please take careful note of the last two lines of the report. Should you not
like the statistics at this point, you should edit the parameters, either
removing or changing components, or adding measurements.
Exercise 3
Run an Analysis
4. Select Petrophysics > Multimin > Run Analysis to display the tp_multimin
Launcher screen.
9. Enter 0.06 for RWS; the formation water resistivity is known to be 0.06 ohmm
at 77 degF (25 degC).
12. Accept the default value of 0 for the KCl content in the mud (KCL_MUD).
13. Double click in the MM_MODEL field, click on the Dropdown List button
and select your model from the displayed list (in this example, halfway_1
is selected).
14. Verify that FTEMP, FPRESS and RMF logs have been selected at the bottom
of the screen.
We do not have any alternate models to run so we are ready to run the
analysis.
15. Click Start to run the module and the numerical optimization process.
The Message area provides you with information about the depths and
intervals that Multimin processed at run time. When the processing is
complete, the results are displayed in the multimin_recon layout (see
Figure 13).
Exercise 4
3. What about the quality indicators? The quality curve indicates a reasonable
result, except for short intervals where it exceeds one.
4. Do the measured and predicted curves match each other sufficiently closely?
What would be the effect on the analysis if the different predicted curves were
forced to match their measurement curves?
Exercise 5
Completion Decision
Your Exploration Manager has just walked in the door; he announces that the
project geologist, reservoir engineer and drilling engineer cannot agree on
anything about this well. He is looking forward to your completion decision, and
would appreciate having it within the next fifteen minutes. He may order a drill
stem test to verify your advice.
What is your completion decision?
_______________________________________________________________
_______________________________________________________________
Procedure
This step will explore, review, challenge, and change the Multimin analysis results of the first
example completed in "First Example - Introduction" on Page 5.
This step will assume you are aware of conventions and procedures, and will start to
explore the geological and petrophysical "what-ifs".
You will:
Exercise 1
1. Go back to the model editor and, if necessary, open the model you built in
"First Example - Introduction".
4. Take a look at the parameters for X Zone Special Fluid. Note that the default
values are for barites. The sonic travel time is set equal to the default for
water, reflecting that barites is part of the fluid filled porosity, and not firmly
attached as part of the rock matrix.
6. Double the uncertainty from 0.32 to 0.64 in the Unc Value/Log field for the
Photoelectric absorption. As we are using a least squares format, this will
mean this curve will now only have a quarter of its previous effect on the
analysis.
7. Select Verify > Display Statistics. The message indicates the model has
been modified, so you must then select Verify > Verify Model.
A new set of descriptive statistics is displayed. The values are more adverse
than those from the previous analysis, but are still quite acceptable for
attempting an analysis.
The reason that most of the minerals have been tagged is the very low
predicted error for X-ZONE SPECIAL FLUID, the barites component. Any
mineral that has a predicted error greater than twenty times the minimum
predicted error is tagged. In reality, the absolute values of the predicted errors
may be acceptable, and we do not know until we run the optimization. See the
Multimin Technical Reference online help documentation for a discussion of
these statistics.
) We are now ready to try this new model. You might want to save as a new
name so you can easily go back to your first model, if necessary.
8. Select Model > Save As, to save the model with a new name.
1. Run your second model with "Allow expansion of claybound water" turned off
(select Method > Parameters).
2. Run your original model (without the special fluid) without using the
photoelectric absorption (select Equations > Wireline).
3. Attempt to add the three minor minerals at the same time to the current
analysis. You will not succeed, but this will illustrate a safety feature of
Multimin. WHY were you not allowed this many volumes?
Exercise 3
Exercise 4
1. Select Petrophysics > Multimin > Plot to display the plotter module which is
partially reproduced in Figure 15.
Exercise 5
Procedure
This section initially focuses on a review of the Model Editor and then introduces an example
that is designed to cover much more complex procedures. Variable uncertainties, model
switching and depth reference conventions are also introduced. You will:
Exercise 1
) There are many subtleties to this analysis, so do not be surprised if you find
some apparently unusual parameters!
Fluid properties:
Oil 45.4 DAPI Nearby well DST
Gas 0.613 SG Nearby well DST
Rw 0.5 ohmm @ 77 degF (25 degC) DST on Heavy
Temperature of target zone = 120 degF (49 degC)
KCl = 3%
5. Select Petrophysics > Parameter Picking > RHOMA-UMA plot to open the
RHOMA-UMA crossplot (see Figure 16).
This plot shows the position of end member minerals. From this plot, what is
the simplest mineral assemblage that would adequately explain the data?
8. Iterate through your model a few times to obtain the best possible result.
9. Once you are happy with your result, open the multimin_recon layout that
shows the spectral gamma data. In particular, have a look at the interval
7513.12-7529.53ft (2290-2295m). What do you think about your prediction?
10. Add the POTA log from the spectral gamma ray to your analysis. What is the
result?
So what does this mean? . . . It probably means that we have either picked a
mineral incorrectly or that a mineral is missing from our model.
Exercise 3
1. Have a look at the Plug surrounds descriptions (Table 5 on page 46) to see if
any clues are given.
2. Also look at the following crossplots to see if there is anything unusual about
the data:
— thor_pota_uran.xplot
— rho_nphi_gr.xplot (try changing the Color Expression to URAN_COR,
THOR_COR and U, as well as GR)
) You may need to adjust the xplot (e.g. limits and color bar) and Datum in
order to pick up the subtleties (see Figure 18).
Having decided that an opaque or heavy mineral may be required, and that this
mineral needs to be high in thorium and uranium, and low in potassium, the first
step would be to go to the rock element catalogs and try to find a suitable
mineral (see "Appendix B - THOR Minerals from Rock Element Catalogue’ on
Page 80).
4. Format the datum such that the well is being processed from the range
WELL_TOP to WELL_BOTTOM.
From the statistics calculated by the module, the relationship for GR from the
individual spectral components can be obtained and it is expected that this is
consistent, as it is a function of tool calibration.
8. Now build a model incorporating the spectral gamma data and a heavy mineral.
10. Also try running the model without one of the spectral gamma logs, using the
GR instead.
Exercise 4
Variable Uncertainties
Uncertainties are values that reflect the precision of measurements and vary as a
result of tool physics, tool calibration, borehole conditions and other factors.
Uncertainty on each wireline log measurement is a function of several factors,
usually similar for most types of measurements. For example, the uncertainty on a
gamma ray log is a function of:
• Changes in hole size and mud weight from nominal conditions.
• Uncertainty in depth recording - does the measurement represent the value at
the exact depth, or part way to the frame of data immediately above or below?
• Precision of detector physics.
• Accuracy of tool calibration.
Exercise 5
Input Sets
There are several different sources of measurements that can be used in Multimin.
We have loaded two kinds of data, with different data structures:
• Wireline data in a set called WIRE. This set is a regularly sampled point data
set. Geolog sampling concepts are explained in the Multimin online help
documentation.
• Core data in a set called CORE. This data is irregularly sampled point data. It
will not match depth for depth, measurement for measurement, the regularly
sampled wireline data.
Other data could have been loaded into two other kinds of sets called (by Multimin
convention):
• XRD_IR contains XRD/IR (X-Ray Diffraction / Infra Red spectroscopy)
measurements of mineral weight fractions.
• ROCK petrography descriptions of cuttings. Their sampling could be either
limited or extensive, and with variable regularity.
Output Sets
Just as we can nominate a default input set, we can also nominate a different
output set. This can be very handy for exploring "what-ifs" and comparing results
from different models. It is also an efficient way of knowing which volume came
from which pass—if some minerals or fluid volumes are included in some models
but not in others then, if all results are sent to the same default set, it is difficult to
know which result log version goes with which pass.
Multimin features logic for model switching between a primary model and up to
10 alternate models.
2. Rerun Multimin with your new model, including the CORE logs selected by
the expression "PHIT_CORE <> missing".
3. What would we need to do to make the model with the CORE logs the
primary model and our original model the alternate?
Procedure
This section is designed to show that the technology behind the Multimin program really
works. To do this, you will:
Overview
After studying the theory of probabilistic petrophysics, and tackling the
mathematics behind response equations and optimization schemes, we have
moved on to perform several practical examples. Complex and involved formulae
have been replaced by an easy to use program that hopefully makes probabilistic
petrophysics fun.
If Multimin has achieved this, there will appear to be little connection between the
theoretical and practical aspects. You may therefore be wondering—how can one
verify that Multimin is really doing what it claims to do?
This process will be performed for an example with only linear response
equations, and then repeated for an example with non-linear conductivity and
neutron equations.
Exercise 2
3. Select Tools > Compile to compile the program to run on the version of
Geolog currently being used for the course.
We are going to create seven distinct lithologies of varying mineralogy and fluids.
Within each lithology, the relative proportion of each mineral remains the same, as
does the relative proportion of each fluid. However, the porosity is increased from
0 to 0.5 over 100 depth frames.
Each of the seven lithology and fluids models is progressively more complex. The
first six are reasonable models, with the volumes of minerals and flushed zone
fluids summing to 1. The seventh model is deliberately unrealistic, as these
components are contrived to sum to 1.1 rather than 1.
Artificial measurements are created using response equations as detailed in
"Response Equations" in the Multimin Technical Reference online help
documentation. The loglan also creates the following artificial log curves:
— RHOB formation density
— TNPH neutron linear
— DT sonic transit time
— U photoelectric absorption cross-section, rather than PEF
— GR total gamma
— CT unflushed zone conductivity, and the corresponding RT
— CXO flushed zone conductivity, and the corresponding RXO
For the TNPH, CT, and CXO linear responses, corresponding non-linear
responses called TNPH_NL, CT_NL and CXO_NL are created.
) Note that the creation of CT, CXO, and U obviates the need to perform
preprocessing with Precalc. In fact, RT, RXO, and PEF are created, so
Precalc could still be used.
7. The default values can be used initially (you might want to experiment with
different parameters and models later).
Exercise 3
5. Set Conductivities Unc Method to Interval and to Value for all other logs
(accept the default values).
7. Select Verify > Edit Parameters and toggle "Calculate fluid responses" to NO.
8. Go back to a Responses screen and edit the orange cells to have the same
values as used in mm_check.model.
13. Enter the values used in mm_check.model for the RMF and FTEMP logs.
FPRESS can take any value in this case—WHY?.
The analysis can be run over the complete range or an interval set could be
created to experiment with different lithologies.
14. Compare the Multimin computed volumes to the known volumes and the
predicted logs to the artificial logs.
15. You may wish to expand this kind of verification example by modifying the
loglan routine to incorporate any components, or response equations, that
interest you.
Exercise 4
A Non-Linear Example
1. Create a new model using non-linear equations for neutron, unflushed and
flushed conductivity.
The curve NFUN shows the number of iterations required during non-linear
optimization. Note that, whereas the linear method gives exact solutions for
component volumes, the non-linear gives close to exact solutions. This is because
the non-linear optimization scheme used minimizes the objective function to a
small, but acceptable, minimum error (see the Multimin Technical Reference
online help documentation).
Procedure
This section presents an example that is designed to examine choice of clay types and their
associated parameters, with particular emphasis on CEC. You will:
Exercise 1
The log header advised that the PEF has been affected by barites
(144.8 lb/bbl of barite) and that the mud was contaminated with diesel. The
following parameter values were taken from the log header:
Fluid Properties
Exercise 2
Handling Coal
At approximately 9908.14 ft (3020 m), a washout zone is present and the logs
would tend to indicate that this may be due to the presence of a coaly horizon. As
this represents a significant change in lithology, we should attempt to use a
different model over this interval.
2. Call the log "coal" and set the values of one (1) in areas of coal and zero (0)
elsewhere by choosing an appropriate DT cutoff.
3. Create a coal model which will only consist of those components that we
should expect to find in a coaly horizon.
Independent Analysis
1. Build a simple model using minerals quartz, illite and kaolinite with gas and
free water as fluids.
2. Run the analysis—run it first without the Coal model and then including the
coal model.
5. What about Rw? How well are the conductivity logs reconstructed?
9. Duplicate the WIRE set and rename the duplicated set to WIRE_SHALE.
12. In the first wireline (Lithology) track, replace the following logs:
VOL_ILLITE with VOL_WCS1
VOL_KAOLINITE with VOL_WCS2
14. In this reservoir, there are probably two shale types—one in the upper section
and one in the lower section. First, run with a separate model in each, and
then run with one model containing two wet shales, and compare the results.
Which one do you think is better?
15. Now try changing the conductivity equation to Juhasz . . . what is the
difference?
DT = 92
GR = 105
φ N = 0.32
ρ b = 2.605
Figure 23: Parameter picks for shales in intervals A and B in Silt well
RT = 4.5
Procedure
This section introduces the Numerical Inversion Modelling By Linear Equations (NIMBLE)
module (tp-multimin_inverse) which solves the inverse problem for a single model over an
extended depth range. The inverse, or calibration, problem is when response parameters for
single or multiple mineral and/or fluid volumes on single or multiple tools are calculated over
an interval. This technique might prove useful when, for example, volumes are well known
from some independent source (like petrography or XRD) but the tool responses are not.
You will:
Exercise 1
1. Using Connect, delete your existing HALFWAY_TEST well and then reload it.
Modelling a "Pseudo-Mineral"
Often there will be heavy, possibly radioactive accessory minerals associated with
clays (or non-clay minerals). In some wells, there might be sufficient logs and
auxiliary data to be able to model these secondary minerals. However, there might
be other wells in the same field where less data is available and the extra volumes
required for the secondary mineral(s) could not be accommodated. NIMBLE could
possibly be used to good effect to model a composite "pseudo-mineral" which
would have some combined response of the clay(s) (or non-clay minerals) and the
accessory mineral(s).
3. Run General > Evaluate to create a VOL_NEW log that is the sum of the
VOL_QUARTZ and VOL_DOLOM logs.
6. Edit the model that was used to produce the current evaluation to remove
VOL_DOLOM. Leave VOL_QUARTZ in, as this now represents the
composite mineral whose response parameters we wish to optimize.
8. You will need to nominate the model that was used to produce the current
evaluation for parameter "MODEL_INVERSE" (i.e., HALFWAY_NIM).
Exercise 3
1. Use the XRD (or other data) in a Multimin model with normal uncertainties set
on their response equations.
4. Optimize the model with normal Multimin passes, producing mineral volumes
in the wireline set.
5. Edit the model to remove the XRD equations and to include all other
response equations with their normal uncertainties.
At the time of writing of this course, a suitable example of this application was not
available.
Procedure
This is an example that is designed to examine how to model the effects of shallow invasion
on tool response, and how to model secondary porosity.
Exercise 1
Independent Analysis
4. Review the section of the Geolog online help which covers modelling of vugs
and fractures, and how to model shallow invasion.
What is the most likely situation here?
Build models for each case and check the results against core again.
Overview
This appendix is intended to review the philosophy and application of aliases within Multimin.
The concept of aliases is discussed, and how aliases are structured and used within a Multimin
framework.
Alias Concepts
Aliases allow the mapping of one or more alternative log mnemonics to a desired log
mnemonic. It is also possible to alias alternative constant names to a desired constant name.
Multimin and other Geolog programs use as inputs particular log or constant names that are
expected to be present in the database. If, on processing, these are absent from the database,
the program will use an alias list to find the next most desirable curve or constant. An example
of the practical use of aliases would be the choice that must be made in using a gamma ray
log. The usual gamma ray curve is GR (a simple total gamma ray log), but it would be
appropriate to use SGR (the spectral total gamma ray log), or CGR (the spectral computed
gamma ray log). To accomplish this, it would be necessary to alias first SGR, then CGR to the
desired GR.
The Geolog6 project specs directory contains the Paradigm supplied alias.alias file.
This file can be viewed at the Unix level or by using an NT text file editor. The user can make
changes or additions to this file but these changes should be stored in the specs directory
at site or project level (otherwise, new releases of Geolog will overwrite your changes). The
alias.alias files are cumulative—the file at the project level is read first, then the site,
then the Geolog6 project level. Only previously undefined aliases are added at each
consecutive level. It is, therefore, only necessary to include changes in your project or site
alias.alias files.
• The alias list for GR allows you to name gamma ray curves with names that
reflect the tool that they were measured with or from different service
companies. Later aliases pick up curves measured with a spectral gamma ray
sonde.
• RHO has the alias NRHO representing a bulk density curve of enhanced vertical
resolution.
• RT has aliases that represent the various resistivity and induction tools, and
enhanced resolution processing.
The alias.alias file shows an interpretation of the plethora of naming schemes available
for the calipers measured with various sondes. For practical purposes, Geolog only uses one
caliper for environmentally correcting porosity curves, and one for correcting resistivity curves.
It is assumed that depth-shifted/edited curves will have the higher version number, and these
will be chosen by Geolog unless an alternative choice is forced by the user.
Most layouts, etc. call for environmentally corrected curves. If the curve is missing, the
priorities for aliases are:
The same priorities are assigned to predicted curves with the only difference being the
addition of the predicted suffix _PRED. Note the following exception to the general table:
Geolog 6.6.1 - Multimin Tutorial Appendix B - THOR Minerals from Rock Element Catalogue 80
Geolog 6.6.1 - Multimin Tutorial Appendix B - THOR Minerals from Rock Element Catalogue 81
Geolog 6.6.1 - Multimin Tutorial Appendix B - THOR Minerals from Rock Element Catalogue 82
Geolog 6.6.1 - Multimin Tutorial Appendix B - THOR Minerals from Rock Element Catalogue 83
Geolog 6.6.1 - Multimin Tutorial Appendix B - THOR Minerals from Rock Element Catalogue 84
Geolog 6.6.1 - Multimin Tutorial Appendix B - THOR Minerals from Rock Element Catalogue 85