您好，您可以使用面填挖方的方法来实现您的功能，并不需要二次开发：
1.打开您的dem数据源，并将dem数据集加载到窗口地图
2.选择 分析>栅格分析>表面分析>面填挖方 ,跳出以下窗口：
3.在源数据中选择您的dem数据集，参考对象中选择您的水库的矢量面数据，当然如果没有矢量面的话可以手动绘制，输入附加高程，并设定缓冲区类型。
4.点击确定，在下方的输出窗口即可看见您的填充体积：
