物探数据处理实验报告
- 格式:docx
- 大小:333.31 KB
- 文档页数:26
物探认识实习实习报告
实习时间,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. 地质预测与勘探方案最后阶段,我们结合实际工程项目,制定了地质预测与勘探方案。
根据前期的勘探与数据分析结果,我们提出了合理的地下构造预测,解释了地质问题可能带来的影响,并提出了优化建议,为后续工程设计与施工提供科学依据。
地球物理勘探中的分析与数据处理研究地球物理勘探是研究地球物理现象和规律及其应用的一门学科,广泛应用于矿产勘探、地埋物探测、水文地质勘探等领域。
而地球物理勘探中的分析与数据处理则是整个勘探过程中最为重要的环节之一。
本文将从地球物理勘探的基本概念入手,探讨分析与数据处理的方法和技术,以及在实际勘探中可能面临的问题和挑战。
一、地球物理勘探的基本概念地球物理勘探是用地球物理学的原理和方法探测地下结构和性质的科学技术,包括测量、分析、解释和应用等一系列过程。
常用的地球物理勘探方法有地震勘探、重力勘探、电磁法勘探、磁法勘探等。
地震勘探是利用声波在不同介质间传播速度差异,研究地下结构和岩石物性的勘探方法。
它通过激发地下震源,记录地面上的震动波形,然后对勘探区域进行三维反演,获得地下结构和物性信息。
重力勘探则是通过测量地球表面上物体引力产生的重力场变化,研究地壳或地球内部结构和密度分布的勘探方法。
它通过精密重力计测量地球表面重力加速度的分布,从而推算出地下的密度分布情况。
电磁法勘探则是利用地球表面上交变电场在地下的电性介质中传播和反射,研究地下介质电性结构的一种勘探方法。
它通过在地面放置线圈建立交变电场,记录地下回波信号,从而推算出地下介质的电性参数。
磁法勘探则是利用地球磁场和磁性物质的相互作用,研究地下磁性物质分布和地下构造的一种勘探方法。
它通过在地面上测量地磁场变化,记录地下磁性物质的异常值,从而推算出地下磁性物质的分布。
二、地球物理勘探中的数据处理技术地球物理勘探所涉及到的数据往往是海量的、复杂的、多维的。
数据的处理和分析是整个勘探过程中至关重要的环节。
地球物理勘探中常用的数据处理技术包括数据重采样、数据滤波、数据校正、数据勘误、反演方法等。
数据重采样是将输入数据用不同采样率的方式重新采样,可以加速数据处理和减小数据存储容量。
数据滤波则是通过将传感器记录到的原始数据进行卷积处理,去除干扰和噪声,提高数据的准确性。
物探综合实习报告一、实习目的和意义本次物探综合实习是地球物理专业重要的教学内容和环节,是理论联系实际、强化技能和能力培养的重要途径。
通过对野外物探工作的综合考察和调查,能够获得典型的物探相关知识,拓展视野,陶冶情操,加深对课堂教学内容的理解,弥补课堂教学的不足,并初步掌握物探野外考察、调查方法和步骤。
二、实习时间、地点和线路2.1实习时间:XX年XX月至XX年XX日2.2实习地点:XX省XX市XX物探队2.3考察路线:野外数据采集——数据处理——成果解释三、实习区域概况3.1地理位置:实习区域位于XX省XX市,地处东经XX度XX分至XX度XX分,北纬XX度XX分至XX度XX分。
3.2地质条件:实习区域地质构造复杂,涵盖多个地质年代,岩性多样,具有良好的物性条件,是进行物探实习的理想场所。
四、实习内容与过程4.1野外数据采集野外数据采集是物探实习的基础工作,主要包括地震数据采集、重力数据采集、磁法数据采集等。
在野外数据采集过程中,我们学会了使用地震仪、重力仪、磁法仪等设备,掌握了数据采集的基本方法和技术要求。
4.2数据处理数据处理是物探实习的重要环节,主要包括数据整理、预处理、正演、反演等。
在数据处理过程中,我们学习了数据处理软件的使用,掌握了数据处理的基本流程和方法。
4.3成果解释成果解释是物探实习的关键环节,主要包括地震解释、重力解释、磁法解释等。
在成果解释过程中,我们学会了如何根据物性特征和地质条件,解释地震、重力、磁法等数据,得出地质构造和矿产资源分布规律。
五、实习收获与体会5.1专业技能提升通过本次实习,我们掌握了物探野外数据采集、数据处理和成果解释等基本方法和技能,提升了专业素质和能力。
5.2团队合作与沟通在实习过程中,我们学会了团队合作和沟通,提高了团队协作能力和沟通技巧。
5.3对物探工作的认识通过实习,我们深刻认识到物探工作的重要性和艰辛性,增强了从事物探工作的信心和决心。
六、实习总结本次物探综合实习使我们受益匪浅,不仅提升了专业技能,还培养了团队合作和沟通能力。
勘查地球物理学实验报告姓名:班级:学号:_______________指导老师:磁法勘探数据采集和资料处理一、实验目的1. 学习磁法勘探数据采集工作方法;2. 了解数据处理的基本流程。
二、磁法勘探的原理磁法勘探是以岩矿石间的磁性差异为基础,通过接收和研究地质体(构造或矿体等)在地表及其周围空间产生的地球磁场的变化和特征来推断地质体存在状态(产状、埋深、规模等)的一种物探方法。
三、实验内容及步骤(一)实验内容本实验在室外使用高精度磁力仪做剖面观测,学习磁法勘探的野外工作过程和仪器操作,对观测的数据进行整理,编写实验报告。
(二)仪器G856质子磁力仪,探头及相关的仪器配件。
(三)实验的主要步骤(1)布置测线、测点;(2)将磁力仪与探头连接;(3)测线测量时通常2-3 人一组,由一人拿探头,一人兼做记录,或单独由一人记录;(4)打开仪器,设置日期和时间、设置线号、设置点号和调谐场等参数,具体操作见仪器使用说明书;(5)逐个测点进行磁场观测,并记录观测值,完成剖面上所有测点的观测;(6)对观测的数据进行整理,绘制磁场变化的剖面图,分析剖面上的磁场变化特征,编写实验报告。
(四)磁法观测的注意事项(1)质子磁力仪操作员应按常规磁法仪器操作规范要求,在进行高精度磁测,在地质勘探方面使用,则需用一台仪器在测区附近做日变观测;(2)操作人员身上一定要清理掉一切铁磁性物体,观测时罗盘应远离仪器5m 以上,要远离汽车30m 以上,其它人员也要与仪器保持距离;(3)对于高精度测量的测区应是电磁干扰小,磁场梯度较小,遇到强磁性干扰时(铁路、高压线等),须合理移动点位,并记录在案;(4)要注意仪器安全,防止碰撞。
四、实验成果及分析(一)实验成果图图1 磁场变化剖面图图2 测区磁场变化平剖图图3 磁场等值线图(二)实验成果分析(1)从剖面图中可以看出磁场变化有一个最高点和最底点,在第6测点处磁场强度最大为47669.0,而在第19个测点处达到最低值为47192.0。
物探认识实习实习报告
实习时间,2021年7月1日-2021年8月1日。
实习地点,某地质勘探公司。
实习内容:
在本次实习中,我主要负责参与物探勘探项目的实地调研和数
据采集工作。
在实地调研中,我学习了如何使用地质仪器进行地质
勘探,掌握了地质勘探的基本原理和方法。
同时,我还参与了地质
勘探现场的勘探工作,学习了如何进行地质勘探的现场操作,包括
地质样品的采集、标注和保存等工作。
在数据采集方面,我学习了如何使用地质勘探软件进行数据采
集和处理,掌握了地质勘探数据的处理方法和技巧。
通过实际操作,我深入了解了地质勘探数据的重要性和用途,对地质勘探工作有了
更深入的认识。
在实习期间,我还参与了公司举办的地质勘探技术培训,学习
了地质勘探技术的最新发展动态和未来发展趋势,对地质勘探行业
的发展有了更深入的了解。
实习收获:
通过本次实习,我对地质勘探行业有了更深入的了解,掌握了
地质勘探的基本理论和技术,提高了自己的实际操作能力和数据处
理能力。
同时,我还学到了与同事合作的重要性和团队合作的技巧,提高了自己的沟通能力和团队协作能力。
总结:
本次实习让我对地质勘探行业有了更深入的了解,提高了自己
的实际操作能力和数据处理能力,对未来的职业发展有了更清晰的
规划。
在未来的工作中,我将继续努力学习,不断提高自己的专业
能力,为地质勘探行业的发展贡献自己的力量。
目录1 实习目的及要求31.1 实习目的31.2 实习要求31.3实习容及安排32 直流电阻率法42.1 直接反演法42.1.1 根本原理42.1.2 数据处理42.1.3 实测数据分析解释62.2 最小二乘法反演62.2.1 根本原理62.2.2 数据处理72.2.3 实测数据分析解释92.3 对称四极测深92.3.1 根本原理92.3.2 数据处理92.3.3 实测数据分析解释102.4 高密度电阻率法102.4.1 根本原理102.4.2 数据处理102.4.3 实测数据分析解释123 电磁法133.1电磁测深Bostic反演133.1.1 根本原理133.1.2 数据处理133.1.3 实测数据分析解释163.2 EH4电磁测深163.2.1 根本原理163.2.2 数据处理163.2.3 实测数据分析解释174 地震勘探184.1 折射波法184.1.1 根本原理184.1.2 数据处理184.1.3 实测数据分析解释214.2 反射波屡次覆盖方法214.2.1 根本原理214.2.2 数据处理214.2.3 实测数据分析解释244.3 面波法244.3.1根本理论244.3.2 数据处理244.3.3 面波地质绘图264.3.4 实测数据分析解释274.4 地震映像法274.4.1 根本原理274.4.2 数据处理274.4.3 实测数据分析解释285 磁法勘探285.1 根本原理295.2 数据处理295.3 实测数据分析解释296 实习总结错误!未定义书签。
6.1物探数据处理的认识错误!未定义书签。
1 实习目的及要求1.1实习目的我们这次的数据处理实习,是为了让我们更加了解物探的工作要求,了解怎么处理数据过程,更加懂得在以后的工作中怎么样处理实测数据,通过每次的数据处理,更加的了解每个方法的运用问题。
1.2 实习要求要求我们每个人都要自己好好的去处理,这样更加懂得怎么样去了解实测数据的修正。
工程与环境物探实验报告(电法实验)班级:空间信息与数字技术二班一,采集系统:1,多功能数字直流激电仪、多路电极转换器、干电池、Res2dinv数据处理软件、BTRC2004接收与格式转换软件2,采用三电位电极系测量装置。
二,处理过程:I,数据采集:1,安装装置:将功能数字直流激电仪,多路电极转换仪连起来,然后连上干电池,选择36V电压。
然后连接的60个电极如下图然后检查装置是否安装妥当。
2,参数设置:以及设置测量的排列类型,分别为α,β,γ排列;3,仪器进入测量状态,在测量过程中仪器自动存储视电阻率值,测量完毕后显示数据总数。
4,然后复位,继续重新设置参数及测量的类型,到测完γ为止。
5,导出数据,进行相应处理。
II,数据处理1,准备相应软件:在电脑上安装好BTRC2004接收与格式转换软件和Res2dinv软件。
2,使用BTRC2004转换数据格式,打开软件导入要处理的数据,然后将负值改为正值然后选择转换格式再保存,则数据转换完成。
重复上述操作,分别将相应实验α,β,γ的值做相应处理。
3,将转换后的数据再用Res2dinv做进一步处理。
打开软件后,选择要处理的数据,再进行反演操作,如下所示,选择最小二乘反演4,接着就会出现相应的反演成果根据使用要求更改数据,以调整误差到要求的范围内,操作如下点击要删除的点,然后保存,。
再重复上述反演操作,直到误差符合规定为止。
5,重复第3步,将每个实验所转换的数据的α,β,γ都做反演变换,保存全部图片,用以后续分析。
三,实验结果及分析:实验1:(高阻大球)α排列结果图像分析:在三电位电极系测量装置42cm到 66cm处下方有一高阻带,电阻率约为53.6~ 67.9Ω·m,高阻体顶部深埋3.75cm处,底部深埋11.5cm处。
在α排列方式下,对于高低阻体分界面,视电阻率形成向低阻方向倾斜的分界面,而反演得到的电阻率断面中,由高阻向低阻的方向向两侧逐渐过渡。
一、实习背景随着我国经济的快速发展,地球物理勘探技术在石油、天然气、矿产、水文地质等领域得到了广泛应用。
为了更好地了解和掌握地球物理勘探的基本原理、技术方法和实际应用,我于2023年在XX地质勘探公司进行了为期一个月的实习。
二、实习目的1. 熟悉地球物理勘探的基本原理、技术方法和实际应用;2. 提高野外实际操作能力,掌握勘探数据的采集、处理和分析;3. 培养团队合作精神,提高沟通协调能力。
三、实习内容1. 地震勘探:学习了地震勘探的基本原理,掌握了地震仪器的操作方法,参与了地震数据的采集、处理和分析工作。
2. 重力勘探:了解了重力勘探的基本原理,掌握了重力仪器的操作方法,参与了重力数据的采集、处理和分析工作。
3. 地磁勘探:学习了地磁勘探的基本原理,掌握了地磁仪器的操作方法,参与了地磁数据的采集、处理和分析工作。
4. 电磁勘探:了解了电磁勘探的基本原理,掌握了电磁仪器的操作方法,参与了电磁数据的采集、处理和分析工作。
5. 实验室工作:参与了勘探数据的预处理、处理和分析工作,学习了数据处理软件的使用。
四、实习收获1. 理论知识与实践相结合:通过实习,我将所学理论知识与实际工作相结合,提高了自己的实际操作能力。
2. 培养了团队合作精神:在实习过程中,我学会了与团队成员沟通、协调,共同完成勘探任务。
3. 拓宽了知识面:实习使我了解了地球物理勘探领域的最新技术和发展趋势,拓宽了自己的知识面。
4. 提高了综合素质:实习过程中,我学会了如何处理突发事件,提高了自己的应变能力和心理素质。
五、实习体会1. 实践是检验真理的唯一标准:通过实习,我深刻体会到理论知识与实践相结合的重要性。
2. 团队合作精神至关重要:在勘探工作中,团队成员之间的协作至关重要,只有团结一致,才能完成各项任务。
3. 严谨的工作态度:地球物理勘探工作要求严谨细致,只有认真对待每一个环节,才能确保数据的准确性。
4. 持续学习:地球物理勘探领域不断发展,我们要不断学习新知识、新技术,以适应行业发展的需求。
物探实习报告一、前言我于2023进行了为期两周的物探实习,通过这次实习,我对物探技术有了更深入的了解,同时也锻炼了自己的实践能力。
在实习过程中,我参与了地震资料的采集、处理和解释等工作,现就实习过程中的所见所闻和收获作如下报告。
二、实习内容1. 实习单位简介我实习的单位是XXX物探公司,成立于1990年,是一家专业从事地球物理勘探的公司。
公司主要业务包括地震勘探、测井、钻井等,服务于石油、天然气、煤炭等矿产资源的开发。
2. 实习任务在实习期间,我主要负责地震资料的采集、处理和解释工作。
具体工作内容包括:(1)参与地震野外数据采集,了解地震仪器的使用和摆放要求;(2)学习地震数据处理技术,包括去噪、速度分析、震源定位等;(3)进行地震资料解释,识别地震反射界面,推断地下构造;(4)撰写实习报告,总结实习收获。
三、实习过程1. 地震数据采集在地震数据采集阶段,我了解了地震仪器的组成、使用方法和摆放要求。
通过实际操作,掌握了地震测线的设计、地震仪器的布置、激发和接收等环节。
此外,我还学会了如何根据地形、地质条件等因素优化测线布局,提高数据采集质量。
2. 地震数据处理在地震数据处理阶段,我学习了去噪、速度分析、震源定位等处理技术。
通过实际操作,掌握了这些技术在地震数据处理中的应用方法和注意事项。
此外,我还学会了如何评估处理结果的质量,以便为后续的解释工作提供可靠的数据基础。
3. 地震资料解释在地震资料解释阶段,我学会了如何识别地震反射界面,推断地下构造。
通过实际操作,掌握了地震剖面图的阅读方法,能够根据地震反射特征判断地层的岩性、厚度和分布规律。
此外,我还学会了如何利用地震资料进行断裂发育程度、油气藏评价等方面的研究。
四、实习收获通过这次实习,我收获颇丰。
首先,我掌握了物探基本原理和地震数据采集、处理、解释等关键技术;其次,我提高了自己的动手能力和实际操作能力,为今后从事物探工作打下了基础;最后,我加深了对物探行业的认识,了解了物探技术在矿产资源开发中的应用价值。
本科生实验报告实验课程物探数据处理学院名称地球物理学院专业名称勘查技术与工程学生姓名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}。