<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">OJG</journal-id><journal-title-group><journal-title>Open Journal of Geology</journal-title></journal-title-group><issn pub-type="epub">2161-7570</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ojg.2018.83015</article-id><article-id pub-id-type="publisher-id">OJG-83286</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Earth&amp;Environmental Sciences</subject></subj-group></article-categories><title-group><article-title>
 
 
  Anomalous Geomagnetic Field Modeling for Crustal Structure of a Graben in West Siberia
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Nikolay</surname><given-names>Pavlovich Kostrov</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Kirill</surname><given-names>Svjatoslavich Ivanov</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Institute of Geology and Geochemistry, Urals Branch of Russian Academy of Sciences, Yekaterinburg, Russia</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>kostrov65@yandex.ru(NPK)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>12</day><month>03</month><year>2018</year></pub-date><volume>08</volume><issue>03</issue><fpage>232</fpage><lpage>246</lpage><history><date date-type="received"><day>21,</day>	<month>February</month>	<year>2018</year></date><date date-type="rev-recd"><day>23,</day>	<month>March</month>	<year>2018</year>	</date><date date-type="accepted"><day>26,</day>	<month>March</month>	<year>2018</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  An intensive magnetic anomaly within the limits of West Siberia Danilov graben-rift indicates magnetic rocks while numerous wells encountered only weakly magnetized Triassic basalts and liparites in the basement covered by thick loose Jurassic and younger sediments. The wells penetrated only the first tens meters of the basement and could not tell us about the liparites structure at depth where supposedly they may form a big single body and magnetic rocks may be situated deeper. Geological ideas on a graben-rift structure may be proved (or rejected) by a computer modeling of its magnetic properties. For the anomalous geomagnetic field interpretation, a method of volume integral equations taking into account demagnetization effect was employed. To fit a model a trial-and-error procedure was utilized. The results show that 1) at the depth some rocks are magnetized in opposite direction to the present field; 2) highly magnetized rocks (magnetic susceptibility 0.06 - 0.1SI) coming up continuously from the bottom of the model and situated under the graben; 3) the studied structure is not a graben but the rift because the continental light crust is absent.
 
</p></abstract><kwd-group><kwd>Graben-Rift</kwd><kwd> Geomagnetic Anomaly</kwd><kwd> Modeling</kwd><kwd> Integral Equation</kwd><kwd> Interpretation</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>The pioneer researchers of West Siberia empirically found that the main oil deposits of the region are related to the system of Triassic rifts ( [<xref ref-type="bibr" rid="scirp.83286-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.83286-ref2">2</xref>] and other) and located predominately in such places in the basement, where Precambrian complexes had been developed. The conducted analysis has shown that two groups of reasons cause this: 1) Rifts destroy continental crust and are regional penetrable zones draining the upper mantle and it is in these zones that deep fluids easy penetrate in the sedimentary cover; 2) Deformations and structures of the sedimentary cover are defined by the basement. Thickness of the basement of West Siberia plate in average 14 times exceeds thickness of the cover. Taking into account average strength of basement rocks (153 &#177; 10 MPa) and other factors one may conclude that its strength of two orders of magnitude is greater than the strength of the sedimentary cover. Sialic blocks (considered earlier as Precambrian) are composed of middle-low crust rocks, predominantly granitoids and metamorphic schists, with density deficiency 0.07 - 0.1 g/cm<sup>3</sup> (relatively to the mean earth crust density 2.67 g/cm<sup>3</sup>). These blocks come to the level of the upper crust in time of Triassic rifting and continue isostatically “floating” after Triassic time. As a result, antiforms in sedimentary cover were created above them. Later the antiforms were filled with oil.</p><p>A modeling area is the central part of Danilov graben-rift that situated in Khanty-Mansi autonomous district (KhMAD) of the central part of the West Siberia oil and gas basin. KhMAD is the main oil region of the country―about half of Russian oil is produced here. The modeling area may be considered as a standard for the entire West Siberia taking into consideration that Danilov and North-Danilov oilfields are here and there are representative geophysical and wellsite data. Besides, vast majority of wells encountered only the uppermost parts of the basement not only in the modeling area but also throughout West Siberia that make it difficult to study graben-rift structures and as well as to detect and classify faults. Therefore, it is of importance to clarify the structure of the standard area utilizing potential field information and geophysical approach in computer modeling. It should be said that such modeling was done for the first time in the study area as scientists worldwide constantly address the potential field interpretation for graben and rift studies [<xref ref-type="bibr" rid="scirp.83286-ref3">3</xref>] - [<xref ref-type="bibr" rid="scirp.83286-ref13">13</xref>] .</p><p>Magnetic anomalies allow us to trace earth crust faults, borders of blocks and structures [<xref ref-type="bibr" rid="scirp.83286-ref14">14</xref>] [<xref ref-type="bibr" rid="scirp.83286-ref15">15</xref>] . In particular, the relatively magnetic serpentinites are able to make chain anomalies along graben borders. On the other hand, the fact that an intense magnetic anomaly is observed in the limits of Danilov graben-rift of West Siberia but numerous wells encountered only weakly magnetized Triassic graben basalts and liparites, leads us to conclusion that magnetic rocks are situated at great depth. In addition, we suppose that the liparites joint together in a single body. In fact, the wells penetrated only the first tens meters of the basement that covered everywhere by loose sediments of Jurassic and younger age with thickness of 1.5 - 2.5 km.</p><p>Besides, geomagnetic data themselves do not give a guaranty of single-valued interpretation because of potential character of the field. Therefore, fitting of magnetic characteristic having in mind reasonable geological constraints, allows not only testing geological hypotheses but also refining existing structural maps. Thus, validity of geological conceptions on graben-rift structure can be proved or denied with the help of computer modeling.</p><p>The purpose of the present study is to build magnetic model of the central part of Triassic Danilov graben situated within Khanty-Mansi autonomous district (KhMAD) of central part of the West Siberia oil and gas basin.</p></sec><sec id="s2"><title>2. Geological Background</title><p>West Siberia basement has been studied for more than fifty years [<xref ref-type="bibr" rid="scirp.83286-ref16">16</xref>] [<xref ref-type="bibr" rid="scirp.83286-ref17">17</xref>] [<xref ref-type="bibr" rid="scirp.83286-ref18">18</xref>] - [<xref ref-type="bibr" rid="scirp.83286-ref24">24</xref>] . The efforts have resulted in detecting many grabens originated in rift formation process that supposedly are responsible for forming oil and gas deposits [<xref ref-type="bibr" rid="scirp.83286-ref25">25</xref>] [<xref ref-type="bibr" rid="scirp.83286-ref26">26</xref>] . It is generally recognized that all West Siberia was subjected to rifting that leaves traces in the form of branching graben system covered by Mesozoic and Cenozoic sediments [<xref ref-type="bibr" rid="scirp.83286-ref27">27</xref>] . In the central part of West Siberia, in the limits of junction zone of the eastern slope of the Urals and Tyumen region of KhMAD, two Triassic graben-rifts were singled out [<xref ref-type="bibr" rid="scirp.83286-ref28">28</xref>] [<xref ref-type="bibr" rid="scirp.83286-ref29">29</xref>] : northern Sosva (at the east slope of the Urals) and Danilov (in the western part of KhMAD). Danilov graben-rift geology (<xref ref-type="fig" rid="fig1">Figure 1</xref>) can be considered as a standard for such structures.</p></sec><sec id="s3"><title>3. Method</title><p>3D magnetic model was built with a trial-and-error procedure when repeatedly solving a direct problem. The method was realized in software program complex “3D Interpretation of potential fields” [<xref ref-type="bibr" rid="scirp.83286-ref30">30</xref>] [<xref ref-type="bibr" rid="scirp.83286-ref31">31</xref>] . Our fitting methodology takes into account the a priori geological information [<xref ref-type="bibr" rid="scirp.83286-ref32">32</xref>] [<xref ref-type="bibr" rid="scirp.83286-ref33">33</xref>] [<xref ref-type="bibr" rid="scirp.83286-ref34">34</xref>] as well as ideas on graben-rift formation processes. Thus, in the given method an interpreter specifies initial conditions that include geometry of the basement, magnetizing field, magnetic susceptibility and remanent magnetization distributions and obtains at the end total force anomalous magnetic field (AMF), which is compared with observed AMF. Then model is changed and procedure is repeated unless satisfactory fitting is obtained. In doing so in the fitting, many times the following equation for magnetic field H ( r ) is solved [<xref ref-type="bibr" rid="scirp.83286-ref35">35</xref>] :</p><p>H ( r ) = H 0 ( r ) − 1 4 π ∇ ∫ V 2 ( κ ( r ′ ) H ( r ′ ) + R ( r ′ ) ) r − r ′ | r − r ′ | 3 d V 2 , (1)</p><p>where | r − r ′ | = ( x − x ′ ) 2 + ( y − y ′ ) 2 + ( z − z ′ ) 2 , r ′ ―the radius-vector of a source point P ′ ( x ′ , y ′ , z ′ ) , r―the radius vector of an observation point P ( x , y , z ) , in which we calculate the field, κ ( r ′ ) ―magnetic susceptibility. An observation point may be inside or outside of a source body. The equation can be obtained from anomalous field definition [<xref ref-type="bibr" rid="scirp.83286-ref36">36</xref>] or from Maxwell equations by standard procedure involving Green’s theorem application [<xref ref-type="bibr" rid="scirp.83286-ref37">37</xref>] [<xref ref-type="bibr" rid="scirp.83286-ref38">38</xref>] [<xref ref-type="bibr" rid="scirp.83286-ref39">39</xref>] [<xref ref-type="bibr" rid="scirp.83286-ref40">40</xref>] . The equation is solved by a numerical method based on [<xref ref-type="bibr" rid="scirp.83286-ref41">41</xref>] [<xref ref-type="bibr" rid="scirp.83286-ref42">42</xref>] .</p><p>The magnetic model was based on our previously constructed density model [<xref ref-type="bibr" rid="scirp.83286-ref43">43</xref>] . The density model of the Danilov graben has shown that rocks with a significantly excess density of 0.3 g/cm<sup>3</sup> and higher, which probably represent the Triassic mantle diapir, should be located under the axial part of the graben.</p><p>Rhyolites, forming small irregular bodies on the geological map of the basement, are combined into the common lopolith-like subvolcanic body at depths. The graben is flattened toward the north, which explains the previously unclear presence of Carboniferous shales among the Triassic rocks in the northern part of the structure. In the course of modeling we achieved good correspondence of magnetic and density geometry. The correctness of the gravitation model was checked by the magnetic model of the graben-rift. Actually, we have a single graben-rift model, in which we fitted different parameters (density, magnetic susceptibility, natural remanent magnetization) and geometry of the bodies.</p></sec><sec id="s4"><title>4. Technical Characteristic of the Model</title><p>In modeling, we used a rectangular system of coordinates with the vertical axis pointed upwards, the abscissa directed eastward and the ordinate directed to the north. Gauss-Kr&#252;ger projection for the 11 zone was used to transfer geological formation contours to the “Complex…”. The final model contains 222,301 elementary cells located in 22 horizontal layers. In horizontal plane, an elementary cell has a dimension of 1.5 &#215; 1.5 km. A layer contains 10,100 elementary cells. Each elementary cell is assigned magnetic susceptibility and remanent magnetization. Basement relief begins with a depth of 1.25 km and continues to a depth of 2.5 km. In this range, a layer thickness of 0.25 km was used. Deeper, from a depth from 2.5 to 12.5 km, we define thickness of a layer equal to 1 km. The grids of the basement and observed AMF coincide and have a dimension of 74 &#215; 76 nodes. The source AMF is presented for epoch of 1965 (<xref ref-type="fig" rid="fig2">Figure 2</xref>, the upper part), which is determined calculation of the magnetizing field in IGRF model [<xref ref-type="bibr" rid="scirp.83286-ref44">44</xref>] . The magnetizing field was calculated at each node of the observed AMF grid and considered constant up to half distance to a next node. The data on magnetic susceptibility were taken typical for the rocks of the area.</p></sec><sec id="s5"><title>5. Accuracy of Fitting</title><p>An accuracy is characterized by residual field (<xref ref-type="fig" rid="fig2">Figure 2</xref>, bottom part) obtained as difference between observed and model AMF. As one can see from the bottom part of <xref ref-type="fig" rid="fig2">Figure 2</xref>, the residual field isolines have mosaic structure with fine details in a range of &#177;30 nT that is about 5% of the range of the observed AMF (changing from −242 to 1053 nT) excluding boundaries of the model area.</p></sec><sec id="s6"><title>6. Results and Discussion</title><p>Our initial supposition that AMF in the area is due to only induced magnetization failed because we cannot model satisfactorily some local negative anomalies above almost nonmagnetic bodies. Another assumption that AMP above the study area contains significant regional component does not confirm. Any filtering of the observed AMF with a purpose to separate regional and local components leads to wide negative anomalies in the local field (taking more than 90% of the area) that requires remanent magnetization opposite to present</p><p>geomagnetic field for all bodies. In the long run, we decide to fit source observed AMF and compensate negative anomalies after taking into account induced magnetization.</p><p>As one can see on <xref ref-type="fig" rid="fig3">Figure 3</xref>, blocks with opposite to the present geomagnetic field magnetization were inserted at depths from 8.5 to 12.5 km. Other model cells have all components of remanent magnetization equal to zero. We assign values only to z-component of remanent magnetization. Reverse magnetization rocks can be explained by tectonic movement, metasomatosis, or inversion of the main geomagnetic field when rock to form. An analysis of the initial fitting reveals big variety of magnetic susceptibility values that can be grouped in several classes in the ranges of 0 - 0.00072, 0.005 - 0.007, 0.02 - 0.04, and 0.07 - 0.1 SI units. In some cases, it was useful to get average value for a class and alter body geometry that rather improve quality of the fit. The result is shown on <xref ref-type="fig" rid="fig4">Figure 4</xref>, where magnetic susceptibility structure. On the <xref ref-type="fig" rid="fig4">Figure 4</xref> magnetic susceptibility in SI units are shown in every horizontal layer of the model where magnetic susceptibility is not equal to zero. The magnetic susceptibility scale is shown to the right. Subscriptions under each graph is a layer range in vertical direction (z-coordinates) beginning from a depth of 1.25 km to the bottom of the model (12.5 km). In the central part of the area, one can trace liparite magnetic susceptibility change with a depth (shown in light blue colors) as well as magnetic</p><p>susceptibility of graben-rift basalts (shown in orange colors). At a depth range 1.75 - 2.5 km the liparites have worm-like structure that encountered by wells, below a depth of 2.5 km the liparites form lopolithlike body. Highly magnetic (having big magnetic susceptibility) areas (shown in red colors) are at the bottom of the model.</p><p>Let us consider vertical cross-sections of the magnetic susceptibility along profiles (AB, CD, MN, RS and OP) presented on <xref ref-type="fig" rid="fig1">Figure 1</xref>. The cross-sections are shown on <xref ref-type="fig" rid="fig5">Figure 5</xref>, vertical and horizontal scales coincide.</p><p>The profile MN from the southwest to the northeast (center of <xref ref-type="fig" rid="fig5">Figure 5</xref>) begins in a serpentinite zone (<xref ref-type="fig" rid="fig1">Figure 1</xref>) bounding the graben-rift from the west. Rocks model the periphery of serpentinites have magnetic susceptibility 0.02 - 0.04 SI units and reach a depth of 5.5 km. There are practically nonmagnetic blocks underneath. Further, the profile intersect effusive rocks with tuffs, a</p><p>negative anomaly under which partly taken into account by bodies with reversed magnetization (<xref ref-type="fig" rid="fig3">Figure 3</xref>, zone 5). Further on the cross-section clear seen basalts of the graben-rift with magnetic susceptibility 0.025 SI units and less magnetic liparites (0.007 SI units). At a mark of 40 km, the profile intersect intense magnetic anomaly (<xref ref-type="fig" rid="fig2">Figure 2</xref>, the upper part) above weakly magnetized graben-rift basalts and liparites (<xref ref-type="fig" rid="fig1">Figure 1</xref>). The given anomaly can be explained by very magnetic rocks under the graben-rift. On the given cross-section this rocks come from a depth of 12.5 km to a depth of 5.5 km. Further along the cross-section one can observe graben-rift bottom rising from the southwest to the northeast. Gabbro-dolerites massif (<xref ref-type="fig" rid="fig1">Figure 1</xref>) are marked by intense magnetic anomaly (about 600 nT in the maximum) (<xref ref-type="fig" rid="fig2">Figure 2</xref>, bottom part), which can be seen at 100 km of the profile. An irregular body (<xref ref-type="fig" rid="fig4">Figure 4</xref>), which cross-section along the profile reaches a depth of 11.5 km, models the anomaly. Further, there is a long low intensity anomaly (170 nT at the maximum). However, according to the geological map (<xref ref-type="fig" rid="fig1">Figure 1</xref>) here graben basalts are situated and to model the anomaly we entered magnetic body at a depth of 12.5 - 6.5 km (<xref ref-type="fig" rid="fig4">Figure 4</xref> and <xref ref-type="fig" rid="fig5">Figure 5</xref>, profile MN). Further, effusive rocks with tuffs have near zero magnetic susceptibility.</p><p>The profile RS (<xref ref-type="fig" rid="fig5">Figure 5</xref>, cross-section RS), coming from the northwest to the southeast (<xref ref-type="fig" rid="fig1">Figure 1</xref>, <xref ref-type="fig" rid="fig2">Figure 2</xref> upper part), gives us an insight into the transversal graben-rift structure. The graben-rift at this cross-section begins 100 km from the profile beginning (magnetic susceptibility 0.025 SI units). To the west of this point, magnetic bodies model Triassic liparites (<xref ref-type="fig" rid="fig1">Figure 1</xref>) with susceptibility 0.007 SI units. The cross-section comes through a circumference of serpentinites limiting graben-rift from the west. To the northwest Carboniferous terrigenous-schist formation has magnetic susceptibility near zero (0.00005 SI units) and give nothing to the anomaly. Surrounding magnetic rocks compensate only a half of intensity of the negative magnetic anomaly (about −130 nT in the minimum) and other half is given by deep seated (from 12.5 to 9.5 km depth) rocks (<xref ref-type="fig" rid="fig3">Figure 3</xref>, zone 2) with remanent magnetization opposite to the present main geomagnetic field. Above granitoids (from 20 to 40 km along profile), an intense magnetic anomaly (492 nT in the maximum) is modelled by rock underneath.</p><p>Coming back to magnetic structure of the graben-rift (110 km along the profile and further) it should be noted its lopolith-like form manifested in magnetic properties and the bottom at a depth of about 5.5 km.</p><p>Profiles AB and CD (<xref ref-type="fig" rid="fig5">Figure 5</xref>) come from the west to the east. At 40 km from the beginning, the profile AB intersects Late Devonian―Early Carboniferous effusives represented by practically nonmagnetic rocks. The profile CD, situated 30 km to the south, intersects magnetic serpentinites, magnetic susceptibility of which changes from 0.06 to 0.1 SI units in the center and from 0.02 to 0.04 SI units at the periphery. They give an anomaly by order of 1000 nT in the maximum and bound the graben-rift at the west.</p><p>Nonmagnetic effusive rocks with tuffs adjoin to them. A negative anomaly (−260 nT) above them are mainly due to magnetic serpentinites at the west, but requires reversely magnetized rock at a deps from 12.5 to 8.5 km (<xref ref-type="fig" rid="fig3">Figure 3</xref>, zone 5) for full fitting.</p><p>The graben-rift is clearly seen in the plane (<xref ref-type="fig" rid="fig4">Figure 4</xref>) and on these cross-sections as rocks with magnetic susceptibility ranging from 0.02 to 0.04 SI units and deepening in some places to 5 km depth. At the end of AB cross-section, further from 120 km and to the end, lowering of AMF (to −170 nT) is observed. This part of the profiles comes through nonmagnetic Carboniferous terrigenous-schist formation. Negative AMF above it cannot be explained by induction of strongly magnetized rocks and require a deep-seated reversely magnetized body (<xref ref-type="fig" rid="fig3">Figure 3</xref>, zone 5). As can be seen from AB cross-section, under relatively weakly magnetized rocks of the graben-rift (basalts and liparites) there is a strongly magnetized block of rocks with magnetic susceptibility varying from 0.06 to 0.1 SI units at a depth from 8.5 to 12.5 km. Magnetic susceptibility of graben-rift rocks (changing from 0.02 to 0.04 SI units) gives a wide positive AMF.</p><p>Profile OP demonstrates transversal cross-section of magnetic structure of the rock at the northeast (<xref ref-type="fig" rid="fig1">Figure 1</xref>), where the bottom of the graben-rift rises to the erosional truncation of the crystal basement. The profile begins from peripheral fragment of serpentinites, the most part of which are situated beyond the modeling area. This is reflected in magnetic susceptibility structure first 5 km along profile. Further, the profile intersects Late Devonian-Early Carboniferous effusives, basalts of the graben-rift (the thin stripe in the upper part of the cross-section with magnetic susceptibility 0.02 - 0.04 SI units) and gabbro-dolerites (0.06 - 0.1 SI units). Negative AMF around 70 km mark of the profile has engaged our attention because it is formed by strongly magnetized bodies (0.06 - 0.1 SI units) situated beyond its outline at a depth of 7.5 km at the northwest and rising up to a depth of 4.5 km at the southeast.</p><p>About magnetic susceptibility structure of granitoids and bordering sialic gneiss and schists laying to the west from the graben-rift on the basis of the RS and AB cross-sections (<xref ref-type="fig" rid="fig5">Figure 5</xref>) and layer by layer structure (<xref ref-type="fig" rid="fig4">Figure 4</xref>) it is possible to say that their magnetic susceptibility (0.005 - 0.007 SI units first 40 km along RS and about zero the first 20 km along AB) is not big enough to generate the observed anomaly. To model the anomaly we enter highly magnetized rocks (0.06 - 0.1 SI units) underneath.</p></sec><sec id="s7"><title>7. Conclusions</title><p>At modeling, we supposed that observed AMF does not have significant regional component above the modeling area.</p><p>The presented magnetic modeling confirms the results of gravimetric modeling that, in the first place, the studied structure is not a graben but a rift because continental light crust is absent in the structure (it pulled apart).</p><p>The modeling of magnetic susceptibility and natural remanent magnetization of the rocks in the limits of Danilov graben-rift of West Siberia reveals the presence of continuous magnetic rocks coming from a depth of 12.5 km with magnetic susceptibility of 0.06 - 0.1 SI units, which can be interpreted as manifestation of a mantle diapir.</p><p>The modeling shows that some local zones of nonmagnetic (according to the geological map) rocks require entering at a depth of 7.5 - 12.5 km reversely magnetized bodies that can be explained if their formation took place in a period of geomagnetic field reversals or recursion. However, the hypotheses about obtaining reverse remanent magnetization in the course of metasomatosis and tectonic movement cannot be excluded.</p><p>The graben-rift structure traced by magnetic susceptibility points to that 1) the bottom of the graben-rift rises from the southwest to the northeast; 2) from the northwest to the southeast the graben-rift has lopolith-like structure with eastern fall of western borders marked by serpentinites and western fall of eastern borders.</p><p>In total, the modeling results allow confirming and refining geological assumptions on structure of Danilov graben-rift of West Siberia.</p><p>Potential field modeling limitation is well known ambiguity: infinitely many models can fit the observed field. The limitation may be overcome by taking into account geological considerations, drilling data and knowledge of genesis of the area. Further development of graben rock property modeling may lay in direction of creating a complex geophysical model incorporating density, magnetic and heat flow parameters.</p></sec><sec id="s8"><title>Funding</title><p>The work was supported by Russian Science Foundation (Project 16-17-10201).</p></sec><sec id="s9"><title>Cite this paper</title><p>Kostrov, N.P. and Ivanov, K.S. (2018) Anomalous Geomagnetic Field Modeling for Crustal Structure of a Graben in West Siberia. Open Journal of Geology, 8, 232-246. https://doi.org/10.4236/ojg.2018.83015</p></sec></body><back><ref-list><title>References</title><ref id="scirp.83286-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Ivanov, K.S., Kucherov, V.G. and Fedorov, Yu.N. (2008) To the Question about Deep Origin of Oil. Proceedings ZapSibNIIGG Conference Condition, Tendencies and Problems Oil and Gas Potential of West Siberia Development, Tyumen, 160-173. (In Russian)</mixed-citation></ref><ref id="scirp.83286-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Nezhdanov, A.A., Ponomarev, V.A., Turenkov, N.A. and Gorbunov, S.A. (2000) Geology and Oil and Gas Bearing of Achmovskaya Thickness of West Siberia. Publisher of Academy of Rock Sciences, Moscow. (In Russian)</mixed-citation></ref><ref id="scirp.83286-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Miller, S.P. and Hey, R.N. (1986) Three-Dimensional Magnetic Modeling of a Propagating Rift, Galapagos 95&amp;#176;30'W. Journal of Geophysical Research, 91, 3395-3406.  
https://doi.org/10.1029/JB091iB03p03395</mixed-citation></ref><ref id="scirp.83286-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Mariano, J. and Hinze, W. (1994) Gravity and Magnetic Models of the Midcontinent Rift in Eastern Lake Superior. Canadian Journal of Earth Sciences, 31, 661-674.  
https://doi.org/10.1139/e94-059</mixed-citation></ref><ref id="scirp.83286-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Mugisha, E., Ebinger, C.J, Strecker, M. and Pope, D. (1997) Two-Stage Rifting in the Kenya Rift: Implications for Half-Graben Models. Tectonophysics, 278, 63-81.  
https://doi.org/10.1016/S0040-1951(97)00095-4</mixed-citation></ref><ref id="scirp.83286-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">MacLean, B.C. and Miles, W. (2002) Potential-Field Modeling of a Proterozoic Half-Graben near Blackwater Lake, Northwest Territories, Canada, and Its Implications for the Fort Simpson Magnetic Anomaly. Canadian Journal of Earth Sciences, 39, 169-187. https://doi.org/10.1139/e01-061</mixed-citation></ref><ref id="scirp.83286-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Bhattacharji, S., Sharma, R. and Chatterjee, N. (2004) Two- and Three-Dimensional Gravity Modeling along Western Continental Margin and Intraplate Narmada-Tapti Rifts: Its Relevance to Deccan Flood Basalt Volcanism. Journal of Earth System Science, 113, 771-784. https://doi.org/10.1007/BF02704036</mixed-citation></ref><ref id="scirp.83286-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Ebbing, J., Skilbrei, J.R. and Olesen, O. (2007) Insights into the Magmatic Architecture of the Oslo Graben by Petrophysically Constrained Analysis of the Gravity and Magnetic Field. Journal of Geophysical Research, 112, Article ID: B04404, 17 p.</mixed-citation></ref><ref id="scirp.83286-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Bridges, D.L. (2011) A Geological and Geophysical Study of the Tendaho Graben in the Afar Depression, Ethiopia: Insights into Transitional Continental Rifting. Doctoral Dissertations, Missouri University of Science and Technology, Rolla, Missouri.  
http://scholarsmine.mst.edu/doctoral_dissertations/2020</mixed-citation></ref><ref id="scirp.83286-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Alatorre-Zamora, M.A., Campos-Enríquez, J.O., Rosas-Elguera, J.G., Pena-García, L., Maciel-Flores, R. and Fregoso-Becerra, E. (2015) Chapala Half-Graben Structure Inferred. A Magnetometric Study. Geofísica Internacional, 54, 323-342.  
http://www.redalyc.org/articulo.oa?id=56842730003</mixed-citation></ref><ref id="scirp.83286-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Lopes de Castro, D. (2011) Gravity and Magnetic Joint Modeling of the Potiguar Rift Basin (NE Brazil): Basement Control during Neocomian Extension and Deformation. Journal of South American Earth Sciences, 31, 186-198.  
https://doi.org/10.1016/j.jsames.2011.01.005</mixed-citation></ref><ref id="scirp.83286-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Langenheim, V.E., Jachens, R.C., Muffler, L.J.P. and Clynne, M.A. (2016) Implications for the Structure of the Hat Creek Fault and Transfer of Right-Lateral Shear from the Walker Lane North of Lassen Peak, Northern California, from Gravity and Magnetic Data. Geosphere, 12, 790-808. https://doi.org/10.1130/GES01253.1</mixed-citation></ref><ref id="scirp.83286-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Obi, D.A., Ilozobhie, A.J., Lebo, S.E. and Zoogbara, E. (2017) Modeling Magnetic Basement in Relationship to Hydrocarbon Habitats in the Central Niger Delta, Nigeria. Journal of Geography, Environment and Earth Science International, 10, 1-13. https://doi.org/10.9734/JGEESI/2017/31888</mixed-citation></ref><ref id="scirp.83286-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Ivanov, K.S., Kormiltsev, V.V., Fedorov, Yu.N., et al. (2003) The Main Characters of Pre-Jurassic Basement Structure in Shaim Oil and Gas Bearing Region. The Way of Realization of Oil and Gas Potential of KhMAD. Proceedings of Viscientific-Pracical Conference, V.1. Khany-Mansiysk, 102-113. (In Russian)</mixed-citation></ref><ref id="scirp.83286-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Rybakov, M.B., Goldshtein, V.I., Holl, D., Ben-Avraam, Z. and Lazar, M. (2011) New Ideas on Sources of Magnetic Anomalies in the East Mediterranean: A Key to Ophyolites and Gabbro Intrusion Localiation. Geologiyaigeofizika, 4, 487-511. (In Russian)</mixed-citation></ref><ref id="scirp.83286-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Rostovtsev, N.N. (1964) Geology of the USSR. West-Siberian lowland. Geological Description, Nedra, Moscow, Vol. 44, 1-550. (In Russian)</mixed-citation></ref><ref id="scirp.83286-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Ivanov, K.S., Puchkov, V.N., Fyodorov, Y.N., Erokhin, Y.V. and Pogromskaya, O.E. (2013) Tectonics of the Urals and Adjacent Part of the West-Siberian Platform Basement: Main Features of Geology and Development. Journal of Asian Earth Sciences, 72, 12-24. https://doi.org/10.1016/j.jseaes.2013.02.029</mixed-citation></ref><ref id="scirp.83286-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Grachev, A.F. (1977) Rift Zones of Earth. Nedra, Leningrad. (In Russian)</mixed-citation></ref><ref id="scirp.83286-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Kucheruk, E.V. and Ushakov, S.A. (1985) Plate Tectonics, Oil, and Gas Bearing (Geophysical Analysis). VINITI, Moscow. (In Russian)</mixed-citation></ref><ref id="scirp.83286-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Milankovskiy, E.E. (1976) Rift Zones of Continents. Nedra, Moscow. (In Russian)</mixed-citation></ref><ref id="scirp.83286-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Elkin, E.A., Krasnnov, V.I., Bakharev, N.K., Belova, E.V., Dubatolov, V.N., Izokh, N.G., Klets, A.G., Kontorovich, A.E., Peregoedov, L.G., Sennikov, V.V., Timokhina, I.G. and Khromykh, V.G. (2001) Stratigraphy of Oil and Gas Bearing Basins of Siberia. Paleozoic of West Siberia. Publisher SO RAS, Subsidiary Geo, Novosibirsk. (In Russian)</mixed-citation></ref><ref id="scirp.83286-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">Bogush, O.I., Bochkarev, V.S. and Yuferev, O.V. (1975) Paleozoic of the South of West Siberia Plane. Novosibirsk, Nauka. (In Russian)</mixed-citation></ref><ref id="scirp.83286-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Klets, A.G., Kontorovich, V.A., Ivanov, K.S., Kazanenkov, V.A., Saraev, S.V., Simonov, V.A. and Fomin, A.N. (2007) Geodynamical Model of Pre-Jurassic Basement—The Basis of Oil and Gas Bearing Zoning of Upper PreCambrian—Low Triassic Structural Stage of West Siberia Oil and Gas Bearing Province. Proceeding of Conference “The Way of Realization of Oil and Gas and Ore Potential of KhMAD”, Khanty-Mansiysk, 13-17 November 2007, Vol. 1, 79-90.</mixed-citation></ref><ref id="scirp.83286-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">Fedorov, Yu.N., Ivanov, K.S., Koroteev, V.A. and Kormiltsev, V.V. (2004) Tectonic Structure of Transuralian Part of West Siberia Megabasin. Proceeding of Conference “The Way of Realization of Oil and Gas and Ore Potential of KhMAD”, Khanty-Mansiysk, 2-5 December 2004, Vol. 1, 91-102. (In Russian)</mixed-citation></ref><ref id="scirp.83286-ref25"><label>25</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Kontorovich</surname><given-names> V.A. </given-names></name>,<etal>et al</etal>. (<year>2007</year>)<article-title>Seismological Criteria of Oil and Gas Bearing of the Contact Zone of Paleozoic and Mesozoic Sediments of West Siberia</article-title><source> Geologiyai Geofizika</source><volume> 48</volume>,<fpage> 538</fpage>-<lpage>547</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.83286-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">Ivanov, K.S. (2016) How Much Oil Does Russia Need to Produce? Uralian Geological Journal, 6, 3-14.  
http://docs.wixstatic.com/ugd/15eb72_f0c0f560110b4165b2e7221570aef26e.pdf</mixed-citation></ref><ref id="scirp.83286-ref27"><label>27</label><mixed-citation publication-type="other" xlink:type="simple">Rudkevich, M.Ya., Ozeranskaya, L.S., Chistyakova, N.F., et al. (1988) Oil and Gas Bearing Complexes of West Siberia Basin. Nedra, Moscow. (In Russian)</mixed-citation></ref><ref id="scirp.83286-ref28"><label>28</label><mixed-citation publication-type="other" xlink:type="simple">Ivanov, K.S., Koroteev, V.A., Pecherkin, M.F., Fedorov, Yu.N. and Erokhin, Yu.V. (2009) History of Geological Development and Structure of the Basement of Western Part of West Siberia Oil and Gas Megabasin. GeologiyaiGeofizika, 50, 484-501. (In Russian)</mixed-citation></ref><ref id="scirp.83286-ref29"><label>29</label><mixed-citation publication-type="other" xlink:type="simple">Fedorov, Yu.N., Ivanov, K.S., Zakharov, S.G., Kormiltsev, V.V., Erokhin, Yu.V., Pogromskaya, O.E., Knyazeva, I.V., Ronkin, Yu.L., Karetin, Yu.S., Surina, O.V., Purtova, S.I. and Glushko, N.K. (2003) Geological Structure and Stratigraphy of Triassic Sediments of the Northern Sosva Graben. Proceeding of Conference “The Way of Realization of Oil and Gas and Ore Potential of KhMAD”, Khanty-Mansiysk, IzdatNaukaServis Vol. 1, 114-123.</mixed-citation></ref><ref id="scirp.83286-ref30"><label>30</label><mixed-citation publication-type="other" xlink:type="simple">Kostrov, N.P., Kormiltsev, V.V. and Fedorov, Yu.N. (2005) System of 3D Interpretation of Gravitational and Magnetic Observations with the Aim of Geological Mapping of Pre-Jurassic West Siberia Complex. GornyeVedomosti, 1, 57-61. (In Russian)</mixed-citation></ref><ref id="scirp.83286-ref31"><label>31</label><mixed-citation publication-type="other" xlink:type="simple">Kostrov, N.P. and Ivanov, K.S. (2017) Complex for Interpretation of the Potential Fields for Geology. Vestnik of IGG Komi Science Center of UrB RAS, 10, 3-8. (In Russian) https://doi.org/10.19110/2221-1381-2017-10-3-8</mixed-citation></ref><ref id="scirp.83286-ref32"><label>32</label><mixed-citation publication-type="other" xlink:type="simple">Ananieva, E.M., Vinnichuk, N.N., Ivanov, K.S., et al. (2008) On Density of Rocks of East Urals and West Siberia Basement. UrB RAS, Yekaterinburg. (In Russian)</mixed-citation></ref><ref id="scirp.83286-ref33"><label>33</label><mixed-citation publication-type="other" xlink:type="simple">Gusev, B.V., Metallova, V.V. and Fainberg, F.S. (1967) Magnetism of Rocks of Traps Formation of Western Part of Siberia Platform. Nedra, Leningrad. (In Russian)</mixed-citation></ref><ref id="scirp.83286-ref34"><label>34</label><mixed-citation publication-type="other" xlink:type="simple">Chursin, A.V., Prutyan, A.M. and Fedorova, N.V. (2008) Digital Anomalous Magnetic Field Map of the Northern, Middle, Southern Urals and Adjacent Territoties of the East-European and West-Siberian Platforms. Litosfera, 6, 63-72. (In Russian)</mixed-citation></ref><ref id="scirp.83286-ref35"><label>35</label><mixed-citation publication-type="other" xlink:type="simple">Kormiltsev, V.V. and Ratushnyak, A.N. (2000) Modeling of Geophysical Fields with the Help of Volume Vector Integral Equations. UrB RAS, Yekaterinburg.</mixed-citation></ref><ref id="scirp.83286-ref36"><label>36</label><mixed-citation publication-type="other" xlink:type="simple">Eloranta, E.H. (1986) Potential Field of a Stationary Electric Current using Fredholm Integral Equations of the Second Kind. Geophysical Prospecting, 34, 856-872.  
https://doi.org/10.1111/j.1365-2478.1986.tb00498.x</mixed-citation></ref><ref id="scirp.83286-ref37"><label>37</label><mixed-citation publication-type="other" xlink:type="simple">Kormiltsev, V.V. and Kostrov, N.P. (1999) Study of Algorithms Based on a Volume Vector Integral Equation for Magnetic Field Strength. VINITI, Moscow, 83-B99. (In Russian)</mixed-citation></ref><ref id="scirp.83286-ref38"><label>38</label><mixed-citation publication-type="other" xlink:type="simple">Kormiltsev, V.V. and Ratushnyak, A.N. (1995) Vector Integral Equations for a Gradient of Potential Geophysical Fields. Rossiyskiygeofizicheskiyzhurnal, 5-6, 4-10. (In Russian)</mixed-citation></ref><ref id="scirp.83286-ref39"><label>39</label><mixed-citation publication-type="other" xlink:type="simple">Kostrov, N.P. (2007) Calculation of Magnetic Anomaliescaused by 2D Bodies of Arbitrary Shape with Consideration of Demagnetization. Geophysical Prospecting, 55, 91-115. https://doi.org/10.1111/j.1365-2478.2006.00579.x</mixed-citation></ref><ref id="scirp.83286-ref40"><label>40</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Sharma</surname><given-names> P.V. </given-names></name>,<etal>et al</etal>. (<year>1996</year>)<article-title>Rapid Computation of Magnetic Anomalies and Demagnetization Effects Caused by Bodies of Arbitrary Shape</article-title><source> Pure and Applied Geophysics</source><volume> 64</volume>,<fpage> 89</fpage>-<lpage>109</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.83286-ref41"><label>41</label><mixed-citation publication-type="other" xlink:type="simple">Atkinson, K.E. (1997) Thenumerical Solution of Integral Equation Soft the Second Kind. Cambridge University Press, Cambridge.  
https://doi.org/10.1017/CBO9780511626340</mixed-citation></ref><ref id="scirp.83286-ref42"><label>42</label><mixed-citation publication-type="other" xlink:type="simple">Blokh, Yu.I. (1987) Magnetic Field Calculation for Three-Dimensional Anisotropic Geological Objects with Corrections for Demagnetization. Izvestiya, Earth Physics, 23, 1003-1008.</mixed-citation></ref><ref id="scirp.83286-ref43"><label>43</label><mixed-citation publication-type="other" xlink:type="simple">Ivanov, K.S., Kostrov, N.P. and Manushko, E.A. (2013) Deep Structure of the Danilovskii Graben in Western Siberia: A Geological Interpretation of the Density Models. Doklady Earth Sciences, 451, 870-874.  
https://doi.org/10.1134/S1028334X13080102</mixed-citation></ref><ref id="scirp.83286-ref44"><label>44</label><mixed-citation publication-type="other" xlink:type="simple">Thébault, E., et al. (2015) International Geomagnetic Reference Field: The 12th Generation. Earth, Planets and Space, 67, 79. 
https://earth-planets-space.springeropen.com/articles/10.1186/s40623-015-0228-9  
https://doi.org/10.1186/s40623-015-0228-9</mixed-citation></ref></ref-list></back></article>