SOFTWARE RELEASE NOTICE 1 PJL 2. Project Title: 3. SRN Title: GM-SYSQ Version 4.8.45b I Project No. Name Read Onlymead- Write Addition/C hangemele te 4. Originator/Requestor: Brandi L. Winfrey II Date: November 30, 2004 5. Sumrnary of Actions X Release of new software 0 Change of access software 0 Release of modified software: 0 Software Retirement 0 Enhancements made 0 Corrections made 6. Validation Status X Validated 0 Lirnited Validation Not Validated Explain: ALL YES NO 8. Element Manager Approval: Date: 8/ 30 -6 n 9. Remarks: Software is licensed, commercial, acqYlfed software not to be modified. ll CNWRA Form TOP-6 (09101)
33
Embed
PJL Project No. 1 · by comparing the model’s gravity and magnetic response to observed measurements. The GM-SYSB 3D Modeling extension to Oasis montaj enables you to design three
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
SOFTWARE RELEASE NOTICE
1 PJL 2. Project Title:
3. SRN Title: GM-SYSQ Version 4.8.45b I
Project No.
Name Read Onlymead- Write
Addi tion/C hangemele te
4. Originator/Requestor: Brandi L. Winfrey II Date: November 30, 2004
5. Sumrnary of Actions
X Release of new software 0 Change of access software
0 Release of modified software: 0 Software Retirement
0 Enhancements made
0 Correc t ions made
6 . Validation Status
X Validated
0 Lirnited Validation
Not Validated Explain:
ALL YES NO
8. Element Manager Approval: Date: 8 / 30 -6 n
9. Remarks: Software is licensed, commercial, acqYlfed software not to be modified. ll
CNWRA Form TOP-6 (09101)
SOFTWARE SUMMARY FORM
01. Summary Date: 1 113 0/2004
04. Software
02. Summary prepared by Brandi L. Winfrey (210)522-5083
05. Short Title: GM-SYSB
03. Summary Action: NEW
07. Internal Software ID: None I 06. Software Title: GM-SYSB Version 4.8.45b
08. Software Type: 09. Processing Mode: 10. Application Area
Automated Data System
II I www.iiga.com
X Computer Program
0 Subroutinehlodule
~ 11. Submitting Organization and Addrcss:
X Interactive
0 Batch
0 Combination
a. General: X ScientificEngineering X Auxiliary Analyses 0 Total System PA 0 Subsystem PA 0 Other
b. Specific: geophysical and magnetic modeling software
12. Technical Contact(s) and Phone: Northwest G eo ph y s i c a1 A s soc i ate s , In c . 1600 S\V Western Boule\ a d , Suitc 300 PO B o x 1063 Corvallis, OK 97333 USA Phone: (541) 757-723 1 Fax: (541) 757-7331
15. Computer Operating System: Windows
XP recommended, Windows 2000
16. Programming Language( s):
N/A
11 13. Software Application:
18. Computer Memory Requirements: 512 MB or more recommended, 128 MB minimum
NGA’s GM-SYSB Profile Modeling and 3D Modeling software is provided by Geosoft as extensions to Oasis montaj. The GM-SYSB Profile Modeling extension to Oasis montaj enables you to create a geologic model and test its accuracy by comparing the model’s gravity and magnetic response to observed measurements. The GM-SYSB 3D Modeling extension to Oasis montaj enables you to design three dimensional models capable of accurately depicting the variation and irregularity of sub surf ace structures .
I I
19. Tape Drives: N/A 20. Disk Units: 400 MB (depending upon software configuration)
14. Computer Platform Windows
23. Software Availability:
X Available 0 Limited 0 In-House ONLY 24. Documentation Availability: X Available 0 Preliminary 0 In-House ONLY
17. Number of Source Program Statements:
N/A
2 1. Graphics: 24-bit graphics card with 3-D acceleration is recommended and required for full color imaging. Recommend 64MB RAM on card.
The physical labeling of software medium (tapes, disks, etc.) contains: Program Name, Module/Name/Title, Module Revision, File type (ASCII, OBJ, EXE), Recording Date, and Operating System(s)?
Yes: 0 No:O N/A:O
(04/0 1 ) Page 1 of 3
CENTER FOR NUCLEAR WASTE REGULATORY ANALYSES QA VERIFICATION REPORT
FOR +ACQUIRED SOFTWARE NOT TO BE MODIFIED C
User Documentation [TOP-01 8, Section 5.5.71 ~~
Is there a Users’ Manual for the software and is it up-to-date? Yes: 0 No:O N/A:O
User’s Manual Version and Date: Comments:
~
Are there basic instructions f a the installation and use of the software? Yes: 0 No:O N/A:O
Location of Instructions: Comments:
Configuration Control [TOP-01 8, Section 5.7,5.9.3] ~
Is the Software Summary Form (Form TOP-4-1) completed and signed? Yes: E/ NO: o N/A: o
Date of Approval:
Is the list of files attached to the Software Summary Form complete and accurate? Yes: 0 No:O N/A:O
Comments:
Is the source code available x, is the esecucable cme zvaihbie i3 the case of (acquiredkommercial codes)? Yes: 0 No:O N/A:D
Location of Source Code:
Have all the scripdmake files and executable files been submitted to the Software Custodian?
Only the executable files are being submitted.
Yes: 0 No:O N/A:U
Location of executable files: Comments:
(04/01) Page2of 3
CENTER FOR NUCLEAR WASTE REGULATORY ANALYSES QA VERIFICATION REPORT
FOR +ACQUIRED SOFTWARE NOT TO BE MODIFIED C
Software Release [TOP-018, Section 5.91
Upon acceptance of the software as verified above, has a Software Release Notice (SRN), Form TOP-6 been I issued and does the version number of the software match the documentation? I SRNNumber: % Comments: I
Yes: W' N ~ : O N/A:O
~ ~ ~ ~~~
L a r e V a l i d a t i o n [TOP-018, Section 5.101
Has a Software Validation Test Plan (SVTP) been prepared for the range of applicatio f the software? I Yes: d o : 0 N/A:O 6.1 I Version and Date of SVTP: 7 / % d P !
Date Reviewed and Approved via QAP-002: &q/dJ
Comments: L h r b t d 3 ~ T f / G l / r Q
Has a Software Validation Test Report (SVTR) been prepared that documents the results of the validation cases, interpretation of the results, and determination if the software has been validated?
Version and Date of SVTR: .7/2 719
E w Reviewed ar?d Appprod vi2 QAP-002: 127,/~9
Comments. :
I Yes: d NO: o N/A:O &.fa
(04/0 1 ) Page 3 of 3
SOFTWARE VALIDATION TEST PLAN AND REPORT GM-SYS, Version 4.8.45b
ASSUMPTIONS AND CONSTRAINTS .................................................................................. 2
TEST CASES ......................................................................................................................... 3
6.1 Gravity and Magnetic Response Over a Horizontal Cylinder of Infinite Strike ............ 3 Objective .......................................................................................................... 3 6.1 . 1
6.1.2 Test Input ......................................................................................................... 3 6.1.3 Test Procedure ................................................................................................ 3 6.1.4 Test Results ..................................................................................................... 5 Gravity and Magnetic Response Over a Semi-Infinite Horizontal SlabNertical Fault . 6 6.2.1 Objective .......................................................................................................... 6 6.2.2 Test Input ......................................................................................................... 6 6.2.3 Test Procedure ................................................................................................ 6 6.2.4 Test Results ..................................................................................................... 8
6.2
APPENDIX A APPENDIX B
i
Figure
FIGURES
Page
Figure 1 . GM-SYS Infinite Cylinder Model and Corresponding Gravity and Magnetic Response ..................................................................................................................... 4
Figure 2 . Gravity Response Over a Horizontal Cylinder of Infinite Strike .................................... 5 Figure 3 . Magnetic Response Over a Cylinder of Infinite Strike .................................................. 6 Figure 4 . GM-SYS Horizontal SlabNertical Fault Model and Corresponding Gravity and
Magnetic Response ...................................................................................................... 7 Figure 5 . Gravity Response Over a Semi-Infinite Horizontal SlabNertical Fault ......................... 9 Figure 6 . Magnetic Response Over a Semi-Infinite Horizontal SlabNertical Fault .................... 10
TABLES
Table Page
Table I . Buried Cylinder Model Parameters ................................................................................ 4 Table 2 . Horizontal SlabNertical Fault Model Parameters .......................................................... 8
iii
1 SCOPE OF VALIDATION
GM-SYS@, developed by Northwest Geophysical Associates, Inc. (NGA) for the OASIS montaj" software platform, is a program for calculating the gravity and magnetic response from a geologic model. GM-SYS provides an easy-to-use interface for interactively creating and manipulating models to fit observed gravity and/or magnetic data. Calculated responses change instantaneously as the model is changed, speeding the interpretation process and allowing quick testing of alternative solutions. Models and data can be imported from ASCII files, entered by hand, or extracted interactively from an OASIS montaj" database.
The software validation tests compare two-dimensional (2D) gravity and magnetic model responses of simple-shaped bodies buried in a homogeneous sediment obtained by (1) GM-SYS and (2) independent calculations using standard equations from the literature. If the GM-SYS and independent calculation gravity and magnetic responses are similar, then the GM-SYS code shall be considered validated and thus usable for the intended regulatory reviewpurposes.
2 REFERENCES
Hinze, W. J., 1990, The role of gravity and magnetic methods in engineering and environmental studies; in Ward, S. H., Ed., Geotechnical and environmental geophysics, Vol. I: Review and Tutorial; Investigations in Geophysics No. 5, SOC. Explor. Geophys., 75-1 26.
Jiracek, G. R, 1994, Gravity Method, in Hydrogeophysics class notes.
The following is a list of registered trademarks used in this document: GM-SYS@ is a trademark of Northwest Geophysical Associates, Inc. OASIS montajTM is a trademark of Geosoft, Inc. PostScript@ is a trademark of Adobe Systems, Inc. Windows@ 95, Windows@ 98, Windows Nl@, Windows@ 2000and Microsoft@ Excel are registered trademarks of Microsoft Corporation. Solaris@ is a trademark or registered trademark of Sun Microsystems. Linux@ is a trademark of Linus Torvalds. IRIX@ is a registered trademark of Silicon Graphics, Inc. A l p is a trademark of the International Business Machines Corporation.
1
3 ENVIRONMENT
3.1 Software
Although it will run independently, GM-SYS is considered an Add-on Tool to the Oasis montaj software platform and thus is typically installed to run concurrently with Oasis montaj.
GM-SYS is cross-platform compatible, supporting Windows’ 95, Windows@ 98, Windows NT@, Windows’ 2000, Solaris@, Linux@, IRK@, and AIX’. However, the test version is installed on an x86 AT/AT compatible 300 MHz PC with 256 Mb RAM running Microsoft Windows NT@ 4.00.1 381.
GM-SYS installation requires the installer to transfer the ‘soft-key’ license protection from the supplied ‘red disk‘, a 3.5-inch diskette, to the local hard drive.
3.2 Hardware
Mini mum hardware requirements include:
16 Mb RAM, 32 Mb recommended. 16-24 bit graphics card recommended and required for full color; VGA minimum resolution, 8-bit (256 color) also supported. CD-ROM drive and 3.5-inch 1.44 HD diskette drive for software installation and ‘soft-key’ license protection transfer.
4 PREREQUISITES
Platform should have licensed version of GM-SYS installed and available for use.
5 ASSUMPTIONS AND CONSTRAINTS
Assume software is installed correctly and functional.
Assume users have a general knowledge of test platform and basic operation of associated a p pl ica t ions and pe ri p herals .
Assume users have general knowledge of basic forward modeling and potential field geophysics .
Assume users are familiar with or can become familiar with the GM-SYS interface.
2
6 TESTCASES
6.1 Gravity and Magnetic Response Over a Horizontal Cylinder of Infinite Strike
6.1 .I Objective
Enter the simple-body model (a horizontal cylinder of infinite strike buried in homogeneous sediment) into GM-SYS to determine the gravity and magnetic response. Independently calculate the gravity and magnetic response of the same simple-body model using standard equations from the literature. Tabulate and plot the GM-SYS and independent calculation responses to directly compare the two methods.
Requirements:
The GM-SYS and independent calculation gravity and magnetic responses should be similar, if not identical, to validate the GM-SYS code 'correctness'. Any differences in the model response (curves) should be explainable (e.g., model geometry).
6.1.2 Test Input
Use selected parameters to build a model in GM-SYS and for the independent calculations.
6.1.3 Test Procedure
A new model space was created in GM-SYS extending from -30,000,000 m to 30,000,000 m in the X direction and from 0 m to 500,000 m in the Z direction. The large model space was constructed to avoid edge effects or boundary interference on the computed responses. The Earth's magnetic field strength was set to 50,000 gammas with an inclination of 90.0 degrees and declination of 0 degrees to simulate vertical magnetic field incidence (i.e., North magnetic pole) and coincidence of the magnetic and geographic poles. The survey profile extended from -50 m to 50 m with gravity and magnetic stations every 1 m, for a total of 101 stations. For simplicity, the surface was assumed flat, thus topography was omitted.
A cylinder of infinite strike (circle), with 10 m radius, was created in the model space. The cylinder was centered at station 0 and buried 10 m below the surface such that the center of the cylinder was 20 m below the surface (Figure 1). The selected air, sediment, and cylinder properties (model parameters) are listed in Table 1. Model parameters were selected to provide known gravity and magnetic responses and allow straightforward, manual independent calculations, yet still provide a realistic validation of the algorithm used in the code.
3
The selected air, sediment, and horizontal slab properties (Le., model parameters) are listed in Table 2. Model parameters were selected to provide known gravity and magnetic responses and allow straightforward, manual independent calculations, yet still provide a realistic validation of the algorithms used in the code.
Density, 17 (gm/cm3) Magnetic Susceptibility , k ( emu/cm3)
Table 2. Horizontal SlabNertical Fault Model Parameters
GM-SYS calculated and plotted the gravity and magnetic response over the buried horizontal slab as shown in Figure 4. To establish a common match point, both response curves were shifted to zero at the left-most station (station -50 m). GM-SYS identifies the shifted response curves by placing a red star at the shift point (Figure 4). The shifted gravity and magnetic response curves were exported in tabular format and subsequently imported into Microsoft@ Excel to allow direct comparison with the independently calculated gravity and magnetic response curves.
0 0 0 0 0 0
Independent gravity response calculations were performed using equations presented in Hinze (1990) and Jiracek (1994) and the magnetic response calculations were performed using the equation presented in Hinze (1 990). Jiracek (1 994) did not have equations for magnetic calculations. All independent calculations were performed using Microsoft@ Excel. The gravity and magnetic equations, unit conversions, and tabulated results for the semi-infinite horizontal slab/vertical fault model are included in Appendix B.
6.2.4 Test Results
GM-SYS and independent calculation results were plotted coincidentally using the scatter plot graphing function in Microsoft@ Excel. The gravity and magnetic response curves are shown in Figures 5 and 6, respectively. The gravity response curves (Figure 5) are effectively coincident for all stations, indicating that the GM-SYS calculations and independently calculated solutions are nearly identical. The magnetic response curves show a slight misfit near the edge of the slab (Figure 6). The slight misfit is likely due to edge effects caused by the GM-SYS approximation of a semi-infinite horizontal slab. However, the overall fit and shape is nearly identical to the independent calculations, thus indicating a comparable response.
The 'gravity and magnetic response over a semi-infinite horizontal slab/vertical fault' test shows that the GM-SYS code is valid and correct, and the test results satisfy the requirements stated in Section 6.2.1. Thus GM-SYS is usable for the intended regulatory review purposes.
8
0.4
-0.05 ! 60 -40 -20 0 20
Position (m)
+GMSYS +Calculated (Jitacek) +Calculated (Hinze)
40
Figure 5. Gravity Response Over a Semi-Infinite Horizontal SlabNertical Fault
9
Appendix A
Infinite Horizontal Cylinder Response Calculations and Results
Appendix A
Infinite Horizontal Cylinder Response Calculations and Results
- CALCULATIONS - Vertical Gravity Component (gz) Response over a Horizontal Cylinder of Infinite Strike
Gravitational constant = G = 6.67E-08 cm3/g/s2
Jiracek: cylinder: radius of cylinder = R =
depth to axis of cylinder = z =
station position = x =
cylinder density, r2 =
earth density, r l =
density contrast, r =
0.01 km 0.02 km
-0.05 to 0.05 km 3 g/cm3
2 g/cm3 1 rr/cm3
k2 = 41.9 1 *((r*R2)/z) rkml
k2 = 0.20955 [km]
Hinze: cylinder:
radius of cylinder = R =
depth to axis of cylinder = d =
station position = x =
cylinder density, r2 =
earth density, r l =
density contrast, r =
1000 cm 2000 cm
-5000 to 5000 cm 3 g/cm3
2 g/cm3 1 g/cm3
gz = (2*o*G*r*d*R2)/(x2+d2) rmGall
GM-SYS Jiracek Hinze Max Amplitude = 0.172434 0.180646552 0.180641425 mGal Min Amplitude = 0 0 0 mGal
- CALCULATIONS - Magnetic Response over a Horizontal Cylinder of Infinite Strike
cylinder calculations:
Earth background field = H = Magnetic Susceptibility = k = 0.001 emu/cm3 = 0.0125664 SI
50000 nT =
Remanent Magnetization = o emu/cm3 = O T = OnT
0.5 oersted
radius of cylinder = R = depth to axis of cylinder = d =
0.01 km = 0.02 km =
1000 cm 2000 cm
useful relations:
1 Y = 1 nT = oersted 1 emu/cm3 = 4n x I O 4 T = I O 3 N m 1 nT = IO-’ T or 1 T = I O 9 nT = I O 9 Y k in cgs has units emu/cm3. k in SI is unitless k (SI) = 4nk (cgs)
Hinze: T = (2*n*k*H*R2*(d2-x2))/(~2+d2)2
T = magnetic anomaly [oersteds] k = magnetic susceptibility contrast [emu/cm31 H = vertical magnetic field [oersteds] R = cylinder radius [cm] d = depth to cylinder axis [cm] x = horizontal distance from observation point to axis [cm]
- CALCULATIONS - Magnetic Response over a Semi-Infinite Horizontal SlabNertical Fault
horizontal slab calculations:
Earth background field = H = Magnetic Susceptibility = k =
Remanent Magnetization = slab thickness = t =
depth to center of slab = d =
50000 nT = 0.5 oersted 0.001 emu/cm3 = 0.012566 SI
0.01 km = 1000 cm 0.01 km = 1500 cm
o emu/cm3 = O T = OnT
useful relations: 1 Y = 1 nT = I O ” oersted 1 emu/cm3 = 47~ x 1 n T = 1 0 ~ 9 T o r 1 T = 1 0 9 n T = 1 0 9 Y k in cgs has units emu/cm3. k in SI is unitless k (SI) = 47Ck (cgs)
T = I O 3 A/m
Hinze: T = (2*k*H*t*x)/(x2+d2)
T = magnetic anomaly [oersteds] k = magnetic susceptibility contrast [emu/cmq t = thickness of slab [cm]
d = depth to center of slab [cm] x = horizontal distance from observation point to axis [cm]
GM-SYS Hinze Maximum Amplitude = 52.2931 7 51.681 96 nT Minimum Amplitude = -15.67101 -14.98471 nT
B-4
- RESULTS - Magnetic Response over a Semi-Infinite Horizontal SlabNertical Fault
dation, x [m] x [cml T [nr], GM-SYS
shifted to zero at -50 m T [oerstedl, Hinze T [nr], Hinze,
shifted to zero at -50 m T [ n l , Hinze -50 -49 -48 -47 -46 -45 -44 -43 -42 -41 -40 -39 -38 -37 -36 -35 -34 -33 -32 -31 -30 -29 -28 -27 -26 -25 -24 -23 -22 -21 -20 -19 -18 -17 -16 -15 -14 -13 -12 -1 1 -10 -9 -8 -7 -6 -5 -4 -3 -2 -1 0