  We calculate the full spectrum of D-wave states in the Upsilon system in
lattice QCD for the first time, using an improved version of NonRelativistic
QCD on coarse and fine "second generation" gluon field configurations from the
MILC collaboration that include the effect of up, down, strange and charm
quarks in the sea. Taking the 2S-1S splitting to set the lattice spacing, we
determine the ^3D_2 - 1S splitting to 2.3%, and find agreement with experiment.
Our prediction of the fine structure relative to the ^3D_2 gives the ^3D_3 at
10.181(5) GeV and the ^3D_1 at 10.147(6) GeV. We also discuss the overlap of
^3D_1 operators with ^3S_1 states.
