/home/bill/OPM/0_OPM-Flow notes.txt ************************** 24Jul2017 09:15 Flow [executables, libraries, includes, shares, tests] +-----+ Executables : $ ls /usr/bin | grep "flow*" flock flow flow_ebos flow_legacy flow_mpi flow_polymer flow_sequential flow_solvent >> Ignore flock $ ls /usr/bin | grep "dune*" dune-am2cmake dune-autogen dunecontrol dunedoxynize dune-git-whitespace-hook dune-mpi-config duneproject +-----+ Includes : /usr/include/dune/[common,geometry] +-----+ Libraries : /usr/lib/dunecontrol/dune-[common,geometry] $ ls /usr/lib/x86_64-linux-gnu/ | grep "dune*" libdunecommon-2.4.1.so libdunecommon.a libdunecommon.so libdunegeometry-2.4.1.so libdunegeometry.a libdunegeometry.so $ ls /usr/lib/x86_64-linux-gnu/ | grep "opm*" libopmcommon.so.2017 libopmcommon.so.2017.04 libopmcore.so.2017 libopmcore.so.2017.04 libopmgrid.so.2017 libopmgrid.so.2017.04 libopmjson.so.2017 libopmjson.so.2017.04 libopmoutput.so.2017 libopmoutput.so.2017.04 libopmparser.so.2017 libopmparser.so.2017.04 libopmsimulators.so.2017 libopmsimulators.so.2017.04 +-----+ Share - cmake setups : $ ls /usr/share | grep "dune*" dune dune-common dune-geometry +-----+ Tests : /home/bill/OPM/opm-data-master ************************** 24Jul2017 09:37 General information on OPM & Flow $ flow >> Information stdout : +-----+ * This is flow_ebos (version 2017.04) * * (eWoms version: 2017.04.0 ("Fabienne")) * * * * Flow is a simulator for fully implicit three-phase black-oil flow, * * and is part of OPM. For more information see: * * http://opm-project.org +-----+ >> in other words - just for black oil simulation!?!?! http://opm-project.org/?page_id=19 OPM Flow is a fully-implicit, black-oil simulator capable of running industry-standard simulation models. The simulator is implemented using automatic differentiation to enable rapid development of new fluid models. web-search "Open Porous Media and Flow and gas well simulation" http://www.uis.no/getfile.php/IOR-senter/Force_15102015.pdf IRIS -> OPM - Open Reservoir Simulator, 15Oct2015 +-----+ OPM Statistics - 24,014 commits made by 92 contributors representing 1,253,466 lines of code - is mostly written in C++ with well-commented source code - has a well-established mature codebase maintained by a very large development team with stable Y-O-Y commits - took an estimated 346 years** of effort (COCOMO model*) starting with its first commit in February, 2008 ending with its most recent commit 11 days ago * Howell - cool, see https://en.wikipedia.org/wiki/COCOMO **Howell : this does NOT include a MUCH larger effort on its pre-requisites, notably BLAS, LAPACK, SuiteSparse, Eigen, dune, etc, etc Norne comparison > 44431 active cells › 10 years of historical injection and production rates › 8 injectors and 28 producers › End-point scaling › Hysteresis >> Howell : great graphs of Norne results.. match pretty exactly Eclipse (also Flow-polymer & Flow-solvent) Norne comparison: Summary › Results match Eclipse › Currently OPM-FLOW is 3 times slower than Eclipse on Norne Sponsors & Contributors : Total, Statoil, Research Council of norway, Universitas Bergensis, CIPR, SINTEF, CLIMIT, IOR Norway +-----+ https://www.sintef.no/en/software/open-porous-media-opm2/ SINTEF is the largest independent research organisation in Scandinavia. 2,000 employees - electrical power, oil etc Norweigian, I think Atgeirr Flø Rasmussen Senior Research Scientist +47 977 48 679 Atgeirr.Rasmussen@sintef.no ************************** 24Jul2017 10:08 Where is the documentation for flow? http://opm-project.org/apidoc/latest/opm-core/html/ very general information Most opm modules have nothing ... http://resinsight.org/docs/home/ Fair documentation? but not what I need now >> I think I'll have to use the full OPM system build for documentation? ************************** 24Jul2017 10:08 web-search "gas wells and Open Porous Media" Images https://www.google.ca/imgres?imgurl=https%3A%2F%2Fimage.slidesharecdn.com%2Fs87brcads7ofqxhwnldm-signature-720151399eac53501c9292693be910dab2566b144a141bf5038f1635ae5eecc3-poli-150323070829-conversion-gate01%2F95%2Fcasing-design-1-638.jpg%3Fcb%3D1427112612&imgrefurl=https%3A%2F%2Fwww.slideshare.net%2Fnarendrakumard%2Fcasing-design&docid=n3dVFVa07FlxjM&tbnid=MfywI_tL2w2gmM%3A&vet=10ahUKEwj1sojRqaLVAhWJez4KHa7BBTkQMwhLKBwwHA..i&w=638&h=826&hl=en&client=browser-ubuntu&bih=957&biw=917&q=gas%20wells%20and%20Open%20Porous%20Media&ved=0ahUKEwj1sojRqaLVAhWJez4KHa7BBTkQMwhLKBwwHA&iact=mrc&uact=8 >> reasons for casing ... https://www.comsol.com/cfd-module CFD Module - Simulation Software for All Fluid Flow Applications https://www.cs.auckland.ac.nz/~john-rugis/pdf/JR_NZGeo2013.pdf J Rugis, P Leary, P Malin, J Pogacnik "GEOCRITICAL RESERVOIR FLOW SIMULATION AND DISPLAY USING OPEN POROUS MEDIUM (OPM) CODE" New Zealand Geothermal Workshop 2013 Proceedings, 17-20 November 2013, Rotorua, New Zealand Institute of Earth Science & Engineering, University of Auckland, New Zealand j.rugis@auckland.ac.nz >> Interesting paper... http://fortune.com/2015/09/08/big-oil-schlumberger/ >> neat article... WAIT!! : https://www.sintef.no/contentassets/8af8db2e42614f7fb94fb0c68f5bc256/mrst-book-2016.pdf Knut-Andreas Lie 13Dec2016 "An Introduction to Reservoir Simulation Using MATLAB : User Guide for the Matlab Reservoir Simulation Toolbox (MRST)" /home/bill/OPM/References/Lie 13Dec2016 An Introduction to Reservoir Simulation Using MATLAB.pdf ... For more computationally challenging cases, we recommend using the open-source Flow simulator from the Open Porous Media (OPM) initiative. +-----+ Part II Single-Phase Flow ... 7.2 A simulator based on automatic differentiation We will now present step-by-step how you can use the AD class in MRST to implement an implicit solver for the compressible, single-phase continu- ity equation (7.1). ... 7.2.3 Well model The production well will appear as a source term in the pressure equation. We therefore need to define an expression for flow rate in all cells in which the well is connected to the reservoir (which we will refer to as well connections). Inside the well, we assuming instantaneous flow so that the pressure drop is always hydrostatic. For a horizontal well, the hydrostatic term is zero and could obviously be disregarded, but we include it for completeness and as a robust precaution, in case we later want to reuse the code with a different well path. ... p226 - example results Can I get their Matlab code for the example? >> probably in : In the modules, you will find many tutorial examples that explain and showcase how the toolbox can be used to make general or fit-for-purpose simulators and workflow tools. +-----+ https://link.springer.com/chapter/10.1007%2F978-3-540-68783-2_9 Jørg E. Aarnes, Tore Gimse, Knut-Andreas Lie "An Introduction to the Numerics of Flow in Porous Media using Matlab", in Editors: Geir Hasle, Knut-Andreas Lie, Ewald Quak 2007 "Geometric Modelling, Numerical Simulation, and Optimization : Applied Mathematics at SINTEF", Springer pp 265-306 >> Seems similar to the user guide above... https://www.sintef.no/projectweb/mrst/ The MATLAB Reservoir Simulation Toolbox (MRST) is developed by the Computational Geosciences group in the Department of Applied Mathematics at SINTEF ICT. Version 2017a was released on the 15th of May 2017, and can be downloaded under the terms of the GNU General Public License (GPL). https://www.mathworks.com/pricing-licensing.html?prodcode=ML Matlab - cost for the standard individual user system : USD 2,150 Neural Network package 1,000 $US Home use - Is not for government, academic, commercial, or other organizational use MATLAB USD 149, Add-On Products USD 45 .. http://www.dumux.com/documents/dumux_awrpaper.pdf B. Flemisch, M. Darcis, K. Erbertseder, B. Faigle, A. Lauser, K. Mosthaf, S. M ̈uthing, P. Nuske, A. Tatomir, M. Wolff, R. Helmig 02Aug2010 "DuMux : DUNE for Multi-{Phase, Component, Scale,Physics, ...} Flow and Transport in Porous Media" Universiẗat Stuttgart, Preprint submitted to Advances in Water Resources >> Good DuMux description ************************** 24Jul2017 15:48 Descriptions of OPM standard tests equilibrium flow_diagnostic_test msw_2d_h norne polymer_simple2D polymer_test_suite README.md reduced_norne solvent_test_suite spe1 spe10model1 spe10model2 spe3 spe5 spe9 wells_test_suite http://petrowiki.org/Reservoir_simulation Awesome description and historical context, solution [types,approaches] http://petrowiki.org/images/3/35/Vol5_Page_1400_Image_0001.png Table 17.1 - SPE comparative solution project problems SPE1 Three-phase black oil, 10*10*03, 300-block grid, 3,650-day depletion with gas injection SPE2 Three-phase black oil, 10*01*15, 150-block r-z grid, 900-day single-well coning depletion SPE3 Nine-component retrograde gas, 09*09*04, 324-block grid, 5,480-day cycling and blowdown SPE4 Cyclic steam injection and steam displacement of heavy oils SPE5 Six-component volatile oil, 07*07*03, 147-block grid, 20-year WAG injection SPE6 Three-phase black oil, Single-block and cross-sectional dual porosity with drainage and gas and water injection cases SPE7 Three-phase black oil, 09*09*06, 486-block grid, with horizontal wells, Eight 1,500-day injection-production cases SPE8 Two -phase black oil, 10*10*04, 400-block grid, 2,500-day compared to 20-block unstructured and locally refined grid results SPE9 Three-phase black oil, 24*25*15, 9,000-block grid, 900-day depletion 25-well with geostatistical description SPE10 Model1 Two -phase gas-oil, 100*01*20, 2,000-block grid, 2,000-day and gas injection to 2,000 days SPE10 Model2 Two -phase water-oil, 60*220*85, 1.12M-block grid, 2,500-day and water injection to 2,000 days Both SPE10 models have geostatistical descriptions ... http://opm-project.org/?page_id=287 19Dec2016 The Norne case is an invaluable benchmark model of a real oil field in the Norwegian Sea, that is available under an open license. A new tutorial shows how to run that case with Flow, and also explains the output and visualization of the results. Norne comparison > 44431 active cells › 10 years of historical injection and production rates › 8 injectors and 28 producers › End-point scaling › Hysteresis ... msw_2d_h -- geophysics and fluid properties are copied from Norne >> Howell - not familiar enough yet with data files to deduce problem specs ... polymer_simple2D ************************** 25Jul2017 12:08 DuMux part of Dune, http://www.dumux.org/ NOT mentioned in : https://en.wikipedia.org/wiki/Reservoir_simulation http://wiki.aapg.org/Reservoir_simulation_study Black oil fluid descriptions are used to describe most oil and gas fields. Primary depletion, waterflooding, and gas injection can all be simulated with black oil models. Volatile oil reservoirs or gas condensate reservoirs generally require compositional models. These models may exhibit such complexities as a fluid whose density is linearly proportional to depth or whose phase switches repeatedly between oil and gas. Thermal models, used to simulate steam injection, may use either black oil or compositional fluid descriptions. Black oil thermal models describe fluid properties as a function of temperature as well as pressure. The importance of oil volatilization in thermal recovery often dictates that compositional models are used to simulate thermal recovery processes. http://www-m2.ma.tum.de/download-publications/cachan_flemisch.pdf Bernd Flemisch, Jochen Fritz, Rainer Helmig, Jennifer Niessner, Barbara Wohlmuth "DuMux: a multi-scale multi-physics toolbox for flow and transport processes in porous media" ECCOMAS Thematic Conference on Multi-scale Computational Methods for Solids and Fluids, A. Ibrahimbegovic, F. Dias, H. Matthies, P. Wriggers (eds.) Cachan, France, November 28-30, 2007 Flemisch etal 28Nov2007 DuMux, a multi-scale multi-physics toolbox for flow and transport processes in porous media Jumbled : https://www.g2crowd.com/categories/oil-and-gas-simulation-and-modeling Best Oil and Gas Simulation and Modeling Software ************************** 25Jul2017 18:07 OPM *.DATA format = Eclipse input format Øystein Pettersen , Dept. of Mathematics, Univ. of Bergen, 2006 "Basics of Reservoir Simulation With the Eclipse Reservoir Simulator" Lecture Notes, http://folk.uib.no/fciop/index_htm_files/ResSimNotes.pdf Pettersen 2006 Basics of Reservoir Simulation With the Eclipse Reservoir Simulator >> Awesome!! +-----+ Well Completions (COMPDAT keyword) Wnm* I_C J_C K_C,TOP K_C,BTM Open/ShutFlag SatTblNbr N/A Dw Kh Skin N/A Dir r0 Wnm* Well name as defined in WELSPECS, or well name wild card. Eclipse will interpret a trailing * as a wildcard, e.g. OP-1* will mean all wells with well names starting with “OP-1”. (Will match OP-1A, OP-11, OP-100B, OP-12,..., but not OP10). Note that the * can only be used as a wildcard at the end of Wnm*. Well name wildcards are so useful that most experienced users choose the well names such that wildcards can be used in the most efficient manner. . I_C J_C I and J indices for connection cells. If defaulted (by setting I, J to 0) the values defined in WELSPECS will be used. . K_C,TOP K_C,BTM K-index (layer) for topmost and bottom cell for which completions (perforations) are defined in this record. All cells in the range (K_C,TOP, K_C,BTM) will be open to flow. E.g. I_C J_C K_C,TOP K_C,BTM = 5 20 6 10 means the cells, (5, 20, 6), (5, 20, 7), (5, 20, 8), (5, 20, 9), and (5, 20, 10) will be perforated. Observe that the syntax is tailored for vertical wells. . Open/ShutFlag Specifies if the connection(s) specified in this record are currently open or closed to flow. Choices are OPEN or SHUT. (There is also a third option, AUTO. Ref. Eclipse manual) . SatTblNbr By default Eclipse will use the relative permeability curves defined for the connection cell for well in/outflow calculations. So far we haven’t even covered use of multiple relative permeability regions, but still can appreciate that there may be good reasons for defining special relative permeabilities for well flow. We contend ourselves with a simple example. Assume a well is completed down to 1960 m, with OWC at 1970 m. Further assume that in the simulation grid, both OWC and the well are in the same cell. In that case, the cell will contain mobile water such that the well will produce water immediately. Obviously in reality, the well will not produce water before the OWC has risen to the well completion level. We could model this behaviour by defining an artificially high S_wc in the well cell, such that the water would be seen as immobile until saturation reached a level which corresponds to the OWC having risen to the perforation. . D_w Wellbore diameter at the connection. This property is defined here rather than in WELSPECS, to allow for the diameter to vary along the wellbore. . Kh Permeability thickness, i.e. the product K*h. This value is a measure for the production (or injection) capacity of the well. The default value is to take K as the cell permeability and h as the connection length, which is the cell thickness for a vertical well. If the well is perforated across the entire cell it is normally OK to use the default value. But often the well is only perforated in part of the cell thickness, and also the average cell permeability may not be representative for the permeability in the interval that is perforated. In such cases it is better to specify the Kh explicitly. . Skin The skin factor enters the well in/outflow calculations (see below), and can be specified here (default value is zero). . Dir The direction the well penetrates the cell. For a vertical well this is the Z-direction, for a horizontal well parallel to the x-axis X. The choices are X, Y, or Z. . r0 Pressure equivalent radius. Generally it is recommended to default this value, in which case the Peaceman equivalent radius will be used (see below). To use an alternative value, specify r0 > 0. +-----+ http://rgn.hr/~dvulin/PRLU/data/tut1/tutorial_1.pdf >> Eclipse tutorial - mostly GUI-based, so not useful for me http://petrofaq.org/wiki/Eclipse_Input_Data >> useful, but no horizontal well description ************************** 28Jul2017 10:08 Local Grid Refinements, Horizontal wells and friction see p79 of : Øystein Pettersen , Dept. of Mathematics, Univ. of Bergen, 2006 "Basics of Reservoir Simulation With the Eclipse Reservoir Simulator" Lecture Notes, http://folk.uib.no/fciop/index_htm_files/ResSimNotes.pdf Pettersen 2006 Basics of Reservoir Simulation With the Eclipse Reservoir Simulator +-----+ -- LGR : -- MAXLGR - Max nbr. of LGRs in model -- MAXCLLGR - Max nbr of cells in one LGR -- MAXAMLGC - Max nbr. of amalg. coarse cells -- MAXAMLGF - Max nbr. of LGR amalgamations -- MAXLGRAM - Max nbr. of LGRs in one amalgamation -- LSTACK - NSTACK for LGRs (NSTACK is defined on page 100) -- INTERP? - Pressure interpolation? (INTERP or NOINTERP (default)) LGR -- MAXLGR MAXCLLGR MAXAMLGC MAXAMLGF MAXLGRAM LSTACK INTERP? 1 500 1* 1* 1* 1* 1* / WELSPECS -- Wname Gname IWH JWH Z_BHP Phase N/A N/A ShutIn XFlow WH1 G1 3 2 1840.0 OIL 1& 1* STOP YES / / WELLSPECSL COMPDAT -- Wnm* I_C J_C K_C,TOP K_C,BTM Open/Shut SatTbl N/A Dw Kh Skin N/A Dir r0 WH1 3 2 2 2 OPEN 0 1* 0.3 1* 0 1* X 1* / WH1 4 2 2 2 OPEN 0 1* 0.3 1* 0 1* X 1* / WH1 5 2 2 2 OPEN 0 1* 0.3 1* 0 1* X 1* / WH1 6 2 2 2 OPEN 0 1* 0.3 1* 0 1* X 1* / / COMPDATL -- Wnm* WH1 WH1 WH1 WH1 / . WCONPROD ... Note that all references to cell indices are relative to the local numbering scheme in the LGR. Note: For complex well paths ECLIPSE sometimes gets the well path all wrong. The user can override the default well-path-construction by the keyword COMPORD, whereby the completion order can be set exactly as defined by the user. +-----+ Also, p85 15.4 Horizontal wells and friction ************************** 28Jul2017 11:32 Horizontal wells in OPM-Flow see "/home/bill/OPM/Howell sims/gas_horizontal/gas_horiz_10_10_10_CASE1.DATA" From my first attempt, it does NOT appear that OPM-Flow can handle horizontal wells. -- Howell 26Jul2017 19:49 bash script error "Program threw an exception: Malformed floating point number ''X''" -- Wnm I J K_hi K_lo Open/Shut SatTblNbr Skin diam_wellbore Kh Skin Dir r0 -- 'PROD' 10 10 3 3 'OPEN' 1* 1* 0.5 1* 0 'X' 1* / -- 'INJ' 1 1 1 1 'OPEN' 1* 1* 0.5 1* 0 'X' 1* / Try again, using : COMPDAT -- Wnm I J K_hi K_lo Open/Shut SatTblNbr Skin diam_wellbore Kh Skin Dir r0 'PROD' 10 10 3 3 'OPEN' 1* 1* 0.5 1* 0 X 1* / 'INJ' 1 1 1 1 'OPEN' 1* 1* 0.5 1* 0 X 1* / / Using "/home/bill/OPM/Howell sims/0_run tests.sh" : $ bash "/home/bill/OPM data/Howell sims/0_run tests.sh" >> This time it ran OK (no "bash script error "Program threw an exception: Malformed floating point number ''X''"" $ cd "/home/bill/OPM/Howell sims/gas_horizontal/" $ python summaryplot WBHP:INJ WWIR:INJ WBHP:PROD WOPR:PROD WGPR:PROD gas_horiz_10_10_10_CASE1.DATA OOPS - bad setup - lost my changes for horizontal!?? $ bash "/home/bill/OPM data/Howell sims/0_run tests.sh" $ python summaryplot WBHP:INJ WWIR:INJ WBHP:PROD WOPR:PROD WGPR:PROD gas_horiz_10_10_10_CASE1.DATA >> All flat-lined at zero!?? >> OOPS - COMPDAT was commented out! $ bash "/home/bill/OPM data/Howell sims/0_run tests.sh" >> In "/home/bill/OPM data/Howell sims/gas_horizontal/0_gas_horiz_10_10_10_CASE1 - variables & run messages 170728 11h54m.txt" : Program threw an exception: Malformed floating point number 'X' $ python summaryplot WBHP:INJ WWIR:INJ WBHP:PROD WOPR:PROD WGPR:PROD gas_horiz_10_10_10_CASE1.DATA >> All flat-lined at zero! CONCLUSION : OPM-FLow won't do horizontal wells !! However - there is a difference between Petterson's variable listing and OPM-Flow examples, so try finding a modern Eclipse *.DATA file for horizontal wells and check >> Can't find +-----+ Try adding a new N/A column as per Petterson example, which DIFFERS from his listing of variables! : COMPDAT -- Wnm* I_C J_C K_C,TOP K_C,BTM Open/Shut SatTbl N/A Dw Kh Skin N/A Dir r0 WH1 3 2 2 2 OPEN 0 1* 0.3 1* 0 1* X 1* / WH1 4 2 2 2 OPEN 0 1* 0.3 1* 0 1* X 1* / WH1 5 2 2 2 OPEN 0 1* 0.3 1* 0 1* X 1* / WH1 6 2 2 2 OPEN 0 1* 0.3 1* 0 1* X 1* / / . I used : COMPDAT -- Wnm I J K_hi K_lo Opn/Sht SatTbl N/A D_w Kh Skin N/A Dir r0 PROD 10 10 3 3 'OPEN' 1* 1* 0.5 1* 0 1* X 1* / INJ 1 1 1 1 'OPEN' 1* 1* 0.5 1* 0 1* X 1* / / 28Jul2017 13:00 $ cd "/home/bill/OPM data/Howell sims/gas_horizontal/" $ bash "/home/bill/OPM data/Howell sims/0_run tests.sh" >> Results (special notes) : seems OK - at least there is no error this time & it ran $ python summaryplot WBHP:INJ WWIR:INJ WBHP:PROD WOPR:PROD WGPR:PROD gas_horiz_10_10_10_CASE1.DATA >> Results (special notes) : Seems sort-of-OK, WOPR:PROD is very unstable (like slugs of oil) $ bash "/home/bill/OPM data/ResInsight bash.sh" >> Results (special notes) : NUTS! wells are vertical, only go down to third level of 10 Fix horizontal : COMPDAT -- Wnm I J K_hi K_lo Opn/Sht SatTbl N/A D_w Kh Skin N/A Dir r0 PROD 10 10 3 3 OPEN 1* 1* 0.5 1* 0 1* X 1* / PROD 9 10 3 3 OPEN 1* 1* 0.5 1* 0 1* X 1* / PROD 8 10 3 3 OPEN 1* 1* 0.5 1* 0 1* X 1* / PROD 7 10 3 3 OPEN 1* 1* 0.5 1* 0 1* X 1* / PROD 6 10 3 3 OPEN 1* 1* 0.5 1* 0 1* X 1* / PROD 5 10 3 3 OPEN 1* 1* 0.5 1* 0 1* X 1* / PROD 4 10 3 3 OPEN 1* 1* 0.5 1* 0 1* X 1* / PROD 3 10 3 3 OPEN 1* 1* 0.5 1* 0 1* X 1* / PROD 2 10 3 3 OPEN 1* 1* 0.5 1* 0 1* X 1* / PROD 1 10 3 3 OPEN 1* 1* 0.5 1* 0 1* X 1* / INJ 1 1 1 10 OPEN 1* 1* 0.5 1* 0 1* X 1* / / . $ cd "/home/bill/OPM data/Howell sims/gas_horizontal/" $ bash "/home/bill/OPM data/Howell sims/0_run tests.sh" >> Results (special notes) : Program threw an exception: Malformed floating point number 'STOP' >> Oops - ill-fdefined variable list problem : WELSPECS -- Wname group I_WH J_WH Z_BHP pref_phase Auto_shut XFlow DensCalc PROD G1 10 10 8775 GAS STOP YES 1* / INJ G1 1 1 8775 WATER STOP YES 1* / >> add two defaults : -- Wname group I_WH J_WH Z_BHP Phase N/A N/A Auto_shut XFlow DensCalc PROD G1 10 10 8775 GAS 1* 1* STOP YES 1* / INJ G1 1 1 8775 WATER 1* 1* STOP YES 1* / Retry : edit "/home/bill/OPM data/Howell sims/0_run tests.sh" to unmask "NOSIM" - check dataset $ cd "/home/bill/OPM data/Howell sims/gas_horizontal/" $ bash "/home/bill/OPM data/Howell sims/0_run tests.sh" >> Results for dataset consistency (special notes) : Program threw an exception: Malformed integer ''PROD'' >> next Program threw an exception: Malformed integer ''OPEN'' >> oops, wellname disappeared... fixed Program threw an exception: Malformed integer ''PROD'' >> wow! hard to find the problem with this, unless my multiple rows of PROD for COMPDAT is the problem >> remove the extra rows temporarily and see what happens Program threw an exception: Malformed integer ''PROD'' >> that didn't work, so put removed data back and look elsewhere again Program threw an exception: Malformed integer ''PROD'' 28Jul2017 14:53 take a break - can't find offending "PROD"!! 28Jul2017 17:56 >> Reverse the order [1 to 10] for COMPDAT COMPDAT -- Wnm I J K_hi K_lo Opn/Sht SatTbl N/A D_w Kh Skin N/A Dir r0 PROD 1 10 3 3 OPEN 1* 1* 0.5 1* 0 1* X 1* / PROD 2 10 3 3 OPEN 1* 1* 0.5 1* 0 1* X 1* / PROD 3 10 3 3 OPEN 1* 1* 0.5 1* 0 1* X 1* / PROD 4 10 3 3 OPEN 1* 1* 0.5 1* 0 1* X 1* / PROD 5 10 3 3 OPEN 1* 1* 0.5 1* 0 1* X 1* / PROD 6 10 3 3 OPEN 1* 1* 0.5 1* 0 1* X 1* / PROD 7 10 3 3 OPEN 1* 1* 0.5 1* 0 1* X 1* / PROD 8 10 3 3 OPEN 1* 1* 0.5 1* 0 1* X 1* / PROD 9 10 3 3 OPEN 1* 1* 0.5 1* 0 1* X 1* / PROD 10 10 3 3 OPEN 1* 1* 0.5 1* 0 1* X 1* / INJ 1 1 1 10 OPEN 1* 1* 0.5 1* 0 1* X 1* / / >> Also - Dir=X, so is that consistent with changing I? (bad notation in Eclipse format) Petterson p9 "... Eclipse convention: I, J, K are used for indices in X, Y, Z – directions, and the ordering is defined as in “normal reading order”, i.e. I runs from left to right, J runs from top of a page to bottom, and K runs from page to page. In a grid oriented along standard geographical axes, I would increase from west to east, J would increase from north to south, and K downwards. ..." $ cd "/home/bill/OPM data/Howell sims/gas_horizontal/" $ bash "/home/bill/OPM data/Howell sims/0_run tests.sh" >> Results for dataset consistency (special notes) : Program threw an exception: Malformed integer ''PROD'' >> change WELSPECS by removing "DensCalc" values ***** I need to know data input conventions - what are the full inputs, etc (FORTRAN-like) ***** Change : WELSPECS -- Wname group I_WH J_WH Z_BHP Phase N/A N/A Auto_shut XFlow DensCalc PROD G1 10 10 8775 GAS 1* 1* STOP YES 1* / INJ G1 1 1 8775 WATER 1* 1* STOP YES 1* / / To : WELSPECS -- Wname group I_WH J_WH Z_BHP Phase N/A N/A Auto_shut XFlow DensCalc PROD G1 10 10 8775 GAS 1* 1* STOP YES / INJ G1 1 1 8775 WATER 1* 1* STOP YES / / $ cd "/home/bill/OPM data/Howell sims/gas_horizontal/" $ bash "/home/bill/OPM data/Howell sims/0_run tests.sh" >> Results for dataset consistency (special notes) : Program threw an exception: Malformed integer ''PROD'' >> So reverse that last change - done >> Move last "/" of COMPDAT to start of line $ cd "/home/bill/OPM data/Howell sims/gas_horizontal/" $ bash "/home/bill/OPM data/Howell sims/0_run tests.sh" >> Results for dataset consistency (special notes) : Program threw an exception: Malformed integer ''PROD'' 28Jul2017 18:56 STOP here for the day - I'm really stuck! ************************* 29Jul2017 06:07 Bug-search in *.DATA files Try inserting "BREAK-STOP" and silence with "-- BREAK-STOP" Right near START : $ cd "/home/bill/OPM data/Howell sims/gas_horizontal/" $ bash "/home/bill/OPM data/Howell sims/0_run tests.sh" >> Results for dataset consistency (special notes) : Program threw an exception: PARSE_RANDOM_TEXT: String 'BREAK-STOP' not formatted/recognized as valid keyword at: "/home/bill/OPM data/Howell sims/gas_horizontal/gas_horiz_10_10_10_CASE1.DATA":72 BINGO!! Just before COMPDAT $ cd "/home/bill/OPM data/Howell sims/gas_horizontal/" $ bash "/home/bill/OPM data/Howell sims/0_run tests.sh" >> Results for dataset consistency (special notes) : Program threw an exception: PARSE_RANDOM_TEXT: String 'BREAK-STOP' not formatted/recognized as valid keyword at: "/home/bill/OPM data/Howell sims/gas_horizontal/gas_horiz_10_10_10_CASE1.DATA":602 Just after COMPDAT $ cd "/home/bill/OPM data/Howell sims/gas_horizontal/" $ bash "/home/bill/OPM data/Howell sims/0_run tests.sh" >> Results for dataset consistency (special notes) : Program threw an exception: PARSE_RANDOM_TEXT: String 'BREAK-STOP' not formatted/recognized as valid keyword at: "/home/bill/OPM data/Howell sims/gas_horizontal/gas_horiz_10_10_10_CASE1.DATA":617 >> OK - so it's not the horizontal well specs that does it! Just before WCONPROD $ cd "/home/bill/OPM data/Howell sims/gas_horizontal/" $ bash "/home/bill/OPM data/Howell sims/0_run tests.sh" >> Results for dataset consistency (special notes) : Program threw an exception: Malformed integer 'WCONPROD' >> Wow! This doesn't make sense yet... Ah Hah! Just after the COMPDAT declaration, I had : -- +-----+ -- Howell 28Jul2017 13:51 from Petterson example on local grid - perhaos way to define horizontal well? COMPDAT -- Wnm* I_C J_C K_C,TOP K_C,BTM Open/Shut SatTbl N/A Dw Kh Skin N/A Dir r0 >> which I changed to : -- +-----+ -- Howell 28Jul2017 13:51 from Petterson example on local grid - perhaos way to define horizontal well? -- COMPDAT -- Wnm* I_C J_C K_C,TOP K_C,BTM Open/Shut SatTbl N/A Dw Kh Skin N/A Dir r0 >> Try again : $ cd "/home/bill/OPM data/Howell sims/gas_horizontal/" $ bash "/home/bill/OPM data/Howell sims/0_run tests.sh" >> Results for dataset consistency (special notes) : Program threw an exception: PARSE_RANDOM_TEXT: String 'BREAK-STOP' not formatted/recognized as valid keyword at: "/home/bill/OPM data/Howell sims/gas_horizontal/gas_horiz_10_10_10_CASE1.DATA":684 >> So it stops Just before WCONPROD Neutralize -- BREAK-STOP $ cd "/home/bill/OPM data/Howell sims/gas_horizontal/" $ bash "/home/bill/OPM data/Howell sims/0_run tests.sh" >> Results for dataset consistency (special notes) : ***************Saturation Functions Diagnostics*************** System: Black-oil system. Relative permeability input format: Saturation Family I. Number of saturation regions: 1 ================ Simulation turned off =============== >> So no data setup errors flagged edit "/home/bill/OPM data/Howell sims/0_run tests.sh" to mask "-- NOSIM" - do NOT check dataset, normal run $ cd "/home/bill/OPM data/Howell sims/gas_horizontal/" $ bash "/home/bill/OPM data/Howell sims/0_run tests.sh" >> Results for dataset consistency (special notes) : Seems OK!! $ python summaryplot WBHP:INJ WWIR:INJ WBHP:PROD WOPR:PROD WGPR:PROD gas_horiz_10_10_10_CASE1.DATA >> Results (special notes) : Seems OK $ bash "/home/bill/OPM data/ResInsight bash.sh" >> Results (special notes) : HOORAY!! It works - I get a horizontal well!!! # enddoc