RU2683817C1 - Method for determining induced and residual magnetization of rocks according to magnetic exploration data - Google Patents
Method for determining induced and residual magnetization of rocks according to magnetic exploration data Download PDFInfo
- Publication number
- RU2683817C1 RU2683817C1 RU2018112370A RU2018112370A RU2683817C1 RU 2683817 C1 RU2683817 C1 RU 2683817C1 RU 2018112370 A RU2018112370 A RU 2018112370A RU 2018112370 A RU2018112370 A RU 2018112370A RU 2683817 C1 RU2683817 C1 RU 2683817C1
- Authority
- RU
- Russia
- Prior art keywords
- magnetization
- section
- values
- induced
- residual
- Prior art date
Links
- 230000005415 magnetization Effects 0.000 title claims abstract description 167
- 238000000034 method Methods 0.000 title claims abstract description 27
- 239000011435 rock Substances 0.000 title claims description 16
- 239000013598 vector Substances 0.000 claims abstract description 33
- 238000010276 construction Methods 0.000 claims abstract description 4
- 230000000694 effects Effects 0.000 abstract description 4
- 239000000470 constituent Substances 0.000 abstract 1
- 239000000126 substance Substances 0.000 abstract 1
- 241000566515 Nedra Species 0.000 description 2
- 230000002547 anomalous effect Effects 0.000 description 2
- 230000001939 inductive effect Effects 0.000 description 2
- 229910052500 inorganic mineral Inorganic materials 0.000 description 2
- 239000011707 mineral Substances 0.000 description 2
- 230000000877 morphologic effect Effects 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 230000005358 geomagnetic field Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000010187 selection method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/40—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for measuring magnetic field characteristics of the earth
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
Description
Изобретение относится к области магниторазведки и может быть использовано при построении разрезов по аномалиям магнитного поля на стадии решения региональных, поисковых и разведочных задач геологии.The invention relates to the field of magnetic exploration and can be used to construct sections from magnetic field anomalies at the stage of solving regional, prospecting and exploration problems of geology.
Одной из проблем при изучении магнитных свойств геологического разреза, является необходимость разделения индуцированной и остаточной намагниченностей в горных породах, так как интерпретация аномалий, связанных только с индуцированным намагничиванием, дает более уверенную информацию о морфологии аномальных объектов и их магнитной восприимчивости. Сильное искажающее влияние остаточной намагниченности не позволяет по данным измерений постоянного магнитного поля над объектами оценить их магнитные свойства.One of the problems in studying the magnetic properties of a geological section is the need to separate the induced and residual magnetizations in rocks, since the interpretation of anomalies related only to induced magnetization gives more reliable information about the morphology of anomalous objects and their magnetic susceptibility. The strong distorting effect of the remanent magnetization does not allow us to estimate their magnetic properties according to measurements of a constant magnetic field over objects.
Известны различные способы раздельного определения индуцированной и остаточной намагниченностей горных пород в лабораторных условиях на образцах с помощью астатических магнитометров или магнитных весов (см., например, «Физические свойства горных пород и полезных ископаемых» (петрофизика). Справочник геофизика. / Под ред. Н.Б. Дортман, - 2-е изд., перераб. и доп. - М.: Недра, 1984, стр. 102-143). Однако данные способы являются достаточно трудоемкой и затратной задачей и, в целом, дают лишь точечную информацию о поведении индуцированной и остаточной намагниченности в породе.There are various methods for the separate determination of the induced and residual magnetizations of rocks in laboratory conditions on samples using astatic magnetometers or magnetic scales (see, for example, “Physical properties of rocks and minerals” (petrophysics). Reference book geophysics. / Ed. N .B. Dortman, - 2nd ed., Revised and additional - M: Nedra, 1984, p. 102-143). However, these methods are a rather time-consuming and costly task and, in general, provide only point information on the behavior of the induced and residual magnetization in the rock.
Определение индуцированной намагниченности осуществляется также электроразведочными методами на основе искусственного подмагничивания с использованием постоянного или низкочастотного переменного токов. Сущность метода искусственного подмагничивания пород состоит в том, что с помощью расположенного на поверхности земли незаземленного контура постоянным или пульсирующим током малой частоты (около 0,7 Гц) возбуждается магнитное поле. Породы, подмагниченные этим полем, создают дополнительный (вторичный) эффект, величина которого зависит от магнитной восприимчивости, глубины до поверхности руд и их формы. Определению подлежит вторичное магнитное поле, создаваемое подмагниченными породами. Для этого непосредственно измеряют магнитное поле на исследуемом участке при отсутствии в контуре тока, при прохождении в контуре тока, а также нормальное поле контура над немагнитными породами (на контрольном участке). По вычисленному вторичному магнитному полю определяют магнитную восприимчивость пород или руд участка. Определение индуцированной намагниченности при применении низкочастотной индуктивной электроразведки, показано, например, в книге «Теоретические основы комплексной магниторазведки», Ю.И. Блох, © 2012, стр. 4-29. Однако данное техническое решение хотя и решает задачу определения векторов индуцированной и остаточной намагниченностей, но только на небольшие глубины, не более чем на несколько сотен метров, в результате чего не обеспечивается эффективность решения задач магниторазведки при глубинных исследованиях.The induced magnetization is also determined by electrical prospecting methods based on artificial magnetization using direct or low-frequency alternating currents. The essence of the method of artificial magnetization of rocks is that, using an ungrounded circuit located on the earth's surface, a constant or pulsating low-frequency current (about 0.7 Hz) is used to excite a magnetic field. The rocks magnetized by this field create an additional (secondary) effect, the magnitude of which depends on the magnetic susceptibility, depth to the surface of the ores and their shape. The secondary magnetic field created by magnetized rocks is subject to determination. To do this, directly measure the magnetic field in the studied area in the absence of current in the circuit, when passing in the current circuit, as well as the normal field of the circuit above non-magnetic rocks (in the control area). The calculated secondary magnetic field determines the magnetic susceptibility of the rocks or ores of the site. The definition of induced magnetization when using low-frequency inductive electrical exploration is shown, for example, in the book “Theoretical Foundations of Integrated Magnetic Exploration”, Yu.I. Bloch, © 2012, p. 4-29. However, this technical solution, although it solves the problem of determining the vectors of the induced and residual magnetizations, but only at shallow depths, no more than several hundred meters, as a result of which the efficiency of solving problems of magnetic prospecting in depth studies is not ensured.
Известны многочисленные способы определения эффективной намагниченности методами подбора на основе решения обратных задач магниторазведки, разработанные еще в прошлом веке (см., например, «Вычислительные математика и техника в разведочной геофизике». Справочник геофизика. Под ред. В.И.) Дмитриева. М., Недра, 1982. стр. 89-98., прототип).Numerous methods are known for determining the effective magnetization by selection methods based on solving inverse magnetic prospecting problems developed in the last century (see, for example, “Computational mathematics and technology in exploratory geophysics.” Reference Geophysics. Edited by V.I.) Dmitriev. M., Nedra, 1982. p. 89-98., Prototype).
Определенная тем или иным образом эффективная намагниченность изучаемого объекта не соответствует действительной, так как действительная намагниченность представляет сумму векторов индуцированной и остаточной намагниченностей. В целом, вектор индуцированной намагниченности совпадает с направлением вектора нормального магнитного поля Земли, а вектор остаточной намагниченности определяется системой иных механизмов намагничивания, таких как, фактор времени, температура, механические напряжения, химические преобразования, происходящие в присутствии магнитного поля. Направление суммарного вектора будет зависеть от соотношения величин намагниченностей (отношение Кеннингсберга). Если величина остаточной намагниченности будет значительно преобладать над величиной индуцированной намагниченности, то направление намагниченности геологических объектов будет отличаться от направления современного геомагнитного поля, что может в свою очередь привести к ошибочным результатам интерпретации аномальных магнитных полей.The effective magnetization determined in one way or another does not correspond to the real one, since the actual magnetization is the sum of the vectors of the induced and residual magnetizations. In general, the vector of induced magnetization coincides with the direction of the Earth’s normal magnetic field vector, and the residual magnetization vector is determined by a system of other magnetization mechanisms, such as time factor, temperature, mechanical stresses, chemical transformations that occur in the presence of a magnetic field. The direction of the total vector will depend on the ratio of the magnetization values (Kenningsberg ratio). If the magnitude of the residual magnetization will significantly prevail over the magnitude of the induced magnetization, the direction of the magnetization of geological objects will differ from the direction of the modern geomagnetic field, which may in turn lead to erroneous results of the interpretation of anomalous magnetic fields.
Задачей, на решение которой направлено изобретение, является повышение эффективности использования магниторазведочных данных за счет их более углубленной интерпретации.The problem to which the invention is directed, is to increase the efficiency of using magnetic data due to their more in-depth interpretation.
Технический результат изобретения выражается в получении магнитной модели изучаемого объекта в виде разделенных между собой моделей, модели магнитного поля, обусловленного индуцированной намагниченностью Ji, и модели магнитного поля, обусловленного остаточной намагниченностью Jr.The technical result of the invention is expressed in obtaining a magnetic model of the studied object in the form of separated models, a magnetic field model due to the induced magnetization Ji, and a magnetic field model due to the residual magnetization Jr.
Заявленный технический результат достигается за счет того, что способ определения индуцированной и остаточной намагниченностей горных пород по данным магниторазведки, согласно изобретению характеризуется тем, чтоThe claimed technical result is achieved due to the fact that the method for determining the induced and residual magnetizations of rocks according to magnetic exploration according to the invention is characterized in that
- по наблюденному магнитному полю создают цифровую модель разреза, в виде совокупности отдельных цифровых блоков, заполняющих изучаемое пространство и имеющих квазиоднородные магнитные свойства, при этом каждый блок представляют как неделимый элемент модели, который помимо пространственных параметров характеризуется намагниченностью и углом намагничивания,- using the observed magnetic field, create a digital model of the section, in the form of a set of separate digital blocks filling the studied space and having quasihomogeneous magnetic properties, each block being presented as an indivisible element of the model, which in addition to spatial parameters is characterized by magnetization and magnetization angle,
- для каждого указанного блока на основе решения обратной задачи определяют значения_эффективной намагниченности,- for each indicated block, based on the solution of the inverse problem, the effective magnetization values are determined,
- после чего уточняют геометрические параметры каждого составляющего разрез цифрового блока и соответствующие значения эффективной намагниченности (с углами намагничивания, совпадающими с направлением вектора нормального магнитного поля Земли),- then specify the geometric parameters of each digital section component of the section and the corresponding values of the effective magnetization (with magnetization angles coinciding with the direction of the Earth’s normal magnetic field vector),
- полученную цифровую модель разреза копируют, с получением второй идентичной модели разреза, при этом один разрез обозначают как разрез, характеризуемый только индуцированной намагниченностью, а второй - только остаточной намагниченностью,- the obtained digital model of the section is copied to obtain a second identical section model, with one section being designated as a section characterized only by induced magnetization, and the second only by residual magnetization,
- после чего осуществляют процедуры определения значений индуцированной и остаточной намагниченностей для каждого цифрового блока обоих разрезов путем решения обратных задач, последовательно изменяя направление вектора остаточной намагниченности,- then carry out procedures for determining the values of the induced and residual magnetizations for each digital block of both sections by solving inverse problems, sequentially changing the direction of the vector of residual magnetization,
- при каждом указанном изменении угла намагничивания в каждом блоке вычисляют значения индуцированной и остаточной намагниченностей, а также значения целевой функции, характеризующей степень сходимости модельного и наблюденного магнитных полей, полученные значения записывают в отдельный файл,- at each indicated change in the angle of magnetization in each block, the values of the induced and residual magnetizations are calculated, as well as the values of the objective function characterizing the degree of convergence of the model and observed magnetic fields, the obtained values are recorded in a separate file,
- для каждого блока формируют собственный файл, включающий изменения значений индуцированной, остаточной намагниченностей и целевой функции, в зависимости от изменения направления вектора остаточной намагниченности,- for each block form its own file, which includes changes in the values of the induced, residual magnetization and the objective function, depending on the change in direction of the residual magnetization vector,
- выделяют минимальные значения целевой функции для каждого входящего в разрез цифрового блока, которые соответствуют истинному направлению вектора остаточной намагниченности, а по известному направлению вектора - истинные значения индуцированной и остаточной намагниченностей,- select the minimum values of the objective function for each digital block included in the section, which correspond to the true direction of the residual magnetization vector, and the true values of the induced and residual magnetizations in the known direction of the vector,
- изменения указанных углов повторяют за несколько полных итерационных циклов до получения минимального значения целевой функции, при этом полученные значения намагниченностей и углов намагничивания определяют как векторные значения искомых индуцированной и остаточной намагниченностей,- changes in these angles are repeated for several complete iterative cycles until the minimum value of the objective function is obtained, while the obtained values of magnetization and magnetization angles are determined as vector values of the desired induced and residual magnetizations,
- в результате чего формируют цифровые модели изучаемого разреза в виде полей индуцированной и остаточной намагниченности,- as a result of which form digital models of the studied section in the form of fields of induced and residual magnetization,
- полученные цифровые модели преобразовывают в графические построения, которые визуализируют в виде карт изолиний, представляющих собой распределение эффективной, индуцированной и остаточной намагниченностей в плоскости изучаемого разреза,- the resulting digital models are converted into graphical constructions, which are visualized in the form of contour maps representing the distribution of the effective, induced and residual magnetizations in the plane of the studied section,
- полученные графические модели изучаемого разреза интерпретируют.- the resulting graphic models of the studied section are interpreted.
При этом способ согласно изобретению преимущественно реализуется с помощью компьютерного устройства.In this case, the method according to the invention is mainly implemented using a computer device.
Фиг. 1 - фиг. 7 иллюстрируют модельный пример реализации способа согласно изобретению, фиг. 8 иллюстрирует практический пример реализации способа.FIG. 1 - FIG. 7 illustrates a model example implementation of the method according to the invention, FIG. 8 illustrates a practical example of the implementation of the method.
Способ согласно изобретению осуществляется следующим образом.The method according to the invention is as follows.
По наблюденному магнитному полю, с использованием априорной информации, создают цифровую модель разреза (в случае решения двухмерной задачи) или серии разрезов (в случае трехмерной задачи). Разрез представляют в виде совокупности блоков (здесь под блоком понимается фигура, имеющая квазиоднородные магнитные свойства), заполняющих все изучаемое пространство. Первоначально геометрические параметры блоков определяют либо по характерным точкам, либо с использованием других геолого-геофизических данных, либо совместно и записывают их в цифровом виде. Затем, на основе решения прямых и обратных задач магниторазведки при первоначально заданной конфигурации блоков (первый тип задач) определяют значение эффективной намагниченности каждого блока. Далее, при известных значениях намагниченности, при помощи решения прямых и обратных задач (второй тип задач) уточняют геометрические параметры каждого блока. В окончательном варианте, при известных геометрических параметрах блоков, на основе решения обратных задач первого типа уточняют эффективную намагниченность каждого блока.Using the observed magnetic field, using a priori information, create a digital model of the section (in the case of solving a two-dimensional problem) or a series of sections (in the case of a three-dimensional problem). A section is represented as a set of blocks (here, a block is a figure having quasihomogeneous magnetic properties) that fill the entire space under study. Initially, the geometric parameters of the blocks are determined either by characteristic points, or using other geological and geophysical data, or jointly and recorded in digital form. Then, based on the solution of direct and inverse problems of magnetic prospecting with the initially given configuration of blocks (the first type of tasks), the value of the effective magnetization of each block is determined. Further, with known values of magnetization, by solving direct and inverse problems (the second type of problems), the geometric parameters of each block are refined. In the final version, with the known geometric parameters of the blocks, on the basis of solving inverse problems of the first type, the effective magnetization of each block is specified.
После того как определены геометрические параметры всех блоков, составляющих разрез и их эффективная намагниченность с углами намагничивания, совпадающими с направлением нормального магнитного поля Земли, создают копию полученной модели разреза. Так как разрезы идентичны, то и намагниченность блоков в разрезах будет попарно идентична.After the geometric parameters of all the blocks making up the section and their effective magnetization with the magnetization angles coinciding with the direction of the normal magnetic field of the Earth are determined, a copy of the obtained section model is determined. Since the sections are identical, the magnetization of blocks in the sections will be pairwise identical.
Произвольно назначается разрез, в котором будут определяться значения индуцированной намагниченности (исходный разрез). В другом разрезе (копия разреза) будут определяться значения остаточной намагниченности. Так как угол (iix) вектора намагничивания для индуцированной намагниченности (Ji) известен в исходном разрезе, то в другом, согласно изобретению, последовательно для каждого блока в копии разреза задают углы (irx) вектора (Jr) остаточной намагниченности и на основе решения прямых и обратных задач первого типа при каждом изменении угла вычисляют соответствующие значения индуцированной (Ji) и остаточной намагниченности (Jr) в обоих разрезах, в каждом блоке одновременно.A section is arbitrarily assigned, in which the values of the induced magnetization (initial section) will be determined. In another section (a copy of the section), the values of the residual magnetization will be determined. Since the angle (i ix ) of the magnetization vector for the induced magnetization (Ji) is known in the initial section, in the other, according to the invention, the angles (i rx ) of the residual magnetization vector (Jr) and, based on the solution direct and inverse problems of the first type, with each change in the angle, the corresponding values of the induced (Ji) and residual magnetization (Jr) are calculated in both sections, in each block at the same time.
Указанные процедуры осуществляют несколькими итерационными циклами. Первая итерация начинается с задания нулевого угла намагничивания для первого блока копии разреза и решением обратной задачи. В ходе решения этой задачи значение угла намагничивания, полученное значение целевой функции δ (целевая функция δ характеризует степень сходимости модельного и наблюденного магнитных полей) вместе с полученными значениями индуцированной и остаточной намагниченности записывают в отдельный файл. Выбирают шаг изменения угла намагниченности (например, 10°). Далее задается следующее значение угла в этом блоке, то есть 10° и решается обратная задача. При каждой последующей итерации значение угла, полученное значение целевой функции δ с соответствующими значениями индуцированной и остаточной намагниченности дополняют в созданный файл. Итерации осуществляют до тех пор, пока не будет пройден полный круг (360°), причем в нижнем полупространстве (от 0° до 180°) значения намагниченностей, всегда положительны, а в верхнем полупространстве (от 180° до 360°) значения намагниченностей, всегда отрицательны. В собранных в указанном файле данных по первому блоку определяют минимальное значение целевой функции δ, которое и определяет значения векторов индуцированной (Ji) и остаточной (Jr) намагниченностей. Полученное значение угла намагничивания присваивается исследуемому блоку.These procedures are carried out in several iterative cycles. The first iteration begins by setting the magnetization angle of zero for the first block of the section copy and solving the inverse problem. In the course of solving this problem, the value of the magnetization angle, the obtained value of the objective function δ (the objective function δ characterizes the degree of convergence of the model and observed magnetic fields), together with the obtained values of the induced and residual magnetization, are recorded in a separate file. Select the step of changing the angle of magnetization (for example, 10 °). Next, the next angle value in this block is set, that is, 10 ° and the inverse problem is solved. At each subsequent iteration, the angle value, the obtained value of the objective function δ with the corresponding values of the induced and residual magnetization are added to the created file. Iterations are carried out until a full circle (360 °) is reached, and in the lower half-space (from 0 ° to 180 °) the magnetization values are always positive, and in the upper half-space (from 180 ° to 360 °) the magnetization values always negative. In the data collected in the indicated file for the first block, the minimum value of the objective function δ is determined, which determines the values of the induced (Ji) and residual (Jr) magnetization vectors. The obtained value of the magnetization angle is assigned to the unit under investigation.
Для каждого последующего блока повторяют те же циклы до тех пор, пока не будет пройдена вся процедура для каждого из блоков, составляющих исследуемый разрез. Причем, для каждого последующего блока значение целевой функции δ будет либо равно, либо меньше чем для предыдущего. Такой итерационный цикл обозначается как полный.For each subsequent block, the same cycles are repeated until the entire procedure for each of the blocks making up the studied section is completed. Moreover, for each subsequent block, the value of the objective function δ will be either equal to or less than for the previous one. Such an iterative cycle is designated as complete.
Полные циклы повторяют до тех пор, пока не будет достигнуто минимальное значение целевой функции δ или, что тоже самое, возможность изменения углов перестанет существовать. Полученные при этом значения векторов намагниченностей соответствуют векторам искомых индуцированной (Ji) и остаточной (Jr) намагниченностям.Complete cycles are repeated until the minimum value of the objective function δ is reached, or, what is the same, the possibility of changing the angles ceases to exist. The obtained values of the magnetization vectors correspond to the vectors of the desired induced (Ji) and residual (Jr) magnetizations.
Далее, решая прямую задачу, осуществляют построение магнитных полей изучаемого объекта: магнитного поля, обусловленного индуцированной намагниченностью, и магнитного поля, обусловленного остаточной намагниченностью.Further, solving the direct problem, the magnetic fields of the object under study are constructed: the magnetic field due to the induced magnetization, and the magnetic field due to the residual magnetization.
Пример модельной реализации способа.An example of a model implementation of the method.
Модельная реализация способа проведена на базе авторской программы "Geolab", см. "Садур О.Г. Моделирование геологических сред на основе вычисления их плотности и магнитных характеристик в классе сложного распределения масс при решении различных геологических задач // Геология и минерально-сырьевые ресурсы Сибири. - 2012. - №1(9). - С. 96-101) с использованием компьютерного устройсва.A model implementation of the method was carried out on the basis of the author's program "Geolab", see "Sadur OG Modeling of geological media based on the calculation of their density and magnetic characteristics in the class of complex mass distribution when solving various geological problems // Geology and Mineral Resources of Siberia . - 2012. - No. 1 (9). - S. 96-101) using a computer device.
На фиг. 1 показана принятая, исходная модель магнитного объекта, разделенная в данном примере на два разреза, в каждом из которых содержится по три пронумерованных блока: (1и-3и) - разрез с заданной индуктивной намагниченностью Ji и (1о-3о) - разрез с заданной остаточной намагниченностью Jr. Для каждого из указанных блоков заданы значения намагниченностей и углы намагничивания, iix, irx, In FIG. Figure 1 shows the adopted, initial model of a magnetic object, divided in this example into two sections, each of which contains three numbered blocks: (1i-3i) - a section with a given inductive magnetization Ji and (1o-3o) - a section with a given residual magnetization jr. For each of these blocks, magnetization values and magnetization angles, i ix, i rx,
Ji - индуцированная намагниченность в А/м,Ji is the induced magnetization in A / m,
iix - направление вектора индуцированной намагниченности (Ji) в градусах,i ix is the direction of the induced magnetization vector (Ji) in degrees,
Jr - остаточная намагниченность в А/м,Jr is the residual magnetization in A / m,
irx - направление вектора остаточной намагниченности (Jr) в градусах.i rx is the direction of the residual magnetization vector (Jr) in degrees.
Как видно из фиг. 1, все значения остаточной намагниченности (Jr) положительны и находятся в нижнем полупространстве, то есть направления векторов irx остаточной намагниченности (Jr) меняются от 0° до 180°.As can be seen from FIG. 1, all values of the residual magnetization (Jr) are positive and are in the lower half-space, that is, the directions of the vectors i rx of the residual magnetization (Jr) change from 0 ° to 180 °.
От модели, представляющей собой совокупность двух разрезов с индуцированной (Ji) и остаточной (Jr) намагниченностями, посчитана прямая задача и получено поле ΔT (фиг. 1, поз. 1).. Аналогично посчитано поле ΔT от разрезов с индуцированной намагниченностью Ji и с остаточной намагниченностью. Jr. (фиг. 1, поз 2, поз. 3, соответственно)From the model, which is a combination of two sections with induced (Ji) and residual (Jr) magnetizations, we calculated the direct problem and obtained the field ΔT (Fig. 1, item 1). Similarly, we calculated the field ΔT from the sections with induced magnetization Ji and с residual magnetization. Jr. (Fig. 1, pos. 2, pos. 3, respectively)
На фиг. 2 приведены графики, иллюстрирующие результаты первого итерационного цикла при определении векторов намагничивания (Ji, Jr) для блоков (1и и 1о) исходной модели. На графиках представлено изменение параметров намагничивания от изменения углов irx вектора Jr остаточной намагниченности. Показано, что в первом итерационном цикле для блоков (1и) минимум целевой функции δ составил 46.2. При этом значение индуцированной намагниченности (Ji) составило 1.572 А/м (задано 2 А/м), угол (irx) остаточной намагниченности равняется 0°, значение Jr=0,455 А/м (при заданных значениях, соответственно, 10° и 1 А/м).In FIG. Figure 2 shows graphs illustrating the results of the first iterative cycle when determining the magnetization vectors (Ji, Jr) for blocks (1i and 1o) of the original model. The graphs show the change in the magnetization parameters from a change in the angles i rx of the residual magnetization vector Jr. It is shown that in the first iteration cycle for blocks (1i), the minimum of the objective function δ was 46.2. In this case, the value of the induced magnetization (Ji) was 1.572 A / m (2 A / m set), the angle (i rx ) of the residual magnetization is 0 °, the value Jr = 0.455 A / m (for given values, respectively, 10 ° and 1 A / m).
На фиг. 3 приведены графики, иллюстрирующие результаты первого итерационного цикла определения векторов намагниченности для блоков 2и и 2о исходной модели (фиг. 1). Согласно иллюстрации видно, что в первом итерационном цикле искомые значения намагниченностей Ji и Jr не совпадают с заданными, хотя угол (irx) вектора остаточной намагниченности (Jr) уже соответствует модельному значению.In FIG. Figure 3 shows graphs illustrating the results of the first iterative cycle for determining the magnetization vectors for blocks 2i and 2o of the original model (Fig. 1). According to the illustration, it can be seen that in the first iteration cycle, the desired magnetization values Ji and Jr do not coincide with the set values, although the angle (i rx ) of the remanent magnetization vector (Jr) already corresponds to the model value.
Аналогичная ситуация (фиг. 4) сложилась и для блока 3. При значении целевой функции δ=0.56, значения намагниченностей Ji, Jr и направление вектора irx остаточной намагниченности (Jr) не совпадают с модельными значениями.A similar situation (Fig. 4) has developed for
Последующее уменьшение значений целевой функции δ за счет использования второго полного итерационного цикла (фиг. 5 - фиг. 7) позволило получить значения и углы irx намагничивания, практически совпадающие с модельными значениями, что обеспечило разделение эффективной намагниченности на две составляющие, индуцированную намагниченность (Ji), и остаточную (Jr).The subsequent decrease in the values of the objective function δ due to the use of the second full iterative cycle (Fig. 5 - Fig. 7) made it possible to obtain magnetization values and angles i rx that practically coincide with the model values, which ensured the separation of the effective magnetization into two components, the induced magnetization (Ji ), and residual (Jr).
Пример практической реализации способа.An example of a practical implementation of the method.
На фиг. 8 показаны геоплотностая (фиг. 8, А) и полученные при реализации способа согласно изобретению геомагнитные модели (фиг. 8, Б, В, Г) в северной части Предверхоянского прогиба (Якутия) по сейсмическому профилю 140305.In FIG. Fig. 8 shows the density density (Fig. 8, A) and the geomagnetic models (Fig. 8, B, C, D) obtained in the implementation of the method according to the invention in the northern part of the Pre-Verkhoyansk trough (Yakutia) along seismic profile 140305.
По данным МОГТ - разрез двухслойный, имеет блоковое строение. Он состоит из осадочных пород (чехол) пермотриасовского, юрского и мелового возрастов примерно одинакового состава, песчаников, глин, аргиллитов и, ограниченно, карбонатных горных пород. Плотностные характеристики (сверху вниз) увеличиваются от 2.54 г/см3 до 2.66 г/см3, а магнитная восприимчивость меняется от 0 до 2.5 *10-5СИ. Кристаллический фундамент составлен породами позднего архея с наложенными впадинами вулкано-плутоническими поясами раннего протерозоя. На разрезе пробурена скважина, вскрывшая фундамент. По сейсмическому профилю создана конструкция геоплотностной и трех геомагнитных моделей, как это показано на фиг. 8.According to the MOGT, the section is two-layer and has a block structure. It consists of sedimentary rocks (cover) of the Permian-Thassian, Jurassic and Cretaceous ages of approximately the same composition, sandstones, clays, mudstones and, to a limited extent, carbonate rocks. Density characteristics (from top to bottom) increase from 2.54 g / cm 3 to 2.66 g / cm 3 , and the magnetic susceptibility varies from 0 to 2.5 * 10-5SI. The crystalline basement is composed of rocks of the Late Archean with superimposed hollows of the volcanic-plutonic belts of the Early Proterozoic. A well was drilled in the section, which revealed the foundation. According to the seismic profile, the design of the geo-density and three geomagnetic models was created, as shown in FIG. 8.
На геоплотностной модели (фиг. 8, А) верхняя часть (чехол) делится на две части: северо-западная часть, попадающая на Оленекское поднятие (ПК 0-140 км), сравнительно однородная и составляет 2.74 г/см3, и неоднородная юго-восточная часть (ПК 140-280 км), совпадающая по координатам с Предверхоянским прогибом, с изменением плотностей от 2.54 до 2.66 г/см3. Нижняя часть осадочных отложений в юго-восточной части разреза (от 240 км до окончания разреза) представлена увеличенными значениями плотностных характеристик от 2.64 до 2.76 г/см3 с изменением направления изолиний плотности вплоть до горизонтального, что соответствует предполагаемому надвигу, сформированному складчатыми образованиями Верхоянского и Сетте-Дебанского антиклинория. Кристаллический фундамент по структурно-морфологическим особенностям плотностного поля разбит (показано пунктирными линиями) на три блока, из которых первый блок (0-100 км) характеризуется плавным увеличением поля плотности на глубину (2.74 до 2.96 г/см3), второй блок (100-280 км) характеризуется волнообразным изменением увеличением плотности на глубину (2.74 до 2.96 г/см3), третий блок характеризуется также увеличением плотности на глубину, но с меньшими значениями (от 2.72 до 2.76 г/см3).In the density-density model (Fig. 8, A), the upper part (cover) is divided into two parts: the north-western part falling on the Olenek uplift (PK 0-140 km) is relatively homogeneous and is 2.74 g / cm 3 and non-uniform south -the tail part (PK 140-280 km), coinciding in coordinates with the Pre-Verkhoyansk trough, with a density change from 2.54 to 2.66 g / cm 3 . The lower part of sedimentary deposits in the southeastern part of the section (from 240 km to the end of the section) is represented by increased values of the density characteristics from 2.64 to 2.76 g / cm 3 with a change in the direction of the density contours up to horizontal, which corresponds to the estimated thrust formed by the folded formations of Verkhoyansk and Sette Deban Anticlinorium. According to the structural and morphological features of the density field, the crystalline basement is divided (shown by dashed lines) into three blocks, of which the first block (0-100 km) is characterized by a smooth increase in the density field to a depth (2.74 to 2.96 g / cm 3 ), the second block (100 -280 km) is characterized by a wave-like change in density increase to a depth (2.74 to 2.96 g / cm 3 ), the third block is also characterized by an increase in density to a depth, but with lower values (from 2.72 to 2.76 g / cm 3 ).
На построенной геомагнитной модели (фиг. 8, Б) с полученной эффективной намагниченностью в верхней части разреза от 260 км и далее на юго-восток в виде субгоризонтальных изолиний, с повышенными значениями намагниченности, более 0.3 А/м, выделяется надвиговая зона, которая достаточно четко проявляется и на геоплотностной модели (фиг. 8, А). Нижняя часть разреза (кристаллический фундамент, граница уточнена на основе решения обратной задачи) по структурно-морфологическим особенностям поведения изолиний намагниченности разбивается (показано пунктирными линиями) на четыре блока: первый блок (от 0 до 105 км) определяется двумя локальными аномалиями намагниченности: первая - с уменьшением ее значений от периферии к центру до глубины 5 км (от 0.5 до 0,3 А/м), а вторая - с пониженными значениями намагниченности по кровле фундамента (до 0.3 А/м) и с дальнейшими повышениями значений на глубину до 0.7 А/м. Второй блок (от 105 до 220 км) представлен тремя локальными аномалиями с направлением изолиний от субвертикальных и субгоризонтальных в верхней части разреза до субгоризонтальных - в нижней части, с общим понижением значений намагниченности на глубину (от 0.6 А/м до 0.05 А/м). Третий блок (от 220 до 285 км) характеризуется мозаичной аномалией со значениями намагниченностей от 0,05 до 0,6 А/м и с увеличением намагниченности по глубине. Горизонтальные размеры четвертого блока составляют около 20 км. Границы блока на восток выходят за пределы разреза. Отличительной чертой структуры этого поля намагниченностей является их пониженные значения (от 0,3 до 0,5 А/м) с углами направления изолиний эффективной намагниченности около 60° на северо-запад, что, очевидно, связано с влиянием зоны надвига.In the constructed geomagnetic model (Fig. 8, B) with the obtained effective magnetization in the upper part of the section from 260 km and further to the southeast in the form of subhorizontal contours, with increased magnetization values of more than 0.3 A / m, a thrust zone is selected that is sufficient clearly manifests itself in the geo-density model (Fig. 8, A). The lower part of the section (crystalline basement, the boundary was refined based on the solution of the inverse problem) is divided into four blocks according to structural and morphological features of the behavior of the magnetization contours (shown by dashed lines): the first block (from 0 to 105 km) is determined by two local magnetization anomalies: the first - with a decrease in its values from the periphery to the center to a depth of 5 km (from 0.5 to 0.3 A / m), and the second with lower magnetization values along the foundation roof (to 0.3 A / m) and with further increases in depth about 0.7 A / m. The second block (from 105 to 220 km) is represented by three local anomalies with the direction of isolines from subvertical and subhorizontal in the upper part of the section to subhorizontal in the lower part, with a general decrease in the magnetization values to a depth (from 0.6 A / m to 0.05 A / m) . The third block (from 220 to 285 km) is characterized by a mosaic anomaly with magnetization values from 0.05 to 0.6 A / m and with an increase in magnetization in depth. The horizontal dimensions of the fourth block are about 20 km. The boundaries of the block east extend beyond the section. A distinctive feature of the structure of this magnetization field is their lower values (from 0.3 to 0.5 A / m) with angles of direction of the isolines of the effective magnetization of about 60 ° to the northwest, which is obviously due to the influence of the thrust zone.
Модель эффективной намагниченности (фиг. 8, Б) по характерным точкам, с учетом данных сейсморазведки разбита на 130 цифровых блоков. Затем, согласно изобретению, как это описано выше на теоретической модели (фиг. 1-фиг. 7), разрез с полученной эффективной намагниченностью (фиг. 8. Б) копировался. Для обоих полученных разрезов использовался график наблюденного магнитного поля ΔТ (фиг. 8, Б). Назначены разрезы для определения значений индуцированной (Ji), и остаточной (Jr) намагниченностей. Первоначально, как для индуцированной (Ji), так и для остаточной (Jr) намагниченности угол намагничивания совпадал с направлением нормального магнитного поля Земли (80°). Для определения значений намагниченности (Ji, Jr) при заданном угле irx намагничивания в каждом блоке решалась обратная задача. Полученные значения намагниченности в каждом блоке строго соответствовали значениям эффективной намагниченности, уменьшенной вдвое. Затем, в назначенном первом блоке, в разрезе, где определяется остаточная намагниченность (Jr), задавалось первое значение угла намагничивания (irx=0°) и решалась обратная задача. Значения индуцированной Ji, остаточной Jr намагниченностей и целевой функции δ записывалась в отдельный файл. Далее, задавалось следующее значение угла irx и вышеописанная процедура повторялась снова. Данная операция продолжалась до тех пор, пока был не пройден полный круг (от 0° до 360°). Следует особо отметить, что в определяемом нижнем полупространстве (от 0° до 180°) создавались условия для определения только положительных значений намагниченности, а в верхнем полупространстве (от 180° до 360°) - только отрицательных.The model of effective magnetization (Fig. 8, B) at characteristic points, taking into account seismic data, is divided into 130 digital blocks. Then, according to the invention, as described above in the theoretical model (Fig. 1-Fig. 7), the section with the obtained effective magnetization (Fig. 8. B) was copied. For both obtained sections, the graph of the observed magnetic field ΔТ was used (Fig. 8, B). Sections are assigned to determine the values of the induced (Ji) and residual (Jr) magnetizations. Initially, for both the induced (Ji) and the residual (Jr) magnetizations, the magnetization angle coincided with the direction of the Earth's normal magnetic field (80 °). To determine the magnetization values (Ji, Jr) for a given magnetization angle i rx, the inverse problem was solved in each block. The obtained magnetization values in each block strictly corresponded to the values of the effective magnetization, halved. Then, in the assigned first block, in the section where the residual magnetization is determined (Jr), the first value of the magnetization angle (i rx = 0 °) was set and the inverse problem was solved. The values of the induced Ji, the residual magnetization Jr, and the objective function δ are recorded in a separate file. Next, the following value of the angle i rx was set and the above procedure was repeated again. This operation continued until a full circle was completed (from 0 ° to 360 °). It should be especially noted that in the determined lower half-space (from 0 ° to 180 °), conditions were created for determining only positive magnetization values, and in the upper half-space (from 180 ° to 360 °) - only negative.
Окончательно, в файле, куда записывались значения намагниченностей и целевой функции δ, определялось минимальное значение целевой функции δ и соответствующий ей угол (irx) намагничивания первого блока. Полученное значение угла намагничивания (irx) закреплялось за первым блоком.Finally, in the file where the values of the magnetizations and the objective function δ were written, the minimum value of the objective function δ and the corresponding magnetization angle (i rx ) of the first block were determined. The obtained value of the magnetization angle (i rx ) was assigned to the first block.
Аналогичные операции проводились для всех 130 цифровых блоков (полный итерационный цикл). Всего для получения моделей с разделением полей эффективной намагниченности на поля индуцированной и остаточной намагниченностей потребовалось шесть полных итерационных циклов.Similar operations were carried out for all 130 digital blocks (full iterative cycle). In total, six complete iterative cycles were required to obtain models with the separation of the fields of effective magnetization into the fields of the induced and residual magnetizations.
Далее полученные цифровые модели (данные) с помощью известных процедур были преобразованы в графические построения и визуализировались в виде карт изолиний, представляющих собой распределение эффективной, индуцированной и остаточной намагниченностей в плоскости разреза (фиг. 8, Б, В, Г).Further, the obtained digital models (data) using well-known procedures were transformed into graphical constructions and visualized in the form of contour maps representing the distribution of effective, induced and residual magnetizations in the section plane (Fig. 8, B, C, D).
На полученных геомагнитных моделях видно, что наиболее дифференцировано и неоднородно в фундаменте магнитное поле индуцированной намагниченности (фиг. 8, В.) Субвертикальное распределение изолиний индуцированной и остаточной намагниченностей в центральной части разреза (2-й блок, пикет 120-165 км), ограниченны градиентными зонами, свидетельствует о наличие здесь тектонической зоны. Как это видно из фиг. 8, В-фиг. 8, Г, более узкая часть тектонической зоны (150-165 км), ограниченная еще более высокоградиентными зонами намагниченности, обладает пониженными значениями намагниченности (0-0.25 А/м) в поле индуцированной намагниченности (Ji) и более высокими (0,48-0,54 А/м) - в поле остаточной намагниченности (Jr). Сопоставляя эти данные с данными по плотностной модели (фиг. 8, А), видно, что эта тектоническая зона отмечается повышенными плотностными характеристиками от 2.92 до 2.76 г/см3. По-видимому, в настоящее время функционирует только эта часть тектонической зоны, а северо-западная часть тектонической зоны "залечена" геологическими процессами. Второе тектоническое нарушение, не выходящее на поверхность, расположено к юго-востоку от первого, на отметке 205 км. В поле остаточной намагниченности (фиг. 8, Г) это нарушение не наблюдается. Аналогично, в поле индуцированной намагниченности (фиг. 8, В) (чехол, юго-восточная часть разреза) субгоризонтальное расположение изолиний наиболее ярко подтверждает наличие надвиговой системы, в то время как в поле остаточной намагниченности в структуре изолиний влияние надвига не проявляется (фиг. 8, Г). Очевидно, что здесь, как и для второго тектонического нарушения, наблюдается разрушение остаточной намагниченности, за счет деформации и сдвигов горных пород при высоком горизонтальном давлении, (см., например Валеев К.А., Абсалямов С.С. «Остаточная намагниченность магнетита при воздействии высоких давлений и сдвиговых деформаций «// Физика Земли. 2000. №3. С. 59-64).The obtained geomagnetic models show that the magnetic field of the induced magnetization is most differentiated and heterogeneous in the foundation (Fig. 8, B.). The vertical distribution of the contour lines of the induced and residual magnetization in the central part of the section (
Таким образом, способ, согласно изобретению, позволяет повысить информативность магниторазведочных работ за счет нового качества интерпретации магниторазведочных данных, а именно за счет возможности раздельной интерпретации магнитных полей индуцированной и остаточной намагниченности изучаемого объекта, что, в целом, позволяет значительно улучшить качество прогноза при поисках полезных ископаемых.Thus, the method according to the invention improves the informational content of magnetic exploration due to the new quality of interpretation of magnetic exploration data, namely due to the possibility of separate interpretation of the magnetic fields of the induced and residual magnetization of the studied object, which, in general, can significantly improve the quality of the forecast when searching for useful fossils.
Достоверность способа, согласно изобретению, проверена автором на многочисленных моделях.The reliability of the method according to the invention is verified by the author on numerous models.
Claims (14)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2018112370A RU2683817C1 (en) | 2018-04-05 | 2018-04-05 | Method for determining induced and residual magnetization of rocks according to magnetic exploration data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2018112370A RU2683817C1 (en) | 2018-04-05 | 2018-04-05 | Method for determining induced and residual magnetization of rocks according to magnetic exploration data |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| RU2683817C1 true RU2683817C1 (en) | 2019-04-02 |
Family
ID=66089754
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| RU2018112370A RU2683817C1 (en) | 2018-04-05 | 2018-04-05 | Method for determining induced and residual magnetization of rocks according to magnetic exploration data |
Country Status (1)
| Country | Link |
|---|---|
| RU (1) | RU2683817C1 (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117148457A (en) * | 2023-08-29 | 2023-12-01 | 长安大学 | Magnetic layer magnetization modulus calculation method and system |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU1125578A1 (en) * | 1983-03-05 | 1984-11-23 | Южное Отделение Всесоюзного Научно-Исследовательского Института Геофизических Методов Разведки | Method of determination of magnetization of mountain rock in natural bedding thereof |
| US5075625A (en) * | 1989-10-10 | 1991-12-24 | Total Compagnie Francaise Des Petroles | Procedure and device for the detection of inversions of the earth's magnetic field by means of measurement taken in a drill shaft |
| RU2303799C1 (en) * | 2005-10-31 | 2007-07-27 | Институт космофизических исследований и аэрономии СО РАН | Method for remote definition of vertical component of magnetic induction of rocks in ancient epochs particularly to paleomagnetism |
| RU2310891C2 (en) * | 2005-10-31 | 2007-11-20 | Институт космофизических исследований и аэрономии СО РАН | Method for remotely determining horizontal component of magnetic induction of rock in ancient epochs |
-
2018
- 2018-04-05 RU RU2018112370A patent/RU2683817C1/en active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| SU1125578A1 (en) * | 1983-03-05 | 1984-11-23 | Южное Отделение Всесоюзного Научно-Исследовательского Института Геофизических Методов Разведки | Method of determination of magnetization of mountain rock in natural bedding thereof |
| US5075625A (en) * | 1989-10-10 | 1991-12-24 | Total Compagnie Francaise Des Petroles | Procedure and device for the detection of inversions of the earth's magnetic field by means of measurement taken in a drill shaft |
| RU2303799C1 (en) * | 2005-10-31 | 2007-07-27 | Институт космофизических исследований и аэрономии СО РАН | Method for remote definition of vertical component of magnetic induction of rocks in ancient epochs particularly to paleomagnetism |
| RU2310891C2 (en) * | 2005-10-31 | 2007-11-20 | Институт космофизических исследований и аэрономии СО РАН | Method for remotely determining horizontal component of magnetic induction of rock in ancient epochs |
Non-Patent Citations (1)
| Title |
|---|
| Вычислительные математика и техника в разведочной геофизике, Справочник геофизика. под ред. В.И. Дмитриева. М., Недра, 1982. с. 89-98. * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117148457A (en) * | 2023-08-29 | 2023-12-01 | 长安大学 | Magnetic layer magnetization modulus calculation method and system |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Khesin et al. | Interpretation of geophysical fields in complicated environments | |
| Pilkington et al. | Mitigating remanent magnetization effects in magnetic data using the normalized source strength | |
| Liu et al. | Particle swarm optimization inversion of magnetic data: Field examples from iron ore deposits in China | |
| Eppelbaum et al. | Geophysical studies in the Caucasus | |
| CN105510993A (en) | Foreland basin deep buried and compressed type complex gypsum-salt rock identification and distribution prediction method | |
| US10534108B2 (en) | Methods and systems for the inversion of magnetic data from remnant and induced sources in geophysical exploration | |
| Abdelazeem et al. | Magnetometric identification of sub-basins for hydrocarbon potentialities in Qattara Ridge, North Western Desert, Egypt | |
| RU2683817C1 (en) | Method for determining induced and residual magnetization of rocks according to magnetic exploration data | |
| Šimonová et al. | Determination of rock densities in the Carpathian-Pannonian Basin lithosphere: based on the CELEBRATION 2000 experiment | |
| Sarvandani et al. | Investigating subsurface structures of Gachsaran oil field in Iran using 2D inversion of magnetotelluric data | |
| Okiwelu et al. | 3D modelling and basement tectonics of the Niger Delta basin from aeromagnetic data | |
| Ghosh | Interpretation of gravity anomaly and crustal thickness mapping of Narmada-Son lineament in central India | |
| Luo et al. | Study on the distribution characteristics of faults and their control over petroliferous basins in the China seas and its adjacent areas | |
| Azab et al. | A new contribution to the geology of the Egyptian Red Sea shelf using geophysical data | |
| Guo et al. | Apparent magnetization mapping in the presence of strong remanent magnetization: The space-domain inversion approach | |
| Eppelbaum et al. | Some common aspects of magnetic, induced polarization and self-potential anomalies interpretation: implication for ore target localization | |
| CN109707375A (en) | A method of passing through underground magnetic gradient tensor data identification underground magnetosphere information | |
| Gonzalez et al. | Magnetic characterization and signature of the basement of eastern Venezuela: Espino Graben | |
| Akca et al. | Imaging gypsum–anhydrite interface by a structure‐based conceptual resistivity model: A case study from Bala (Turkey) gypsum deposits | |
| Egbelehulu et al. | Interpretation of Aeromagnetic Data of Part of Gwagwalada Abuja Nigeria for Potential Mineral Targets | |
| Shu‐Ling et al. | Research on Pre‐Cenozoic Residual Basin Structure and Basement Property of the Huanghua Depression | |
| Darijani et al. | Synthetic modeling and joint inversion of gravity and seismic refraction data for overburden stripping in the Athabasca Basin, Canada | |
| Totoş | 3D modelling and inversion of airborne gravity gradiometry and aeromagnetic data from Budgell Harbour, North-Central Newfoundland | |
| Seminsky et al. | Geoelectrical images of normal fault zones: tectonophysical interpretation of the shallow-depth electrical resistivity tomography data on the Buguldeika-Chernorud graben in the Western Baikal region | |
| Ahmed et al. | Groundwater Inspection and Aquifer Assessment Using Magnetotellurics and Magnetic Data at the Reclamation Area Around New Sphinx City, Egypt |