重力异常正演资料

布格重力异常布格重力异常（Bouguer gravity anomaly）：重力仪的观测结果(g测)，经过纬度改正(g 纬)、高度改正(g高)，中间层改正(g中)和地形改正(g形)以后，再减去正常重力值（γ）后所得到的重力差(Δg= g测- g高- g中- g纬- g形-γ)称为布格重力异常。

布格重力异常一般用相对重力测最方法获得。

如某点对总基点的实测重力差为Δg'，则布格重力异常为Δg=Δg'+ g高+ g 中+ g纬+ g形（g高、g中、g纬、g形分别为测点与总基点的相对高度、中间层、纬度和地形变化的改正项)。

布格重力异常资料是重力勘探的基础资料。

重力（g）观测首先要从地球内部研究：地壳，地幔，地核。

地壳分为上地壳和下地壳，上地壳主要由密度较小的富含硅铝岩的花岗岩组成，称为花岗岩层，平均密度为2.7g/cm^3;下地壳是由富含硅镁的玄武岩组成，称为玄武层，平均密度 3.1/cm^3。

花岗层与玄武层的分界面称为康德界面，简称”康氏面“。

从地壳底部一直向下倒2900km处称为地幔，地幔也分为上地幔和下地幔，上地幔深度为33-900km，主要为基性岩，超基性岩，物质密度在3.32~4.49g/cm^3，下地幔物质则主要为铁镍等金属氧化物，深度为900~2700km，密度约为4.60~5.68km，地幔与地壳的分界面称为莫霍维奇面，简称”莫霍面“。

从2700km以下一直到地心称为地核，地核又分为外核，过渡带，内核三部分，组成地核的物质目前不清楚，推测主要为铁镍物质，故称为铁镍核，密度约为从9.69g/cm^3曾至12。

17g/cm^3.地球平均密度为5.53g/cm^3。

重力（g）观测可以通过对地壳和地幔的勘测，推断区域布格重力异常和均衡重力异常及地壳深部的构造起伏情况。

球体重力异常正演程序报告球体重力异常正演是地球物理学中的一种重要方法，用于研究地下物质分布和地球内部结构。

本报告将重点介绍球体重力异常正演程序的原理、步骤和应用。

一、原理球体重力异常正演是基于牛顿引力定律和球体模型的数学计算方法。

根据牛顿引力定律，在球体表面上的任意一点，重力加速度可以表示为：g = G * (M / r^2)其中，g为重力加速度，G为引力常数，M为球体的质量，r为球心到该点的距离。

根据球体模型，球体的质量可以表示为：M = (4/3) * π * ρ * R^3其中，ρ为球体的密度，R为球体的半径。

将质量公式代入重力加速度公式，可得到球体表面上的重力加速度公式：g = (4/3) * G * π * ρ * R / r^2二、步骤球体重力异常正演程序的步骤如下：1. 确定观测点的位置和高度，以及球体模型的半径和密度。

2. 计算球体表面上的重力加速度，根据上述公式进行计算。

3. 根据观测点与球心的距离，计算球体表面上的重力加速度的投影值。

4. 重复步骤3，直到计算出所有观测点的重力加速度投影值。

5. 计算观测点的球体重力异常值，即观测点的重力加速度减去球体表面上的重力加速度投影值。

三、应用球体重力异常正演程序在地球物理勘探中有广泛的应用，主要包括以下几个方面：1. 地质勘探：通过球体重力异常正演，可以对地下的岩石密度分布进行推测，从而帮助地质勘探人员确定地质构造和找到潜在的矿产资源。

2. 油气勘探：油气藏通常与地下的密度异常有关，通过球体重力异常正演，可以对潜在的油气藏进行初步判断，指导油气勘探的方向和深度。

3. 地壳构造研究：地球内部的构造和演化与地下岩石的密度分布密切相关，通过球体重力异常正演，可以揭示地壳的变形和演化过程，为地壳构造研究提供重要的参考依据。

4. 火山和地震研究：火山和地震活动通常与地下的岩浆和断层有关，球体重力异常正演可以帮助科学家们理解火山和地震的发生机制，预测可能的灾害风险。

关于布格重力异常计算及资料处理与反演和解释的报告姓名：***班级：061084-27学号：**********指导老师：***日期：2011.4.14目录前言 (2)目的 (2)任务要求 (2)工作过程 (2)成果 (2)工作内容及步骤 (3)§1-布格重力异常计算 (3)§2-布格重力异常处理 (3)1.绘制平面等值线图 (3)2.异常处理（分离区域异常和局部异常） (6)§3-布格重力异常反演——特征点法反演 (11)§4-布格重力异常的解释 (13)评述与结论 (13)评述 (13)结论 (14)关于布格重力异常计算及资料处理与反演和解释的报告 前言目的：熟悉并掌握布格重力异常计算及资料处理与反演和解释 任务要求：根据在一个地区重力测量的结果，计算出布格重力异常，并根据异常进行资料处理和解释，并完成一份工作报告。

工作过程：(1)利用实测的相对重力值、相对高程值和X,Y 坐标值，计算各种校正（地形校正除外），纬度校正用 计算，自由空间（或高度）校正用计算，中间层校正用 计算，已知地表物质密度为2.50g/cm 3，起算点纬度为45°；(2)获得各点处的布格重力异常值后，绘出平面等值线图，等值线距为0.5mGal ；(3)根据异常（平面或剖面）特征，选用适当的方法进行处理（如压制干扰、消除区域场等）进行处理，并对处理效果进行描述； (4)将处理后的异常进行反演；(5)写出全部过程和所采用的处理与反演方法之应用理由。

成果：根据布格重力异常数据计算及资料处理与反演初步结果判断，该异常应由地区下一球体引起，球体埋深98.8m ，剩余质量t 6103.07⨯，球体中心在地面的投影点坐标为（248.8,248.8）m 。

):,()2sin(814.0mkX mGal X g ∆⋅-=ϕδϕ):,(3086.0m h mGal h g f ∆⋅=δ):/:,(0419.03m h cm g mGal h g ρρδσ∆⋅-=工作内容及步骤§1-布格重力异常计算用excel先算出各项校正（除地形校正外），需要注意的是在纬度校正中为测点到总基点间纬向距离，由于测点都位于总基点以北，故取正值；为总基点纬度即45°; 单位要划为km。

简述重力场的正反演问题
重力场的正反演问题涉及重力异常的正演和反演。

正演问题是给定地下某种地质体的形状、产状和剩余密度等，通过理论计算来求得它在地面上产生的异常大小、特征和变化规律，这是正向思维的问题。

反演问题则是依据已获得的异常特征、数值大小、分布情形等并结合物性资料来求解地下地质体的形状和空间位置等，这是逆向思维的问题。

重力正演是指根据地下地质体的形状、大小、密度等物理参数，利用重力场理论计算其在地球表面产生的重力异常。

重力反演则是根据实测的重力异常数据，结合物性资料，推断地下地质体的形状、大小、空间位置等信息。

重力正演是解决正问题的过程，它从地下地质体的物理参数出发，预测其在地球表面产生的重力异常。

重力反演则是解决反问题的过程，它从实测的重力异常数据出发，推断地下地质体的形状、大小、空间位置等信息。

重力场的正反演问题在地球物理学中具有重要的应用价值，例如在矿产资源勘探、地质构造研究、地下水资源调查等领域都有广泛的应用。

通过正反演问题的解决，可以更好地理解地球内部结构和动力学过程，为资源开发和环境保护提供科学依据。

三维盐丘模型重力异常正演与约束反演研究的开题报告一、研究背景地球物理勘探技术常常利用重力资料来探测地下构造与岩性变化，其中三维盐丘模型对勘探盆地油气资源具有重要意义。

盐丘是由地壳下积累的盐岩物质隆起而形成的椭圆形几何体，具有天然的覆盖形式与沉积填充的同位素标志等特点。

在目前的研究中，重力异常探测技术是三维盐丘模型分析的基础之一。

重力异常是由于地球内部质量分布不均造成的引力异常，可以通过建立数学模型来计算及反演。

然而，由于盐体的高密度与强反射特性，其在重力场中呈现出明显的异常表现。

在分析盐体结构时，需要对其重力异常进行定量计算及反演。

二、研究目的本研究旨在建立三维盐丘模型重力异常正演模型，通过对盐体结构特征进行研究，提高盆地油气资源的勘探效率与准确性，并采用约束反演技术对重力异常进行反演，以实现对盐体结构特征的定量分析。

三、研究内容及方法1. 建立三维盐丘模型：利用地震资料提取出盐体结构，在建立三维模型时采用地质解释、测量技术、成像技术及计算机模拟技术等多种手段，确定盐体的几何形态、大小及空间位置等特征。

2. 重力异常正演模型：通过网格化的方式，建立三维盐丘模型重力场分布的数学模型，并计算出盐体结构的重力异常。

3. 约束反演方法：采用先验信息、约束优化等数学方法，对重力异常进行反演，得到盐体结构的定量分析结果。

四、预期成果通过建立三维盐丘模型重力异常正演与约束反演模型，实现对盐体结构特征的定量研究，提高油气资源勘探效率与准确度，为油气勘探事业的发展提供相关技术支持。

五、研究意义本研究具有重要的理论与应用价值，对盐体结构的勘探及油气资源的发现与开采具有实际意义。

同时，通过建立数学模型及采用约束反演技术，对重力异常进行反演，提高分析结果的准确性与可靠性，在盆地油气资源开发中具有重要的应用前景。

地球物理学中的重力异常解释理论地球物理学是研究地球内部构造和物理学现象的学科，其中涉及到很多现代的科技手段。

其中，重力异常解释理论就是其中的一种重要的研究工具。

这种理论可以帮助地球物理学家更好地理解地球的内部结构和地质属性，并应用于地质勘探和自然资源的开发与利用。

一、重力异常的定义和分类重力异常是地球重力场在不同处产生的扰动，通俗地说，它就是地球表面处处的重量，与一定位置处平均重量之差。

重力异常可分为正负两类，正重力异常表示在正常重力场的基础上，某一区域的重力场比平均重力场增强，而负重力异常则相反。

值得注意的是，地球表面的重力异常分布是高度不均匀的，其中包含着丰富的地质信息。

二、形成重力异常的地质现象形成重力异常的地质现象主要有三种。

一种是密度不同的岩石侵占了另一种密度不同的岩石，这就形成了“重力胶合”，导致特定区域的重力场异常。

第二种情况是，在地球的构造运动过程中可能会形成溶洞、洞穴和空隙，这些洞穴等物质的密度较低，因此造成了区域重力场的异常。

第三种情况是地球内部物质组成与分布的异质性所导致，如地球内部自然放射性元素含量的不均匀分布，会形成类似的地球重力异常。

三、重力异常解释理论重力异常解释理论是一种利用重力异常研究区域地质信息的方法。

该理论是基于物理学原理构建的，对于区域的重力场异常进行解释和分析，以确定区域内地质物质的特征和分布。

该理论包括多种科学手段，如重力数据的采集和处理、地质勘探技术、地震勘探技术等。

在实践中，重力异常解释理论被广泛应用于地质勘探、矿产资源的开发和利用、遥感技术等领域。

例如，科学家使用重力异常解释理论研究了海洋地壳形成的过程和机制、大地构造和板块漂移等问题，并预测了地球内部某些物质的分布、某些石油和天然气资源的分布等。

此外，重力异常解释理论还被广泛应用于地震预测和监测等工作中。

四、结语地球物理学作为一门交叉学科，涉及的学科领域广泛，有着广泛的应用前景。

重力异常解释理论是其中的一种重要研究方法，通过研究地球表面区域重力场的异常来分析地质属性和分布等信息。

《应用地球物理学》课程作业基于MATLAB的球体重力正演程序实验报告1一程序简介本程序基于MATLAB软件的GUI模块编写，旨在实现球体重力正演结果的可视化分析。

MATLAB是一个高级的编程语言，其矩阵思想方便了地球物理的编程工作。

随着该语言和相应软件的发展，其内部也集成了许多模块，如该实验用到的GUI模块。

在该模块中，可以通过窗口、按键和赋值框等基本元素的组合，编写出可视化的应用程序，再配合MATLAB强大的作图功能，可以实现正演结果的展示与分析。

该程序应包含以下内容：1.可以自由输入参数，如球体半径，埋深和剩余密度。

2.可以计算出Δg、V ZZ、V XZ和V ZZZ这四种重力异常及其导数的对应值。

3.可以绘制剖面图及平面图两种图像。

二源程序由于GUI程序的头文件均大同小异，这里只列出赋值框及绘图按键的程序代码。

% --- Executes on selection change in popupmenu1.function popupmenu1_Callback(hObject, eventdata, handles)% hObject handle to popupmenu1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Hints: contents = cellstr(get(hObject,'String')) returns popupmenu1contents as cell array% contents{get(hObject,'Value')} returns selected item from popupmenu1s=get(hObject,'value');handles.s = s;guidata(hObject, handles);function edit1_Callback(hObject, eventdata, handles)% hObject handle to edit1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit1 as text% str2double(get(hObject,'String')) returns contents of edit1 as a doubled = str2double(get(hObject,'string'));handles.d = d;guidata(hObject, handles);function edit2_Callback(hObject, eventdata, handles)% hObject handle to edit2 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit2 as text% str2double(get(hObject,'String')) returns contents of edit2 as a doubler = str2double(get(hObject,'string'));handles.r = r;guidata(hObject, handles);function edit3_Callback(hObject, eventdata, handles)% hObject handle to edit3 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)% Hints: get(hObject,'String') returns contents of edit3 as text% str2double(get(hObject,'String')) returns contents of edit3 as a doublero = str2double(get(hObject,'string'));handles.ro = ro;guidata(hObject, handles);% --- Executes on button press in pushbutton1.function pushbutton1_Callback(hObject, eventdata, handles)% hObject handle to pushbutton1 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)ro = handles.ro;r = handles.r;d = handles.d;x=-2*d:2*d;G=6.67e-11;pi=3.14159;switch handles.scase 2z=4*ro*pi*G*r^3*d./(3*(x.^2+d^2).^1.5);plot(x,z.*1e6);xlabel('X/m');ylabel('\Deltag/g.u.');case 3z=4*ro*pi*G*r^3*(2*d^2-x.^2)./(3*(x.^2+d^2).^2.5);plot(x,z.*1e9);xlabel('X/m');ylabel('Vzz/E');case 4z=-4*ro*pi*G*r^3*(d.*x)./(x.^2+d^2).^2.5;plot(x,z.*1e9);xlabel('X/m');ylabel('Vxz/E');case 5z=4*ro*pi*G*r^3*(2*d^2-3.*x.^2)./(x.^2+d^2).^3.5;plot(x,z.*1e9);xlabel('X/m');ylabel('Vzzz/nMKS');end% --- Executes on button press in pushbutton3.function pushbutton3_Callback(hObject, eventdata, handles)% hObject handle to pushbutton3 (see GCBO)% eventdata reserved - to be defined in a future version of MATLAB% handles structure with handles and user data (see GUIDATA)ro = handles.ro;r = handles.r;d = handles.d;x=-2*d:2*d;y=-2*d:2*d;[X,Y]=meshgrid(x,y);G=6.67e-11;pi=3.14159;switch handles.scase 2c=4*ro*pi*G*r^3*d./(3*(X.^2+Y.^2+d^2).^1.5);contour(X,Y,c.*1e6,'showtext','on');colormap(winter);xlabel('X/m');yl abel('Y/m');case 3c=4*ro*pi*G*r^3*(2*d^2-X.^2-Y.^2)./(3*(X.^2+Y.^2+d^2).^2.5);contour(X,Y,c.*1e9,'showtext','on');colormap(winter);xlabel('X/m');yl abel('Y/m');case 4c=-4*ro*pi*G*r^3*(d.*X.*Y)./(X.^2+Y.^2+d^2).^2.5;contour(X,Y,c.*1e9,'showtext','on');colormap(winter);xlabel('X/m');yl abel('Y/m');case 5c=4*ro*pi*G*r^3*(2*d^2-3.*X.^2-Y.^2)./(X.^2+Y.^2+d^2).^3.5;contour(X,Y,c.*1e9,'showtext','on');colormap(winter);xlabel('X/m');yl abel('Y/m');end三程序运行结果假设埋深为100m，球体半径为10m，剩余密度为2kg/m3.Δg的剖面图和平面图如下：V zz的剖面图和平面图如下：V XZ的剖面图和平面图如下：V ZZZ的剖面图和平面图如下：。