海洋模拟源码解析：从Gerstner波到水面渲染全流程

拿到这套海洋模拟项目源码时，我做的第一件事不是跑Demo，而是直接翻目录结构和核心Shader。原因很简单——市面上的海洋模拟开源项目，十有八九要么是数学公式堆得漂亮但工程完整度极差，要么是美术效果唬人但稍微改两个参数就全线崩坏。这套源码能在标题里挂上“模拟项目”而不是“演示Demo”，本身说明作者对工程化是有要求的。看完之后，整体印象是结构干净、选型务实，适合做图形入门进阶或水体相关产品开发的人当作主线参考。

这篇文章我会从源码模块拆解讲起，把Gerstner波建模、顶点位移、法线计算、水面着色到性能优化这条完整管线逐个过一遍，同时穿插我在阅读和改造这套代码时踩过的问题，以及针对特定场景的取舍思路。无论你是第一次接触水体模拟，还是已经在做渲染优化但想看看别人怎么组织这套工程，这篇内容都应该对你有用。

1. 拿到源码先看什么：整体框架与技术选型

1.1 项目模块划分与数据流

好的源码一定在目录结构上就给你划好边界。这套海洋模拟项目大致分成了四大块：波浪核心算法层、渲染着色层、网格管理组件、工具与调试脚本。没有把一堆东西全塞进单文件里硬写，也没有过度抽象出十几层接口让你追代码追到崩溃。层与层之间的调用关系是单向的，数据流也很明确：CPU端根据时间累积和相机位置计算波浪参数，更新网格顶点，写入顶点缓冲；GPU端通过Shader对顶点做二次精修（比如细节波浪、顶点色扰动），再进行光照、反射、折射和雾效合成。

这个数据流设计里我最认可的一点是：CPU端只负责低频的、大尺度的波浪位移，高频细节全部交给GPU。这样做有两个直接好处，一是CPU开销小，同一帧内的Draw Call和同步压力低；二是可以做到顶点级波浪动态变化，同时不牺牲数量级很高的细节表现。很多刚入门的开发者容易犯的错误是试图在CPU端把所有波浪细节全算一遍，结果一扩大水面区域帧率直接跳水，这就是分工意识不足导致的。

1.2 技术路线：为什么是这个方案

海洋模拟的技术路线，通常绕不开三个流派：Gerstner波叠加、基于FFT的统计海浪模型（如Tessendorf的海浪频谱）、以及基于物理的浅水方程或光滑粒子流体动力学（SPH）。

这套源码选了Gerstner波作为主结构，而不是视觉效果上限更高的FFT方案，原因很实在：第一，FFT方案需要设计频谱模型并进行每帧的二维FFT计算，要么依赖GPU Compute Shader，要么在CPU上做离线预计算，对运行平台的兼容性要求很高；第二，FFT的调参难度明显高于Gerstner波，你很难直观地从一个频谱公式系数上判断海面最终会呈现什么样。Gerstner波则不同，每个波都有自己的方向、波长、振幅、速度，改一个参数就能直观看到对应波浪形态的变化，工程调试成本低得多。

当然这不代表Gerstner波没有短板。它的表现力上限在于“波峰”“波谷”的形态相对规律，大面积铺开时容易出现波浪走向重复感。这个项目用两个方式来缓解：一是叠加层数足够多，源码里默认是六层Gerstner波叠加，层间方向随机化；二是在Shader里加了一组高频细节扰动，贴片式地打碎规律感。这套组合方案在中等距离表现力上完全不输给某些只做单层FFT的引擎水体，而性能开销更可控，适合作为移动端或中端PC平台的主流水面方案。

2. 波浪生成器的核心实现：从数学公式到顶点位移

2.1 Gerstner波公式的源码表达

Gerstner波本质是一种周期性的表面行进波，它和普通正弦波最大的区别在于：粒子的运动不是简单的上下振动，而是在垂直方向和水平方向都有分量，形成一种“水质点绕平衡位置做圆轨迹运动”的效果。这套源码在核心计算上严格按照这个物理模型来展开。

以下是代码里核心的波形计算简化逻辑：

csharp复制// 单个Gerstner波对顶点位置的贡献
Vector3 GerstnerWave(float time, Vector3 vertex, WaveParams wave)
{
    float freq = 2.0f * Mathf.PI / wave.wavelength;   // 由波长求频率
    float phase = freq * (Vector3.Dot(wave.dir, vertex) + wave.speed * time);

    float sinVal = Mathf.Sin(phase);
    float cosVal = Mathf.Cos(phase);

    // 振幅（波高）在水平面法向方向的投影
    float steepness = wave.steepness / freq;

    vertex.x += wave.dir.x * steepness * cosVal;
    vertex.z += wave.dir.z * steepness * cosVal;
    vertex.y += (wave.amplitude / 1.0f) * sinVal;
    return vertex;
}

关键在于源码对变量steepness的处理，它将“陡峭度”从纯粹的视觉参数转化成了参与水平位移的量。如果你仔细推敲数学公式会发现，水平位移和垂直位移的比例直接决定了水面的陡峭程度。这个陡峭度不能无限调高，否则水面会出现大面积翻转和自交，形成灾难性的网格穿插。我实测下来，陡峭度控制在0.2到0.5之间是比较安全的范围，超过0.6时波浪底部会出现明显反折，效果从“海浪”变成“卷曲的金属箔片”，视觉和物理上都不可接受。

2.2 多波叠加与方向权重设计

多波叠加是这套源码里让水面“活”起来的核心手段。默认六层波形，每一层都有自己的独立参数空间，包括方向、波长、振幅、速度、陡峭度。方向不是随便给的，而是围绕海浪主风向产生扇形分布。比如主风向朝向东北，那其他波的方向就会在主风向左右各偏转30到60度的范围内随机生成。这样的方向分布能模拟出真实海面上浪向不完全一致、但又具有主传播方向的视觉效果。

有一点源码处理得很聪明：它对每一层波的力量权重做了归一化，避免所有层都用力过猛导致总位移爆炸。代码里是每加一层波形之后，都会将总位移和总振幅除以一个权重和，让最终水面高度保持在一个可预测的范围内。这个细节很多开源项目不做，带来的直接后果就是你往波形列表里塞五层、八层没问题，一旦塞到十几层，水面就会因为位移量累加而整体上下漂移，离开初始平面越来越远。

我在改造这套源码时额外加了一个需求：让不同区域的波高不一样，比如近岸比深海区更高一些。原项目的处理方式比较朴素，是全局统一参数。我通过一个自定义的波高衰减遮罩纹理，在合成总位移时乘以一个高度缩放系数，效果立竿见影，而且完全不影响原代码结构。综合来看，Gerstner波叠加的核心其实就是对“层数”和“参数空间”的把控，调得合理就能在性能和表现力之间取得不错的平衡。

2.3 法线计算的三种方式与源码选型

水面的光照效果完全依赖于法线质量，这个部分几乎决定了一个水面Shader的档次。这套源码在法线计算上提供了两条路径：一条是解析法直接通过波形函数的偏导数求法线；另一条是使用法线贴图扰动。

解析法的核心思路是：既然位移函数我们已经有了，那就直接对这个函数求梯度，从梯度推导出切线向量和副切线向量，再做叉积得到法线。这种方式得到的法线非常精准，尤其在近距离观察波浪轮廓时不会出现贴图法线的模糊感。代价是计算量较大，因为每个顶点都要对所有波形层求一次偏导。

源码的优化策略很值得学习：它把解析法用在了顶点Shader里，但只对最重要的低频主波计算，细节高频波不再参与解析法线计算，而是通过一张高频法线贴图在片元阶段继续扰动。这等于将“低频几何信息走顶点计算，高频细节信息走材质贴图”这一套经典的渲染管线分工贯彻到了水体渲染中。在场景表现上，近处水面波浪结构清晰，远处的碎浪由法线贴图承担，两不误。

glsl复制// 片元阶段法线扰动示例（简化）
float3 normal = normalize(vertexNormal);
float3 detailNormal = UnpackNormal(tex2D(_DetailNormalMap, uv));
normal = normalize(float3(normal + detailNormal));

这种混合法线方案我强烈建议直接照抄，因为它是水面项目中性价比最高的法线处理方式。全解析法太贵，纯贴图法在近距离容易破功，混合方案在绝大多数中等距离水面场景中都能取得良好效果。

3. 水面渲染瞒不住的细节：Shader中的水体着色

3.1 菲涅尔近似与反射折射权重

水看着像水，很大程度靠反射和折射的比例控制。这套源码的Shader里用的是Schlick菲涅尔近似，色值接近水和空气界面的真实光学行为，同时又能实时运算：

hlsl复制float fresnel = pow(1.0 - saturate(dot(viewDir, normal)), 5.0);
fresnel = lerp(_WaterReflectivity, 0.02, fresnel); // 0.02是折射主导时的最小菲涅尔值

这个公式和数据流设计让我觉得作者是真正调过水面渲染的：直接给固定菲涅尔系数是错得离谱的做法。真实水面的反射率随观察角度变化极大——垂直俯瞰水面时，水几乎是透明的，你能看到水底的泥沙；水平贴视水面时，水几乎变成一面镜子，反射出全部的天空和周围物体。Schlick近似把这个效应模拟得相当到位，而且只有一条GPU原生支持的指令，零额外开销。

这里的_WaterReflectivity是一个可调的基准反射率，通常设成0.02到0.05之间即可，配合一个经验性的pow指数来控制急变程度。在水面渲染的参数面板里，这是第一优先级调整的数值，很多新手先调颜色不调这个，结果水面怎么看都不通透。

3.2 深度差与近岸颜色过渡

很多海洋模拟项目把颜色直接写死成深蓝或浅蓝，一旦场景出现沙滩或浅滩区域，色彩的断裂感立刻显现。这套源码为此设计了一个深度采样流程：用带位置的深度纹理（Camera Depth Texture）采样当前位置的水深，然后把水深值映射成三个区间的颜色——浅水区偏青绿，中水区偏蓝，深水区则加入更多的深蓝和微小泡沫。

深度采样常见的坑是：没有把深度转换为线性深度就直接拿来算差分，导致浅水区颜色变化剧烈、深水区完全失去层次。源码在处理这个环节时，严格执行了从非线性深度到线性视空间的转换，再用线性深度做平滑插值。这一点虽然技术含量不算高，但恰恰是很多手写水面Shader容易出Bug的重灾区。我检查过不少项目，报“水面颜色一坨一坨的”问题的，十有八九就是这里出错了。

另外，它还在近岸处叠了一层高光增强。原理是当水深小于某个阈值时，不是直接压低颜色，而是额外增加一层反射高光亮度，用来模拟阳光穿透浅水后照亮近岸沙地的效果。这个细节虽然只有几行代码，但对视觉真实度的提升非常显著。

3.3 天空反射与向量贴图取巧

实时计算精确的天空反射是开销极高的方案，比如屏幕空间反射（SSR）或立方体贴图反射，都对平台性能有较高要求。这套源码选了一个更省的做法：用一张反射贴图配合粗糙度来模拟环境反射，本质上类似于镜面贴图的变体。

具体来说，它在水面片元阶段采样一张预计算的天空环境贴图（Cubemap），将采样结果乘以粗糙度系数来模糊反射内容，再和菲涅尔权重相乘输出。同时，它利用水面法线和视线方向计算反射向量，保证反射内容会随视角变化而变化，产生与真实反射近似的视差效果。

这种做法在开阔海面场景下能取得相当令人信服的效果，而在码头或者都市湖面场景下劣势会暴露——因为反射内容无法包含动态物体（船只、行人、建筑物），容易造成“水面里少了几栋楼”的违和感。所以方案是好方案，但选型时一定要看项目背景。如果要做的是海岛大世界项目，反射贴图完全够用；如果要做的是港口城市水面，就得考虑加一份实时反射相机或屏幕空间反射。

4. 性能与工程量产：LOD、实例化与纹理扰动

4.1 网格LOD切换与顶点密度策略

大面积水面的网格管理是整个项目里最容易被忽视的性能黑洞，没有之一。常见做法是给水面生成一个家用面积大小的网格，顶点均匀分布。这套源码没有那么做，它在网格初始化时根据距离相机的远近来分配顶点密度，近处密集，远处稀疏，实现了半径等差递减的LOD结构。原理并不复杂，关键是LOD级别的切换阈值要设置合理。

我在调试过程中发现一个很典型的问题：LOD切换阈值设太小，远处看不到什么波浪细节，却仍然分配了海量顶点，白白浪费渲染性能；设太大，相机移动或方向改变时，不同密度区域之间会产生明显的网格切换“跳变”。源码里通过给顶点级数设置平滑过渡因子，并在切换瞬间用高度差值做插值，成功掩盖了大部分突变。这是很成熟的工程处理手法，比很多只切换不插值的方案质感高出一截。

如果你要在自己的项目里直接套用这套LOD逻辑，我给一个经过实测的参数参考：近区顶点间距0.5米，覆盖半径25米；中区顶点间距2米，覆盖半径150米；远区顶点间距8到10米，覆盖半径500米以上。在这个配置下，水面场景的三角形数量可以压缩到合理范围，同时近处细节不缺席。

4.2 GPU Instancing与合批优化

水面的渲染批次这一块，稍不留神就会击穿Draw Call预算。这套源码自然避不开这个坑，好在它给每个水面区块都开启了GPU Instancing，配合Unity的SRP Batcher，在同一个材质属性的前提下，渲染多个相邻区块只需要一个Draw Call。这个能力在移动端表现尤为关键，因为移动端的Draw Call预算通常只有几十到一两百个，而水面场景动不动就需要几十个区块。

源码中一个值得借鉴的小技巧是：把水面的材质属性（波浪参数、颜色、法线贴图）统一存进MaterialPropertyBlock，以逐区块而非全局的方式下发。这样一来，每个区块可以局部调整波浪振幅（比如潮汐影响），同时又不破坏合批条件。很多人误以为一旦用MaterialPropertyBlock就一定会打断合批，实际上Unity对内置标准管线和SRP管线里使用属性块的水面渲染处理得相当好，只要不改变材质对象本身的渲染队列，这块开销是可以接受的。

4.3 法线贴图扰动替代全顶点动画

这套源码最值得夸奖的优化，是它在远处海面实现了一种“零顶点消耗”的波浪观感方案。它在片元阶段仅通过高频率法线贴图扰动，加一点点uv流动，就让远距离海面看起来仍然在动。远观察者通常辨别不出细节波形，但会很敏锐地感知“这块水面是不是静止的”。静态水面在视觉上会瞬间暴露“这是贴图而不是水面”的破绽。这个片元级别的扰动，近乎零成本地弥补了远距离的动感缺失。

再配合前面提到的解析法线思路，项目的整体性能模型是这样：近距离，CPU有少量波形位移计算，GPU有强度中等的顶点法和法线贴图混合；中距离，以法线贴图为主，波形计算降频；远距离，几乎不参与波形计算，只做低开销的贴图流动。性能模型阶梯清晰，是完整的量产级水面方案，而不只是跑个效果炫光的Demo。

5. 我踩过的坑和排查思路

5.1 常见问题速查表

下面把我在研究这套源码和改造过程中遇到的典型问题汇总成一张表，方便你排查时直接对照参考。

5.2 现场调试实录：一个法线翻转引发的“诡异高光”

我在把一个定制场景接入这套海洋模拟系统时，遇到了一个很诡异的现象：水面出现了大片跟随视角变化的异常高光大块，像是有水银液体在水面乱滚。初步排查时我首先怀疑是法线贴图方向错了，结果查了半天没发现异常。后来我把顶点输出的法线可视化出来，发现波浪陡峭度过高时，部分顶点的法线发生了翻转——波浪折叠导致法线指向了下方的水面。

这是一个非常微观但破坏力极大的问题。解决办法无非两种：一是限制陡峭度参数，并在Shader端增加一个法线“拉直”操作，把法线的垂直分量钳位到一个正范围内；二是降低多层波的叠加强度，尤其是方向互相冲突的波。我采取了第二种方案，把多变波的方向权重做了重新分布，并手动削弱了冲突层的振幅。调试完再看水面，那些光斑全部消失，波浪轮廓也恢复了立体感。

这个案例值得单独拿出来分享，是因为它说明了一个道理：海洋模拟这类项目，表面上拼的是Shader效果，实际上拼的是信号处理能力和对数值稳定性的把控。参数稍微越界一点，视觉立刻给你脸色看，而且是那种非常难以排查的脸色。

5.3 值得扩展的后续方向

如果你打算把这套源码当作自己项目的基础，而不是只看完就算，我建议从三个方向做扩展。第一，把Gerstner波叠加替换成FFT频谱方案，虽然调参难度提升，但效果上限明显更高，适合制作写实类产品。第二，给水面渲染加入地形交互，也就是在近岸区域让水面高度自动贴合海底地形采样，让浪在浅水区自然爬升和散开。第三，加入并行计算优化，把波形计算全部迁移到Compute Shader或Job System上，彻底释放CPU端压力。

我个人在实际操作中的体会是：这套源码给的架构足够稳，给的优化空间也足够大。波浪核心逻辑和渲染层是解耦的，你替换任何一层都不会牵连到另外一层。这一点在开源项目里反而是稀缺品质，大部分海洋模拟代码一改动Shader参数就得回去重新调顶点数据，改得人想砸键盘。

最后再分享一个小技巧：调试水面渲染时，把相机视角切到水面正上方往下看，最容易暴露颜色和法线的问题；切到水平视角则最容易暴露反射权重和轮廓的问题。两个角度交叉验证，能帮你快速定位八成以上的水面渲染Bug。这也是我调试这套源码时最常用的手段。

这篇内容希望能给正在啃海洋模拟源码的读者一些直接的启发。如果你手头也有这套代码，或者准备引入类似的方案，按我上面的路径走一遍，应该比盲目改参数靠谱得多。