A geothermal field is located in deep-seated basement aquifers in the northeastern part of the North China Plain near Tianjin, China. Carbonate rocks of Ordovician and Middle and Upper Proterozoic age on the Cangxian Uplift are capable of yielding 960 to 4200 [m.sup.3]/d of 57[degrees]C to 96[degrees]C water to wells from a depth of more than 1000 m. A three-dimensional nonisothermal numerical model was used to simulate and predict the spatial and temporal evolution of pressure and temperature in the geothermal system. The density of the geothermal water, which appears in the governing equations, can be expressed as a linear function of pressure, temperature, and total dissolved solids. A term describing the exchange of heat between water and rock is incorporated in the governing heat transport equation. Conductive heat flow from surrounding formations can be considered among the boundary conditions. Recent data of geothermal water production from the system were used for a first calibration of the numerical model. The calibrated model was used to predict the future changes in pressure and temperature of the geothermal water caused by two pumping schemes. The modeling results indicate that both pressure and temperature have a tendency to decrease with time and pumping. The current withdrawal rates and a pumping period of five months followed by a shut-off period of seven months are helpful in minimizing the degradation of the geothermal resource potential in the area.