参考解析中提到了因为D为椭圆形区域，所以使用广义极坐标进行计算。但具体的计算过程和结果并没有给出。需要知道具体的函数表达式、积分区域以及计算步骤才能给出准确的答案。