For a Bose-Einstein condensate (BEG) confined in a double lattice consisting of two weak laser standing waves we find the Melnikov chaotic solution and chaotic region of parameter space by using the direct perturbation method. In the chaotic region, spatial evolutions of the chaotic solution and the corresponding distribution of particle number density are bounded but unpredictable between their superior and inferior limits. It is illustrated that when the relation fci fa k% between the two laser wave vectors is kept, the adjustment from k2 <k1 to k2 ≥ k1 can transform the chaotic region into regular one or the other way round. This suggests a feasible scheme for generating and controlling chaos, which could lead to an experimental observation in the near future.