下载文档原格式
物探认识实习实习报告
实习时间，2021年7月1日至2021年8月31日。
实习地点，某物探公司。
实习内容：
在这两个月的实习中，我有幸加入了某物探公司，参与了多个项目的实地勘察和数据分析工作。
在实习的过程中，我学到了许多关于地质勘探和数据处理方面的知识，也提升了自己的实际操作能力和团队合作能力。
首先，我参与了一些地质勘探项目的实地勘察工作。
在实地勘察中，我学会了如何正确使用地质勘探仪器和设备，如何进行地质剖面的绘制和数据采集等技能。
通过实际操作，我对地质勘探的流程和方法有了更深入的了解，并且在实践中不断提高了自己的操作技巧。
其次，我还参与了一些数据处理和分析的工作。
在这方面，我学到了如何使用地质勘探软件对采集到的数据进行处理和分析，包
括地质剖面的绘制、地层分析、勘探成果的评价等内容。
通过这些工作，我对地质勘探数据的处理和分析方法有了更深入的了解，也提升了自己的数据处理能力和分析能力。
在实习的过程中，我还有幸结识了一些优秀的同事，他们在工作中给予了我许多指导和帮助。
通过与他们的交流和合作，我不仅学到了许多专业知识和技能，也提高了自己的团队合作能力和沟通能力。
总的来说，这两个月的物探认识实习经历让我收获颇丰。
我不仅学到了许多专业知识和技能，也锻炼了自己的实际操作能力和团队合作能力。
我相信这段宝贵的实习经历将对我的未来职业发展产生积极的影响。
感谢公司给予我这次实习的机会，我会继续努力，不断提升自己，为公司的发展贡献自己的力量。
物探实习报告一、前言我于2023进行了为期两周的物探实习，通过这次实习，我对物探技术及其应用有了更深入的了解。
物探，即地球物理勘探，是利用各种物理场和方法探测地下地质结构、矿产资源和其他地下物质的技术。
在这次实习中，我参与了野外数据采集、数据处理和解释等工作，收获颇丰。
二、实习内容1. 野外数据采集野外数据采集是物探工作的基础。
我们使用了地震勘探、电磁法、重力法等多种方法进行数据采集。
在这个过程中，我学会了如何使用地震仪、电磁法设备、重力仪等设备，并了解了各种方法的原理和适用范围。
2. 数据处理数据处理是物探工作的重要环节。
我们使用了地震数据处理软件、电磁法数据处理软件等工具对采集到的数据进行处理。
在这个过程中，我学会了如何进行数据预处理、数据编辑、数据解释等工作，并了解了数据处理的重要性。
3. 数据解释数据解释是物探工作的最终目的。
我们根据处理后的数据，分析了地下地质结构、矿产资源等情况，并提出了相应的建议。
在这个过程中，我学会了如何分析数据、绘制地质图、建立地质模型等技能，并了解了物探成果的应用价值。
三、实习收获1. 知识与技能的提升通过这次实习，我掌握了物探基本原理、数据采集、数据处理和数据解释等方面的知识，提高了自己的物探技能。
2. 团队合作与沟通能力的培养在实习过程中，我与同学们一起完成了各项任务，学会了团队合作和沟通，提高了自己的组织协调能力。
3. 实践经验的积累这次实习让我有了实际操作的经验，对我今后从事物探工作具有很大的帮助。
四、实习总结通过这次物探实习，我对物探技术及其应用有了更深入的了解，收获颇丰。
我将以此为契机，继续努力学习物探知识，为将来从事物探工作打下坚实的基础。
2.3高精度磁法磁法勘探是勘探矿产资源的地球物理方法之一。
它是研究磁力在地表的分布与解决地质任务的方法,以在地磁场中岩石磁性的不同为基础,这些岩石的磁化程度彼此各异,而且造成磁异常。
在实际工作中,磁法勘探用来测量地表或一定高空中的磁场强度,它具有这样的详细程度和精确度,能够确定在普查对象影响下的磁场变化的规律性。
地磁场背景值大概为47000nT，周边无引起明显异常的岩石、矿山和剩磁物质，主要是电力电线、电信电缆等一些人为干扰。
我们在马面测区的6条测线均进行了高精度磁法测量。
我们在马面测区的6条测线均进行了高精度磁法测量,并运用以及自带的wcz-1质子磁力仪软件对测得数据进行日变校正。
绘制了磁异常剖面图，及此法异常平面等值线图2.4放射性勘探利用仪器测定土壤中的天然放射性，将放射性同位数作为示踪剂，用来测定地下水运动状态与有关参数或利用放射性同位数测量岩土密度与含水量的一种地球物理勘探方法。
测量地壳内放射性元素放出的射线强度，以发现放射性元素矿床，探明矿体大小，确定放射性元素含量，并指导铀、钍矿和钾盐矿的开采。
此法还用来寻找与放射性元素共生的其他非放射性矿床，如磷块岩矿，钼土矿、稀土和稀有元素矿床等；解决各类有关问题，如地质填图，寻找同油气、地下水等有关的构造，以及对环境污染的监测等。
放射性勘探的主要缺点是勘探深度不大，必须与其他地质工作配合进行。
此次实习我们在重点马面测区的3、4线进行了放射性勘探（主要是氡气异常探测）。
2.5地震勘探方法2.5.1地震干扰波野外进行地球物理勘探时，首先最好是能够找到工区附近已知的钻孔资料，如果没有钻孔资料必须进行干扰波调查，地震勘探的常见干扰波主要有声波速度在340m/s，工业电流干扰，频率在50-60Hz，机械振动干扰，了解规则和不规则干扰的存在情况，以便选择压制干扰的观测系统和设计采集数据的工作参数。
主要参数表现为视数度、波至时间、振幅和波形宽度的差异。
干扰剖面调查是一种小道间距、小偏移距、多个排列组成的剖面。
一、实习背景为了提高地质专业学生的实践能力和专业技能，我于2021年7月15日至8月15日参加了某物探队的野外实习。
本次实习主要在我国的某地开展，实习期间，我跟随物探队进行了野外数据采集、数据处理、成果解释等工作，对物探工作有了更加深入的了解。
二、实习目的1. 熟悉物探野外工作的基本流程和方法；2. 掌握物探数据采集、处理、解释等方面的技能；3. 提高团队合作能力和沟通能力；4. 增强对地质专业的认识和热爱。
三、实习内容1. 数据采集在实习期间，我参与了物探数据采集工作。
主要采用地震勘探方法，包括地震测线布设、激发、观测等环节。
在采集过程中，我学会了如何正确操作地震仪器，如地震检波器、地震炮等，并掌握了地震数据采集的基本原理。
2. 数据处理数据采集完成后，我参与了物探数据处理工作。
在数据处理过程中，我学习了地震数据预处理、地震成像、解释等方面的知识。
通过实际操作，我掌握了地震数据处理的基本流程和常用软件的使用方法。
3. 成果解释在成果解释阶段，我学习了如何根据物探数据解释地质构造、油气藏等。
在实习过程中，我参与了多个项目的成果解释工作，了解了物探成果在实际生产中的应用。
四、实习收获1. 理论与实践相结合通过本次实习，我将所学理论知识与实际工作相结合，加深了对物探工作的理解。
在野外工作中，我学会了如何将理论知识运用到实际操作中，提高了自己的实践能力。
2. 团队合作与沟通能力在实习过程中，我与团队成员密切配合，共同完成了各项任务。
这使我认识到团队合作的重要性，并提高了自己的沟通能力。
3. 对地质专业的认识通过本次实习，我对地质专业有了更加深入的了解，增强了对地质工作的热爱。
我认识到，地质工作是一项具有挑战性的工作，需要我们不断学习、探索，为我国地质事业贡献自己的力量。
五、实习体会1. 实践是检验真理的唯一标准在实习过程中，我深刻体会到实践的重要性。
只有将理论知识与实际工作相结合，才能不断提高自己的专业技能。
为了更好地了解物探技术在地质勘探中的应用，提高我们的实践操作能力，我们于2023年9月4日至9月10日在河南省登封市进行了为期一周的物探实习。
实习期间，我们通过实地操作、数据采集、成果分析等环节，对物探技术有了更深入的认识。
二、实习内容1. 实地操作实习期间，我们先后在登封市的太室山、少室山、少林水库等地进行了实地操作。
主要操作内容包括：（1）地震勘探：使用地震仪进行地震波采集，分析地层结构、构造和岩性特征。
（2）电磁勘探：利用电磁仪进行电磁波探测，获取地下介质信息。
（3）重力勘探：使用重力仪测量重力异常，分析地质构造和岩性。
2. 数据采集在实习过程中，我们严格按照操作规程进行数据采集，确保数据的准确性。
主要采集内容包括：（1）地震勘探：记录地震波传播时间、振幅等信息。
（2）电磁勘探：记录电磁场强度、相位等信息。
（3）重力勘探：记录重力加速度、梯度等信息。
3. 成果分析采集到的数据经过处理后，我们进行了以下分析：（1）地震勘探：根据地震波传播时间、振幅等信息，绘制地震剖面图，分析地层结构、构造和岩性特征。
（2）电磁勘探：根据电磁场强度、相位等信息，绘制电磁剖面图，分析地下介质信息。
（3）重力勘探：根据重力加速度、梯度等信息，绘制重力剖面图，分析地质构造和岩性。
1. 提高了实践操作能力：通过实地操作，我们掌握了物探仪器的使用方法，熟悉了物探技术的基本流程。
2. 增强了团队协作意识：实习过程中，我们分工合作，共同完成了各项任务，培养了团队协作精神。
3. 深化了对物探技术的认识：通过实习，我们对物探技术在地质勘探中的应用有了更深入的了解，为今后的工作奠定了基础。
4. 培养了严谨的工作态度：在实习过程中，我们严格按照操作规程进行数据采集和处理，培养了严谨的工作态度。
四、实习总结本次登封物探实习让我们受益匪浅，不仅提高了我们的实践操作能力，还加深了对物探技术的认识。
在今后的学习和工作中，我们将继续努力，将所学知识运用到实际工作中，为我国地质事业做出贡献。
物探,实习报告篇一：物探实习报告目录第一章实习目的和要求 (2)一实习目的 (2)二实习要求 (2)第二章工作内容和方法 (4)一工作内容 (4)二工作方法 (4)第三章收获及体会 (8)第四章问题和不足 (9)第五章致谢 (10)第一章实习目的和要求一实习目的作为一名地质学校的在校学生，我们不仅应该认真学习课本上的基础理论知识，掌握野外地质勘查工作的基本工作原理，更应该积极参加野外地质勘查工作，以便更加深入的学习和掌握野外地质勘查工作的工作方法和流程。
而这次顶岗实习不仅是对我们所学习的专业知识的加强和巩固，培养我们的实际动手能力和对野外施工的认识，培养我们实事求是、严谨认真的工作态度和勇于探索、不畏艰苦的工作作风。
而且，通过这次顶岗实习，有利于巩固我们在学校所学习的专业知识。
二实习要求这次顶岗实习中，我们主要的任务是学习在野外施工过程中需要注意的各种问题及工作方法。
学习野外施工过程中的基本工作方法及大概工作流程。
初步掌握地质野外勘察的具体工作方法和技术要求，能熟练的操作工作过程中所遇到的各种专业勘察仪器，切实掌握仪器的操作方法及保证仪器安全的主要措施。
在实习过程中，我们要服从公司的工作安排，积极向老员工学习工作过程中的各种注意事项以及野外施工过程中各个工作岗位上的工作方法及技术。
有问题，及时向老员工请教。
实习工作过程中，要将在学校中学到的各种专业知识付诸于实践当中去，养成严谨认真、不畏艰苦的工作作风。
争取早日将自己培养成一名合格的地质工作者。
第二章工作内容和方法一工作内容此次我们的实习单位是山东中煤物探测量总公司电法队，实习地点是新疆伊宁市霍城县，工作的主要目的任务包括（1）查明测区范围内小窑采空区的分布范围（2）尽量查明测区内大的断裂构造的位置、产状和煤层的埋藏深度、起伏状态（3）查明测区内火烧区的分布范围。
为了更好的完成本次工作的地质任务，我们的主要工作方法是：首先，利用GPS RTK放样进行野外地质测量放样；再在勘察范围内利用瞬变电磁法（PROTEM 67）进行野外地下数据采集；最后，通过对采集数据进行推导和反演，从而得到勘察结果。
一、实习单位及时间实习单位：XX石油物探公司实习时间：2023年7月1日至2023年8月31日二、实习目的通过本次实习，我将所学理论知识与实际工作相结合，深入了解石油物探行业的基本工作流程，掌握物探野外采集、数据处理、解释分析等方面的技能，为今后的工作打下坚实基础。
三、实习内容1. 物探野外采集在实习期间，我参与了地震数据采集工作。
具体内容包括：（1）学习地震仪器操作，了解地震数据采集原理。
（2）跟随经验丰富的工程师进行野外实地作业，掌握地震数据采集流程。
（3）熟悉地震数据处理设备，了解地震数据采集质量评价标准。
2. 数据处理在实习期间，我参与了地震数据处理工作。
具体内容包括：（1）学习地震数据处理软件，了解数据处理流程。
（2）参与地震数据预处理，包括去噪、增益、滤波等。
（3）进行地震数据反演，提取地质信息。
3. 解释分析在实习期间，我参与了地震解释分析工作。
具体内容包括：（1）学习地震解释分析方法，了解地质构造、储层等地质信息。
（2）参与地震剖面解释，绘制地质构造图、断层图等。
（3）对地震解释结果进行验证，提高解释准确性。
四、实习收获1. 提升了专业技能：通过本次实习，我对物探野外采集、数据处理、解释分析等方面的技能有了深入了解，为今后的工作打下了坚实基础。
2. 增强了团队协作能力：在实习过程中，我学会了与团队成员沟通协作，共同完成工作任务。
3. 拓宽了视野：实习让我对石油物探行业有了更全面的认识，为我今后的发展提供了更多可能性。
4. 培养了敬业精神：在实习期间，我深刻体会到石油物探行业的工作严谨性和责任感，为今后从事相关工作奠定了基础。
五、实习体会通过本次实习，我深刻认识到理论知识与实际工作相结合的重要性。
在今后的学习和工作中，我将努力将所学知识运用到实际工作中，不断提升自己的专业素养，为我国石油物探事业贡献自己的力量。
总结：本次实习让我受益匪浅，不仅提升了我的专业技能，还让我对石油物探行业有了更深入的了解。
物探单位实习报告一、实习目的与意义作为一名地质学校的学生，理论知识的学习和野外实践经验的积累同等重要。
实习是连接学校教育与实际工作的桥梁，它不仅能够巩固我们所学的专业知识，还能培养我们的实际动手能力和对野外施工的认识。
通过这次在物探单位的顶岗实习，我们有望培养出实事求是、严谨认真的工作态度和勇于探索、不畏艰苦的工作作风。
二、实习单位与实习内容本次实习单位是我国某知名物探公司，实习期间，我们主要参与了野外物探数据的采集、处理和解释工作。
实习内容涵盖了地震勘探、地质调查、地球物理勘探等多个方面，使我们全面了解了物探工作的整个流程。
三、实习过程与收获1. 实习过程在实习过程中，我们跟随经验丰富的工程师和技师，参与了野外数据采集、地震资料处理和解释等工作。
我们学会了使用各种物探仪器，掌握了数据采集的基本技巧，了解了地震资料处理的方法和流程，并学会了如何根据资料解释地质情况。
2. 实习收获（1）专业技能的提升通过实习，我们学会了使用各种物探仪器，掌握了数据采集、处理和解释的基本方法，使我们的专业技能得到了很大提升。
（2）团队协作能力的培养在实习过程中，我们学会了如何与团队成员有效沟通、协作，提高了我们的团队协作能力。
（3）严谨工作作风的养成实习过程中，我们深刻体会到严谨认真、一丝不苟的工作态度的重要性，使我们养成了严谨的工作作风。
（4）勇于探索的精神在实习过程中，我们面对未知的地质情况，学会了勇于探索、积极思考，培养了我们解决问题的能力。
四、实习总结通过这次实习，我们不仅巩固了所学的专业知识，还培养了实际动手能力和团队协作精神。
我们认识到，理论知识与实践经验是相辅相成的，只有将二者结合起来，才能更好地应对实际工作中的挑战。
在今后的工作中，我们将继续努力学习，不断提高自己，为我国物探事业的发展贡献自己的力量。
最后，感谢实习单位给予我们的关心与帮助，使我们能够在短时间内掌握物探工作的基本技能。
也感谢学校的组织与安排，让我们有机会参加这次实习，提升自己的专业素养。
物探的实习报告一、实习背景介绍物探是地球科学中的一项重要研究领域，通过使用各种物理探测方法，探测地球内部的物质性质和结构特征，对地质资源勘探、地下水资源评价、地震预测等方面具有重要意义。
本次实习是我作为一名地学专业的学生，在公司X参与的物探实习项目。
二、实习目标与任务本次实习的主要目标是学习和掌握物探方法的基本原理和实践操作技能，通过实地勘探、数据处理和分析，了解地下地质构造情况，并结合实际工程项目，制定合理的地质预测与勘探方案。
在实习的过程中，我主要承担了以下任务：1.参与实地勘探工作，包括地质地貌调查、岩芯采集与分析等。
2.学习和掌握地球物理仪器的使用方法，包括地震仪、电磁测深仪等。
3.参与数据处理和分析工作，包括地震波形解译、测深数据处理等。
4.结合实际工程项目，制定地质预测与勘探方案，提出优化建议。
三、实习过程与心得体会1. 实地勘探在实习的初期，我与同组的同学一起进行了实地勘探工作。
我们首先对勘探区域进行了地质地貌的调查，了解了区域的地形、岩性分布等信息。
随后，我们进行了岩芯的采集与分析工作，通过对采集到的岩芯样品进行观察与测试，了解了地下岩石的构造和特征。
2. 地球物理仪器的使用在实习的过程中，我们还学习并掌握了地球物理仪器的使用方法。
其中，地震仪是我们使用频率较高的一种仪器。
我们通过设置合理的参数，记录地震波形数据，并进行数据的处理和解译。
电磁测深仪是另一种常用的仪器，通过测量地下电磁场的变化，了解地下地质结构和含水层的情况。
3. 数据处理与分析在实习的过程中，我们还进行了大量的数据处理与分析工作。
通过使用专业软件，对采集到的地震波形数据进行解译与解析，了解地下地质构造和岩性分布情况。
同时，我们还对测深数据进行处理，计算出地下水位的位置与深度，为后续工程预测提供参考依据。
4. 地质预测与勘探方案最后阶段，我们结合实际工程项目，制定了地质预测与勘探方案。
根据前期的勘探与数据分析结果，我们提出了合理的地下构造预测，解释了地质问题可能带来的影响，并提出了优化建议，为后续工程设计与施工提供科学依据。
地球物理勘探中的分析与数据处理研究地球物理勘探是研究地球物理现象和规律及其应用的一门学科，广泛应用于矿产勘探、地埋物探测、水文地质勘探等领域。
而地球物理勘探中的分析与数据处理则是整个勘探过程中最为重要的环节之一。
本文将从地球物理勘探的基本概念入手，探讨分析与数据处理的方法和技术，以及在实际勘探中可能面临的问题和挑战。
一、地球物理勘探的基本概念地球物理勘探是用地球物理学的原理和方法探测地下结构和性质的科学技术，包括测量、分析、解释和应用等一系列过程。
常用的地球物理勘探方法有地震勘探、重力勘探、电磁法勘探、磁法勘探等。
地震勘探是利用声波在不同介质间传播速度差异，研究地下结构和岩石物性的勘探方法。
它通过激发地下震源，记录地面上的震动波形，然后对勘探区域进行三维反演，获得地下结构和物性信息。
重力勘探则是通过测量地球表面上物体引力产生的重力场变化，研究地壳或地球内部结构和密度分布的勘探方法。
它通过精密重力计测量地球表面重力加速度的分布，从而推算出地下的密度分布情况。
电磁法勘探则是利用地球表面上交变电场在地下的电性介质中传播和反射，研究地下介质电性结构的一种勘探方法。
它通过在地面放置线圈建立交变电场，记录地下回波信号，从而推算出地下介质的电性参数。
磁法勘探则是利用地球磁场和磁性物质的相互作用，研究地下磁性物质分布和地下构造的一种勘探方法。
它通过在地面上测量地磁场变化，记录地下磁性物质的异常值，从而推算出地下磁性物质的分布。
二、地球物理勘探中的数据处理技术地球物理勘探所涉及到的数据往往是海量的、复杂的、多维的。
数据的处理和分析是整个勘探过程中至关重要的环节。
地球物理勘探中常用的数据处理技术包括数据重采样、数据滤波、数据校正、数据勘误、反演方法等。
数据重采样是将输入数据用不同采样率的方式重新采样，可以加速数据处理和减小数据存储容量。
数据滤波则是通过将传感器记录到的原始数据进行卷积处理，去除干扰和噪声，提高数据的准确性。
本科生实验报告实验课程物探数据处理学院名称地球物理学院专业名称勘查技术与工程学生姓名00000学生学号*************指导教师李勇实验地点5417实验成绩二〇一七年九月二〇一七年十月物探数据处理实验报告实验一1.1实验目的本次实验的目的是将课本上给出的地质体情况，和计算参量的代数公式用c 语言通过程序表现出来。
在得知已知参量，例如形状大小、空间位置的情况下将其带入编写好的程序，从而得出地下板状体的物性值数据，并将所得数据其通过成图数据成图后在用反演原理与提前已知的异常体情况对比，判断所得结果的正误，并分析最终的实验结果。
1.2实验原理下图1-1为实验原理示意图。
图1-1 二维板模型利用如下公式进行c语言程序的实现：1.3实验参数中点横坐标x0(m)：1000 中点深度坐标z0(m)：1000 板高度(m)：400 板长度(m)：200 磁化强度：2000 二维板角度：45、90 磁化角度：90 测点k：100}1.4 实验结果图示图1-2 磁异常X分量图1-3 磁异常Z分量图1-4 重力异常1.5 实验结果分析根据实验数据结果图可知，通过观察x方向的磁分量可以看出：图中存在着明显的两个相互对称的正负异常，可知在地质体的上部积累了一定的负磁荷，而在地质体的下部积累了同样大小的正磁荷，实验中磁化角度取90度，二维板倾角取不同的角度。
通过观察z方向上的磁异常可知：地质体呈现出中部正异常的对称图形，总体呈现先增后减的趋势。
通过观察重力异常的图像可知：得出的结果基本与已知相符，所以通过重力异常、磁场x方向异常、磁场z方向异常数据曲线，就可对所求数据做相应的分析。
实验二2.1 实验目的通过给定的已知相关物性参数，建立截面为多边形的水平柱体模型，带入所编写好的c 语言程序中得到对应的正演结果。
主要考查了c 语言程序进编程，和对所得结果的理解，分析，判断能力。
2.2 实验原理图2-1为实验原理图图2-1 多边形截面水平柱体截面形状根据下列公式进行编程得到实验数据，作图分析。
；]1111[)1(2)1(2)11)(1(221212ln )1(2)1(2)11)(1(21{12x i zi tg x i z i tg x i x i z i z i z i x i z i x i x i x i z i x i z i x i x i x i z i z i z i x i z i x i z i z i n i G g ++---⋅-++-++-+-++++++⋅-++-++-+-+⋅=∑=∆σ；)(21Q M z P M x X +=∆π ；)(21P M z Q M x Z -=∆π式中：；z ix i z i x i x i x i z i z i z i z i x i z i tg x i z i tg x i x i z i z i x i x i z i z i n i Q 221212ln )1(2)1(2)1(221]1111[)1(2)1(2)1)(1({1++++⋅+-+-+-+⋅-++---⋅+-+-++--+=∑=；z ix i z i x i x i x i z i z i x i x i z i z i x i z i tg x i z i tg x i x i z i z i z i z i n i P 221212ln )1(2)1(2)1()1(21]1111[)1(2)1(2)1(2{1++++⋅+-+-++--+⋅+++---⋅+-+-+-+=∑=2.3 实验参数多边形半径r:200m ；磁化强度M:2000 拟合边数:10 磁化角度ct:90,120 测点k:100 中点深度z0(m):1000 中点横坐标x0(m):10002.3.1 讨论abtg 1-的值ab实际为由坐标原点对多边形第i 边的夹角，在具体计算中应当进行讨论： 当0.115->a 时，如果a b >0.0；当a>0.0,b>0.0,则θ=a btg 1-；当a<0.0，b<0.0，则θ=—π+abtg 1-；如果a b <0.0，则θ=π+ab tg 1-；当0.115-≤a 时，若b<0.0，则θ=—2π；若b>=0,则θ=2π2.4 实验结果图2-2 重力异常图2-3 磁异常X分量图2-4 磁异常Z分量2.5 实验结果分析通过重力异常图可知，通过观察x方向的磁分量可以看出：图中仍然存在着明显的两个相互对称的正负异常，故我们仍然可以认为地质体的上部积累了一定的负磁荷，而在地质体的下部积累了同样大小的正磁荷。
所以左端出现了较大的正异常，右端出现了较大的正异常。
通过观察z方向上的磁异常可知：地质体在中线（x=50）呈现出中部正异常的对称图形，总体呈现先增后减的趋势。
通过观察重力异常的图像可知：得出的结果基本与已知相符，出现了一个重力异常值曲线，最大值点对应的测线坐标与地质体实际位置一致。
实验三3.1 实验目的本次实验的目的是通过给定的地下异常体参量，对地下的正长方体进行正演计算，然后将计算参量的代数公式用c 语言通过程序表现出来。
在得知已知参量，例如形状大小、空间位置的情况下将其带入编写好的程序，并通过积分来实现正演。
从而得出地下板状体的物性值数据，并将所得数据其通过成图数据成图后在用反演原理与提前已知的异常体情况对比，判断所得结果的正误，并分析最终的实验结果。
将其运用的主要方法有通过划分单元来进行体积分，即8个体积元叠加即为立方体所对应的积分值，偶数为正，奇数为负，则可以计算出相应重磁异常。
3.2 实验原理利用如下公式进行编程：；2020202020201ln ln ({c z c z b y b y a x ax }zk)R(z yk)xk)(y (x tg zk)(z R]xk)[(x yk)(y R]yk)[(y xk)x G Δg -+-+-+----⋅-++--++---=σ；2020********yk)]}-(y ln[R M z zk)]-(z ln[R M y x k)R-(x zk)-yk)(z -(y tg 1-M x {-41c z cz b y b y a x a x Δx -+-+-+++++=π；2020********x k)]}-(x ln[R M z ]yk)R-(y zk)-x k)(z -(x tg 1-M y -zk)]-(z ln[R M x {41c z c z b y b y a x a x Δy -+-+-+++++=π；2020********}zk)R-(z yk)-x k)(y -(x tg 1-M z x k)]-(x ln[R M y -yk)]-(y ln[R M x {41c z c z b y b y a x a x Δz -+-+-+-++=π式中：])(2)(2)(2[2/1zk z yk y xk x R -+-+-=；3.3 实验参数列数:101 行数: 101 磁化强度:2000 磁倾角:45 磁偏角:45 立方体上坐标:（300,300,10） 立方体下坐标:（700,700,410） 密度:1 起终坐标:（0,1000）3.4 实验结果图示图3-1 重力异常平面图 图3-2 磁异常x 分量平面图3-3 磁异常y 分量平面 图3-4 磁异常z 分量平面图3.5实验结果分析根据实验数据结果图可知，通过观察x方向的磁分量可以看出：图中存在着两个对称的正负异常，且异常呈现出相似性形态幅值基本一致。
通过观察z方向上的磁异常可知：地质体呈现出中部正异常的对称图形，切呈现出对角线对称的形态上方为正下方为负的负异常。
通过观察重力异常的图像可知：得出的结果基本与已知相符，得到一个明显的类似同心圆的重力异常，边界较为清晰。
实验程序1#include<stdio.h>#include<math.h>#define PI 3.1415926#define G 6.67e-6void main(){FILE *fp1,*fp2,*fp3;fp1=fopen("g.txt","w");fp2=fopen("deltax.txt","w");fp3=fopen("deltaz.txt","w");double r1,r2,r3,r4; //定义板截面角点A,B,C,D到P点之间的距离double f1,f2,f3,f4; //定义r1.r2.r3.r4与X轴正向夹角int k,i; //定义P点的坐标xk,zk;double xk,zk;double cgm=2.67e+3;double x0=1000.0; //中点横坐标double z0=1000.0;double l=400.0; //板高度double b=200.0; //板长度double I,M=2000.0;double z1,z2,g,deltax,deltaz;double J=90.0*PI/180.0; //计算二维板角度I=90.0*PI/180.0; //计算磁化角度double a[100]={0.0};double b1[100]={0.0};double ct[100]={0.0};for(k=0;k<101;k++){xk=k*20.0;zk=0.0;b1[0]=z0-zk-l*sin(J);a[0]=xk-x0+b+l*cos(J);b1[1]=z0-zk+l*sin(J);a[1]=xk-x0+b-l*cos(J);b1[2]=z0-zk-l*sin(J);a[2]=xk-x0-b+l*cos(J);b1[3]=z0-zk+l*sin(J);a[3]=xk-x0-b-l*cos(J);r1=sqrt(pow(b1[0],2.0)+pow(a[0],2.0));r2=sqrt(pow(b1[1],2.0)+pow(a[1],2.0));r3=sqrt(pow(b1[2],2.0)+pow(a[2],2.0));r4=sqrt(pow(b1[3],2.0)+pow(a[3],2.0));for(i=0;i<4;i++){if(fabs(a[i])>1.0e-15){if((b1[i]/a[i])>0.0){if(a[i]>0.0&&b1[i]>0.0){ct[i]=atan(b1[i]/a[i]);}else if(a[i]<0.0&&b1[i]<0.0){ct[i]=atan(b1[i]/a[i])-PI;}}else if((b1[i]/a[i])<0.0){ct[i]=PI+atan(b1[i]/a[i]);}}else if(fabs(a[i])<=1.0e-15){if(b1[i]<0.0){ct[i]=-PI/2.0;}else{ct[i]=PI/2.0;}}f1=PI-ct[0];f2=PI-ct[1];f3=PI-ct[2];f4=PI-ct[3];z1=z0-l*sin(J);z2=z1+l*sin(J);g=2.0*G*cgm*((z2*(f2-f4)-z1*(f1-f3))+xk*(sin(J)*sin(J)*log((r2*r3)/(r 1*r4))+cos(J)*sin(J)*(f1-f2-f3+f4))+2.0*b*(sin(J)*sin(J)*log(r4/r3)+c os(J)*sin(J)*(f3-f4)));deltax=((M*sin(J))/(2.0*PI))*(log((r2*r3)/(r1*r4))*cos(J-I)-(sin(J-I) )*(f1-f2-f3+f4));deltaz=((M*sin(J))/(2.0*PI))*((sin(J-I))*log((r2*r3)/(r1*r4))+(cos(J-I))*(f1-f2-f3+f4));fprintf(fp1,"%f\n",g); //将所得数据写入文件fprintf(fp2,"%f\n",deltax);fprintf(fp3,"%f\n",deltaz);}}fclose(fp1);fclose(fp2);fclose(fp3);实验程序2include"stdio.h"#include"math.h"#define PI 3.1415926 //圆周率#define G 6.67e-11 // 万有引力常数#define r 200.0 //多边形半径#define M 2000.0 //磁化强度#define thegama 2.67e+3void main(){double ct11(double a,double b); //声明子函数FILE *fp1,*fp2,*fp3;fp1=fopen("deltag.txt","w");fp2=fopen("deltax1.txt","w");fp3=fopen("deltaz1.txt","w");double x0=1000.0;double z0=1000.0;int n=10,N=100;double R=200.0;double M=2000.0;double cgm=2.67e+3;double ct=PI/4.0;double Mx=M*cos(ct);double Mz=M*sin(ct);double x[10000]={0.0},x1[1000]={0.0};double z[10000]={0.0},z1[1000]={0.0};double xk[10000]={0.0},zk[1000]={0.0};double sum1=0.0;double Q=0.0;double P=0.0;doublea1,a2,a3,a4,a5,a6,a7,aa,bb,ct1,ct2,deltag,deltax1,deltaz1,a,b;int i,k;for(i=0;i<n+1;i++)/////////////////////////∠坐标{x1[i]=x0+R*cos(i*2*PI/n);z1[i]=z0-R*sin(i*2*PI/n);}for(k=0;k<N;k++){xk[k]=20.0*k;///////////////////平面坐标zk[k]=0.0;}for(k=0;k<N;k++){sum1=0.0;Q=0.0;P=0.0;for(i=0;i<n+1;i++){x[i]=xk[k]-x1[i];z[i]=zk[k]-z1[i];}for(i=0;i<n;i++){a1=(z[i+1]-z[i])*(x[i]*z[i+1]-x[i+1]*z[i])/(pow(z[i+1]-z[i],2.0)+ pow(x[i+1]-x[i],2.0));a2=(x[i+1]-x[i])*(x[i]*z[i+1]-x[i+1]*z[i])/(pow(z[i+1]-z[i],2.0)+pow( x[i+1]-x[i],2.0));a3=(z[i+1]-z[i])*(x[i]-x[i+1])/(pow(z[i+1]-z[i],2.0)+pow(x[i]-x[i +1],2.0));a4=pow(z[i+1]-z[i],2.0)/(pow(z[i+1]-z[i],2.0)+pow(x[i]-x[i+1],2.0));aa=pow(x[i+1],2.0)+pow(z[i+1],2.0);bb=pow(x[i],2.0)+pow(z[i],2.0);ct1=ct11(z[i],x[i]); //调用子函数ct2=ct11(z[i+1],x[i+1]);a5=(a1/2.0)*log(aa/bb)+a2*(atan(ct1)-atan(ct2));sum1+=a5;a6=a3*(atan(ct1)-atan(ct2))-(a4/2.0)*log(aa/bb);Q+=a6;a7=a4*(atan(ct1)-atan(ct2))+(a3/2.0)*log(aa/bb);P+=a7;}deltag=2.0*G*cgm*sum1;deltax1=(Mx*P+Mz*Q)/(2.0*PI);deltaz1=(Mx*Q-Mz*P)/(2.0*PI);fprintf(fp1,"%lf\n",deltag);fprintf(fp2,"%lf\n",deltax1);fprintf(fp3,"%lf\n",deltaz1);}fclose(fp1);fclose(fp2);fclose(fp3);}///////////////////////////////////////////double ct11(double a,double b) //判断子函数{double ct;if(fabs(a)>1.0e-15){if((b/a)>0.0){if(a>0.0&&b>0.0)ct=atan(b/a);elsect=atan(b/a)-PI;}if((b/a)<0.0){ct=PI+atan(b/a);}}else{if(b<0.0){ct=-PI/2.0;}else{ct=PI/2.0;}}return(ct);}实验程序3#include <stdlib.h>#include <stdio.h>#include <math.h>#define l 101 // 列数#define h 101 // 行数#define M 2000.0 //磁化强度#define II 60 // 磁倾角#define aa 60 //磁偏角#define PI 3.1415926 //圆周率#define G 6.674E-11 // 万有引力常数#define SI 1e6 // 将m/s^2换算为g.u.的比例因子#define DENSITY 1E3 // 将g/cm^3换算为kg/m^3的比例因子//计算atan值double batan(double m,double n,double x,double r);//Kernel函数,用途：为Cuboid函数准备double Kernel(double dXk, double dYk, double dZk,//dXk, dYk, dZk: 观测点坐标double dXs, double dYs, double dZs);//dYs, dZs, dZs: 立方体坐标，s取1,2 //Cuboid函数,用途: 计算立方体重力异常,参数说明double Cuboid(double dXk, double dYk, double dZk,//dXk, dYk, dZk: 观测点坐标double dX1, double dY1, double dZ1,// dX1, dY1, dZ1: 立方体坐标1double dX2, double dY2, double dZ2,//dX2, dY2, dZ2: 立方体坐标2double dDens); //dDens: 立方体密度// Cix函数double Cix(double dXk, double dYk, double dZk,double dXs, double dYs, double dZs,double Mx, double My, double Mz);// Cixxx函数: 用于计算立方体△x异常double Cixxx(double dXk, double dYk, double dZk,double dX1, double dY1, double dZ1,double dX2, double dY2, double dZ2,double Mx, double My, double Mz );// Ciy函数double Ciy(double dXk, double dYk, double dZk,double dXs, double dYs, double dZs,double Mx, double My, double Mz);// Ciyyy函数: 用于计算立方体△Y磁异常double Ciyyy(double dXk, double dYk, double dZk,double dX1, double dY1, double dZ1,double dX2, double dY2, double dZ2,double Mx, double My, double Mz );// Ciz函数double Ciz(double dXk, double dYk, double dZk,double dXs, double dYs, double dZs,double Mx, double My, double Mz);// Cizzz函数: 用于计算立方体△Z磁异常double Cizzz(double dXk, double dYk, double dZk,double dX1, double dY1, double dZ1,double dX2, double dY2, double dZ2,double Mx, double My, double Mz );void main(){FILE *f1,*f2,*f3,*f4; //定义文件指针f1 = fopen("重力异常.txt", "w"); //打开文件"重力异常.txt"f2 = fopen("磁异常X分量.txt", "w"); //打开文件"磁异常X分量.txt"f3 = fopen("磁异常Y分量.txt", "w"); //打开文件"磁异常Y分量.txt" f4 = fopen("磁异常Z分量.txt", "w"); //打开文件"磁异常Z分量.txt"int i, j;// 网格参数double dMinX = 0.0, dMaxX = 1000.0; // x坐标范围double dMinY = 0.0, dMaxY = 1000.0; // y坐标范围double dDX = (dMaxX - dMinX) / (l-1); // x方向上的网格间距double dDY = (dMaxY - dMinY) / (h-1); // y方向上的网格间距double px[101]={0.0},py[101]={0.0}; // x,y坐标构成的向量double pz[101][101]; // z坐标构成的矩阵double dg[101][101]={0.0}; // 定义重力异常数组double dx[101][101]={0.0},dy[101][101]={0.0},dz[101][101]={0.0}; //定义磁异常X分量、磁异常Y分量以及磁异常Z分量数组double I=II*PI/180; //角度转换double a=aa*PI/180; //角度转换// 指定坐标范围for (i = 0; i < h; i++){py[i] = dMinY + i*dDY;for (j = 0; j < l; j++){pz[i][j] = 0.0;}}for(i=0; i<l; i++){px[i] = dMinX + i*dDX;}// 定义立方体参数double dX1 = 300, dY1 = 300, dZ1 = 10; // 坐标1double dX2 = 700, dY2 = 700, dZ2 = 410; // 坐标2double dDens = 1.0;// 密度double Mx = M*cos(I)*cos(a);double My = M*cos(I)*sin(a);double Mz = M*sin(I); //磁分量// 生成重力与磁分量异常for (i = 0; i < h; i++){for (j = 0; j < l; j++){dg[i][j] = Cuboid(px[i], py[j], pz[i][j],dX1, dY1, dZ1,dX2, dY2, dZ2, dDens);dx[i][j] = Cixxx(px[i], py[j], pz[i][j],dX1, dY1, dZ1,dX2, dY2, dZ2, Mx, My, Mz);dy[i][j] = Ciyyy(px[i], py[j], pz[i][j],dX1, dY1, dZ1,dX2, dY2, dZ2, Mx, My, Mz);dz[i][j] = Cizzz(px[i], py[j], pz[i][j],dX1, dY1, dZ1,dX2, dY2, dZ2, Mx, My, Mz);fprintf(f1,"%d\t%d\t%f\n",i*10,j*10,dg[i][j]);//将所得到的重力异常数据写入文件fprintf(f2,"%d\t%d\t%f\n",i*10,j*10,dx[i][j]);//将所得到的磁异常X分量数据写入文件 fprintf(f3,"%d\t%d\t%f\n",i*10,j*10,dy[i][j]);//将所得到的磁异常Y分量数据写入文件 fprintf(f4,"%d\t%d\t%f\n",i*10,j*10,dz[i][j]);//将所得到的磁异常Z分量数据写入文件}}}//arctan的子程序double batan(double m,double n,double x,double r){double q;if(fabs(x*r)>1.0e-15){if((m*n/(x*r))>0.0){if((m*n)>0.0&&(x*r)>0.0) q=atan(m*n/(x*r));else q=-PI+atan(m*n/(x*r));}else q=PI+atan(m*n/(x*r));}else{if((m*n)<0.0) q=-PI/2;else q=PI/2;}return q;}// Kernel函数double Kernel(double dXk, double dYk, double dZk,double dXs, double dYs, double dZs){double dX=dXs-dXk, dY=dYs-dYk, dZ=dZs-dZk;double dR=sqrt(dX*dX+dY*dY+dZ*dZ);if(dZ==0){return -dX*log(dR+dY)-dY*log(dR+dX);}else{return -dX*log(dR+dY)-dY*log(dR+dX)+dZ*atan(dX*dY/(dZ*dR));}}// Cuboid函数: 用于计算立方体重力异常double Cuboid(double dXk, double dYk, double dZk,double dX1, double dY1, double dZ1,double dX2, double dY2, double dZ2,double dDens){return SI*DENSITY*G*dDens*(-Kernel(dXk, dYk , dZk, dX1, dY1, dZ1) //SI取数据为毫伽，见常数定义+Kernel(dXk, dYk , dZk, dX1, dY1, dZ2) //8个体积叠加即为立方体所对应的积分值+Kernel(dXk, dYk , dZk, dX1, dY2, dZ1) //偶数为正，奇数为负-Kernel(dXk, dYk , dZk, dX1, dY2, dZ2)+Kernel(dXk, dYk , dZk, dX2, dY1, dZ1)-Kernel(dXk, dYk , dZk, dX2, dY1, dZ2)-Kernel(dXk, dYk , dZk, dX2, dY2, dZ1)+Kernel(dXk, dYk , dZk, dX2, dY2, dZ2));}// Cix函数double Cix(double dXk, double dYk, double dZk,double dXs, double dYs, double dZs,double Mx, double My, double Mz){double dX=dXs-dXk, dY=dYs-dYk, dZ=dZs-dZk;double dR=sqrt(dX*dX+dY*dY+dZ*dZ);return My*log(dR+dZ)+Mz*log(dR+dY)-Mx*batan(dY,dZ,dX,dR);}// Cixxx函数: 用于计算立方体△X磁异常double Cixxx(double dXk, double dYk, double dZk,double dX1, double dY1, double dZ1,double dX2, double dY2, double dZ2,double Mx, double My, double Mz ){return 1.0/(4*PI)*(-Cix(dXk, dYk , dZk, dX1, dY1, dZ1, Mx, My,Mz) +Cix(dXk, dYk , dZk, dX1, dY1, dZ2, Mx, My,Mz)+Cix(dXk, dYk , dZk, dX1, dY2, dZ1, Mx, My,Mz) -Cix(dXk, dYk , dZk, dX1, dY2, dZ2, Mx, My,Mz)+Cix(dXk, dYk , dZk, dX2, dY1, dZ1, Mx, My,Mz)-Cix(dXk, dYk , dZk, dX2, dY1, dZ2, Mx, My,Mz)-Cix(dXk, dYk , dZk, dX2, dY2, dZ1, Mx, My,Mz)+Cix(dXk, dYk , dZk, dX2, dY2, dZ2, Mx, My,Mz));}// Ciy函数double Ciy(double dXk, double dYk, double dZk,double dXs, double dYs, double dZs,double Mx, double My, double Mz){double dX=dXs-dXk, dY=dYs-dYk, dZ=dZs-dZk;double dR=sqrt(dX*dX+dY*dY+dZ*dZ);return Mx*log(dR+dZ)+Mz*log(dR+dX)-My*batan(dX,dZ,dY,dR);}// Ciyyy函数: 用于计算立方体△Y磁异常double Ciyyy(double dXk, double dYk, double dZk,double dX1, double dY1, double dZ1,double dX2, double dY2, double dZ2,double Mx, double My, double Mz ){return 1.0/(4*PI)*(-Ciy(dXk, dYk , dZk, dX1, dY1, dZ1, Mx, My,Mz) +Ciy(dXk, dYk , dZk, dX1, dY1, dZ2, Mx, My,Mz)+Ciy(dXk, dYk , dZk, dX1, dY2, dZ1, Mx, My,Mz) -Ciy(dXk, dYk , dZk, dX1, dY2, dZ2, Mx, My,Mz)+Ciy(dXk, dYk , dZk, dX2, dY1, dZ1, Mx, My,Mz)-Ciy(dXk, dYk , dZk, dX2, dY1, dZ2, Mx, My,Mz)-Ciy(dXk, dYk , dZk, dX2, dY2, dZ1, Mx, My,Mz)+Ciy(dXk, dYk , dZk, dX2, dY2, dZ2, Mx, My,Mz));}// Ciz函数double Ciz(double dXk, double dYk, double dZk,double dXs, double dYs, double dZs,double Mx, double My, double Mz){double dX=dXs-dXk, dY=dYs-dYk, dZ=dZs-dZk;double dR=sqrt(dX*dX+dY*dY+dZ*dZ);return Mx*log(dR+dY)-My*log(dR+dX)-Mz*batan(dX,dY,dZ,dR);}// Cizzz函数: 用于计算立方体△Z磁异常double Cizzz(double dXk, double dYk, double dZk,double dX1, double dY1, double dZ1,double dX2, double dY2, double dZ2,double Mx, double My, double Mz ){return 1.0/(4*PI)*(-Ciz(dXk, dYk , dZk, dX1, dY1, dZ1, Mx, My,Mz) +Ciz(dXk, dYk , dZk, dX1, dY1, dZ2, Mx, My,Mz)+Ciz(dXk, dYk , dZk, dX1, dY2, dZ1, Mx, My,Mz) -Ciz(dXk, dYk , dZk, dX1, dY2, dZ2, Mx, My,Mz)+Ciz(dXk, dYk , dZk, dX2, dY1, dZ1, Mx, My,Mz)-Ciz(dXk, dYk , dZk, dX2, dY1, dZ2, Mx, My,Mz)-Ciz(dXk, dYk , dZk, dX2, dY2, dZ1, Mx, My,Mz)+Ciz(dXk, dYk , dZk, dX2, dY2, dZ2, Mx, My,Mz}。
