数值模拟实验报告
中尺度数值模拟实习报告
实习报告实习单位：XX大学气象研究所实习时间：20XX年XX月XX日—20XX年XX月XX日实习内容：中尺度数值模拟一、实习背景及目的随着科学技术的不断发展，数值模拟技术在气象领域得到了广泛的应用。
中尺度数值模拟作为气象研究中的一种重要方法，可以有效地模拟和研究气象现象及其演变过程。
本次实习旨在通过实际操作，掌握中尺度数值模拟的基本原理和方法，提高自己在气象领域的实际工作能力。
二、实习内容及过程1. 实习前的准备工作在实习开始前，我参加了由导师组织的中尺度数值模拟培训课程，学习了中尺度数值模拟的基本原理、方法和实际应用。
此外，我还通过查阅相关文献，对中尺度数值模拟的研究现状和发展趋势有了更深入的了解。
2. 实习内容实习期间，我主要进行了以下几个方面的工作：（1）了解并熟悉中尺度数值模拟所使用的软件和硬件设备，掌握相关操作技巧。
（2）根据实际气象资料，设置数值模拟的初始条件和边界条件。
（3）运行中尺度数值模拟软件，观察模拟过程中气象参数的变化，并对模拟结果进行分析。
（4）根据模拟结果，撰写实习报告，总结中尺度数值模拟的特点和不足。
3. 实习过程在实习过程中，我首先选择了某一地区的一次降水过程作为研究对象。
根据实测气象资料，我设置了数值模拟的初始条件和边界条件，并使用中尺度数值模拟软件进行模拟。
在模拟过程中，我密切关注气象参数的变化，并对模拟结果进行了详细的分析。
通过实习，我了解到中尺度数值模拟在气象研究中的应用广泛，可以模拟各种气象现象，如降水、台风、沙尘暴等。
同时，中尺度数值模拟还能够为气象预报提供重要的参考依据。
三、实习收获及体会通过本次实习，我对中尺度数值模拟的基本原理和方法有了更加深入的了解，提高了自己的实际操作能力。
同时，实习过程中我对气象研究的重要性有了更深刻的认识，为自己的未来研究方向和方法的选择提供了宝贵的经验。
实习期间，我也发现中尺度数值模拟存在一些不足之处，如模拟结果受初始条件和边界条件的影响较大，有时难以准确反映实际气象状况。
中尺度气象学数值模拟实习报告
中尺度气象学数值模拟实习报告一、实习背景及目的近年来，随着计算机技术的快速发展，中尺度气象学数值模拟在气象预报、气候研究以及自然灾害预警等方面发挥着越来越重要的作用。
本次实习旨在通过实际操作，掌握中尺度气象学数值模拟的基本原理和方法，提高对气象现象的认识和预报能力。
二、实习内容与过程1. 实习前的准备在实习开始前，我们对中尺度气象学的基本概念、数值模拟的原理和方法进行了系统的学习。
通过查阅相关文献和资料，了解了中尺度气象学数值模拟在气象预报中的应用及其重要性。
2. 实习内容实习主要分为两个部分：一是中尺度气象学数值模拟的基本操作；二是对模拟结果进行分析与讨论。
（1）中尺度气象学数值模拟的基本操作我们使用了我国自主研发的中尺度气象数值模拟系统（MM5）进行实习。
实习过程中，我们根据实际气象资料，设置模拟的初始条件和边界条件，运行数值模拟系统，得到模拟气象场的分布。
（2）对模拟结果进行分析与讨论通过对模拟结果的对比分析，我们评估了数值模拟的准确性。
同时，针对模拟的气象场，分析了其中的中尺度特征，探讨了气象现象的产生和发展机制。
三、实习成果与反思1. 实习成果通过本次实习，我们掌握了中尺度气象学数值模拟的基本操作，了解了数值模拟在气象预报中的应用。
同时，通过对模拟结果的分析，提高了我们对气象现象的认识，为今后的气象研究工作奠定了基础。
2. 实习反思在实习过程中，我们深刻体会到中尺度气象学数值模拟在气象预报中的重要地位。
然而，由于气象现象的复杂性和不确定性，数值模拟仍存在一定的局限性。
因此，在实际应用中，我们需要不断优化模拟参数和模型，提高数值模拟的准确性，并结合其他气象观测资料，为气象预报和气候研究提供更为可靠的依据。
四、总结本次实习让我们对中尺度气象学数值模拟有了更深入的了解，提高了我们的实际操作能力和气象分析能力。
通过实习，我们认识到中尺度气象学数值模拟在气象预报和气候研究中的重要作用，也为今后的气象研究工作奠定了基础。
铸造过程数值模拟综合实验报告书
大连理工大学本科生实验报告课程名称：材料成型过程计算机应用实验学院（系）：材料科学与工程学院专业：材料成型与控制工程班级：学号：学生姓名：年月日铸造过程数值模拟综合实验报告要求铸造过程数值模拟综合实验课程共包含三个实验内容：激光点光源加热过程的ANSYS模拟、小方坯连铸过程数值模拟实验、铸造应力的PROCAST数值模拟实验。
实验前通过实验资料进行预习，实验结束得到的计算结果也上传到自己的邮箱里，课后将下载的适当图片粘贴到实验报告相应位置，并打印实验报告。
具体要求如下：1、六个实验全为必做实验，并完成实验报告（每个实验中的选做部分可根据个人条件选择完成，不做硬性规定）。
2、进行实验之前，学生需要进行课程预习，了解相关原理；同时参照实验说明书（下载打印）完成实验报告的前四项内容（手写）——实验目的、实验原理、实验设备和实验步骤。
实验老师对学生的预习情况进行检查。
实验完成后要填写实验报告中实验数据结果分析以及实验思考题中的相应内容。
3、实验过程中按照实验教材的实验步骤要求展开实验，如遇到难题可向实验指导老师或者助教咨询。
4、实验课资料的下载和实验结果的上传：实验过程中：实验中学生需在工作目录下新建“学号_实验名称”的文件夹，将拟上传的实验结果保存到此文件夹中。
允许同学将实验结果发送到自己邮箱中供自己下载。
5、实验报告的提交：根据教材的输出要求将下载的适当图片粘贴到实验报告相应位置，并打印实验报告，手写补充实验报告的其他内容。
在实验课结束后一周时间内，以班级为单位将实验报告统一提交到实验指导老师处。
6、未提交实验报告，或经查实属于抄袭实验报告或者编造原始数据者，实验课成绩以零分处理。
实验名称：激光点光源加热过程的ANSYS模拟实验目的实验原理和内容实验仪器设备（软件）实验步骤与操作方法实验数据记录与处理实验结果分析实验思考题实验名称：小方坯连铸过程数值模拟实验实验目的实验原理和内容实验仪器设备（软件）实验步骤与操作方法实验数据记录与处理实验结果分析实验思考题实验名称：铸造应力的PROCAST数值模拟实验实验目的实验原理和内容实验仪器设备（软件）实验步骤与操作方法实验数据记录与处理实验结果分析实验思考题10。
数值模拟沙盘认知实习报告
实习报告实习时间：2023年7月1日实习地点：XX公司数值模拟沙盘实验室实习组人员：组长、副组长、成员（包括：数据分析师、模型工程师、可视化工程师、项目经理）实习设备：数值模拟沙盘设备、计算机、投影仪、绘图板等一、实习目的1. 了解数值模拟沙盘的基本原理和应用领域；2. 学习数值模拟沙盘在实际项目中的应用流程；3. 提高自己在团队协作中的沟通与协作能力；4. 掌握一定的数值模拟沙盘操作技巧，为今后的工作积累经验。
二、实习内容1. 数值模拟沙盘的基本原理及操作方法；2. 数值模拟沙盘在工程项目、生产运营、市场营销等方面的应用；3. 团队协作完成一个实际项目的数值模拟沙盘演练；4. 针对模拟结果进行数据分析和总结。
三、实习过程1. 实习初期，我们参加了数值模拟沙盘的培训课程，学习了数值模拟沙盘的基本原理、操作方法和应用领域。
通过培训，我们对数值模拟沙盘有了初步的认识，为后续的实际操作打下了基础。
2. 在掌握了基本原理和操作方法后，我们开始在实验室进行实际操作。
首先，我们选择了一个小型的工程项目作为模拟对象，通过搭建沙盘模型、设定参数、运行模拟等步骤，完成了第一个模拟项目。
在此过程中，我们学会了如何运用数值模拟沙盘进行项目分析和管理。
3. 随着实习的深入，我们逐渐将模拟对象扩展到生产运营、市场营销等领域。
在团队合作中，我们充分发挥了各自的专业优势，共同完成了多个模拟项目。
通过这些项目的实践，我们提高了自己在团队协作中的沟通与协作能力，同时也积累了丰富的数值模拟沙盘操作经验。
4. 最后，我们对模拟结果进行了数据分析，总结了数值模拟沙盘在实际项目中的应用规律。
通过分析，我们发现数值模拟沙盘在预测项目风险、优化资源配置、提高运营效率等方面具有显著优势。
四、实习收获1. 掌握了数值模拟沙盘的基本原理和操作方法；2. 学会了如何在实际项目中运用数值模拟沙盘进行分析和管理；3. 提高了团队协作能力，学会了在团队中发挥各自优势，共同完成任务；4. 积累了数值模拟沙盘的实际操作经验，为今后的工作打下了基础。
泥浆沉积与速度关系数值模拟实验报告
泥浆沉积与速度关系数值模拟实验报告泥浆沉积是地质学中的一个重要研究领域，对于理解地质历史和资源勘探具有重要意义。
泥浆沉积速度是泥浆沉积过程中的一个重要参数，它直接影响着沉积物的分布和性质。
为了研究泥浆沉积与速度之间的关系，我们进行了数值模拟实验。
实验采用了计算机模拟的方法，模拟了泥浆在不同速度下的沉积过程。
首先，我们建立了一个沉积槽模型，模拟了泥浆在槽中的流动和沉积过程。
然后，我们通过改变泥浆的流速，观察了泥浆沉积的变化。
实验结果显示，泥浆的沉积速度与沉积物的分布和性质有着密切关系。
当泥浆的流速较慢时，沉积物呈现出均匀分布的趋势，沉积层较薄。
而当泥浆的流速增大时，沉积物呈现出不均匀分布的趋势，沉积层较厚。
同时，随着泥浆流速的增加，沉积物的粒度也会变大，沉积物中的颗粒越大，沉积层越粗。
进一步分析实验结果，发现泥浆沉积速度的增加会导致沉积物中的颗粒沉降速度增加，从而使得沉积物的分选效应更加明显。
而当泥浆流速较慢时，由于沉降速度较小，较小的颗粒更容易沉积，形成了均匀的沉积层。
实验还发现泥浆沉积速度与沉积物的结构和组成有关。
泥浆流速较快时，沉积物中的颗粒容易发生堆积，形成了较大的颗粒团。
而泥浆流速较慢时，沉积物中的颗粒更加分散，形成了较细的颗粒团。
总结以上实验结果，我们可以得出结论：泥浆沉积速度与沉积物的分布、性质和组成有着密切关系。
泥浆流速较慢时，沉积物呈均匀分布，颗粒较细；而泥浆流速较快时，沉积物呈不均匀分布，颗粒较大。
这一结论对于地质学研究和资源勘探具有重要意义，可以帮助我们更好地理解地质历史和开展资源勘探工作。
需要注意的是，本实验是基于数值模拟的方法得出的结果，与实际地质情况可能存在一定差异。
因此，在实际研究中，还需要结合野外观测和实验室分析等方法，进行综合研究，以得出更准确的结论。
通过数值模拟实验，我们研究了泥浆沉积与速度之间的关系。
实验结果表明，泥浆的沉积速度直接影响沉积物的分布、性质和组成。
流体力学数值模拟报告
流体力学数值模拟报告一、背景随着计算机技术的快速发展，数值模拟方法成为流体力学研究和工程设计中的重要工具。
传统的实验室试验和理论计算往往只能针对特定的情况进行研究，而数值模拟可以模拟各种复杂的流体问题，研究范围更广，具有更大的灵活性。
因此，流体力学数值模拟在科学研究、工程设计和产品优化中得到了广泛应用。
二、方法流体力学数值模拟通常采用计算流体力学（CFD）方法，通过将流体问题转化为求解计算网格中的运动方程来模拟流体运动。
其中，最常用的方法是有限体积法和有限元法。
有限体积法将流体区域分割为有限个体积单元，在每个单元上求解运动方程，得到整个流场的数值解。
有限元法则是通过在每个单元内定义一个函数，并在这个函数上求解运动方程，最后将每个单元上的解拼接在一起得到整个流场的数值解。
这些方法通常需要进行网格生成、边界条件处理和数值计算等工作。
三、应用流体力学数值模拟广泛应用于航空航天、汽车工程、建筑设计等领域。
例如，在航空航天领域，数值模拟可以预测飞机在不同速度下的空气动力学性能，优化机翼形状和发动机设计。
在汽车工程中，数值模拟可以模拟车辆行驶时的气动特性和燃烧过程，优化车身外形和发动机参数。
在建筑设计中，数值模拟可以预测建筑结构受风荷载时的动态响应，指导建筑物的抗风设计。
四、重要性流体力学数值模拟在工程领域具有重要的意义。
首先，它可以帮助工程师更好地理解和预测流体的行为，指导产品设计和性能优化。
其次，数值模拟可以大大减少试验成本和周期，提高工程设计的效率。
最后，数值模拟可以模拟各种复杂情况，例如高速气流、湍流和多相流等，为工程师提供更全面、准确的数据，帮助他们做出更好的决策。
五、未来发展随着计算机技术的进一步发展和数值模拟方法的日益成熟，流体力学数值模拟将在未来得到更广泛的应用。
首先，随着高性能计算机的出现，数值模拟可以处理更大规模的流体问题，提供更准确的结果。
其次，数据驱动的流体力学模拟方法将成为一个新的研究方向，通过机器学习和数据挖掘技术，将实验数据和模拟数据相结合，提高模拟结果的准确性和鲁棒性。
油藏数值模拟实验报告
目录1.前言1.1上机实践的目的及要求.1.2主要完成的实践内容...2.油藏特征分析2.1储层物性特征2.2流体物性特征2.3储层岩石物性特征3................................................. XX气藏数值模型建立3.1模型网格的划分3.2模型物性3.3模型流体性质及相渗曲线3.4XX气藏地质储量4.XX气藏方案优选4.1开发方案的优选4.2采速与稳产时间的关系5.结论认识5.1结论5.2对本实践课程的建议1 1 1 1 1 1 1 1 1 2 21.前言1.1上机实践的目的及要求1.掌握油藏数值模拟的上机操作流程；2.掌握ECLIPSE软件的数据录入、编辑和修改方法；3.掌握ECLIPSE软件结果输出及三维可视化方法；4.掌握机理模型研究方案设计的思路及方法1.2主要完成的实践内容1.油藏数值模拟数值整理；2.依据现有数据，应用块中心网络系统建立一个三维油藏数值模拟模型；3.预测单口气藏天然能量开发的最终采收率（20年）（不考虑水体能量）；4.预测多口气井采收率（20年）；5.预测不同稳产年限下，气井的合理产量（稳产5年）；6.水平井开发和直井开发效果对比；2.油藏特征分析2.1储层物性特征表2-1储层物性特征2.2流体物性特征表2-4岩石的PVT参数3.1模型网格的划分表2-2流体密度（地面条件下）表2-3水的PVT参数2.3储层岩石物性特征表3-2网格物性气藏X网格数Y网格数X步长Y步长470 50 80 100 100 X渗透率Y渗透率Z渗透率XYZ净毛比XYZ孔隙度0.95 0.95 0.095 0.1 0.15 3.2模型物性气藏X网格数Y网格数X步长Y步长470 50 80 100 100 3.3模型流体性质及相渗曲线SECONDCLASS_E10QFitEdit V KW Options-日目直砲已s ua3⑧型如sSGWFM(Gai-WalerSaiurationFunction时口口Kjg-w-*—SU闻—D.5O0.10X:■0--Z7SG7dfnen»onlH*.Y:0.43-512dtnenscHni3.4XX气藏地质储量iNa^igartianGirail>.dO 0-75L旳and4.XX气藏方案优选4.1开发方案的优选水平井方案水平井方案4Y4.H4.2采速与稳产时间的关系采油速度越快，稳产时间越短。
数值模拟报告：利用模型和计算预测结果
数值模拟报告：利用模型和计算预测结果引言：数值模拟在现代科学和工程领域中扮演着重要的角色。
借助数学模型和计算方法，数值模拟可以对复杂的现象和过程进行预测和分析。
基于已知的初始条件和边界条件，数值模拟可以得出一系列预测结果，为决策和规划提供参考。
本报告将介绍数值模拟的基本原理和方法，并通过具体案例阐述其在不同领域的应用。
1. 流体力学模拟1.1 模型基础在流体力学模拟中，最常用的模型是Navier-Stokes方程，它描述了流体在不同条件下的运动。
通过离散化和数值解法，我们可以得到流体的速度、压力、密度等关键参数的分布情况，从而预测流体流动的行为。
1.2 应用实例以风洞实验为例，我们可以利用数值模拟来预测空气在不同气流速度下对建筑物或车辆的压力分布，从而为建筑设计和风力发电规划提供有力的支持。
2. 电磁场模拟2.1 模型基础在电磁场模拟中，通过Maxwell方程组描述电磁场的分布和变化。
通过数值方法，我们可以得到电场、磁场、电流、电荷等关键信息的分布情况，进而揭示电磁场的特性。
2.2 应用实例以电子设备设计为例，我们可以利用数值模拟来预测电磁场对电路中信号传输的影响，优化电路布局和材料选择，提高电子设备的性能和可靠性。
3. 结构力学模拟3.1 模型基础结构力学模拟是通过求解弹性力学方程来分析结构的应力和变形情况。
通过数值方法，我们可以得到结构的位移、应力、应变等关键参数的分布情况，从而评估结构的稳定性和安全性。
3.2 应用实例以桥梁设计为例，我们可以利用数值模拟来预测桥梁在不同荷载下的应力分布和变形情况，为优化桥梁的结构和材料选择提供依据。
4. 生物医学模拟4.1 模型基础生物医学模拟是利用数学模型和计算方法对生物系统进行分析和预测。
通过建立生物系统的数学模型和参数化，我们可以模拟生物过程的动力学和变化，如细胞生长、药物传递等。
4.2 应用实例以药物研发为例，我们可以利用数值模拟来预测药物在人体内的分布与代谢，评估药物的治疗效果和安全性，加速药物研发过程。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
一、实验题目地震记录数值模拟的这几模型法 二、实验目的掌握褶积模型基本理论、实现方法与程序编制，由褶积模型初步分析地震 信号的分辨率问题 三、实验原理 1、褶积原理地震勘探的震源往往是带宽很宽的脉冲，在地下传播、反射、绕射到测线，传播经过中 高频衰减，能量被吸收。
吸收过程可以看成滤波的过程，滤波可以用褶积完成。
在滤波中， 反射系数与震源强弱关联，吸收作用与子波关联。
最简单的地震记录数值模拟，可以看成反 射系数与子波的褶积。
通常，反射系数是脉冲，子波取雷克子波。
（1） 雷克子波Wave(t) = (1−2π2f 2t 2)e −2π2f 2t 2（2） 反射系数：rflct(z)={1 z =z 反射界面0 z =others（3） 褶积公式：数值模拟地震记录trace(t): trace(t) =rflct(t)*wave(t)反射系数的参数由z 变成了t ，怎么实现？在简单水平层介质，分垂直和非垂直入射两 种实现，分别如图1 和图2 所示。
1） 垂直入射： t =2hv图一 垂直入射 2） 非垂直入射：t =2√h 2+x 2v图二非垂直入射 2、褶积方法（1） 离散化（数值化） 计算机数值模拟要求首先必须针对连续信号离散化处理。
反射系数在空间模型中存在，不同深度反射系数不同，是深度的 函数。
子波是在时间记录上一延续定时间的信号，是时间的概念。
在离散化时，通过深度采 样完成反射系数的离散化，通过时间采样完成子波的离散化。
如果记录是Trace(t)，则记录 是时间的函数，以时间采样离散化。
时间采样间距以∆t 表示，深度采样间距以∆z 表示。
在 做多道的数值模拟时，还有横向x 的概念，横向采样间隔以∆x 表示。
离散化的实现：t=It×∆t ；x=Ix×∆x ；z=Iz×∆z 或：It=t/∆t; Ix=x/∆x; Iz=z/∆z （2） 离散序列的褶积 trace (It )=∑rflct(Itao) ×wave(It −Itao)∞Itao=−∞ 四、实验内容1、 垂直入射地震记录数值模拟的褶积模型；2、 非垂直入射地震记录数值模拟的褶积模型；3、点绕射的地震记录数值模拟的褶积模型；五、方法路线根据褶积模型的实验原理编写C++程序，完成对于垂直入射波的褶积。
改变子波的长度与主频的大小，关注其对于实验结果的影响。
通过增加一个地层来模拟地下两层界面的反射情况，通过改变界面的高度来说明其对于实验结果的影响，同时改变子波长度与主频。
非垂直入射改变时间t来改变褶积结果显示地面情况。
点绕射模型通过时间的改变，任意位置检波器的自激自搜时间来改变褶积结果，同时改变子波长度与主频来分析影响。
六、实验结果：1.垂直入射：单层界面1>h=1000，fm=10，Nw=80v=2000v=3250 v=4500 v=5750v=70002>h=1000，v=4000，fm=10Nw=80Nw=100Nw=120Nw=140Nw=1603>h=1000，v=4000，Nw=80fm=5fm=10fm=15fm=20fm-464>h=1000，v=4000，Nw=160fm=5fm=10fm=15fm=20fm=465>v=4000，fm=10，Nw=80h=800h=1000h=1200 h=1400h=1600双层界面1>H12=20h=1000，v=4000，Nw=80fm=5fm=10fm=15fm=20fm=46h=1000，v=4000，Nw=160fm=5fm=10fm=15fm=20fm=462>H12=40h=1000，v=4000，Nw=80fm=5fm=10fm=15fm=20fm=46h=1000，v=4000，Nw=160fm=5fm=10fm=15fm=20fm=463>H12=60h=1000，v=4000，Nw=80fm=5fm=10fm=15fm=20fm=46h=1000，v=4000，Nw=160 fm=5fm=10fm=15fm=20fm=46H12=80h=1000，v=4000，Nw=80 fm=5fm=10fm=15fm=20fm=46h=1000，v=4000，Nw=160fm=5fm=10fm=15fm=20fm=462.非垂直入射1> h=1000，fm=10，Nw=80v=3000v=4000 v=5000 v=6000 v=7000 2> h=1000，v=4000，fm=10 Nw=80Nw=100Nw=120Nw=140Nw=1604>h=1000，v=4000，Nw=80fm=5fm=10fm=15fm=20fm=55 h=1000，v=4000，Nw=160 fm=5fm=10fm=15fm=20fm=55v=4000，fm=10，Nw=80 h=800h=1000 h=1200 h=1400 h=16003.点绕射卷积模型图件Nw=80 fm=15相关程序：#include<math.h>#include<stdio.h>#include<string.h>int Cnltn(float,float);int Rflct(float,float,float);int Wave(float,float);#define Nx 128#define Nt 256#define Nw 32void main(){floatdt=0.004,dx=20.0,fm=25.0,h=1000.0,v=3000.0;int iflag_Co,iflag_Re,iflag_Wv;if(iflag_Wv=Wave(fm,dt)!=1)printf("Waveiserror"); if(iflag_Re=Rflct(dt,h,v)!=1)printf("Reflectioniserror ");if(iflag_Co=Cnltn(dt,dx)!=1)printf("Convosioniserror ");}//Wave Formaing functionint Wave(float fm,float dt){FILE*fpw;int It;float Wa[Nw],t;double pai=3.1415926;if((fpw=fopen("wave.dat","wb"))==NULL) printf("Connotopenfile""wave""");for(It=0;It<Nw;It++){t=(float)(It+1)*dt;Wa[It]=(float)(1.-2*pai*pai*fm*fm*t*t)*(float) exp(-2*pai*pai*fm*fm*t); //雷克子波的计算？fwrite(&Wa[It],sizeof(Wa[It]),1,fpw);}fclose(fpw);return(1);}//Reflect Formaing functionint Rflct(float dt,float h,float v){FILE*fpr;int It,Ix,Ltdpth;float t;float Re[Nt];if((fpr=fopen("Reflect.dat","wb"))==NULL)printf("C onnot open file""Reflect""");for(Ix=0;Ix<Nx;Ix++){for(It=0;It<Nt;It++){Re[It]=0.;}t=(float)2.*h/v;Ltdpth=(int)(t/dt);Re[Ltdpth]=1.;for(It=0;It<Nt;It++){fwrite(&Re[It],sizeof(Re[It]),1,fpr);}}fclose(fpr);return(1);}//Convolution functionint Cnltn(float dt,float dx){FILE*fpc,*fpw,*fpr;int It,Ix,Itao;floatWa1[Nw],Wa[Nw],Re[Nt+Nw+Nw],Re1[Nt];float Con[Nt+Nw];if((fpc=fopen("Convosion.dat","wb"))==NULL) printf("Connotopenfile""Convosion""");if((fpw=fopen("wave.dat","rb"))==NULL)printf("Con notopenfile""wave""");if((fpr=fopen("Reflect.dat","rb"))==NULL)printf("Co nnotopenfile""Reflect""");for(Ix=1;Ix<2;Ix++){for(It=0;It<Nw;It++){fread(&Wa1[It],sizeof(Wa1[It]),1,fpw);Wa[Nw-1-It]=Wa1[It];}//子波数据反褶？}fclose(fpw);for(Ix=0;Ix<Nx;Ix++){for(It=0;It<Nt;It++){fread(&Re1[It],sizeof(&Re1[It]),1,fpr);}for(It=0;It<Nt+2*Nw;It++){Re[It]=0.;}for(It=0;It<Nt;It++){Re[It]=Re1[It+Nw];//反射系数数据移动:"0到Nt"移为"Nw到Nw+Nt"}for(It=0;It<Nt+Nw;It++){Con[It]=0;for(Itao=0;Itao<Nw;Itao++){Con[It]=Con[It]+Wa[Itao]*Re[-Itao+It];//褶积运算}}for(It=Nw/2;It<Nt+Nw/2;It++){fwrite(&Con[It],sizeof(Con[It]),1,fpc); }}fclose(fpw); fclose(fpr); fclose(fpc); return(1);}对比分析：1.垂直入射的结果是一条直线，反映了地下为水平界面；非垂直入射的结果是半条的双曲线，其反映了地下的非水平界面。